A network theory approach for a better understanding of overland flow connectivity

Hydrological connectivity describes the physical coupling (linkages) of different elements within a landscape regarding (sub‐) surface flows. A firm understanding of hydrological connectivity is important for catchment management applications, for example, habitat and species protection, and for flood resistance and resilience improvement. Thinking about (geomorphological) systems as networks can lead to new insights, which has also been recognized within the scientific community, seeing the recent increase in the use of network (graph) theory within the geosciences. Network theory supports the analysis and understanding of complex systems by providing data structures for modelling objects and their linkages, and a versatile toolbox to quantitatively appraise network structure and properties. The objective of this study was to characterize and quantify overland flow connectivity dynamics on hillslopes in a humid sub‐Mediterranean environment by using a combination of high‐resolution digital‐terrain models, overland flow sensors and a network approach. Results showed that there are significant differences between overland flow connectivity on agricultural areas and semi‐natural shrubs areas. Significant positive correlations between connectivity and precipitation characteristics were found. Significant negative correlations between connectivity and soil moisture were found, most likely because of soil water repellency and/or soil surface crusting. The combination of structural networks and dynamic networks for determining potential connectivity and actual connectivity proved a powerful tool for analysing overland flow connectivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Considering large earthquake clustering in seismic risk analysis

Earthquakes cluster in space and time resulting in nonlinear damage effects. We compute earthquake interactions using the Coulomb stress transfer theory and dynamic vulnerability from the concept of ductility capacity reduction. We combine both processes in the generic multi-risk framework where risk scenarios are simulated using a variant of the Markov chain Monte Carlo method. We apply the proposed approach to the thrust fault system of northern Italy, considering earthquakes with characteristic magnitudes in the range ~[6, 6.5], different levels of tectonic loading \(\dot{\tau }\) = {10⁻⁴, 10⁻³, 10⁻²} bar/year and a generic stock of fictitious low-rise buildings with different ductility capacities μ _Δ = {2, 4, 6}. We describe the process’ stochasticity by non-stationary Poisson earthquake probabilities and by binomial damage state probabilities. We find that earthquake clustering yields a tail fattening of the seismic risk curve, the effect of which is amplified by damage-dependent fragility due to clustering. The impact of clustering alone is in average more important than dynamic vulnerability, the spatial extent of the former phenomenon being greater than of the latter one.

     

Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.     

An Alternative Protocol for Single Zircon Dissolution with Application to (U‐Th‐Sm)/He Thermochronometry

