Relation between mean stress,thermodynamic, and lithostatic pressure

Pressure is one of the most important parameters to be quantified in geological problems. However, in metamorphic systems the pressure is usually calculated with two different approaches. One pressure calculation is based on petrological phase equilibria and this pressure is often termed thermodynamic pressure. The other calculation is based on continuum mechanics, which provides a mean stress that is commonly used to estimate the thermodynamic pressure. Both thermodynamic pressure calculations can be justified by the accuracy and applicability of the results. Here, we consider systems with low‐differential stress (<1 kbar) and no irreversible volumetric deformation, and refer to them as conventional systems. We investigate the relationship between mean stress and thermodynamic pressure. We discuss the meaning of thermodynamic pressure and its calculation for irreversible processes such as viscous deformation and heat conduction, which exhibit entropy production. Moreover, it is demonstrated that the mean stress for incompressible viscous deformation is essentially equal to the mean stress for the corresponding viscous deformation with elastic compressibility, if the characteristic time of deformation is five times longer than the Maxwell viscoelastic relaxation time that is equal to the ratio of shear viscosity to bulk modulus. For typical lithospheric rocks, this Maxwell time is smaller than c. 10,000 years. Therefore, numerical simulations of long‐term (>10 kyr) geodynamic processes, employing incompressible deformation, provide mean stress values that are close to the mean‐stress value associated with elastic compressibility. Finally, we show that for conventional systems the mean stress is essentially equal to the thermodynamic pressure. However, mean stress and, hence, thermodynamic pressure can be significantly different from the lithostatic pressure.     

Timing and mechanisms of Central Himalayan exhumation: discriminating between tectonic and erosion processes

The ability to deduce exhumation mechanisms from thermochronological data is hampered by the fact that assumptions on the thermal state of the lithosphere have to be made. Additional argumentation is generally required to discriminate between erosion-controlled and tectonically induced exhumation. This problem can be overcome by studying the spatial distribution of zircon and apatite (U-Th)/He and fission track data. In this work the variation of four different low temperature isotopic systems generating age trends along a sampling line is used to infer mechanisms of Quaternary exhumation in the Central High Himalayan Metamorphic Belt. Observed zircon age trends with southwards increasing cooling ages (from 0.5 to 1.7 Ma) are attributed to tectonically induced exhumation. The uniform apatite cooling ages clustered c. 0.5 Ma are attributed to erosion.     

Sediment related impacts of climate change and reservoir development in the Lower Mekong River Basin: a case study of the Nam Ou Basin,Lao PDR

This study applies the soil and water assessment tool (SWAT), with climate (precipitation and temperature) outputs from four general circulation models (GCMs) and a regional circulation model (PRECIS), to evaluate (1) the impacts of climate change on reservoir sedimentation and (2) the impacts of climate change and reservoir development on sediment outflow in the Nam Ou River Basin located in northern Laos. Three reservoir–density scenarios, namely one reservoir (1R), three reservoirs in series (3R), and five reservoirs in series (5R), were evaluated for both no climate change and climate change conditions. The results show that under no climate change conditions, by 2070, around 17, 14, and 15% of the existing reservoir storage volume in the basin will be lost for 1R, 3R, and 5R scenarios, respectively. Notably, under climate change scenario with highest changes in erosion and sediment outflux from the basin, the additional reduction in reservoir storage capacity due to sedimentation is estimated to be nearly 26% for 1R, 21% for 3R, and 23% for 5R. Climate change alone is projected to change annual sediment outflux from the basin by ?20 to 151%. In contrast, the development of reservoirs in the basin will reduce the annual sediment outflux from the basin varying from 44 to 80% for 1R, 44–81% for 3R, and 66–89% for 5R, considering climate change. In conclusion, climate change is expected to increase the sediment yield of the Nam Ou Basin, resulting in faster reduction of the reservoir’s storage capacity. Sediment yield from the Nam Ou River Basin is likely to decrease significantly due to the trapping of sediment by planned reservoirs. The impact of reservoirs is much more significant than the impact of climate change on the sediment outflow of the basin. Hence, it is necessary to investigate appropriate reservoir sediment management strategies.     

Application of Numerical Modeling for Groundwater Flow System Analysis in the Akaki Catchment,Central Ethiopia

