Anthropogenic impacts on nitrogen fixation rates between restored and natural Mediterranean salt marshes

To test the effects of site and successional stage on nitrogen fixation rates in salt marshes of the Venice Lagoon, Italy, acetylene reduction assays were performed with Salicornia veneta‐ and Spartina townsendii‐vegetated sediments from three restored (6–14 years) and two natural marshes. Average nitrogen fixation (acetylene reduction) rates ranged from 31 to 343 μmol C₂H₄·m^?2·h^?1 among all marshes, with the greatest average rates being from one natural marsh (Tezze Fonde). These high rates are up to six times greater than those reported from Southern California Spartina marshes of similar Mediterranean climate, but substantially lower than those found in moister climates of the Atlantic US coast. Nitrogen fixation rates did not consistently vary between natural and restored marshes within a site (Fossei Est, Tezze Fonde, Cenesa) but were negatively related to assayed plant biomass within the acetylene reduction samples collected among all marshes. Highest nitrogen fixation rates were found at Tezze Fonde, the location closest to the city of Venice, in both natural and restored marshes, suggesting possible site‐specific impacts of anthropogenic stress on marsh succession.     

Hydrothermal Fluid Sources of the Fengjia Barite–fluorite Deposit in Southeast Sichuan,China: Evidence from Fluid Inclusions and Hydrogen and Oxygen Isotopes

The Fengjia barite–fluorite deposit in southeast Sichuan is a stratabound ore deposit which occurs mainly in Lower Ordovician carbonate rocks. Here we present results from fluid inclusion and oxygen and hydrogen isotope studies to determine the nature and origin of the hydrothermal fluids that generated the deposit. The temperature of the ore‐forming fluid shows a range of 86 to 302 °C. Our detailed microthermometric data show that the temperature during mineralization of the fluorite and barite in the early ore‐forming stage was higher than that during the formation of the calcite in the late ore‐forming stage. The salinity varied substantially from 0.18% to 21.19% NaCl eqv., whereas the density was around 1.00 g/cm³. The fluid composition was mainly H₂O (>91.33%), followed by CO₂, CH₄ and traces of C₂H₆, CO, Ar, and H₂S. The dominant cation was Na⁺ and the dominant anion Cl^‐, followed by Ca²⁺, SO₄^2‐, K⁺, and Mg²⁺, indicating a mid–low‐temperature, mid‐low‐salinity, low‐density NaCl–H₂O system. Our results demonstrate that the temperature decreased during the ore‐forming process and the fluid system changed from a closed reducing environment to an open oxidizing environment. The hydrogen and oxygen isotope data demonstrate that the hydrothermal fluids in the study area had multiple sources, primarily formation water, as well as meteoric water and metamorphic water. Combined with the geological setting and mineralization features we infer that the stratabound barite–fluorite deposits originated from mid–low‐temperature hydrothermal fluids and formed vein filling in the fault zone.     

Spatial and temporal constraints on regional-scale groundwater flow in the Pampa del Tamarugal Basin,Atacama Desert,Chile

Aquifers within the Pampa del Tamarugal Basin (Atacama Desert, northern Chile) are the sole source of water for the coastal city of Iquique and the economically important mining industry. Despite this, the regional groundwater system remains poorly understood. Although it is widely accepted that aquifer recharge originates as precipitation in the Altiplano and Andean Cordillera to the east, there remains debate on whether recharge is driven primarily by near-surface groundwater flow in response to periodic flood events or by basal groundwater flux through deep-seated basin fractures. In addressing this debate, the present study quantifies spatial and temporal variability in regional-scale groundwater flow paths at 20.5°S latitude by combining a two-dimensional model of groundwater and heat flow with field observations and δ¹⁸O isotope values in surface water and groundwater. Results suggest that both previously proposed aquifer recharge mechanisms are likely influencing aquifers within the Pampa del Tamarugal Basin; however, each mechanism is operating on different spatial and temporal scales. Storm-driven flood events in the Altiplano readily transmit groundwater to the eastern Pampa del Tamarugal Basin through near-surface groundwater flow on short time scales, e.g., 10⁰–10¹ years, but these effects are likely isolated to aquifers in the eastern third of the basin. In addition, this study illustrates a physical mechanism for groundwater originating in the eastern highlands to recharge aquifers and salars in the western Pampa del Tamarugal Basin over timescales of 10⁴–10⁵ years.     

Molybdenum Mass Fractions and Isotopic Compositions of International Geological Reference Materials

The double‐spike method with multi‐collector inductively coupled plasma‐mass spectrometry was used to measure the Mo mass fractions and isotopic compositions of a set of geological reference materials including the mineral molybdenite, seawater, coral, as well as igneous and sedimentary rocks. The long‐term reproducibility of the Mo isotopic measurements, based on two‐year analyses of NIST SRM 3134 reference solutions and seawater samples, was ≤ 0.07‰ (two standard deviations, 2s, n = 167) for δ^98/95Mo. Accuracy was evaluated by analyses of Atlantic seawater, which yielded a mean δ^98/95Mo of 2.03 ± 0.06‰ (2s, n = 30, relative to NIST SRM 3134 = 0‰) and mass fraction of 0.0104 ± 0.0006 μg g^?1 (2s, n = 30), which is indistinguishable from seawater samples taken world‐wide and measured in other laboratories. The comprehensive data set presented in this study serves as a reference for quality assurance and interlaboratory comparison of high‐precision Mo mass fractions and isotopic compositions.     

