Comparative hydrogeology – reference analysis of groundwater dynamics from neighbouring observation wells

ABSTRACT

Groundwater level fluctuations are caused by spatial and temporal superposition of processes within and outside the aquifer system. Most of the subsurface processes are usually observed on a small scale. Upscaling to the regional scale, as required for future climate change scenarios, is difficult due to data scarcity and increasing complexity. In contrast to the limited availability of system characteristics, high-resolution data records of groundwater hydrographs are more generally available. Exploiting the information contained in these records should thus be a priority for analysis of the chronical lack of data describing groundwater system characteristics. This study analyses the applicability of 63 indices derived from daily hydrographs to quantify different dynamics of groundwater levels in unconfined gravel aquifers from three groundwater regions (Bavaria, Germany). Based on the results of two different skill tests, the study aids index selection for different dynamic components of groundwater hydrographs.     

Exploring large wood retention and deposition in contrasting river morphologies linking numerical modelling and field observations

Large wood tends to be deposited in specific geomorphic units within rivers. Nevertheless, predicting the spatial distribution of wood deposits once wood enters a river is still difficult because of the inherent complexity of its dynamics. In addition, the lack of long‐term observations or monitored sites has usually resulted in a rather incomplete understanding of the main factors controlling wood deposition under natural conditions. In this study, the deposition of large wood was investigated in the Czarny Dunajec River, Polish Carpathians, by linking numerical modelling and field observations so as to identify the main factors influencing wood retention in rivers. Results show that wood retention capacity is higher in unmanaged multi‐thread channels than in channelized, single‐thread reaches. We also identify preferential sites for wood deposition based on the probability of deposition under different flood scenarios, and observe different deposition patterns depending on the geomorphic configuration of the study reach. In addition, results indicate that wood is not always deposited in the geomorphic units with the highest roughness, except for low‐magnitude floods. We conclude that wood deposition is controlled by flood magnitude and the elevation of flooded surfaces in relation to the low‐flow water surface. In that sense, the elevation at which wood is deposited in rivers will differ between floods of different magnitude. Therefore, together with the morphology, flood magnitude represents the most significant control on wood deposition in mountain rivers wider than the height of riparian trees. Copyright © 2015 John Wiley & Sons, Ltd.     

Diagnostic analysis of future climate scenarios applied to urban flooding in the Chicago metropolitan area

Past heavy precipitation events in the Chicago metropolitan area have caused significant flood-related economic and environmental damages. A key component in flood management policies and actions is determining flood magnitudes for specified return periods. This is a particularly difficult task in areas with a complex and changing climate and land-use, such as the Chicago metropolitan area. The standard design storm methodology based on the NOAA Atlas 14 and ISWS Bulletin 70 has been used in the past to estimate flood hydrographs with variable return periods in this region. In a changing climate, however, these publications may not be accurate. This study presents and illustrates a methodology for diagnostic analysis of future climate scenarios in the framework of urban flooding, and assesses the corresponding uncertainties. First, the design storms are calculated using data downscaled by a regional climate model (RCM) at 30-km spacing for the present and 2050s under the IPCC A1Fi (high) and B1 (low) emissions scenarios. Next, the corresponding flood discharges at six watersheds in suburban Chicago are estimated using a hydrologic event model. The resulting scenarios in flood frequency were first assessed through a set of diagnostic tests for precipitation timing and frequency. The study did not reveal any significant changes in the 2050s in the average timing of heavy storms, but their regularity decreased. The average timing did not exhibit any significant spatial variability throughout the region. The precipitation frequency analysis revealed distinct differences between the northern and southeastern subregions of the Chicago metropolitan area. The quantiles in the northern subregion averaged for 2-year, 5-year, and 10-year return periods exhibited a 20% and 16% increase in daily precipitation for scenarios B1 and A1Fi, respectively. The southeastern subregion, however, exhibited a decrease of 12% for scenario B1 and a minor increase of 3% for scenario A1Fi. The hydrologic effects of changing precipitation on the flood quantiles were illustrated using six small watersheds in the region. The relative increases or decreases in precipitation translated into even larger relative increases or decreases in flood peaks, due to the nonlinear nature of the rainfall-runoff process. Simulations using multiple climate models, for longer periods, finer spatio-temporal resolution, and larger areal coverage could be used to more accurately account for numerous uncertainties in the precipitation and flood projections.     