Zircon (U‐Th‐Sm)/He (ZHe) thermochronometry is a powerful tool that has been widely used in geology to constrain the exhumation histories of orogens. In this study, we present an alternative protocol for dissolving zircon grains for determination of parent nuclides. This new alkali fusion procedure developed at the SARM (Service d'Analyse des Roches et des Minéraux) in Nancy, France, is fast (requiring only 2 d, including cleaning steps) and offers several advantages over conventional methods by avoiding: (i) use of HF pressure dissolution and (ii) complete removing of grains from the metal microvials. After dissolution, U, Th and Sm were measured using an ICP‐MS. We tested the new procedure on two different ZHe reference materials, the Fish Canyon Tuff and Buluk Tuff; these provided precision values for ZHe‐age estimations of 9 and 6% (1s), respectively. In addition, using this method, zircons from the Buluk Tuff are shown to be chemically more homogenous and more suitable for assessing the uncertainty of the entire integrated procedure.     

Implementing landslide path dependency in landslide susceptibility modelling

Landslide susceptibility modelling—a crucial step towards the assessment of landslide hazard and risk—has hitherto not included the local, transient effects of previous landslides on susceptibility. In this contribution, we implement such transient effects, which we term “landslide path dependency”, for the first time. Two landslide path dependency variables are used to characterise transient effects: a variable reflecting how likely it is that an earlier landslide will have a follow-up landslide and a variable reflecting the decay of transient effects over time. These two landslide path dependency variables are considered in addition to a large set of conditioning attributes conventionally used in landslide susceptibility. Three logistic regression models were trained and tested fitted to landslide occurrence data from a multi-temporal landslide inventory: (1) a model with only conventional variables, (2) a model with conventional plus landslide path dependency variables, and (3) a model with only landslide path dependency variables. We compare the model performances, differences in the number, coefficient and significance of the selected variables, and the differences in the resulting susceptibility maps. Although the landslide path dependency variables are highly significant and have impacts on the importance of other variables, the performance of the models and the susceptibility maps do not substantially differ between conventional and conventional plus path dependent models. The path dependent landslide susceptibility model, with only two explanatory variables, has lower model performance, and differently patterned susceptibility map than the two other models. A simple landslide susceptibility model using only DEM-derived variables and landslide path dependency variables performs better than the path dependent landslide susceptibility model, and almost as well as the model with conventional plus landslide path dependency variables—while avoiding the need for hard-to-measure variables such as land use or lithology. Although the predictive power of landslide path dependency variables is lower than those of the most important conventional variables, our findings provide a clear incentive to further explore landslide path dependency effects and their potential role in landslide susceptibility modelling.     

The impact of overfishing and El Niño on the condition factor and reproductive success of Peruvian hake, Merluccius gayi peruanus

In this study we examine physiological responses of Peruvian hake (Merluccius gayi peruanus) to changes in their population structure and physical environment during the period 1971–2004. In particular, we assess the relative investment of energy in growth and/or reproduction of small (20–35 cm) and large (35–50 cm) hake. We calculated the (i) condition (Fulton’s K), (ii) gonad and (iii) gut fullness indices for 42,761 female hakes sampled from commercial landings; these indices indicate fish somatic, reproductive and feeding condition, respectively. Using Generalized Additive Models we then examined potential relationships between these indices and sea surface temperature anomalies and date. Drastic energy exhaustion and a decrease in female hake fecundity were observed during El Niño events. The long-term trend showed a general increase in condition factor and a decrease in gonad index for large hake between 1971 and 2004. Small hake exhibited a different trend with an increase in reproductive activity, which was accompanied by an earlier maturation. We hypothesise that the observed low investment of energy in reproduction by large female hake might be related to the lack of large males, due to a sex-selective fishery and the impact of El Niño. We suggest that fishing diminished hake reproductive capacity, modified the sex ratio in favour of females and increased population vulnerability to environmental stress, in particular to the El Niño. The impact of multidecadal variability and predators like the squid, Dosidicus gigas, remain unresolved until longer time series become available.     

Revisiting Peruvian anchovy (Engraulis ringens) trophodynamics provides a new vision of the Humboldt Current system

The Peruvian anchovy or anchoveta (Engraulis ringens) forages on plankton and is a main prey for marine mammals, seabirds, fish, and fishers, and is therefore a key element of the food web in the Humboldt Current system (HCS). Here, we present results from the analysis of 21,203 anchoveta stomach contents sampled during 23 acoustic surveys over the period 1996–2003. Prey items were identified to the genus level, and the relative dietary importance of different prey was assessed by determination of their carbon content. Variability in stomach fullness was examined relative to the diel cycle, the distance from the coast, sea surface temperature, and latitude, using generalized additive models (GAMs). Whereas phytoplankton largely dominated anchoveta diets in terms of numerical abundance and comprised >99% of ingested prey items, the carbon content of prey items indicated that zooplankton was by far the most important dietary component, with euphausiids contributing 67.5% of dietary carbon followed by copepods (26.3%). Stomach fullness data showed that anchoveta feed mainly during daytime between 07h00 and 18h00, although night-time feeding also made a substantial contribution to total food consumption. Stomach fullness also varied with latitude, distance from the coast, and temperature, but with substantial variability indicating a high degree of plasticity in anchoveta feeding behaviour. The results suggest an ecological role for anchoveta that challenges current understanding of its position in the foodweb, the functioning of the HCS, and trophic models of the HCS.     

Trophic significance of the kelp Laminaria digitata (Lamour.) for the associated food web: a between-sites comparison

This study aimed at establishing the trophic significance of the kelp Laminaria digitata for consumers inhabiting two rocky shores of Northern Brittany (France), displaying contrasted ecological conditions. The general trophic structure did not vary between these two sites, with a wide diversity of filter-feeders and predators, and only 14% of the species sampled belonging to the grazers' trophic group. The diversity of food sources fueling the food web appeared also similar. The food webs comprised four trophic levels and the prevalence of omnivory appeared relatively low compared to previous studies in the same area. Conversely, to the food web structure, which did not differ, the biochemical composition of L. digitata differed between the two sites, and was correlated to a larger diversity of grazers feeding on this kelp in sheltered conditions. This indicated that the spatial variability occurring in the nutritive value of L. digitata is likely to deeply affect the functioning of kelp-associated food webs. The contribution of L. digitata-derived organic matter to the diet of filter-feeders inhabiting these two environments was assessed using the mixing model Isosource, which showed the higher contribution of kelp matter in sheltered conditions. These results highlight the spatial variability that may occur in the functioning of kelp-associated food webs. Moreover, this suggests that hydrodynamics is likely to control the availability of kelp-derived organic matter to local filter-feeders, probably through an increase of detritus export in exposed areas.     

Palaeoflood activity and climate change over the last 1400 years recorded by lake sediments in the north‐west European Alps

A high‐resolution sedimentological and geochemical study of a high‐altitude proglacial lake (Lake Blanc, Aiguilles Rouges, 2352 m a.s.l.) revealed 195 turbidites, 190 of which are related to flood events over the last 1400 years. We used the coarsest sediment fraction of each turbidite as a proxy for the intensity of each flood event. Because most flood events at this locality are triggered by localized summer convective precipitation events, the reconstructed sedimentary record reveals changes in the frequency and intensity of such events over the last millennium. Comparisons with other temperature, palaeohydrological and glacier reconstructions in the region suggest that the most intense events occurred during the warmest periods, i.e. during the Medieval Climate Anomaly (AD 800–1300) and the current period of global warming. On a multi‐decadal time scale, almost all the flood frequency peaks seem to correspond to warmer periods, whereas multi‐centennial variations in flood frequency appear to follow the regional precipitation pattern. Consequently, this new Alpine flood record provides further evidence of a link between climate warming and an increase in the frequency and intensity of flooding on a multi‐decadal time scale, whereas the centennial variability in flood frequencies is related to regional precipitation patterns. Copyright © 2013 John Wiley & Sons, Ltd.