Simulating rainfall time-series: how to account for statistical variability at multiple scales?

Daily rainfall is a complex signal exhibiting alternation of dry and wet states, seasonal fluctuations and an irregular behavior at multiple scales that cannot be preserved by stationary stochastic simulation models. In this paper, we try to investigate some of the strategies devoted to preserve these features by comparing two recent algorithms for stochastic rainfall simulation: the first one is the modified Markov model, belonging to the family of Markov-chain based techniques, which introduces non-stationarity in the chain parameters to preserve the long-term behavior of rainfall. The second technique is direct sampling, based on multiple-point statistics, which aims at simulating a complex statistical structure by reproducing the same data patterns found in a training data set. The two techniques are compared by first simulating a synthetic daily rainfall time-series showing a highly irregular alternation of two regimes and then a real rainfall data set. This comparison allows analyzing the efficiency of different elements characterizing the two techniques, such as the application of a variable time dependence, the adaptive kernel smoothing or the use of low-frequency rainfall covariates. The results suggest, under different data availability scenarios, which of these elements are more appropriate to represent the rainfall amount probability distribution at different scales, the annual seasonality, the dry-wet temporal pattern, and the persistence of the rainfall events.     

Record of enhanced late Pliocene tectono‐volcanic activity in the Lesser Antilles forearc basin: Insights into the long‐term behaviour of the Lesser Antilles subduction zone

In this study, we investigate the Lesser Antilles forearc basin, focusing on the late Pliocene to Pleistocene sedimentary archives in order to track the occurrence of extreme events triggered by enhanced subduction‐related tectono‐volcanic activity. We identify late Piacenzian deposits covering a major regional erosional surface, displaying sedimentary dykes and large marine boulders embedded in a mixed continental–marine matrix, characteristic of tsunamites. We interpret this episode of platform emersion and the successive cataclysmic deposits as resulting from enhanced tectonic activity at the interface of the subduction zone, synchronous with the initiation of the Lesser Antilles volcanic arc. We then discuss the implications in terms of the mechanical behaviour of the Lesser Antilles subduction zone.     

Simulation of site-scale water fluxes in desert and natural oasis ecosystems of the arid region in Northwest China

The study of water fluxes is important to better understand hydrological cycles in arid regions. Data-driven machine learning models have been recently applied to water flux simulation. Previous studies have built site-scale simulation models of water fluxes for individual sites separately, requiring a large amount of data from each site and significant computation time. For arid areas, there is no consensus as to the optimal model and variable selection method to simulate water fluxes. Using data from seven flux observation sites in the arid region of Northwest China, this study compared the performance of random forest (RF), support vector machine (SVM), back propagation neural network (BPNN), and multiple linear regression (MLR) models in simulating water fluxes. Additionally, the study investigated inter-annual and seasonal variation in water fluxes and the dominant drivers of this variation at different sites. A universal simulation model for water flux was built using the RF approach and key variables as determined by MLR, incorporating data from all sites. Model performance of the SVM algorithm (R² = 0.25–0.90) was slightly worse than that of the RF algorithm (R² = 0.41–0.91); the BPNN algorithm performed poorly in most cases (R² = 0.15–0.88). Similarly, the MLR results were limited and unreliable (R² = 0.00–0.66). Using the universal RF model, annual water fluxes were found to be much higher than the precipitation received at each site, and natural oases showed higher fluxes than desert ecosystems. Water fluxes were highest during the growing season (May–September) and lowest during the non-growing season (October–April). Furthermore, the dominant drivers of water flux variation were various among different sites, but the normalized difference vegetation index (NDVI), soil moisture and soil temperature were important at most sites. This study provides useful insights for simulating water fluxes in desert and oasis ecosystems, understanding patterns of variation and the underlying mechanisms. Besides, these results can make a contribution as the decision-making basis to the water management in desert and oasis ecosystems.     

Predicted present-day evolution patterns of ice thickness and bedrock elevation over Greenland and Antarctica

