Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

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.     

Optimizing interpolation of shoot density data from a Posidonia oceanica seagrass bed

A case study on the optimization of Posidonia oceanica density interpolation, using a data set from a large meadow at Porto Conte Bay (NW Sardinia, Italy), is presented. Ordinary point kriging, cokriging and a weighted average based on inverse square distance were used to interpolate density data measured in 36 sampling stations. The results obtained from different methods were then compared by means of a leave‐one‐out cross‐validation procedure. The scale at which interpolation was carried out was defined on the basis of the Hausdorff dimension of the variogram. Optimizing spatial scale and data points search strategy allowed obtaining more accurate density estimates independently of the interpolation method.     

Mt. Melbourne volcano,antarctica: Evidence of seismicity related to volcanic activity

Characterization of the microseismic activity (M _L <2.0) has been performed at Mt. Melbourne since 1990. We recorded a group of low frequency events with common morphological characteristics, i.e., an emerging onset, an unclear second phase and a sharply dropping coda. Spectral analysis of events recorded at more than one station indicates that the seismogram characteristics and spectral content are largely due to source effects. A polarization filter applied to a set of three component data revealed a first phase made up ofP waves followed (after about 0.9–1.4 sec.) by a second phase probably composed ofSH-type waves. Particle motion analysis detected a seismic ray angle direction mainly between N70°E and N110°E and apparent angle of incidence between 35° and 48° for the first phase. The studied seismicity was localized in an area on the eastern slope of Mt. Melbourne Volcano which presents a surface temperature anomaly (Mazzarini andSalvini, 1994). We formulate two hypotheses for the type of earthquakes recorded: 1) long-period events involving active presence of magmatic fluids in the source processes; 2) or the result of fracturing processes (shear?) in a medium characterized by transition between brittle and plastic behaviors. In the latter hypothesis the superficial thermal anomaly may be a symptom of this behavior at depth and is confirmed by the lown values observed for the exponential fit in the codaQ analysis.     

A calderas field at the Marifil Formation, new volcanogenic interpretation, Norpatagonian Massif, Argentina

Major volcanogenic structures within the Marifil Formation allowed to determine the location of a claderas field that spreads along 400 kilometers. This field has at least three large calderas with diameters near to 100 kilometers, recognizable by regional distribution of associated dikes and rhyolite lava flows present at the caldera border that separates a monotonous ignimbritic plateau from a depressed (150 meters) inner volcanic collapsed caldera. The Marifil ignimbritic plateau cover more than 50,000 square kilometers with thicknesses that reach 800 meters. Associated with these plateau ignimbrites there is an important fluorite mineralization.     

Slab–mantle flow interaction: influence on subduction dynamics and duration

We investigate the influence of mantle flow relative to the lithosphere on subduction dynamics. We use 2D thermo‐mechanical models assuming incompressible non‐Newtonian fluid rheology. Different mantle flow velocities consistent with absolute plate motion models are tested, as well as both directions of flow, either sustaining or opposing slab dip. The effects of different inflow/outflow velocity profiles, slab strengths and upper–lower mantle viscosity contrasts are also evaluated. Slab dip deviations between models with opposite mantle flow directions range from 37° for relatively strong slabs (η_max = 10²⁵ Pa s) to 50° for weaker slabs (η_max = 10²⁴ Pa s), accounting for a significant amount of natural slab dip variability. For imposed mantle flow supporting the slab, the initial stage of slab steepening is followed by a stage of continuous slab dip decrease. This slab shallowing eventually leads to mantle wedge closure, subduction cessation and slab break‐off, possibly driving subduction flips.     

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.