Wind-Driven Dissolved Organic Matter Dynamics in a Chesapeake Bay Tidal Marsh-Estuary System

Controls on organic matter cycling across the tidal wetland-estuary interface have proved elusive, but high-resolution observations coupled with process-based modeling can be a powerful methodology to address shortcomings in either methodology alone. In this study, detailed observations and three-dimensional hydrodynamic modeling are used to examine biogeochemical exchanges in the marsh-estuary system of the Rhode River, MD, USA. Analysis of observations near the marsh in 2015 reveals a strong relationship between marsh creek salinity and dissolved organic matter fluorescence (fDOM), with wind velocity indirectly driving large amplitude variation of both salinity and fDOM at certain times of the year. Three-dimensional model results from the Finite Volume Community Ocean Model implemented for the wetland system with a new marsh grass drag module are consistent with observations, simulating sub-tidal variability of marsh creek salinity. The model results exhibit an interaction between wind-driven variation in surface elevation and flow velocity at the marsh creek, with northerly winds driving increased freshwater signal and discharge out of the modeled wetland during precipitation events. Wind setup of a water surface elevation gradient axially along the estuary drives the modeled local sub-tidal flow and thus salinity variability. On sub-tidal time scales (>36 h, <1 week), wind is important in mediating dissolved organic matter releases from the Kirkpatrick Marsh into the Rhode River.     

Invasions in Marine Communities: Contrasting Species Richness and Community Composition Across Habitats and Salinity

While many studies of non-native species have examined either soft-bottom or hard-bottom marine communities, including artificial structures at docks and marinas, formal comparisons across these habitat types are rare. The number of non-indigenous species (NIS) may differ among habitats, due to differences in species delivery (trade history) and susceptibility to invasions. In this study, we quantitatively compared NIS to native species richness and distribution and examined community similarity across hard-bottom and soft-sediment habitats in San Francisco Bay, California (USA). Benthic invertebrates were sampled using settlement panels (hard-bottom habitats) and sediment grabs (soft-bottom habitats) in 13 paired sites, including eight in higher salinity areas and five in lower salinity areas during 2 years. Mean NIS richness was greatest in hard-bottom habitat at high salinity, being significantly higher than each (a) native species at high salinity and (b) NIS richness at low salinity. In contrast, mean NIS richness in soft-bottom communities was not significantly different from native species richness in either high- or low-salinity waters, nor was there a difference in NIS richness between salinities. For hard-bottom communities, NIS represented an average of 79% of total species richness per sample at high salinity and 78% at low salinity, whereas the comparable values for soft bottom were 46 and 60%, respectively. On average, NIS occurred at a significantly higher frequency (percent of samples) than native species for hard-bottom habitats at both salinities, but this was not the case for soft-bottom habitats. Finally, NIS contributed significantly to the existing community structure (dissimilarity) across habitat types and salinities. Our results show that NIS richness and occurrence frequency is highest in hard-bottom and high-salinity habitat for this Bay but also that NIS contribute strongly to species richness and community structure across each habitat evaluated.     

Site scale modeling of slow-moving landslides,a 3D viscoplastic finite element modeling approach

This paper presents an advanced 3D numerical methodology to reproduce the kinematics of slow active landslides, more precisely, to reproduce the nearly constant strain rate (secondary creep) and the acceleration/deceleration of the moving mass due to hydrological changes. For this purpose, finite element analyses are performed in a large area covering a long time-span (12 years), in order to exhibit different interacting slope movements. First, we perform a stability analysis using the shear strength reduction (SSR) technique with a Mohr-Coulomb failure criteria. It is done in order to compute factors of safety (FS) and to identify two different scenarios, the first one being stable (FS > 1) and the second one being unstable (FS < 1). In the studied test case, the Portalet landslide (Central Spanish Pyrenees), the first scenario corresponds to an initial stable configuration of the slope and the second one to an unstable excavated configuration. Second, taking the first scenario as an initial condition, a time-dependent analysis is performed using a coupled formulation to model solid skeleton and pore fluids interaction, and a simplified ground water model that takes into account daily rainfall intensity. In this case, a viscoplastic constitutive model based on Perzyna’s theory is applied to reproduce soil viscous behavior and the delayed creep deformation due to the excavation. The fluidity parameter is calibrated to reproduce displacements measured by the monitoring systems. Our results demonstrate that 3D analyses are preferable to 2D ones for reproducing in a more realistic way the slide behavior. After calibration, the proposed model is able to simulate successfully short- and medium-term predictions during stages of primary and secondary creep.     