This paper discusses predicted evolution patterns of present-day changes of ice thickness, surface elevation, and bedrock elevation over the Greenland and Antarctic continents. These were obtained from calculations with dynamic 3-D ice sheet models which were coupled to a visco-elastic solid Earth model. The experiments were initialized over the last two glacial cycles and subsequently averaged over the last 200 years to obtain the current evolution. The calculations indicate that the Antarctic Ice Sheet is still adjusting to the last glacial-interglacial transition yielding a decreasing ice volume and a rising bedrock elevation of the order of several centimetres per year. The Greenland Ice Sheet was found to be close to a stationary state with a mean thickness change of only a few millimetres per year, but the calculations revealed large spatial differences. Predicted patterns over Greenland are characterized by a small thickening over the ice sheet interior and a general thinning of the ablation area. In Antarctica, almost all of the predicted changes are concentrated in the West Antarctic Ice Sheet, which is still retreating at both the Weddell and Ross Sea margins. Over most of both ice sheets, the model indicates that the surface elevation trend is dominated by ice thickness changes rather than by bedrock elevation changes.     

Climate Warming Mitigation from Nationally Determined Contributions

Individual countries are requested to submit nationally determined contributions(NDCs) to alleviate global warming in the Paris Agreement. However, the global climate effects and regional contributions are not explicitly considered in the countries’ decision-making process. In this study, we evaluate the global temperature slowdown of the NDC scenario(?T =0.6°C) and attribute the global temperature slowdown to certain regions of the world with a compact earth system model.Considering reductions ...     

Controlling factors on the morphology and spatial distribution of methane-related tubular concretions – Case study of an Early Eocene seep system

Up to 10 m in length and >1 m in diameter tubular, calcite-cemented sandstone concretions are hosted by the faulted Dikilitash unconsolidated sands and sandstones. These structures document shallow subsurface pathways of Early Eocene methane seepage in the Balkan Mountains foreland (NE Bulgaria). Their exceptional exposure allowed a unique study of the factors governing the morphology and spatial distribution of such fossilized fluid conduits. The large dimensions and subvertical, cylindrical shape of the most common tube type primarily reflects the buoyancy-driven, vertical path of an ascending gas-bearing fluid through permeable, mainly unconsolidated sandy host sediments. Tube morphology was also influenced by local stratigraphic anisotropies and might as well document differences in former seepage conditions. Mapping of >800 tubular concretions showed the NNW–NNE elongation and alignment of tube clusters and massive cemented sandstone structures. This suggests that Paleogene fault systems played a major role in directing the movement of fluids. However, within a single tube cluster, tubes are preferentially aligned, over distances up to 50 m along directions at an angle between 10° and 36° with respect to the inferred NNW–NNE, cluster parallel fault traces. In addition, cylindrical tubes of analogue dimensions are aligned over distances >100 m along N15° to N25°-oriented directions. It is hypothesized that this spatial geometry of tubular concretions reflects the complex geometry of deformations structures in fault damage zones along which fluids were preferentially channelled.     

Modeling the impact of microbial activity on redox dynamics in porous media

The present study investigates the interaction between microbial growth and activity and the redox dynamics in natural porous media. The impact the transport regime has on this interaction is also addressed. Expressions for microbial growth are incorporated into a geochemical reaction network linking redox reaction rates to the activity of the microorganisms. A flexible simulation environment, the Biogeochemical Reaction Network Simulator (BRNS) is used for this purpose. Two reactive transport applications relevant to fields of contaminant hydrology and early diagenesis are simulated with the BRNS. Model results are evaluated based on a comparison with comprehensive datasets on the biodegradation of lactate in a sand column experiment and on the distribution of redox-sensitive chemical species in marine sediments of the Skagerrak, Denmark. It is shown that, despite quite different transport regimes, the geomicrobiological model performs equally well in the reproduction of measured chemical species distribution for both applications. This result emphasizes the broad applicability of the proposed approach. Our simulations support that the competitive behavior between various microbial groups is a process controlling the development of redox stratified environments. Furthermore, it is also shown that the transport regime is a key controlling factor for the degree of spatial correlation between microbial biomass distributions and redox reaction rates. Although all our simulations yield a pronounced stratification of the redox processes in the system, the biomass distribution is related to the associated reaction rates only in case of the advection controlled column experiment. In the early diagenetic application, mixing due to bioturbation is the dominant transport process for particulate matter, hence leading to fairly homogeneous distribution of bacterial biomasses which are unrelated to the spatial distribution of redox reaction rates. This homogeneous biomass distribution combined with the 1G carbon degradation model approach might explain why the steady state concentration profiles in such systems can be reproduced by diagenetic models without explicit representation of microbial growth.     