A three dimensional steady-state finite difference groundwater flow model is used to quantify the groundwater fluxes and analyze the subsurface hydrodynamics in the Akaki catchment by giving particular emphasis to the well field that supplies water to the city of Addis Ababa. The area is characterized by Tertiary volcanics covered with thick residual and alluvial soils. The model is calibrated using head observations from 131 wells. The simulation is made in a two layer unconfined aquifer with spatially variable recharge and hydraulic conductivities under well-defined boundary conditions. The calibrated model is used to forecast groundwater flow pattern, the interaction of groundwater and surface water, and the effect of pumping on the well field under different scenarios. The result indicates that the groundwater flows regionally to the south converging to the major well field. Reservoirs and rivers play an important role in recharging the aquifer. Simulations made under different pumping rate indicate that an increase in pumping rate results in substantial regional groundwater level decline, which will lead to the drying of springs and shallow hand dug wells. Also, it has implications of reversal of flow from contaminated rivers into productive aquifers close to main river courses. The scenario analysis shows that the groundwater potential is not enough to sustain the ever-growing water demand of the city of Addis Ababa. The sensitivity and scenario analysis provided important information on the data gaps and the specific sites to be selected for monitoring, and may be of great help for transient model development. This study has laid the foundation for developing detailed predictive groundwater model, which can be readily used for groundwater management practices.     

A space-time discontinuous Galerkin method applied to single-phase flow in porous media

A space-time discontinuous Galerkin finite element method is proposed and applied to a convection-dominant single-phase flow problem in porous media. The numerical scheme is based on a coupled space-time finite element discretization allowing for discontinuous approximations in space and in time. The continuities on the element interfaces are weakly enforced by the flux treatments, so that no extra penalty factor has to be determined. The resulting space-time formulation possesses the advantage of capturing the steep concentration front with sharp gradients efficiently. The stability and reliability of the proposed approach is demonstrated by numerical experiments. The author is grateful to the DFG (German Science Foundation—Deutsche Forschungsgemeinschaft) for the financial support under the grant number Di 430/4-2.     

New insights into the crustal structure of the North German Basin from reprocessing of seismic reflection data using the Common Reflection Surface stack

The influence of deep crustal processes on basin formation and evolution and its relation to current morphology is not well understood yet. A key feature to unravel these issues is a detailed seismic image of the crust. A part of the data recorded by the hydrocarbon industry in the late 1970s and 1980s in the North German Basin were released to the public recently. The seismic reflection data were recorded down to 15 s two-way travel time. The mean Common Midpoint fold of about 20 is relatively low compared to contemporary seismic acquisitions. The processing of the 1980s focussed on the sedimentary structures to explore the hydrocarbon potential of this area. We applied the Common Reflection Surface stack technique to the data sets, which is well suited for low-fold data. The reprocessing was focussed on the imaging of the subsedimentary crustal range. The reprocessed images show enhanced reflections, especially in the mid and lower crustal part. Also, the image of the salt structures in the graben area was improved. Furthermore, the reprocessed images indicate an almost flat Moho topography in the area of the Glückstadt Graben and an additional lower crustal structure, which can be correlated with a high-density body found in recent gravity modeling studies.     

Corrosion-induced gas generation in a nuclear waste repository: Reactive geochemistry and multiphase flow effects

Corrosion of steel canisters, stored in a repository for spent fuel and high-level nuclear wastes, leads to the generation and accumulation of H₂ gas in the backfilled emplacement tunnels, which may significantly affect long-term repository safety. Previous studies have used H₂ generation rates based on the volume of the waste or canister material and the stoichiometry of the corrosion reaction. However, Fe corrosion and H₂ generation rates vary with time, depending on factors such as amount of Fe, water availability, water contact area and aqueous and solid chemistry. To account for these factors and feedback mechanisms, a chemistry model was developed related to Fe corrosion, coupled with two-phase (liquid and gas) flow phenomena that are driven by gas-pressure buildup associated with H₂ generation and water consumption. Results indicate that by dynamically calculating H₂ generation rates based on a simple model of corrosion chemistry, and by coupling this corrosion reaction with two-phase flow processes, the degree and extent of gas-pressure buildup could be much smaller compared to a model that neglects the coupling between flow and reactive transport mechanisms. By considering the feedback of corrosion chemistry, the gas pressure increases initially at the canister, but later decreases and eventually returns to a stabilized pressure that is slightly higher than the background pressure. The current study focuses on corrosion under anaerobic conditions for which the coupled hydrogeochemical model was used to examine the role of selected physical parameters on H₂ gas generation and corresponding pressure buildup in a nuclear waste repository. The developed model can be applied to evaluate the effect of water and mineral chemistry of the buffer and host rock on the corrosion reaction for future site-specific studies.