An analysis of failure mechanism constraints on pre-failure rock block deformation using TLS and roto-translation methods

Rockfall hazards increase the risk of train derailment along railway corridors in western Canada. In this study, repeated terrestrial laser scanning (TLS) datasets were collected every 2–3 months at three different sites along the Thompson and Fraser River corridors in British Columbia, referred to as the Goldpan, White Canyon, and Mile 109 sites. A total of 207 rockfall events occurring across all three sites between November 11, 2014 and October 18, 2016 were recorded in a database. For each of these rockfalls, pre-failure deformation was measured using a method of three-dimensional roto-translation block tracking. Each rockfall was classified by its deformation behaviour and further categorised based on failure mechanism, volume, lithology, and the roughness condition of the failure plane. Results reveal that detectable levels of deformation were measured in 33% of the total number of rockfall events using the present methods. Rotation deformation was most commonly observed in toppling failures with relatively steep joint orientations. Conversely, planar sliding blocks generally exhibited the least measurable deformation, with the majority not showing any precursory translation or rotation. It is postulated that overhanging rockfall configurations may suppress the expression of deformation in rockfall source blocks, though additional research is required to confirm this.     

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.     

Were merely storm-landslides driven by the 2015-2016 Niño in the Mendoza River valley?

Since Holocene time, above-mean precipitations recorded during the El Niño warm ENSO phase have been linked to the occurrence of severe debris flows in the arid Central Andes. The 2015–2016 El Niño, for its unusual strength, began driving huge and dangerous landslides in the Central Andes (32°) in the recent South Hemisphere summer. The resulting damages negatively impacted the regional economy. Despite this, causes of these dangerous events were ambiguously reported. For this reason, a multidisciplinary study was carried out in the Mendoza River valley. Firstly, a geomorphological analysis of affected basins was conducted, estimating morphometric parameters of recorded events such as velocity, stream flow, and volume. Atmospheric conditions during such events were analyzed, considering precipitations, snow cover, temperature range, and the elevation of the zero isotherm. Based on our findings, the role of El Niño on the slope instability in the Central Andes is more complex in the climate change scenario. Even though some events were effectively triggered by intense summer rainstorm following expectations, the most dangerous events were caused by the progressive uplifting of the zero isotherm in smaller basins where headwaters are occupied by debris rock glaciers. Our research findings give light to the dynamic coupled system ENSO–climate change–landslides (ECCL) at least in this particular case study of the Mendoza River valley. Landslide activity in this Andean region is driven by wetter conditions linked to the ENSO warm phase, but also to progressive warming since the twentieth century in the region. This fact emphasizes the future impact of the natural hazards on Andean mountain communities.     

Study of landslides at the archaeological site of Abakainon necropolis (NE Sicily) by geomorphological and geophysical investigations

