Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

The distribution of chlorophyll-a and dominant planktonic components in the coastal transition zone off Concepción, central Chile, during different oceanographic conditions

The oceanographic setting and the planktonic distribution in the coastal transition zone off Concepción (∼35-38°S, ∼73-77°W), an area characterized by its high biological production, were assessed during two different seasons: austral spring with equatorward upwelling favorable winds and austral winter with predominately northerly winds. Oceanographic and biological data (total chlorophyll-a, particulate organic carbon, microplankton, large mesozooplankton >500 μm as potential consumers of microplankton) were obtained during two cruises (October 1998, July 1999) together with satellite imagery for wind stress, geostrophic flow, surface temperature, and chlorophyll-a data. The physical environment during the spring sampling was typical of the upwelling period in this region, with a well-defined density front in the shelf-break area and high concentrations of surface chlorophyll-a (>5 mg m⁻³) on the shelf over the Itata terrace. During the winter sampling, highly variable though weakly upwelling-favorable winds were observed along with lower surface chlorophyll-a values (<2 mg m⁻³) on the shelf. In the oceanic area (>100 km from the coast), cyclonic and anti-cyclonic eddies were evident in the flow field during both periods, the former coinciding with higher chlorophyll-a contents (∼1 mg m⁻³) than in the surrounding waters. Also, a cold, chlorophyll-a rich filament was well defined during the spring sampling, extending from the shelf out to 350-400 km offshore. Along a cross-shelf transect, the micro- and meso-planktonic assemblages displayed higher coastal abundances during the spring cruise but secondary peaks appeared in the oceanic area during the winter cruise, coinciding with the distribution of the eddies. These results suggest that the mesoscale features in this region, in combination with upwelling, play a role in potentially increasing the biological productivity of the coastal transition zone off Concepción.     

Groundwater flow systems in turbidites of the Northern Apennines (Italy): natural discharge and high speed railway tunnel drainage

Turbidites crop out extensively in the Northern Apennine mountains (Italy). The huge amounts of groundwater drained by tunnels, built for the high speed railway connection between Bologna and Florence, demonstrate the aquifer-like behaviour of these units, up to now considered as aquitards. A conceptual model of groundwater flow systems (GFS) in fractured aquifers of turbidites is proposed, taking into account both system natural state and the perturbation induced by tunnel drainage. Analysis of hydrological data (springs, streams and tunnel discharge), collected over 10 years, was integrated with analysis of hydrochemical and isotopic data and a stream-tunnel tracer test. Hydrologic recession analysis of undisturbed conditions is a key tool in studying turbiditic aquifer hydrogeology, permitting the discrimination of GFS, the estimation of recharge relative to the upstream reach portion and the identification of springs most vulnerable to tunnel drainage impacts. The groundwater budgeting analysis provides evidence that the natural aquifer discharge was stream-focused through GFS, developed downslope or connected to main extensional tectonic lineaments intersecting stream beds; now tunnels drain mainly active recharge groundwater and so cause a relevant stream baseflow deplenishment (approximately two-thirds of the natural value), possibly resulting in adverse effects on local ecosystems.     

Salinity and periodic inundation controls on the soil‐plant‐atmosphere continuum of gray mangroves

Salinity and periodic inundation are both known to have a major role in shaping the ecohydrology of mangroves through their controls on water uptake, photosynthesis, stomatal conductance, gas exchanges, and nutrient availability. Salinity, in particular, can be considered one of the main abiotic regulating factors for halophytes and salt‐tolerant species, due to its influence on water use patterns and growth rate. Ecohydrological literature has rarely focused on the effects of salinity on plant transpiration, based on the fact that the terrestrial plants mostly thrive in low‐saline, unsaturated soils where the role of osmotic potential can be considered negligible. However, the effect of salinity cannot be neglected in the case of tidal species like mangroves, which have to cope with hyperosmotic conditions and waterlogging. We introduce here a first‐order ecohydrological model of the soil/plant‐atmosphere continuum of Avicennia marina—also known as gray mangrove—a highly salt‐tolerant pioneer species able to adapt to hyperarid intertidal zones and characterized by unique morphological and ecophysiological traits. The A. marina's soil‐plant‐atmosphere continuum takes explicitly into account the role of water head, osmotic water potential, and water salinity in governing plant water fluxes. A. marina's transpiration is thus modeled as a function of salinity based on a simple parameterization of salt exclusion mechanisms at the root level and a modified Jarvis' expression accounting for the effects of salinity on stomatal conductance. Consistently with previous studies investigating the physiology of mangroves in response to different environmental drivers, our results highlight the major influence of salinity on mangrove transpiration when contrasted with other potential stressors such as waterlogging and water stress.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Short‐ and long‐term evapotranspiration rates at ecological restoration sites along a large river receiving rare flow events

