Blocky regularization schemes for Full‐Waveform Inversion★

With ill‐posed inverse problems such as Full‐Waveform Inversion, regularization schemes are needed to constrain the solution. Whereas many regularization schemes end up smoothing the model, an undesirable effect with FWI where high‐resolution maps are sought, blocky regularization does not: it identifies and preserves strong velocity contrasts leading to step‐like functions. These models might be needed for imaging with wave‐equation based techniques such as Reverse Time Migration or for reservoir characterization. Enforcing blockiness in the model space amounts to enforcing a sparse representation of discontinuities in the model. Sparseness can be obtained using the ?¹ norm or Cauchy function which are related to long‐tailed probability density functions. Detecting these discontinuities with vertical and horizontal gradient operators helps constraining the model in both directions. Blocky regularization can also help recovering higher wavenumbers that the data used for inversion would allow, thus helping controlling the cost of FWI. While the Cauchy function yields blockier models, both ?¹ and Cauchy attenuate illumination and inversion artifacts.     

Diagenesis and the reconstruction of paleoenvironments: A method to restore original δO values of carbonate and phosphate from fossil tooth enamel

Intra-tooth δ¹⁸O variations within the carbonate (δ¹⁸Oc) and phosphate (δ¹⁸Op) components of tooth apatite were measured for Miocene and Pliocene hypsodont mammals from Afghanistan, Greece and Chad in order to evaluate the resistance of enamel to diagenetic alteration. Application of water-apatite interaction models suggest that the different kinetic behaviours of the phosphate-water and carbonate-water systems can be used to detect subtle oxygen isotope disequilibria in fossil enamel when intra-individual variations are considered. Selective alteration of the oxygen isotope composition from the carbonate component of Afghan and Greek enamels suggests inorganic isotopic exchange processes. Microbially-induced isotopic exchange for phosphate is demonstrated for the first time in enamel samples from Chad, in association with extensive recrystallization. In Chad, δ¹⁸Op values were derived from partial isotopic exchange with fossil groundwater during early diagenesis. Mass balance calculations using average carbonate content in enamel as a proxy for recrystallization, and the lowest δ¹⁸Op value of dentine as a proxy for the isotopic composition of the diagenetic fluid, indicate that diagenesis can alter δ¹⁸Op by as much as 3‰ in some enamel samples. This diagenetic alteration is also responsible for a decrease in intra-individual variations of up to 1‰ in affected specimens. The effects of diagenesis on δ¹⁸Op values of fossil enamel are not systematic, however, and can only be estimated if sequential δ¹⁸Op and δ¹⁸Oc analyses are performed on fossil enamel and dentine. Reconstruction of large temporal- or spatial-scale paleoclimates based on δ¹⁸Op analyses from mammalian teeth cannot be considered valid if enamel has been affected by bacterial activity or if the data cannot be corrected for diagenetic effects.     

An improved process-based representation of stream solute transport in the soil and water assessment tools

Hydrological models have long been used to study the interactions between land, surface and groundwater systems, and to predict and manage water quantity and quality. The soil and water assessment tool (SWAT), a widely used hydrological model, can simulate various ecohydrological processes on land and subsequently route the water quality constituents through surface and subsurface waters. So far, in-stream solute transport algorithms of the SWAT model have only been minimally revised, even though it has been acknowledged that an improvement of in-stream process representation can contribute to better model performance with respect to water quality. In this study, we aim to incorporate a new and improved solute transport model into the SWAT model framework. The new process-based model was developed using in-stream process equations from two well established models—the One-dimensional Transport with Inflow and Storage model and the Enhanced Stream Water Quality Model. The modified SWAT model (Mir-SWAT) was tested for water quality predictions in a study watershed in Germany. Compared to the standard SWAT model, Mir-SWAT improved dissolved oxygen (DO) predictions by removing extreme low values of DO (<6 mg/L) simulated by SWAT. Phosphate concentration peaks were reduced during high flows and a better match of daily predicted and measured values was attained using the Mir-SWAT model (R² = 0.17, NSE = −0.65, RSR = 1.29 with SWAT; R² = 0.28, NSE = −0.04, RSR = 1.02 with Mir-SWAT). In addition, Mir-SWAT performed better than the SWAT model in terms of Chlorophyll-a content particularly during winter months, improving the NSE and RSR for monthly average Chl-a by 74 and 42%, respectively. With the new model improvements, we aim to increase confidence in the stream solute transport component of the model, improve the understanding of nutrient dynamics in the stream, and to extend the applicability of SWAT for reach-scale analysis and management.     

Fluid flow through the sedimentary cover in northern Switzerland recorded by calcite–celestite veins (Oftringen borehole,Olten)

Abundant veins filled by calcite, celestite and pyrite were found in the core of a 719 m deep borehole drilled in Oftringen near Olten, located in the north-western Molasse basin, close to the thrust of the Folded Jura. Host rocks are calcareous marl, argillaceous limestone and limestone of the Dogger and Malm. The δ¹⁸O values of vein calcite are lower than in host rock carbonate and, together with microthermometric data from fluid inclusions in vein calcite, indicate precipitation from a seawater-dominated fluid at average temperatures of 56–68°C. Such temperatures were reached at the time of maximum burial of the sedimentary pile in the late Miocene. The depth profile of δ¹³C and ⁸⁷Sr/⁸⁶Sr values and Sr content of both whole-rock carbonate and vein calcite show marked trends towards negative δ¹³C, high ⁸⁷Sr/⁸⁶Sr, and low Sr content in the uppermost 50–150 m of the Jurassic profile (upper Oxfordian). The ⁸⁷Sr/⁸⁶Sr of vein minerals is generally higher than that of host rock carbonate, up to very high values corresponding to Burdigalian seawater (Upper Marine Molasse, Miocene), which represents the last marine incursion in the region. No evidence for internally derived radiogenic Sr (clay minerals) has been found and so an external source is required. S and O isotope composition of vein celestite and pyrite can be explained by bacterial reduction of Miocene seawater sulphate. The available data set suggests the vein mineralization precipitated from descending Burdigalian seawater and not from a fluid originating in the underlying Triassic evaporites.     

