Numerical solution to the oblique derivative boundary value problem on non-uniform grids above the Earth topography

The paper presents a numerical solution of the oblique derivative boundary value problem on and above the Earth’s topography using the finite volume method (FVM). It introduces a novel method for constructing non-uniform hexahedron 3D grids above the Earth’s surface. It is based on an evolution of a surface, which approximates the Earth’s topography, by mean curvature. To obtain optimal shapes of non-uniform 3D grid, the proposed evolution is accompanied by a tangential redistribution of grid nodes. Afterwards, the Laplace equation is discretized using FVM developed for such a non-uniform grid. The oblique derivative boundary condition is treated as a stationary advection equation, and we derive a new upwind type discretization suitable for non-uniform 3D grids. The discretization of the Laplace equation together with the discretization of the oblique derivative boundary condition leads to a linear system of equations. The solution of this system gives the disturbing potential in the whole computational domain including the Earth’s surface. Numerical experiments aim to show properties and demonstrate efficiency of the developed FVM approach. The first experiments study an experimental order of convergence of the method. Then, a reconstruction of the harmonic function on the Earth’s topography, which is generated from the EGM2008 or EIGEN-6C4 global geopotential model, is presented. The obtained FVM solutions show that refining of the computational grid leads to more precise results. The last experiment deals with local gravity field modelling in Slovakia using terrestrial gravity data. The GNSS-levelling test shows accuracy of the obtained local quasigeoid model.     

Multi-scale aquifer characterization and groundwater flow model parameterization using direct push technologies

Direct push (DP) technologies are typically used for cost-effective geotechnical characterization of unconsolidated soils and sediments. In more recent developments, DP technologies have been used for efficient hydraulic conductivity (K) characterization along vertical profiles with sampling resolutions of up to a few centimetres. Until date, however, only a limited number of studies document high-resolution in situ DP data for three-dimensional conceptual hydrogeological model development and groundwater flow model parameterization. This study demonstrates how DP technologies improve building of a conceptual hydrogeological model. We further evaluate the degree to which the DP-derived hydrogeological parameter K, measured across different spatial scales, improves performance of a regional groundwater flow model. The study area covers an area of ~60 km² with two overlying, mainly unconsolidated sand aquifers separated by a 5–7 m thick highly heterogeneous clay layer (in north-eastern Belgium). The hydrostratigraphy was obtained from an analysis of cored boreholes and about 265 cone penetration tests (CPTs). The hydrogeological parameter K was derived from a combined analysis of core and CPT data and also from hydraulic direct push tests. A total of 50 three-dimensional realizations of K were generated using a non-stationary multivariate geostatistical approach. To preserve the measured K values in the stochastic realizations, the groundwater model K realizations were conditioned on the borehole and direct push data. Optimization was performed to select the best performing model parameterization out of the 50 realizations. This model outperformed a previously developed reference model with homogeneous K fields for all hydrogeological layers. Comparison of particle tracking simulations, based either on the optimal heterogeneous or reference homogeneous groundwater model flow fields, demonstrate the impact DP-derived subsurface heterogeneity in K can have on groundwater flow and solute transport. We demonstrated that DP technologies, especially when calibrated with site-specific data, provide high-resolution 3D subsurface data for building more reliable conceptual models and increasing groundwater flow model performance.     

Anthropogenic influence on trace element geochemistry of healing mud (peloid) from Makirina Cove (Croatia)

Due to their balneotherapeutic features, the organic-rich sediments in Makirina Cove are an important source of healing mud. An environmental geochemistry approach using normalization techniques was applied to evaluate the anthropogenic contribution of trace metals to sediments used as healing mud. Sediment geochemistry was found to be associated with land-use change and storm events, as well as with proximity of a road with heavy traffic in the summer months. Local valley topography preferentially channels lithogenic and pollutant transport to the cove. Concentrations and distribution of trace metals indicate lithogenic (Ni, Cr, Co) and anthropogenic (Pb, Cu, Zn and Se) contributions to the sediments. The calculation of enrichment factors indicates a moderate (EFs between 2–3.5) input of anthropogenic Cu and Pb in surficial sediments to a depth of 10 cm. Patients using the Makirina Cove sediments as healing mud could be to some extent exposed to enhanced uptake of metals from anthropogenic sources via dermal contact.     

Effects of demand estimates on the evaluation and optimality of service centre locations

Public service systems, such as emergency health care, police or fire brigades, are critical for day-to-day functioning of the society. To design and operate these systems efficiently much data needs to be collected and properly utilised. Here, we use the OpenStreetMap (OSM) data to model the demand points (DPs), which approximate the geographical location of customers, and the road network, which is used to access or distribute services. We consider all inhabitants as customers, and therefore to estimate the demand, we use the available population grids. People are changing their location in the course of the day and thus the demand for services is changing accordingly. In this paper, we investigate how the used demand estimate affects the optimal design of a public service system. We calculate and compare efficient designs corresponding to two demand models, a night-time demand model when the majority of inhabitants rest at home and the demand model derived from the 24-hour average of the population density. We propose a simple measure to quantify the differences between population grids and we estimate how the size of differences affects the optimal structure of a public service system. Our analyses reveal that the efficiency of the service system is not only dependent on the placement strategy, but an inappropriate demand model has significant effects when designing a system as well as when evaluating its efficiency.     

Seismic demand in medium- and long-period structures

Based on the results of an extensive parametric study of elastic and inelastic response of SDOF systems, in which the most important structural parameters were varied and ground motions of very different characteristics were taken into account, simple formulae for determining the seismic demand in SDOF systems with natural periods in the medium- and long-period range are proposed. Seismic demand is expressed in terms of the mean values of maximum relative displacements and maximum input energy. These results can be used to provide rough estimates of structural behaviour when different damage models are applied. As well as this, the proposed formulae can be used to construct design spectra of the Newmark-Hall type.     