A groundwater salinity hotspot and its connection to an intermittent stream identified by environmental tracers (Mt Lofty Ranges,South Australia)

High and variable levels of salinity were investigated in an intermittent stream in a high-rainfall area (～800 mm/year) of the Mt. Lofty Ranges of South Australia. The groundwater system was found to have a local, upslope saline lens, referred to here as a groundwater salinity ‘hotspot’. Environmental tracer analyses (δ¹⁸O, δ²H, ^87/86Sr, and major elements) of water from the intermittent stream, a nearby permanent stream, shallow and deep groundwater, and soil-water/runoff demonstrate seasonal groundwater input of very saline composition into the intermittent stream. This input results in large salinity increases of the stream water because the winter wet-season stream flow decreases during spring in this Mediterranean climate. Furthermore, strontium and water isotope analyses demonstrate: (1) the upslope-saline-groundwater zone (hotspot) mixes with the dominant groundwater system, (2) the intermittent-stream water is a mixture of soil-water/runoff and the upslope saline groundwater, and (3) the upslope-saline-groundwater zone results from the flushing of unsaturated-zone salts from the thick clayey regolith and soil which overlie the metamorphosed shale bedrock. The preferred theory on the origin of the upslope-saline-groundwater hotspot is land clearing of native deep-rooted woodland, followed by flushing of accumulated salts from the unsaturated zone due to increased recharge. This cause of elevated groundwater and surface-water salinity, if correct, could be widespread in Mt. Lofty Ranges areas, as well as other climatically and geologically similar areas with comparable hydrogeologic conditions.     

First VLTI observations of Mira stars

We present first observations of Mira stars obtained at the Very Large Telescope Interferometer (VLTI) at 2 microns, equipped with the VINCI instrument, using different baselines with both the test siderostats and the 8.2 m unit telescopes. These observations, collected in the course of the VLTI commissining program, have targeted so far about sixty cool giant stars. In this talk, we list and describe the measurements obtained for 14 Mira stars, many of them representing first-time determinations of the visibility. In particular, we devote special attention to a study of Mira itself, for which very accurate determinations of the visibility at several baseline lenghts and orientations could be obtained. We find that a two-component model is more consistent with the data than a single stellar disc. Further observations are needed for a better understanding of this source. The VLTI will constitute an ideal instrument for such studies in the future. In particular it will permit detailed investigations of southern AGB stars, such as accurate measurements of surface structureparameters (diameters, diameter variations, asymmetries, center-to-limbvariations, special features like hot spots) and of circumstellar envelopes. This revised version was published online in July 2006 with corrections to the Cover Date.     

Nuclear and incommensurate magnetic structure of NaFeGe2O6 between 5?K and 298?K and new data on multiferroic NaFeSi2O6

The compound NaFeGe₂O₆ was grown synthetically as polycrystalline powder and as large single crystals suitable for X-ray and neutron-diffraction experiments to clarify the low temperature evolution of secondary structural parameters and to determine the low temperature magnetic spins structure. NaFeGe₂O₆ is isotypic to the clinopyroxene-type compound aegirine and adopts the typical HT-C2/c clinopyroxene structure down to 2.5?K. The Na-bearing M2 polyhedra were identified to show the largest volume expansion between 2.5?K and room temperature, while the GeO₄ tetrahedra behave as stiff units. Magnetic susceptibility measurements show a broad maximum around 33?K, which marks the onset of low-dimensional magnetic ordering. Below 12?K NaFeGe₂O₆ transforms to an incommensurately modulated magnetic spin state, with k?=?[0.323, 1.0, 0.080] and a helical order of spins within the M1-chains of FeO₆ octahedra. This is determined by neutron-diffraction experiments on a single crystal. Comparison of NaFeGe₂O₆ with NaFeSi₂O₆ is given and it is shown that the magnetic ordering in the latter compound, aegirine, also is complex and is best described by two different spin states, a commensurate one with C2??/c?? symmetry and an incommensurate one, best being described by a spin density wave, oriented within the (1 0 1) plane.     

Reconstructing spatio-temporal patterns of debris-flow activity using dendrogeomorphological methods

