Analytical hydrologic models and the design of policy instruments for groundwater-quality management

This paper explores how analytical hydrologic models can inform the effective design and choice of policy instruments to manage groundwater quality by coupling a social-planner’s problem of optimal groundwater-quality management with analytical solutions from the hydrology literature. A theoretical analysis is performed in order to characterize the properties of an optimal emissions policy. The model is then applied in a numerical analysis of groundwater contamination by chloride from highway deicers, demonstrating the relevance of the theoretical results to practical management settings. This analytical approach can help determine which policy instruments are likely to be effective in controlling groundwater pollution, especially if costly numerical groundwater models are not available for the aquifer in question. Unlike previous economic studies of groundwater contamination, this approach defines optimal emissions policies as a function of geophysical parameters employed by hydrologists such as distance between source and sink, groundwater velocity, and aquifer dispersivity. The theoretical section of this paper also demonstrates that the relationship between geophysical parameters and optimal emissions levels may be ambiguous.     

Mesoscale Reproducibility in Regional Ocean Modelling with a Three-Dimensional Stratification Estimate Based on Aviso-Argo Data

ABSTRACT

For dynamically consistent, high-resolution, yet cost-effective regional oceanic downscaling modelling, an empirical three-dimensional (3D) density estimate based on publicly available datasets is utilized for the Regional Oceanic Modeling System (ROMS) with simple data assimilation (i.e., TS nudging, where TS stands for temperature and salinity). We rely on a method built upon the two-layer model to reconstruct a mesoscale 3D temperature and salinity field, referred to as Tokyo University of Marine Science and Technology (TUMSAT)-TS, using near real-time altimeter-derived dynamic height along with Argo float profiling data. The TUMSAT-TS is first validated using in situ hydrographic data, then is implemented in the Japan Coastal Ocean Predictability Experiment (JCOPE2)-ROMS downscaling system for the Kuroshio region off Japan. We explore the usability of TUMSAT-TS by carrying out three comparative simulations with temperature and salinity nudging towards the (i) TUMSAT-TS and (ii) JCOPE2-TS fields, and (iii) without the nudging. Whereas the unassimilated case fails to properly account for the Kuroshio, both datasets individually are found to help reproduce the mesoscale variability of the Kuroshio, as well as its transient paths, volume transport, associated kinetic energy (KE) and eddy KE, and seasonally varying stratification.     

Morphodynamic effects of vegetation life stage on experimental meandering channels

Previous studies have demonstrated that riparian vegetation leads to channel transformation from a multi-bar to a single-thread channel planform. However, it still remains unclear how the presence of pioneer and mature vegetation affects the morphodynamics of single-thread meandering rivers. In this study, we therefore investigated the effects of vegetation strength on the morphodynamic evolution of an experimental meandering channel. Three physical laboratory experiments were conducted using alfalfa sprouts in different life stages – no vegetation, immature vegetation, and mature vegetation – to simulate different floodplain vegetation strengths. Our results demonstrate that vegetation plays a key role in mediating bank erosion and point-bar accretion, and that this is reflected in both the evolution of the channel bed as well as the sediment flux. The presence of mature vegetation maintained a deep, single-thread channel by reducing bank erosion, thereby limiting both channel widening and sediment storage capacity. Conversely, an unvegetated floodplain led to channel widening and high sediment storage capacity. Channel evolution in the unvegetated scenario showed that the active sediment supply from outer bank erosion led to slightly delayed point-bar accretion on the inner banks due to helical flow, deflecting the surface flow toward the outer banks and causing further erosion. In contrast, in the immature vegetation scenario, the outer banks were also initially eroded, but point-bar accretion did not clearly progress. This led to a greater width-to-depth ratio, resulting in a transition from a single- to a multi-thread channel with minor flow paths on the floodplain. The experimental results suggest that the eco-morphodynamic effects of young (low-strength) and mature (high-strength) vegetation are different. Notably, low-strength, early-stage vegetation increases channel complexity by accelerating both channel widening and branching, and therefore might promote the coexistence of multi-bars and pioneer vegetation.     

Irreversible thermodynamic and viscoelastic model for power-law relaxation and attenuation of rocks

