Review: From multi-scale conceptualization to a classification system for inland groundwater-dependent ecosystems

Aquifers provide water, nutrients and energy with various patterns for many aquatic and terrestrial ecosystems. Groundwater-dependent ecosystems (GDEs) are increasingly recognized for their ecological and socio-economic values. The current knowledge of the processes governing the ecohydrological functioning of inland GDEs is reviewed, in order to assess the key drivers constraining their viability. These processes occur both at the watershed and emergence scale. Recharge patterns, geomorphology, internal geometry and geochemistry of aquifers control water availability and nutritive status of groundwater. The interface structure between the groundwater system and the biocenoses may modify the groundwater features by physicochemical or biological processes, for which biocenoses need to adapt. Four major types of aquifer-GDE interface have been described: springs, surface waters, peatlands and terrestrial ecosystems. The ecological roles of groundwater are conditioned by morphological characteristics for spring GDEs, by the hyporheic zone structure for surface waters, by the organic soil structure and volume for peatland GDEs, and by water-table fluctuation and surface floods in terrestrial GDEs. Based on these considerations, an ecohydrological classification system for GDEs is proposed and applied to Central and Western-Central Europe, as a basis for modeling approaches for GDEs and as a tool for groundwater and landscape management.     

CAT-scan analysis of sedimentary sequences: An ultrahigh-resolution paleoclimatic tool

Paleoclimate research is essential to determine the natural variability of climate and to place the current climate change into its natural context. The current need is to generate the highest temporal resolution paleoclimatic reconstructions possible in order to assess the natural variability of the climate system, but also to test the ability of numerical models to simulate conditions different from the ones observed with the relatively short instrumental records. In this paper, we show that CAT-scan analysis of sedimentary sequences, with its 1 mm downcore resolution, can be used to identify millennial to seasonal cycles in sedimentary sequences. In examples from the St. Lawrence Estuary, Eastern Canada, spectral analysis of the CAT-scan data from Holocene postglacial sediments revealed millennial- to centennial-scale oscillations possibly associated with either solar variability, changes in relative sea-level or tidal amplitude. Similarly, spectral analysis of Holocene and Sangamonian glaciomarine sequences revealed decadal- to annual-scale oscillations with periods close to the one previously associated with the North Atlantic Oscillation (NAO), whereas spectral analysis of the CAT-scan data from the Sangamonian rhythmites possibly revealed seasonal cycles.     

Influence of the Laurentide Ice Sheet and relative sea‐level changes on sediment dynamics in the Estuary and Gulf of St. Lawrence since the last deglaciation

Physical properties, grain size, bulk mineralogy, elemental geochemistry and magnetic parameters of three sediment piston cores recovered in the Laurentian Channel from its head to its mouth were investigated to reconstruct changes in detrital sediment provenance and transport related to climate variability since the last deglaciation. The comparison of the detrital proxies indicates the succession of two sedimentary regimes in the Estuary and Gulf of St. Lawrence (EGSL) during the Holocene, which are associated with the melting history of the Laurentide Ice Sheet (LIS) and relative sea‐level changes. During the early Holocene (10–8.5 cal. ka BP), high sedimentation rates together with mineralogical, geochemical and magnetic signatures indicate that sedimentation in the EGSL was mainly controlled by meltwater discharges from the local retreat of the southeastern margin of the LIS on the Canadian Shield. At this time, sediment‐laden meltwater plumes caused the accumulation of fine‐grained sediments in the ice‐distal zones. Since the mid‐Holocene, postglacial movements of the continental crust, related to the withdrawal of the LIS (c. 6 cal. ka BP), have triggered significant variations in relative sea level (RSL) in the EGSL. The significant correlation between the RSL curves and the mineralogical, geochemical, magnetic and grain‐size data suggest that the RSL was the dominant force acting on the sedimentary dynamics of the EGSL during the mid‐to‐late Holocene. Beyond 6 cal. ka BP, characteristic mineralogical, geochemical, magnetic signatures and diffuse spectral reflectance data suggest that the Canadian Maritime Provinces and western Newfoundland coast are the primary sources for detrital sediments in the Gulf of St. Lawrence, with the Canadian Shield acting as a secondary source. Conversely, in the lower St. Lawrence Estuary, detrital sediments are mainly supplied by the Canadian Shield province. Finally, our results suggest that the modern sedimentation regime in the EGSL was established during the mid‐Holocene.     