Many large rivers around the world no longer flow to their deltas, due to ever greater water withdrawals and diversions for human needs. However, the importance of riparian ecosystems is drawing increasing recognition, leading to the allocation of environmental flows to restore river processes. Accurate estimates of riparian plant evapotranspiration (ET) are needed to understand how the riverine system responds to these rare events and achieve the goals of environmental flows. In 2014, historic environmental flows were released into the Lower Colorado River at Morelos Dam (Mexico); this once perennial but now dry reach is the final stretch to the mighty Colorado River Delta. One of the primary goals was to supply native vegetation restoration sites along the reach with water to help seedlings establish and boost groundwater levels to foster the planted saplings. Patterns in ET before, during, and after the flows are useful for evaluating whether this goal was met and understanding the role that ET plays in this now ephemeral river system. Here, diurnal fluctuations in groundwater levels and Moderate Resolution Imaging Spectroradiometer (MODIS) data were used to compare estimates of ET specifically at 3 native vegetation restoration sites during 2014 planned flow events, and MODIS data were used to evaluate long‐term (2002–2016) ET responses to restoration efforts at these sites. Overall, ET was generally 0–10 mm d^?1 across sites, and although daily ET values from groundwater data were highly variable, weekly averaged estimates were highly correlated with MODIS‐derived estimates at most sites. The influence of the 2014 flow events was not immediately apparent in the results, although the process of clearing vegetation and planting native vegetation at the restoration sites was clearly visible in the results.     

Carbon fluxes through major phytoplankton groups during the spring bloom and post-bloom in the Northwestern Mediterranean Sea

The carbon flux through major phytoplankton groups, defined by their pigment markers, was estimated in two contrasting conditions of the Northwestern Mediterranean open ocean ecosystem: the spring bloom and post-bloom situations (hereafter Bloom and Post-bloom, respectively). During Bloom, surface chlorophyll a (Chl a) concentration was higher and dominated by diatoms (53% of Chl a), while during Post-bloom Synechococcus (42%) and Prymnesiophyceae (29%) became dominant. The seawater dilution technique, coupled to high pressure liquid chromatography (HPLC) analysis of pigments and flow cytometry (FCM), was used to estimate growth and grazing rates of major phytoplankton groups in surface waters. Estimated growth rates were corrected for photoacclimation based on FCM-detected changes in red fluorescence per cell. Given the 30% average decrease in the pigment content per cell between the beginning and the end of the incubations, overlooking photoacclimation would have resulted in a 0.40 d^?1 underestimation of phytoplankton growth rates. Corrected average growth rates (μ_o) were 0.90±0.20 (SD) and 0.40±0.14 d^?1 for Bloom and Post-bloom phytoplankton, respectively. Diatoms, Cryptophyceae and Synechococcus were identified as fast-growing groups and Prymnesiophyceae and Prasinophyceae as slow-growing groups across Bloom and Post-bloom conditions. The higher growth rate during Bloom was due to dominance of phytoplankton groups with higher growth rates than those dominating in Post-bloom. Average grazing rates (m) were 0.58±0.20 d^?1 (SD) and 0.31±0.07 d^?1. The proportion of phytoplankton growth consumed by microzooplankton grazing (m/μ_o) tended to be lower in Bloom (0.69±0.34) than in Post-bloom (0.80±0.08). The intensity of nutrient limitation experienced by phytoplankton indicated by μ_o/μ_n (where μ_n is the nutrient-amended growth rate), was similar during Bloom (0.78) and Post-bloom (0.73). Primary production from surface water (PP) was estimated with ¹⁴C incubations. A combination of PP and Chl a synthesis rate yielded C/Chl a ratios of 34±21 and 168±75 (g:g) for Bloom and Post-bloom, respectively. Transformation of group-specific Chl a fluxes into carbon equivalents confirmed the dominant role of diatoms during Bloom and Synechococcus and Prymnesiophyceae during Post-bloom.     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.