A constitutive law for viscoelastic behaviour of rocks is derived from irreversible thermodynamics. To this model, two specific parameters are introduced; one is an internal state variable which is a variable concerning the microstructures such as defects in crystals or microcracks, and the other is a temperature reduced time obtained by normalizing the various temperature behaviours. A large number of internal state variables have the respective relaxation times and show the respective time evolutions, while a set of the time evolutions generates temporal power-law behaviour of rocks. The time evolutions of internal states are regarded as dynamics of elements of the generalized Maxwell model, and the stress–strain relation is represented by a response function following a temporal power-law in terms of linear system theory. This relation is inversely formulated to investigate the source field from output data. This model enables us to explain experimentally-based constitutive laws for transient and steady-state behaviour of rocks (e.g., lherzolite) following a temporal power-law and for attenuation behaviour of polycrystals (e.g., olivine) represented by a relation between the quality factor and frequency. Both laws show power-laws on deformation time or frequency depending on the fractal structure in polycrystals or rocks, and the experimental high-temperature behaviours can be extrapolated to long deformation time or high frequency behaviour.     

Coupled uranium isotope and sea-level variations in the oceans

Globally, rivers supply uranium to the oceans with excess ²³⁴U relative to secular equilibrium and ²³⁴U taken-up by corals can be used for dating. In addition, the ²³⁴U abundance in sea water, at the time the coral was growing, can be measured independently. The veracity of U-series ages used in determining past sea-level variations is dependent on selecting pristine corals free from diagenetic alteration. A quantitative test for alteration assumes invariant ²³⁴U abundances in the oceans for at least the past half a million years and results from samples outside of a narrow range in modern ocean ²³⁴U abundance are excluded from data sets. Here, we have used previously published data to show that ²³⁴U in the oceans appears to be variable and directly related to changes in sea-level, not only over long glacial-interglacial timescales but also at very short, centennial timescales. Most of the previously discarded data can be used to provide valuable additional sea-level information. The process permits a unique insight into the interplay between sources and sinks of uranium in the oceans mediated by sea-level changes at rates far faster than previously thought possible. Similar, rapid sea-level, forcing of other trace element abundances in the oceans is likely.     

Holocene Refreshening and Reoxygenation of a Bothnian Sea Estuary Led to Enhanced Phosphorus Burial

Salinity variations in restricted basins like the Baltic Sea can alter their vulnerability to hypoxia (i.e., bottom water oxygen concentrations <2 mg/l) and can affect the burial of phosphorus (P), a key nutrient for marine organisms. We combine porewater and solid-phase geochemistry, micro-analysis of sieved sediments (including XRD and synchrotron-based XAS), and foraminiferal δ¹⁸O and δ¹³C analyses to reconstruct the bottom water salinity, redox conditions, and P burial in the Ångermanälven estuary, Bothnian Sea. Our sediment records were retrieved during the Integrated Ocean Drilling Program (IODP) Baltic Sea Paleoenvironment Expedition 347 in 2013. We demonstrate that bottom waters in the Ångermanälven estuary became anoxic upon the intrusion of seawater in the early Holocene, like in the central Bothnian Sea. The subsequent refreshening and reoxygenation, which was caused by gradual isostatic uplift, promoted P burial in the sediment in the form of Mn-rich vivianite. Vivianite authigenesis in the surface sediments of the more isolated part of the estuary ultimately ceased, likely due to continued refreshening and an associated decline in productivity and P supply to the sediment. The observed shifts in environmental conditions also created conditions for post-depositional formation of authigenic vivianite, and possibly apatite formation, at ～8 m composite depth. These salinity-related changes in redox conditions and P burial are highly relevant in light of current climate change. The results specifically highlight that increased freshwater input linked to global warming may enhance coastal P retention, thereby contributing to oligotrophication in both coastal and adjacent open waters.     

Crustal anorthosite formation by deep-seated hydrothermal circulation beneath fast-spreading axis: Constraints from chronological approach,Sr isotope,and fluid–chromite inclusion investigation