Soil charcoal stability over the Holocene across boreal northeastern North America

The analysis of macroscopic wood charcoal fragments extracted from soils is frequently used as a palaeoecological tool for reconstructing stand-scale forest composition and fire history. Here we explored the putative loss of palaeoecological information due to charcoal degradation through time and in different biogeographical settings. We compared the relationship between charcoal mass and abundance for soil samples from five biogeographical regions of boreal northeastern North America spanning most of the Holocene period. We verified whether charcoal (Ø ≥ 2 mm) conservation differed as a consequence of different taphonomical processes between organic and mineral soil types. We also assessed the mass/abundance relationship as a function of charcoal residence time in soil. Overall, the slope of the regression between charcoal particles mass (g) and abundance (number of particles) was 0.0042. The slope was not significantly different in samples from organic and mineral soil, and all biogeographical regions had similar slope values except one (higher charcoal fragmentation, probably due to high colluvial activity). Charcoal conservation also did not vary according to residence time in soil. This study shows that macroscopic soil charcoal particles resist fragmentation over millennia in different biogeographical settings and under the influence of various taphonomical processes.     

Bentho-planktonic evidence from the Austrian Alps for a decline in sea-surface carbonate production at the end of the Triassic

A high-resolution micropalaeontological study, combined with geochemical and sedimentological analyses was performed on the Tiefengraben, Schlossgraben and Eiberg sections (Austrian Alps) in order to characterize sea-surface carbonate production during the end-Triassic crisis. At the end-Rhaetian, the dominant calcareous nannofossil Prinsiosphaera triassica shows a decrease in abundance and size and this is correlated with a increase in δ¹⁸O and a gradual decline in δ¹³C_carb values. Simultaneously, benthic foraminiferal assemblages show a decrease in diversity and abundance of calcareous taxa and a dominance of infaunal agglutinated taxa. The smaller size of calcareous nannofossils disturbed the vertical export balance of the biological carbon pump towards the sea-bottom, resulting in changes in feeding strategies within the benthic foraminiferal assemblages from deposit feeders to detritus feeders and bacterial scavengers. These micropalaeontological data combined with geochemical proxies suggest that changes in seawater chemistry and/or cooling episodes might have occurred in the latest Triassic, leading to a marked decrease of carbonate production. This in turn culminated in the quasi-absence of calcareous nannofossils and benthic foraminifers in the latest Triassic. The aftermath (latest Triassic earliest Jurassic) was characterised by abundance peaks of “disaster” epifaunal agglutinated foraminifera Trochammina on the sea-floor. Central Atlantic Magmatic Province (CAMP) paroxysmal activity, superimposed on a major worldwide regressive phase, is assumed to be responsible for a deterioration in marine palaeoenvironments. CAMP sulfuric emissions might have been the trigger for cooling episodes and seawater acidification leading to disturbance of the surface carbonate production at the very end-Triassic.     

The prehistoric Mt Wilberg rock avalanche, Westland, New Zealand

The Mt Wilberg rock avalanche in Westland, New Zealand occurred before 1300 AD and may have occurred as a consequence of an Alpine fault earthquake in ca. 1220 AD or earlier. Its ∼40 × 10⁶ m³ deposit may have briefly obstructed the Wanganui River, but only about 25% of its surface morphology still survives, on terraces isolated from river erosion. The landslide appears to have moved initially as a block, in a direction controlled by a strong rock mass at the base of the source area, before disintegrating and spreading across terraces, fans, and floodplains. Rock avalanche deposits in Westland have relatively short expected lifetimes in the rugged terrain and high rainfall of the area; hence, the hazard from such events is under-represented by their current remnants.     

The Geoid Slope Validation Survey 2014 and GRAV-D airborne gravity enhanced geoid comparison results in Iowa

