Solving the nonstationary Richards equation with adaptive hp-FEM

This paper examines the potential of the adaptive hp-FEM method for the numerical solution of time-dependent variably saturated Darcian flow problems described by the Richards equation. The method is illustrated on three model problems: a benchmark with known exact solution, groundwater seepage into a dry lysimeter box with time-dependent boundary conditions, and capillary barrier behavior under an intense infiltration. In the second part of the paper we present the weak formulation of the Richards equation for the Newton’s and Picard’s methods, give a brief overview of adaptive hp-FEM with emphasis on aspects that are usually not discussed in the literature, and we briefly introduce the open source adaptive hp-FEM library HERMES that was used to generate numerical results for this paper. All computations that we present are easily reproducible.     

Bi-normalized response spectra and seismic intensity in Bucharest for 1986 and 1990 Vrancea seismic events

The Vrancea subcrustal earthquakes of August 30,1986 and May 30,1990 are the two most recent seismic events that have occurred in Romania with moment magnitudes M W ≥ 6.9.The spectral analysis of the strong ground motions recorded in Bucharest reveals that despite small differences in magnitude between the 1986 and 1990 earthquakes,their frequency contents are very different,sometimes even opposing.The main focus of this study is to conduct a comparative analysis of the response spectra in terms of the bi-normalized response spectra(BNRS) proposed by Xu and Xie(2004 and 2007) for strong ground motions recorded in Bucharest during these two seismic events.The mean absolute acceleration and relative velocity response spectra for the two earthquakes are discussed and compared.Furthermore,the mean bi-normalized absolute acceleration and normalized relative velocity response spectra with respect to the control period T C are computed for the ground motions recorded in Bucharest in 1986 and 1990.The predominant period T P is also used in this study for the normalization of the spectral period axis.Subsequently,the methodology proposed by Martinez-Perreira and Bommer(1998) is applied in order to estimate the seismic intensity of the two events.The results are discussed and several conclusions regarding the possibility of using the bi-normalized response spectra(BNRS) are given.     

Thermal precursors in satellite images of the 1999 eruption of Shishaldin Volcano

Shishaldin Volcano, Unimak Island Alaska, began showing signs of thermal unrest in satellite images on 9 February 1999. A thermal anomaly and small steam plume were detected at the summit of the volcano in short-wave thermal infrared AVHRR (advanced very high resolution radiometer) satellite data. This was followed by over 2 months of changes in the observed thermal character of the volcano. Initially, the thermal anomaly was only visible when the satellite passed nearly directly over the volcano, suggesting a hot source deep in the central crater obscured from more oblique satellite passes. The "zenith angle" needed to see the anomaly increased with time, presumably as the thermal source rose within the conduit. Based on this change, an ascent rate of ca. 14 m per day for the thermal source was estimated, until it reached the summit on around 21 March. It is thought that Strombolian activity began around this time. The precursory activity culminated in a sub-Plinian eruption on 19 April, ejecting ash to over 45,000 ft. (13,700 m). The thermal energy output through the precursory period was calculated based on geometric constraints unique to Shishaldin. These calculations show fluctuations that can be tied to changes in the eruptive character inferred from seismic records and later geologic studies. The remote location of this volcano made satellite images a necessary observation tool for this eruption. To date, this is the longest thermal precursory activity preceding a sub-Plinian eruption recorded by satellite images in the region. This type of thermal monitoring of remote volcanoes is central in the efforts of the Alaska Volcano Observatory to provide timely warnings of volcanic eruption, and mitigate their associated hazards to air-traffic and local residents.     

Modelling of the sensing system of a vibrational magnetometer by digital computer

Summary An analytico-numerical computation of the transfer function of the sensing system of a vibrational magnetometer is presented. In discussing the results, the properties of the method are pointed out and a proposal for the sensing system is made with regard to the required sensitivity of the instrument.     

Time-series of δ26Mg values in a headwater catchment reveal decreasing magnesium isotope variability from precipitation to runoff