Hydrothermal circulation beneath the spreading axis plays a significant role in the exchange of energy and mass between the solid Earth and the oceans. Deep-seated hydrothermal circulation down to the crust/mantle boundary in the fast-spreading axis has been introduced by a number of studies regarding geological investigations and numerical models. In order to assess a reaction between hydrothermal fluid and host rock around the crust/mantle boundary, we conducted bulk trace element and Sr isotope analyses with a series of in situ investigations for crustal anorthosite, a reaction product between hydrothermal fluid and gabbro in the lowermost crustal section along Wadi Fizh, northern Oman ophiolite. In addition, we conducted titanite U–Pb isotope analyses to evaluate timing of the crustal anorthosite formation in the framework of the evolutional process of the Oman ophiolite. We estimated the formation age of the crustal anorthosite at 97.5 Ma ± 5.0 Ma, overlapping with the timing of the crust formation in the paleo spreading axis. The crustal anorthosite shows high-Th/U ratio (~2.5) and high-initial ⁸⁷Sr/⁸⁶Sr ratio (0.7050) due to seawater-derived hydrothermal fluid ingress into the precursor gabbro. With using analytical technique of micro-excavation at cryo-temperature, we detected Cl from a few micrometer-sized inclusion of aqueous fluid and chromite grains. The solubility of Cr was enhanced by complexation reactions with Cl in the hydrothermal fluid. Regarding reconstructed three-dimensional mass distribution of the inclusion and chromite composition, maximum Cr content of parental fluid was estimated at ~69 000 μg/g. The exceptionally high-Cr content was achieved locally by leaking of fluid and synchronous chromite crystallization during fluid entrapment. Presence of the deep-seated hydrothermal circulation could be assigned to the segment end, where cold seawater penetrates into the lowermost crust and extract heat along widely spaced network-like fluid channel.     

Comparative Performance of Steel Drainage Pipes Against Flood-Induced Deformation in River Levee

Geotechnical and Geological Engineering - Steel drainage pipes that can provide both drainage and reinforcement functions are expected to give better performance in levee protection against...     

Historical Auroras in the 990s: Evidence of Great Magnetic Storms

A significant carbon-14 enhancement has recently been found in tree rings for the year 994, suggesting an extremely strong and brief cosmic ray flux event. The origin of this particular cosmic ray event has not been confirmed, but one possibility is that it might be of solar origin. Contemporary historical records of low-latitude auroras can be used as supporting evidence of intense solar activity around that time. We investigate previously reported as well as new records that have been found in contemporary observations from the 990s to determine potential auroras. Records of potential red auroras in late 992 and early 993 were found around the world, i.e. in the Korean Peninsula, Saxonian cities in modern Germany, and the Island of Ireland, suggesting the occurrence of an intense geomagnetic storm driven by solar activity.     

Influence of the Atmospheric Model on Hanle Diagnostics

We clarify the uncertainty in the inferred magnetic field vector via the Hanle diagnostics of the hydrogen Lyman-\(\upalpha\) line when the stratification of the underlying atmosphere is unknown. We calculate the anisotropy of the radiation field with plane-parallel semi-empirical models under the nonlocal thermal equilibrium condition and derive linear polarization signals for all possible parameters of magnetic field vectors based on an analytical solution of the atomic polarization and Hanle effect. We find that the semi-empirical models of the inter-network region (FAL-A) and network region (FAL-F) show similar degrees of anisotropy in the radiation field, and this similarity results in an acceptable inversion error (e.g., \({\sim}\, 40~\mbox{G}\) instead of 50 G in field strength and \({\sim}\,100^{\circ}\) instead of \(90^{\circ}\) in inclination) when FAL-A and FAL-F are swapped. However, the semi-empirical models of FAL-C (averaged quiet-Sun model including both inter-network and network regions) and FAL-P (plage regions) yield an atomic polarization that deviates from all other models, which makes it difficult to precisely determine the magnetic field vector if the correct atmospheric model is not known (e.g., the inversion error is much larger than 40% of the field strength; \({>}\,70~\mbox{G}\) instead of \(50~\mbox{G}\)). These results clearly demonstrate that the choice of model atmosphere is important for Hanle diagnostics. As is well known, one way to constrain the average atmospheric stratification is to measure the center-to-limb variation of the linear polarization signals. The dependence of the center-to-limb variations on the atmospheric model is also presented in this paper.