Double trouble: subsidence and CO<Subscript>2</Subscript> respiration due to 1,000 years of Dutch coastal peatlands cultivation

Coastal plains are amongst the most densely populated areas in the world. Many coastal peatlands are drained to create arable land. This is not without consequences; physical compaction of peat and its degradation by oxidation lead to subsidence, and oxidation also leads to emissions of carbon dioxide (CO₂). This study complements existing studies by quantifying total land subsidence and associated CO₂ respiration over the past millennium in the Dutch coastal peatlands, to gain insight into the consequences of cultivating coastal peatlands over longer timescales. Results show that the peat volume loss was 19.8 km³, which lowered the Dutch coastal plain by 1.9 m on average, bringing most of it below sea level. At least 66 % of the volume reduction is the result of drainage, and 34 % was caused by the excavation and subsequent combustion of peat. The associated CO₂ respiration is equivalent to a global atmospheric CO₂ concentration increase of ~0.39 ppmv. Cultivation of coastal peatlands can turn a carbon sink into a carbon source. If the path taken by the Dutch would be followed worldwide, there will be double trouble: globally significant carbon emissions and increased flood risk in a globally important human habitat. The effects would be larger than the historic ones because most of the cumulative Dutch subsidence and peat loss was accomplished with much less efficient techniques than those available now.     

Late Quaternary changes in surface productivity and oxygen minimum zone (OMZ) in the northwestern Arabian Sea: Micropaleontologic and sedimentary record at ODP site 728A

Changes in the abundance of selected planktic foraminiferal species and some sedimentological parameters at ODP site 728A were examined to understand the fluctuations in the surface productivity and deep sea oxygenation in the NW Arabian Sea during last ∼540 kyr. The increased relative abundances of high fertility taxa, i.e., Globigerinita glutinata and Globigerina bulloides mainly during interglacial intervals indicate intense upwelling. Strong SW summer monsoon probably increased the upwelling in the western Arabian Sea during interglacial intervals and caused high surface productivities due to the lateral transport of eutrophic waters. Most of the glacial periods (i.e., MIS 2, 4, 6, 8 and 12) are characterized by higher relative abundances of Neogloboquadrina pachyderma and Neogloboquadrina dutertrei associated with Globigerinoides ruber. The more stratified condition and deep mixed layer due to increased NE winter monsoon are mainly responsible for the higher relative abundances of N. pachyderma during glacial periods. Some of the glacial intervals (i.e., MIS 6 and 8) are also characterized by pteropod spikes reflecting deepening of aragonite compensation depth (ACD) and relatively less intense oxygen minimum zone (OMZ) in this region due to deep sea mixing and thermocline ventilation, and relatively less intense surface productivity during winter monsoon. The interglacial periods are largely devoid of pteropod shells indicating more aragonite dissolution due to increased intensity of OMZ in the northwestern Arabian Sea.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Determining the optimum cell size of digital elevation model for hydrologic application

Scale is one of the most important but unsolved issues in various scientific disciplines that deal with spatial data. The arbitrary choice of grid cell size for contour interpolated digital elevation models (DEM) is one of the major sources of uncertainty in the hydrologic modelling process. In this paper, an attempt was made to identify methods for determining an optimum cell size for a contour interpolated DEM in prior to hydrologic modelling. Twenty-meter interval contour lines were used to generate DEMs of five different resolutions, viz., 30, 45, 60, 75, and 90 m using TOPOGRID algorithm. The obtained DEMs were explored for their intrinsic quality using four different methods, i.e., sink analysis, fractal dimension of derived stream network, entropy measurement and semivariogram modelling. These methods were applied to determine the level artifacts (interpolation error) in DEM surface as well as derived stream network, spatial information content and spatial variability respectively. The results indicated that a 90 m cell size is sufficient to capture the terrain variability for subsequent hydrologic modelling in the study area. The significance of this research work is that it provides methods which DEM users can apply to select an appropriate DEM cell size in prior to detailed hydrologic modelling.     

Analytical model for screening potential CO<Subscript>2</Subscript> repositories

Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small.     

Crystal bending,subgrain boundary development,and recrystallization in orthopyroxene during granulite-facies deformation

