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.     

A first Holocene leaf wax isotope-based paleoclimate record from the semi-humid to semi-arid south-eastern Caucasian lowlands

The Holocene paleoclimate of the Caucasus region is rather complex and not yet well understood: while existing studies are mainly based on pollen records from high-altitude and humid lowland regions, no records are available from the semi-humid to semi-arid south-eastern Caucasian lowlands. Therefore, this study investigated compound-specific δ²H and δ¹³C isotopes of leaf wax biomarkers from Holocene floodplain soils in eastern Georgia. Our results show that the leaf wax δ²H signal from the paleosols mostly reflects changes in the moisture source and its isotopic composition. Depleted δ²H values before ~8 cal ka bp change towards enriched values after ~5 cal ka bp and become again depleted after ~1.6 cal ka bp. This trend could be caused by Holocene changes of the isotopic compositions of the Black and eastern Mediterranean Sea, and/or by varying contribution of both moisture sources linked with the North Atlantic Oscillation. The leaf wax δ¹³C signal from the paleosols directly indicates varying local water availability and drought stress. Depleted δ¹³C values before ~8 and after ~5 cal ka bp indicate wetter local conditions with higher water availability, whereas more enriched values during the middle Holocene (~8 until at least 5 cal ka bp ) indicate drier conditions with increased drought stress.     

Conditional heavy-tail behavior with applications to precipitation and river flow extremes