Data on temporal variability in Mg isotope ratios of atmospheric deposition and runoff are critical for decreasing the uncertainty associated with construction of isotope mass balances in headwater catchments, and statistical evaluation of isotope differences among Mg pools and fluxes. Such evaluations, in turn, are needed to distinguish between biotic and abiotic contributions to Mg²⁺ in catchment runoff. We report the first annual time-series of δ²⁶Mg values simultaneously determined for rainfall, canopy throughfall, soil water and runoff. The studied 55-ha catchment, situated in western Czech Republic, is underlain by Mg-rich amphibolite and covered by mature spruce stands. Between 1970 and 1996, the site received extremely high amounts of acid deposition and fly ash form nearby coal-burning power plants. The δ²⁶Mg values of open-area precipitation (median of −0.79‰) at our study site were statistically indistinguishable from the δ²⁶Mg values of throughfall (−0.73‰), but significantly different from the δ²⁶Mg values of soil water (−0.55‰) and runoff (−0.55‰). The range of δ²⁶Mg values during the observation period decreased in the order: open-area precipitation (0.57‰) > throughfall (0.27‰) > runoff (0.21‰) > soil water (0.16‰). The decreasing variability in δ²⁶Mg values of Mg²⁺ from precipitation to soil water and runoff reflected an increasing homogenization of atmospheric Mg in the catchment and its mixing with geogenic Mg. In addition to atmospheric Mg, runoff also contained Mg mobilized from the three major solid Mg pools, bedrock (δ²⁶Mg of −0.32‰), soil (−0.28‰), and vegetation (−0.31‰). The drought of summer 2019 did not affect the nearly constant δ²⁶Mg value of runoff. Collectively, our data show that within-catchment processes buffer the Mg isotope variability of the atmospheric input.     

Use of Groundwater Level Fluctuations near an Operating Water Supply Well to Estimate Aquifer Transmissivity

We developed a method to estimate aquifer transmissivity from the hydraulic-head data associated with the normal cyclic operation of a water supply well thus avoiding the need for interrupting the water supply associated with a traditional aquifer test. The method is based on an analytical solution that relates the aquifer's transmissivity to the standard deviation of the hydraulic-head fluctuations in one or more observation wells that are due to the periodic pumping of the production well. We analyzed the resulting analytical solution and demonstrated that when the observation wells are located near the pumping well, the solution has a simple, Dupuit like form. Numerical analysis demonstrates that the analytical solution can also be used for a quasi-periodic pumping of the supply well. Simulation of cyclic pumping in a statistically heterogeneous medium confirms that the method is suitable for analyzing the transmissivity of weakly or moderately heterogeneous aquifers. If only one observation well is available, and the shift in the phase of hydraulic-head oscillations between the pumping well and the observation well is not identifiable. Prior knowledge of aquifer's hydraulic diffusivity is required to obtain the value of the aquifer transmissivity.     

Global streamflow and thermal habitats of freshwater fishes under climate change

Climate change will affect future flow and thermal regimes of rivers. This will directly affect freshwater habitats and ecosystem health. In particular fish species, which are strongly adapted to a certain level of flow variability will be sensitive to future changes in flow regime. In addition, all freshwater fish species are exotherms, and increasing water temperatures will therefore directly affect fishes’ biochemical reaction rates and physiology. To assess climate change impacts on large-scale freshwater fish habitats we used a physically-based hydrological and water temperature modelling framework forced with an ensemble of climate model output. Future projections on global river flow and water temperature were used in combination with current spatial distributions of several fish species and their maximum thermal tolerances to explore impacts on fish habitats in different regions around the world. Results indicate that climate change will affect seasonal flow amplitudes, magnitude and timing of high and low flow events for large fractions of the global land surface area. Also, significant increases in both the frequency and magnitude of exceeding maximum temperature tolerances for selected fish species are found. Although the adaptive capacity of fish species to changing hydrologic regimes and rising water temperatures could be variable, our global results show that fish habitats are likely to change in the near future, and this is expected to affect species distributions.     

Marine transient electromagnetic sounding of deep buried hydrocarbon reservoirs: principles,methodologies and limitations

The possibility of a time‐domain electromagnetic sounding method using excitation and measurement of vertical electric fields to search for and identify deeply buried reservoirs of hydrocarbons offshore is investigated. The method operates on source–receiver offsets, which are several times less than the depth of the reservoir. Geoelectric information is obtained from the transient responses recorded in the pauses between the pulses of electric current in the absence of the source field. The basics of the method, as well as its sensitivity, resolution, and the highest accessible depth of soundings for various geological conditions in a wide range of sea depths, are analyzed. For the analysis, 1D and 3D geoelectric models of hydrocarbon reservoirs are used. It is shown that under existing technologies of excitation and measurement of vertical electric fields, the highest accessible depth of soundings can be up to 4 km. Technology for the inversion and interpretation of transient responses is demonstrated on experimental data.