A prominent feature of a granulite-facies shear zone from the Hidaka Main Zone (Japan) is the folding of orthopyroxene (opx) porphyroclasts. Dislocation density estimated by transmission electron microscope (TEM) and chemical etching in homogeneously folded domains is too low to account for the amplitude of crystallographic bending, leading us to propose a model similar to “flexural slip” folding, where folded layers are micrometer-wide opx layers between thin planar clinopyroxene (cpx) exsolutions. Extension (compression) in the extrados (intrados) of the folded layer is accommodated by dislocations at the cpx–opx interfaces. Alternatively to distributed deformation, crystal bending also localizes in grain boundaries (GBs), mostly oriented close to the (001) plane and with various misorientation angles but misorientation axes consistently close to the b-axis. For misorientation up to a few degrees, GBs were imaged as tilt walls composed of regularly spaced (100)[001] dislocations. For misorientation angles of 7°, individual dislocations are no longer visible, but high-resolution TEM (HRTEM) observation showed the partial continuity of opx tetrahedral chains through the boundary. For 21° misorientation, the two adjacent crystals are completely separated by an incoherent boundary. In spite of these atomic-scale variations, all GBs share orientation and rotation axis, suggesting a continuous process of misorientation by symmetric incorporation of (100)[001] dislocations. In addition to the dominant GBs perpendicular to the (100) plane, boundaries at low angle with (100) planes are also present, incorporating dislocations with a component of Burgers vector along the a-axis. The two kinds of boundaries combine to delimit subgrains, which progressively rotate with respect to host grains around the b-axis, eventually leading to recrystallization of large porphyroclasts.     

Evaluation of four supervised learning methods for groundwater spring potential mapping in Khalkhal region (Iran) using GIS-based features

One important tool for water resources management in arid and semi-arid areas is groundwater potential mapping. In this study, four data-mining models including K-nearest neighbor (KNN), linear discriminant analysis (LDA), multivariate adaptive regression splines (MARS), and quadric discriminant analysis (QDA) were used for groundwater potential mapping to get better and more accurate groundwater potential maps (GPMs). For this purpose, 14 groundwater influence factors were considered, such as altitude, slope angle, slope aspect, plan curvature, profile curvature, slope length, topographic wetness index (TWI), stream power index, distance from rivers, river density, distance from faults, fault density, land use, and lithology. From 842 springs in the study area, in the Khalkhal region of Iran, 70 % (589 springs) were considered for training and 30 % (253 springs) were used as a validation dataset. Then, KNN, LDA, MARS, and QDA models were applied in the R statistical software and the results were mapped as GPMs. Finally, the receiver operating characteristics (ROC) curve was implemented to evaluate the performance of the models. According to the results, the area under the curve of ROCs were calculated as 81.4, 80.5, 79.6, and 79.2 % for MARS, QDA, KNN, and LDA, respectively. So, it can be concluded that the performances of KNN and LDA were acceptable and the performances of MARS and QDA were excellent. Also, the results depicted high contribution of altitude, TWI, slope angle, and fault density, while plan curvature and land use were seen to be the least important factors.     

Evaluation of bank filtration as a pretreatment method for the provision of hygienically safe drinking water in Norway: results from monitoring at two full-scale sites

Two case studies were carried out in central Norway in order to assess the performance of bank filtration systems in cold-climate fluvial aquifers relying on recharge from humic-rich surface waters with moderate microbial contamination. Three municipal wells and two surface-water sources at operative bank filtration systems were monitored for naturally occurring bacteriophages, fecal indicators, natural organic matter (NOM) and physico-chemical water quality parameters during a 4-month period. Aquifer passage effectively reduced the microorganism and NOM concentrations at both study sites. Bacteriophages were detected in 13 of 16 (81%) surface-water samples and in 4 of 24 (17%) well-water samples, and underwent 3 ± 0.3 log₁₀ reduction after 50–80-m filtration and 20–30 days of subsurface passage. NOM reductions (color: 74–97%; dissolved organic carbon: 54–80%; very hydrophobic acids: 70%) were similar to those achieved by conventional water-treatment processes and no further treatment was needed. Both groundwater dilution and sediment filtration contributed to the hygienic water quality improvements, but sediment filtration appeared to be the most important process with regard to microbial and NOM reductions. A strengths-weaknesses-opportunities-threats analysis showed that bank filtration technology has a high potential as a pretreatment method for the provision of hygienically safe drinking water in Norway.     

Investigating the impact of the properties of pilot points on calibration of groundwater models: case study of a karst catchment in Rote Island,Indonesia

A robust configuration of pilot points in the parameterisation step of a model is crucial to accurately obtain a satisfactory model performance. However, the recommendations provided by the majority of recent researchers on pilot-point use are considered somewhat impractical. In this study, a practical approach is proposed for using pilot-point properties (i.e. number, distance and distribution method) in the calibration step of a groundwater model. For the first time, the relative distance–area ratio (d/A) and head-zonation-based (HZB) method are introduced, to assign pilot points into the model domain by incorporating a user-friendly zone ratio. This study provides some insights into the trade-off between maximising and restricting the number of pilot points, and offers a relative basis for selecting the pilot-point properties and distribution method in the development of a physically based groundwater model. The grid-based (GB) method is found to perform comparably better than the HZB method in terms of model performance and computational time. When using the GB method, this study recommends a distance–area ratio of 0.05, a distance–x-grid length ratio (d/X _grid) of 0.10, and a distance–y-grid length ratio (d/Y _grid) of 0.20.