Three Geoid Slope Validation Surveys were planned by the National Geodetic Survey for validating geoid improvement gained by incorporating airborne gravity data collected by the “Gravity for the Redefinition of the American Vertical Datum” (GRAV-D) project in flat, medium and rough topographic areas, respectively. The first survey GSVS11 over a flat topographic area in Texas confirmed that a 1-cm differential accuracy geoid over baseline lengths between 0.4 and 320 km is achievable with GRAV-D data included (Smith et al. in J Geod 87:885–907, 2013). The second survey, Geoid Slope Validation Survey 2014 (GSVS14) took place in Iowa in an area with moderate topography but significant gravity variation. Two sets of geoidal heights were computed from GPS/leveling data and observed astrogeodetic deflections of the vertical at 204 GSVS14 official marks. They agree with each other at a \({\pm }1.2\,\, \hbox {cm}\) level, which attests to the high quality of the GSVS14 data. In total, four geoid models were computed. Three models combined the GOCO03/5S satellite gravity model with terrestrial and GRAV-D gravity with different strategies. The fourth model, called xGEOID15A, had no airborne gravity data and served as the benchmark to quantify the contribution of GRAV-D to the geoid improvement. The comparisons show that each model agrees with the GPS/leveling geoid height by 1.5 cm in mark-by-mark comparisons. In differential comparisons, all geoid models have a predicted accuracy of 1–2 cm at baseline lengths from 1.6 to 247 km. The contribution of GRAV-D is not apparent due to a 9-cm slope in the western 50-km section of the traverse for all gravimetric geoid models, and it was determined that the slopes have been caused by a 5 mGal bias in the terrestrial gravity data. If that western 50-km section of the testing line is excluded in the comparisons, then the improvement with GRAV-D is clearly evident. In that case, 1-cm differential accuracy on baselines of any length is achieved with the GRAV-D-enhanced geoid models and exhibits a clear improvement over the geoid models without GRAV-D data. GSVS14 confirmed that the geoid differential accuracies are in the 1–2 cm range at various baseline lengths. The accuracy increases to 1 cm with GRAV-D gravity when the west 50 km line is not included. The data collected by the surveys have high accuracy and have the potential to be used for validation of other geodetic techniques, e.g., the chronometric leveling. To reach the 1-cm height differences of the GSVS data, a clock with frequency accuracy of \(10^{-18}\) is required. Using the GSVS data, the accuracy of ellipsoidal height differences can also be estimated.     

Voronoi grids conforming to 3D structural features

Flow simulation in a reservoir can be highly impacted by upscaling errors. These errors can be reduced by using simulation grids with cells as homogeneous as possible, hence conformable to horizons and faults. In this paper, the coordinates of 3D Voronoi seeds are optimized so that Voronoi cell facets honor the structural features. These features are modeled by piecewise linear complex (PLC). The optimization minimizes a function made of two parts: (1) a barycentric function, which ensures that the cells will be of good quality by maximizing their compactness; and (2) a conformity function, which allows to minimize the volume of cells that is isolated from the Voronoi seed w.r.t., a structural feature. To determine the isolated volume, a local approximation of the structural feature inside the Voronoi cells is used to cut the cells. It improves the algorithm efficiency and robustness compared to an exact cutting procedure. This method, used jointly with an adaptive gradient solver to minimize the function, allows dealing with complex 3D geological cases. It always produces a Voronoi simulation grid with the desired number of cells.     

Estimate of recharge of a rising water table in semiarid Niger from 3H and 14C modeling

A hydrodynamic survey carried out in semiarid southwest Niger revealed an increase in the unconfined ground water reserves of approximately 10% over the last 50 years due to the clearing of native vegetation. Isotopic samplings (3H, 18O, 2H for water and 14C, 13C for the dissolved inorganic carbon) were performed on about 3500 km2 of this silty aquifer to characterize recharge. Stable isotope analyses confirmed the indirect recharge process that had already been shown by hydrodynamic surveys and suggested the tracers are exclusively of atmospheric origin. An analytical model that takes into account the long-term rise in the water table was used to interpret 3H and 14C contents in ground water. The natural, preclearing median annual renewal rate (i.e., recharge as a fraction of the saturated aquifer volume) lies between 0.04% and 0.06%. For representative characteristics of the aquifer (30 m of saturated thickness, porosity between 10% and 25%), this implies a recharge of between 1 and 5 mm/year, which is much lower than the estimates of 20 to 50 mm/year for recent years, obtained using hydrological and hydrodynamic methods and the same aquifer parameters. Our study, therefore, reveals that land clearing in semiarid Niger increased ground water recharge by about one order of magnitude.