Debris flows are a major threat in many parts of the Alps, where they repeatedly cause severe damage to infrastructure and transportation corridors or even loss of life. Nonetheless, the spatial behavior of past debris-flow activity and the analysis of areas affected during particular events have been widely neglected in reconstructions so far. It was therefore the purpose of this study to reconstruct spatio-temporal patterns of past debris flows on a forested cone in the Swiss Alps (Bruchji torrent, Blatten, Valais). The analysis of past events was based on a detailed geomorphic map (1:1000) of all forms related to debris flows as well as on tree-ring series from 401 heavily affected trees (Larix decidua Mill. and Picea abies (L.) Karst.) growing in or next to deposits. The samples were analyzed and growth disturbances related to debris-flow activity assessed, such as tangential rows of traumatic resin ducts, the onset of reaction wood or abrupt growth suppression or release.In total, 960 growth disturbances were identified in the samples, belonging to 40 different event years between A.D. 1867 and 2005. In addition, the coupling of tree-ring data with the geomorphic map allowed reconstruction of eleven formerly active channels and spatial representation of individual events. Based on our results we believe that before 1935, debris flows preferentially used those channels located in the western part of the cone, whereas the eastern part of the cone remained widely unaffected. The spatial representation of the 40 events also allowed identification of five different spatial patterns for debris flows at the study site.     

Projected climate change impact on Baltic Sea cyanobacteria

Compared to other phytoplankton groups, nitrogen-fixing cyanobacteria generally prefer high water temperatures for growth and are therefore expected to benefit from global warming. We use a coupled biological-physical model with an advanced cyanobacteria life cycle model to compare the abundance of cyanobacteria in the Baltic Sea during two different time periods (1969–1998; 2069–2098). For the latter, we find prolonged growth and a more than twofold increase in the climatologically (30 years) averaged cyanobacteria biomass and nitrogen fixation. Additional sensitivity experiments indicate that the biological-physical feedback mechanism through light absorption becomes more important with global warming. In general, we find a nonlinear response of cyanobacteria to changes in the atmospheric forcing fields as a result of life-cycle related feedback mechanisms. Overall, the sensitivity of the cyanobacteria-driven system suggests that biological-physical and life-cycle related feedback mechanisms are important and must therefore be included in future projection studies.     

Calcium isotope fractionation in coccoliths of cultured Calcidiscus leptoporus, Helicosphaera carteri, Syracosphaera pulchra and Umbilicosphaera foliosa

Four species of marine calcifying algae, the coccolithophores Calcidiscus leptoporus, Helicosphaera carteri, Syracosphaera pulchra and Umbilicosphaera foliosa were grown in laboratory cultures under temperatures varying between 14 and 23 °C, and one species, C. leptoporus, under varying [CO₃²⁻], ranging from 105 to 219 μmol/kg. Calcium isotope compositions of the coccoliths resemble in both absolute fractionation and temperature sensitivity previous calibrations of marine calcifying species e.g. Emiliania huxleyi (coccolithophores) and Orbulina universa (planktonic foraminifera) as well as inorganically precipitated CaCO₃, but also reveal small species specific differences. In contrast to inorganically precipitated calcite, but similar to E. huxleyi and O. universa, the carbonate ion concentration of the medium has no statistically significant influence on the Ca isotope fractionation of C. leptoporus coccoliths; however, combined data of E. huxleyi and C. leptoporus indicate that the observed trends might be related to changes of the calcite saturation state of the medium. Since coccoliths constitute a significant portion of the global oceanic CaCO₃ export production, the Ca isotope fractionation in these biogenic structures is important for defining the isotopic composition of the Ca sink of the ocean, one of the key parameters for modelling changes to the marine Ca budget over time. For the present ocean our results are in general agreement with the previously postulated and applied mean value of the oceanic Ca sink (Δ_sed) of about − 1.3‰, but the observed inter- and intra-species differences point to possible changes in Δ_sed through earth history, due to changing physico-chemical conditions of the ocean and shifts in floral and faunal assemblages.     

Identification of groundwater mean transit times of precipitation and riverbank infiltration by two‐component lumped parameter models

Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments, studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for 3 to 9 years and analysed for stable isotopes (δ¹⁸O and δ²H), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two‐component lumped parameter models (LPMs) that are solved using δ¹⁸O records. The study illustrates that two‐component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low‐permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration, short mTTs (<40 weeks) were determined for sites closer to the river (<200 m), whereas for the precipitation infiltration, the mTTs were longer (>80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface–groundwater interaction can be conceptualized by exploiting two‐component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.