Quantifying multiple Permian–Recent exhumation events during the break‐up of eastern Gondwana: sonic transit time analysis of the central and southern Perth Basin

The central and southern Perth Basin in southwestern Australia has a geological history involving multiple regional unconformity‐forming events from the Permian to Recent. This study uses sonic transit time analysis to quantify the magnitudes of net and gross exhumation for four stratigraphic periods from 43 wells. Most importantly, we quantify gross exhumation of the Permian–Triassic, Triassic–Jurassic, Valanginian break‐up and post‐Early Cretaceous events. Post‐Early Cretaceous gross exhumation averages 900‐m offshore and 600‐m onshore. Up to 200 m of this exhumation may be attributed to localized fault block rotation during extension in the Late Cretaceous and/or reverse fault re‐activation due to the compressive stresses in Australia in the last 50 Ma. The remainder is attributed to regional exhumation caused by epeirogenic processes either during the Cenozoic or at the Aptian–Albian boundary. Maximum burial depths prior to the Valanginian unconformity‐forming event were less than those reached subsequently, so that the magnitude of Valanginian break‐up exhumation cannot be accurately quantified. Gross exhumation prior to the break‐up of Gondwana was defined by large magnitude differences (up to 2500 m) between adjoining sub‐basins. At the end of Triassic, exhumation is primarily attributed to reverse re‐activation of faults that were driven by short‐wavelength inversion and exhumation at the end Permian is likely caused by uplift of rotated fault blocks during extension. The evidence from quantitative exhumation analysis indicates a switch in regime, from locally heterogeneous before break‐up to more regionally homogeneous after break‐up.     

MORPHOMETRIC ANALYSIS OF KARST DEPRESSIONS ON A MEDITERRANEAN KARST MASSIF

Dolines are the most typical landform in karst landscapes. They play an important role in karst hydrogeology; they are effective sediment traps that provide information about past climates and karst evolution, and can also be indicators of tectonic activity. The morphometic analysis of those depressions can help in their study and interpretation. However, detailed identification and delineation of dolines becomes increasingly difficult when the study area is very large, of limited access or covered by vegetation. Here, we have used a method developed by the authors for automatic identification and delineation of karst depressions based on the use of a pit removal algorithm. Next a complete morphometric analysis of the karst depressions has been carried out. The main advantage of the methodology is that by delineating and identifying the depth of each doline from its rim, many morphometric statistics can be obtained, including the relation of doline depth and area with respect to altitude, the main doline orientations, doline density, the fractal dimension of doline area and frequency, and point field analysis, among others. The results are interpreted from a geomorphological point of view and the mapped dolines can be easily integrated into geomorphological mapping. A case study is presented from the large Sierra Gorda karst massif in southern Spain, where 3100 depressions have been identified. A strong structural control of doline shape and orientation is demonstrated; the main areas of doline density are highlighted and doline size is shown to follow a fractal law.     

Adjoint analysis of Buckley-Leverett and two-phase flow equations

This paper analyzes the adjoint equations and boundary conditions for porous media flow models, specifically the Buckley-Leverett equation, and the compressible two-phase flow equations in mass conservation form. An adjoint analysis of a general scalar hyperbolic conservation law whose primal solutions include a shock jump is initially presented, and the results are later specialized to the Buckley-Leverett equation. The non-convexity of the Buckley-Leverett flux function results in adjoint characteristics that are parallel to the shock front upstream of the shock and emerge from the shock front downstream of the shock. Thus, in contrast to the behavior of Burgers’ equation where the adjoint is continuous at a shock, the Buckley-Leverett adjoint, in general, contains a discontinuous jump across the shock. Discrete adjoint solutions from space-time discontinuous Galerkin finite element approximations of the Buckley-Leverett equation are shown to be consistent with the derived closed-form analytical solutions. Furthermore, a general result relating the adjoint equations for different (though equivalent) primal equations is used to relate the two-phase flow adjoints to the Buckley-Leverett adjoint. Adjoint solutions from space-time discontinuous Galerkin finite element approximations of the two-phase flow equations are observed to obey this relationship.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

Effect of loading duration on uncertainty in creep analysis of clay

Some studies suggest that creep parameters should be determined using a greater quantity of creep test data to provide more reliable prediction regarding the deformation of soft soils. This study aims to investigate the effect of loading duration on model updating. One‐dimensional consolidation data of intact Vanttila clay under different loading durations collected from the literature is used for demonstration. The Bayesian probabilistic method is used to identify all unknown parameters based on the consolidation data during the entire consolidation process, and their uncertainty can be quantified through the obtained posterior probability density functions. Additionally, the optimal models are also determined from among 9 model candidates. The analyses indicate that the optimal models can describe the creep behavior of intact soft soils under different loading durations, and the adopted method can evaluate the effect of loading duration on uncertainty in the creep analysis. The uncertainty of a specific model and its model parameters decreases as more creep data are involved in the updating process, and the updated models that use more creep data can better capture the deformation behavior of an intact sample. The proposed method can provide quantified uncertainty in the process of model updating and assist engineers to decide whether the creep test data are sufficient for the creep analysis.