This article deals with the right-tail behavior of a response distribution \(F_Y\) conditional on a regressor vector \({\mathbf {X}}={\mathbf {x}}\) restricted to the heavy-tailed case of Pareto-type conditional distributions \(F_Y(y|\ {\mathbf {x}})=P(Y\le y|\ {\mathbf {X}}={\mathbf {x}})\), with heaviness of the right tail characterized by the conditional extreme value index \(\gamma ({\mathbf {x}})>0\). We particularly focus on testing the hypothesis \({\mathscr {H}}_{0,tail}:\ \gamma ({\mathbf {x}})=\gamma _0\) of constant tail behavior for some \(\gamma _0>0\) and all possible \({\mathbf {x}}\). When considering \({\mathbf {x}}\) as a time index, the term trend analysis is commonly used. In the recent past several such trend analyses in extreme value data have been published, mostly focusing on time-varying modeling of location or scale parameters of the response distribution. In many such environmental studies a simple test against trend based on Kendall’s tau statistic is applied. This test is powerful when the center of the conditional distribution \(F_Y(y|{\mathbf {x}})\) changes monotonically in \({\mathbf {x}}\), for instance, in a simple location model \(\mu ({\mathbf {x}})=\mu _0+x\cdot \mu _1\), \({\mathbf {x}}=(1,x)'\), but the test is rather insensitive against monotonic tail behavior, say, \(\gamma ({\mathbf {x}})=\eta _0+x\cdot \eta _1\). This has to be considered, since for many environmental applications the main interest is on the tail rather than the center of a distribution. Our work is motivated by this problem and it is our goal to demonstrate the opportunities and the limits of detecting and estimating non-constant conditional heavy-tail behavior with regard to applications from hydrology. We present and compare four different procedures by simulations and illustrate our findings on real data from hydrology: weekly maxima of hourly precipitation from France and monthly maximal river flows from Germany.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Development and validation of a finite element morphological model for shallow water basins

Recognising the importance of understanding sediment dynamics to evaluate the status of a coastal lagoon environment, this work has been focused on the investigation of the hydrodynamic and sediment transport processes occurring in such basins. In order to describe the lagoon system, a modelling approach combining hydrodynamics, waves and sediment dynamics has been developed. The framework of the numerical model consists of a finite element hydrodynamic model, a third generation finite element spectral wave model and a sediment transport and morphodynamic model for both cohesive and non-cohesive sediments. The model adopts the finite element technique for spatial integration, which has the advantage to describe more accurately complicated bathymetry and irregular boundaries for shallow water areas. The developed model has been applied to test cases and to a very shallow tidal lagoon, the Venice Lagoon, Italy. Numerical results show good agreement with water level, waves and turbidity measurements collected in several monitoring stations inside the Lagoon of Venice. Such a model represents an indispensable tool in analysing coastal problems and assessing morphological impacts of human interference.     

Zircon evidence for a ~200 k.y. supereruption-related thermal flare-up in the Miocene southern Black Mountains,western Arizona,USA

The Silver Creek caldera (southern Black Mountains, western Arizona) is the source of the 18.8 Ma, >700 km³ Peach Spring Tuff (PST) supereruption, the largest eruption generated in the Colorado River Extensional Corridor (CREC) of the southwestern United States. Within and immediately surrounding the caldera is a sequence of volcanics and intrusions ranging in age from ~19 to 17 Ma. These units offer a record of magmatic processes prior to, during, and immediately following the PST eruption. To investigate the thermal evolution of the magmatic center that produced the PST, we applied a combination of Ti-in-zircon thermometry, zircon saturation thermometry, and high-precision U–Pb CA–TIMS zircon dating to representative pre- and post-supereruption volcanic and intrusive units from the caldera and its environs. Similar to intracaldera PST zircons, zircons from a pre-PST trachytic lava (19 Ma) and a post-PST caldera intrusion (18.8 Ma) yield exceptionally high-Ti concentrations (most >20 ppm, some up to nearly 60 ppm), corresponding to calculated temperatures that exceed 900 °C. In these units, Ti-in-zircon temperatures typically surpass zircon saturation temperatures (ZSTs), suggesting the entrainment of zircon that had grown in hotter environments within the magmatic system. Titanium concentrations in younger volcanic and intrusive units (~18.7–17.5 Ma) decline through time, corresponding to an average cooling rate of 10^?3.5 °C/year. The ~200 k.y. thermal peak evident at Silver Creek caldera is spatially limited: elsewhere in the Miocene record of the northern CREC, Ti-in-zircon concentrations and ZSTs are much lower, suggesting that felsic magmas were generally substantially cooler.     

Vertical column abundances of HCN deduced from ground-based infrared solar spectra: Long-term trend and variability

A set of high-resolution IR solar spectra recorded at the International Scientific Station of the Jungfraujoch, Switzerland, from 84/06 to 93/06, and at the National Solar Observatory McMath-Pierce solar telescope facility on Kitt Peak, Arizona, U.S.A. from 78/05 to 92/07 have been analyzed to determine the vertical column abundances of hydrogen cyanide, HCN, above the two stations. The analysis was based on least-squares curve fitting of calculated spectra to the observations encompassing the P4 and the P8 lines of HCN respectively located at 3299.5273 and 3287.2483 cm^–1. The results obtained for the two stations indicate that no significant long-term trend affects either of the two databases; however, this analysis reveals variable increases during springtime of up to a factor of 2 in the HCN total column above the Jungfraujoch and even up to 3 above Kitt Peak. The calculated mean vertical column abundances, excluding the spring observations, are equal to (2.55±0.30)×10¹⁵ molec./cm² (S.D.) and (2.75±0.30)×10¹⁵ molec./cm² respectively above the Jungfraujoch and the Kitt Peak observatories. Based on a realistic volume mixing ratio profile, these columns translate into mean volume mixing ratios equal to 190×10^–12 ppv at the respective altitudes of the stations.     

Creep in rock as a stochastic process

The heterogeneous character of rock suggests the use of statistical mechanics for the formulation of models of plasticity and creep. A stochastic creep model which allows for large variations in strength and size of the involved rock mass is outlined. The model, which is analogous to a previously developed concept for clastic materials, is shown to account rather well for typical creep behaviour of rock at moderate stresses and temperatures. Thus, it yields the often observed linear relation between the creep strain and the logarithm of time.     

The energy balance experiment EBEX-2000. Part II: Intercomparison of eddy-covariance sensors and post-field data processing methods

The eddy-covariance method is the primary way of measuring turbulent fluxes directly. Many investigators have found that these flux measurements often do not satisfy a fundamental criterion—closure of the surface energy balance. This study investigates to what extent the eddy-covariance measurement technology can be made responsible for this deficiency, in particular the effects of instrumentation or of the post-field data processing. Therefore, current eddy-covariance sensors and several post-field data processing methods were compared. The differences in methodology resulted in deviations of 10% for the sensible heat flux and of 15% for the latent heat flux for an averaging time of 30 min. These disparities were mostly due to different sensor separation corrections and a linear detrending of the data. The impact of different instrumentation on the resulting heat flux estimates was significantly higher. Large deviations from the reference system of up to 50% were found for some sensor combinations. However, very good measurement quality was found for a CSAT3 sonic together with a KH20 krypton hygrometer and also for a UW sonic together with a KH20. If these systems are well calibrated and maintained, an accuracy of better than 5% can be achieved for 30-min values of sensible and latent heat flux measurements. The results from the sonic anemometers Gill Solent-HS, ATI-K, Metek USA-1, and R.M. Young 81000 showed more or less larger deviations from the reference system. The LI-COR LI-7500 open-path H₂O/CO₂ gas analyser in the test was one of the first serial numbers of this sensor type and had technical problems regarding direct solar radiation sensitivity and signal delay. These problems are known by the manufacturer and improvements of the sensor have since been made. The National Center for Atmospheric Research is supported by the National Science Foundation.