Lithofacies,depositional environments,regional biostratigraphy and age of the Chitarwata Formation in the Bugti Hills,Balochistan, Pakistan

The Oligocene-early Miocene Chitarwata Formation records a critical interval of terrestrial sedimentation that predates the Siwalik deposits on the Potwar Plateau of north-central Pakistan. This Oligocene-early Miocene time interval has long been considered as lacking in the entire Indo-Pakistan region. The Chitarwata Formation is widely exposed in the Sulaiman Range, but has never been described in detail in the Sulaiman Lobe, where the famous fossiliferous strata called ‘Bugti Bone Beds’ have been known for over a century and half. The Chitarwata Formation represents coastal-delta at the base, and plain and fluvial environments at the top. Lithofacies and sedimentary structures of the Chitarwata Formation in the Bugti area are described in detail, and show a clearly distinct lithologic pattern, different from that reported from the Zinda Pir area. The Chitarwata Formation also records an important transition in the evolution of the drainage systems in the area during the late Paleogene and early Neogene. This transition from the west-flowing paleo-Indus fluvial system to the development of the ancestral Indus drainage system may explain the numerous hiatuses that characterize the Chitarwata Formation. The abundance of fossil mammals from the Chitarwata and overlying Vihowa formation in the Bugti Hills provides critical biochronologic information that sheds new light on biostratigraphic correlation with the Zinda Pir area and for the entire Sulaiman Range.     

The impact of Dansgaard–Oeschger cycles on the loessic environment and malacofauna of Nussloch (Germany) during the Upper Weichselian

A loess sequence has been sampled continuously at high resolution in Nussloch (Rhine Valley, Germany) for malacological and grain-size analyses between ca. 34 and 20 ka. Molluscan abundance and richness, percentage in hygrophilous species and grain-size index show cyclical variations related to the lithological loess–gley alternation. Major molluscan abundance maxima were triggered by temperature increases through an enhancement of the reproduction cycle, whereas cyclical richness fluctuations and percentage in hygrophilous species reflect variations in local humidity and changes in the environmental mosaic. Malacological parameters allow the distinction of four environmental phases organised in cyclical successions correlated with most of the loess–gley doublets. The correlation of the grain-size index of the Nussloch loess sequence with the dust content of the GRIP ice core demonstrates the synchronicity of major molluscan abundance maxima and δ¹⁸O increases characterising temperature increases during Dansgaard–Oeschger interstades. A schematic model is proposed to link the North Atlantic Dansgaard–Oeschger climatic oscillations with local environmental changes indicated by both malacofauna and pedostratigraphy. This malacological study of the Nussloch loess sequence thus provides new information about the response of terrestrial loessic palaeoenvironments to millennial-timescale climatic fluctuations during the Upper Weichselian ( marine isotope stage 2 (MIS 2) and end of MIS 3).     

On the sensitivity of Antarctic sea ice model biases to atmospheric forcing uncertainties

Climate Dynamics - Although atmospheric reanalyses are an extremely valuable tool to study the climate of polar regions, they suffer from large uncertainties in these data-poor areas. In this work,...     

Geostatistical approach for the assessment of the water reservoir capacity in arid regions: a case study of the Akdarya reservoir,Uzbekistan

The paper presents the results of a new geostatistical approach to generate bathymetric surface models from point measurement converted into continuous contour surfaces of reservoir bottoms in Uzbekistan. Sedimentation of reservoirs raises engineering, environmental and economical issues for the communities around the world in areas affected by a strong water deficit. Because of Uzbekistan’s arid climatic conditions, and uneven spatial and temporal water resources distribution, responsive and innovative water availability assessment surveys of all major water reservoirs are required. Bathymetric surveying is a traditional method that is carried out for the estimation of reservoir volumes and surface areas of the corresponding reservoir stages in order to assess the water availability. Volume and surface area differences derived from multiple surveys of a reservoir provide storage loss estimates over time due to sedimentation. However, two main factors, such as intensive field data measurement and post data-processing, often limit the frequency of these surveys. Alternatively, innovative depth measurement technologies coupled with contouring and surface mapping programs provide automated reservoir volume and surface area calculations. This significantly reduces time, workload and financial burdens for reservoir sedimentation projects. This research proposes the use of geostatistical approach to assess the reservoir sedimentation in the Akdarya reservoir of Uzbekistan. The geostatistical approach includes (semi-) variogram analysis and interpolation (kriging and simulations—turning bands) techniques predicting values at unsampled locations for generating digital bathymetric surface models of reservoir bottom conditions in order to calculate the volume and surface area at a given water elevation. The simulation enables to have range of reservoir volumes and surface areas with the same probability, in comparison to the kriging and traditional methods. This gives a real estimation of the resource availability for water operators to manage natural resources and hydraulic infrastructure in a sustainable manner.