This paper reports on the investigation of the old landslide body on which the necropolis of Abakainon, an Italian archaeological site, was built about 2400 years ago. This site is located in a geologically complex area of northeastern Sicily, where active tectonics has strongly modified the morphological landscape, as testified by the meandering path of watercourses, widespread presence of paleosurfaces, and deviations of crest line. Previous studies have hypothesized that the necropolis, suddenly destroyed by a catastrophic event in historical times, could have been built on an old landslide body. This represents the starting point of this paper, focused on geomorphological analysis and geophysical surveys, which allowed discovering peculiar elements in support of the landslide theory. In particular, the elaboration of a Digital Terrain Model allowed locating the main morphostructural elements, thus highlighting how the occurrence of landslides in the study area is strongly linked to recent tectonics. Passive seismic surveys, carried out through the acquisition of ambient noise, allowed detecting the presence of impedance contrasts, which were associated to peculiar features, such as the landslide body itself. Moreover, local strong directional effects testify the presence of faults, affecting even the necropolis area. Such results, further corroborated by the outcomes of an electrical resistivity survey, lay the foundations for interesting consideration on the possible cause of the destruction of this site, probably related to the occurrence of a historical earthquake in the first century AD.     

Activity and kinematic behaviour of deep-seated landslides from PS-InSAR displacement rate measurements

Large landslides and deep-seated gravitational slope deformations (DSGSD) represent an important geo-hazard in relation to the deformation of large structures and infrastructures and to the associated secondary landslides. DSGSD movements, although slow (from a few millimetres to several centimetres per year), can continue for very long periods, producing large cumulative displacements and undergoing partial or complete reactivation. Therefore, it is important to map the activity of such phenomena at a regional scale. Ground surface displacements at DSGSD typically range close to the detection limit of monitoring equipment but are suitable for synthetic aperture radar (SAR) interferometry. In this paper, permanent scatterers (PSInSAR?) and SqueeSAR? techniques are used to analyse the activity of 133 DSGSD, in the Central Italian Alps. Statistical indicators for assigning a degree of activity to slope movements from displacement rates are discussed together with methods for analysing the movement and activity distribution within each landslide. In order to assess if a landslide is active or not, with a certain degree of reliability, three indicators are considered as optimal: the mean displacement rate, the activity index (ratio of active PS, displacement rate larger than standard deviation, overall PS) and the nearest neighbor ratio, which allows to describe the degree of clustering of the PS data. According to these criteria, 66% of the phenomena are classified as active in the monitored period 1992–2009. Finally, a new methodology for the use of SAR interferometry data to attain a classification of landslide kinematic behaviour is presented. This methodology is based on the interpretation of longitudinal ground surface displacement rate profiles in the light of numerical simulations of simplified failure geometries. The most common kinematic behaviour is rotational, amounting to 41 DSGSDs, corresponding to the 62.1% of the active phenomena.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Relationships between natural terrain landslide magnitudes and triggering rainfall based on a large landslide inventory in Hong Kong

Rain-induced landslides are recognized as one of the most catastrophic hazards on hilly terrains. To develop strategies for landslide risk assessment and management, it is necessary to estimate not only the rainfall threshold for the initiation of landslides, but also the likely magnitudes of landslides triggered by a storm of a given intensity. In this study, the frequency distributions of both open hillside landslides and channelized debris flows in Hong Kong are established on the basis of the Enhanced Natural Terrain Landslide Inventory (ENTLI) with 19,763 records in Hong Kong up to 2013. The landslide magnitudes are measured in terms of the number, scar area, volume, or density of landslides. The mean values of the scar areas and volumes are 55.2 m² and 102.0 m³, respectively, for the open hillside landslides and 91.3 m² and 166.5 m³, respectively, for the channelized debris flows. Empirical correlations between the numbers, scar areas, and volumes of hillside landslides or channelized debris flows and the maximum rolling rainfall intensities of different periods have been derived. The maximum rolling 4- to 24-h rainfall amounts provide better predictions compared with those with the maximum rolling 1-h rainfall. Maximum rolling rainfall intensity-duration thresholds identifying the likely rainfall conditions that yield natural terrain landslides or debris flows of different magnitudes are also proposed. The initiation rainfall thresholds are identified as 75, 90, 100, 120, 150, 180, and 200 mm for the maximum rolling 1-, 2-, 4-, 6-, 8-, 12-, and 24-h rainfall, respectively.