Effects of magnetic interactions in anisotropy of magnetic susceptibility: Models, experiments and implications for igneous rock fabrics quantification

Besides granites of the ilmenite series, in which the anisotropy of magnetic susceptibility (AMS) is mainly controlled by paramagnetic minerals, the AMS of igneous rocks is commonly interpreted as the result of the shape-preferred orientation of unequant ferromagnetic grains. In a few instances, the anisotropy due to the distribution of ferromagnetic grains, irrespective of their shape, has also been proposed as an important AMS source. Former analytical models that consider infinite geometry of identical and uniformly magnetized and coaxial particles confirm that shape fabric may be overcome by dipolar contributions if neighboring grains are close enough to each other to magnetically interact. On these bases we present and experimentally validate a two-grain macroscopic numerical model in which each grain carries its own magnetic anisotropy, volume, orientation and location in space. Compared with analytical predictions and available experiments, our results allow to list and quantify the factors that affect the effects of magnetic interactions. In particular, we discuss the effects of (i) the infinite geometry used in the analytical models, (ii) the intrinsic shape anisotropy of the grains, (iii) the relative orientation in space of the grains, and (iv) the spatial distribution of grains with a particular focus on the inter-grain distance distribution. Using documented case studies, these findings are summarized and discussed in the framework of the generalized total AMS tensor recently introduced by Cañon-Tapia (Cañon-Tapia, E., 2001. Factors affecting the relative importance of shape and distribution anisotropy in rocks: theory and experiments. Tectonophysics, 340, 117–131.). The most important result of our work is that analytical models far overestimate the role of magnetic interaction in rock fabric quantification. Considering natural rocks as an assemblage of interacting and non-interacting grains, and that the effects of interaction are reduced by (i) the finite geometry of the interacting clusters, (ii) the relative orientation between interacting grains, (iii) their heterogeneity in orientation, shape and bulk susceptibility, and (iv) their inter-distance distribution, we reconcile analytical models and experiments with real case studies that minimize the role of magnetic interaction onto the measured AMS. Limitations of our results are discussed and guidelines are provided for the use of AMS in geological interpretation of igneous rock fabrics where magnetic interactions are likely to occur.     

Early contamination of European flounder (Platichthys flesus) by PCDD/Fs and dioxin-like PCBs in European waters

Contamination levels and profiles of 7 polychlorinated-p-dioxins, 10 polychlorinated furans (PCDD/Fs) and 12 dioxin-like polychlorinated biphenyls (dl-PCBs) were investigated in juvenile European flounder (Platichthys flesus) captured in different nursery areas in the northeastern Atlantic coast across its geographical distribution range. The toxic equivalent concentrations (WHO-TEQ_fish) were also determined in order to evaluate which P. flesus population was more exposed to dioxin-like toxicity. Juveniles caught in the Sørfjord (Norway) showed the lowest WHO-TEQ_fish concentration (0.052 pg WHO-TEQ_fish g⁻¹ wet weight) whereas the highest value was observed in fish from the Wadden Sea (The Netherlands; 0.291 pg WHO-TEQ_fish g⁻¹ ww), mainly due to the greater contribution of 2,3,7,8-tetrachlorodibenzo-p-dioxin, the most toxic congener. Nonetheless, when comparing the results with existent tissue residue-based toxicity benchmarks, no adverse effects resulting from PCDD/Fs and dl-PCBs are expected to occur in flounder from the studied systems.