Towards an operational use of space imagery for oil pollution monitoring in the Mediterranean basin: a demonstration in the Adriatic Sea

Studies of operational pollution carried out by European commission - Joint Research Centre in the Mediterranean Sea for the years 1999-2004 are briefly introduced. The specific analysis of the Adriatic Sea for the same period demonstrates that this area has been characterized by a relevant number of illegal discharges from ships. After setting the historical background of the project AESOP (aerial and satellite surveillance of operational pollution in the Adriatic Sea), the content, partners and aim of the project are presented. Finally, the results of the first phase of the AESOP project are presented. The results seem very encouraging. For the first time in the Adriatic, real time detection of oil spills in satellite images and an immediate verification by the Coast Guard has been undertaken. An exploratory activity has also been carried out in collaboration with the University of Ljubljana to use automatic information system (AIS) to identify the ships detected in the satellite images.     

Five different types of lowered snow line in alpine valleys

The snow line for the same precipitation event in alpine valleys can vary greatly from valley to valley,even when the valleys are close together.Often this is n...     

Fabric transpositions in granite plutons — An insight from non-scaled analogue modelling

Investigations on a set of experimental models of highly viscous intrusions were carried out in order to study the internal strain pattern during vertical ascent and emplacement of granite intrusions. The strain pattern was determined by means of anisotropy of magnetic susceptibility (AMS) resulting from the orientation of magnetite particles in a liquid plaster medium. The modelled intrusions show distinct fabrics reflecting the flow of a rheologically complex, non-Newtonian material. During the vertical growth of the intrusion, constrictional vertical fabrics are transposed into flattening fabrics, and along with the development of low-intensity fabric domains are passively transported upwards. Vertical growth takes place along subvertical thrust shear zones that satisfactorily explain the discordant magmatic fabrics in granites along intrusion sides. The resulting complex fabric patterns suggest that the vertical movement of material in ascending intrusions is accommodated by various flow mechanisms operating simultaneously.     

Large wood load and transport in a flood-free period within an inter-dam reach: a decade of monitoring the Dyje River,Czech Republic

This study examined the large wood (LW) load and transport during the non-flood period (2009–2018) following major floods that occurred in 2002 and 2006 within the inter-dam reach of the Dyje River (Czech–Austrian border). The LW load was examined in 36 river corridor segments scattered within the reach in the 2009–2018 period. Two whole-reach surveys (2011 and 2019) on LW frequency and distribution were conducted, and the export of LW to the downstream reservoir was analysed between June 2013 and December 2018. In the period of non-flood discharges, the recruitment and depletion of LW were highly variable processes in space and time, leading to a considerable change in the total LW quantity. Whereas the total number of LW pieces decreased, the total LW volume increased because of the increasing dimensions of newly recruited pieces. The annual variability in the quantity of newly recruited pieces was better explained by the variation in the maximum annual discharges (y = 41.043ln(x) + 3.2737, R² = 0.5352) than by the variability in the number of days with wind gusts >17.2 m/s (y = 1.5004x + 82.096, R² = 0.118). The land use change with the abandonment of human settlements after World War II and the progressive expansion of forest was the major historical factor driving the increased recruitment of LW to the river corridor. While the 2006 (>100-year RI) flood brought approximately 1,250 LW pieces to the reservoir, the 2013 (1.5-year RI) flood delivered 45 pieces. The long-term average monthly input of LW to the reservoir was 7.7 pieces. The exceptional low-magnitude flood of 2013, which occurred at the beginning of the monitoring period, was shown to be a threshold above which the number of LW pieces that floated to the dam significantly increased. © 2020 John Wiley & Sons, Ltd.     

Structure and chemistry of (111) twin boundaries in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel crystals from Mogok

The atomic scale structure and chemistry of (111) twins in MgAl₂O₄ spinel crystals from the Pinpyit locality near Mogok (Myanmar, formerly Burma) were analysed using complementary methods of transmission electron microscopy (TEM). To obtain a three-dimensional information on the atomic structure, the twin boundaries were investigated in crystallographic projections and Using conventional electron diffraction and high-resolution TEM (HRTEM) analysis we have shown that (111) twins in spinel can be crystallographically described by 180° rotation of the oxygen sublattice normal to the twin composition plane. This operation generates a local hcp stacking in otherwise ccp lattice and maintains a regular sequence of kagome and mixed layers. In addition to rotation, no other translations are present in (111) twins in these spinel crystals. Chemical analysis of the twin boundary was performed by energy-dispersive X-ray spectroscopy (EDS) using a variable beam diameter (VBD) technique, which is perfectly suited for analysing chemical composition of twin boundaries on a sub-nm scale. The VBD/EDS measurements indicated that (111) twin boundary in spinel is Mg-deficient. Quantitative analyses of HRTEM (phase contrast) and HAADF-STEM (Z-contrast) images of (111) twin boundary have confirmed that Mg²⁺ ions are replaced with Be²⁺ ions in boundary tetrahedral sites. The Be-rich twin boundary structure is closely related to BeAl₂O₄ (chrysoberyl) and BeMg₃Al₈O₁₆ (taaffeite) group of intermediate polysomatic minerals. Based on these results, we conclude that the formation of (111) twins in spinel is a preparatory stage of polytype/polysome formation (taaffeite) and is a result of thermodynamically favourable formation of hcp stacking due to Be incorporation on the {111} planes of the spinel structure in the nucleation stage of crystal growth. The twin structure grows as long as the surrounding geochemical conditions allow its formation. The incorporation of Be induces a 2D-anisotropy and exaggerated growth of the crystal along the (111) twin boundary.