Using Simple Dilution Models to Predict New Zealand Estuarine Water Quality

A tool based on simple dilution models is developed to predict potential nutrient concentrations and flushing times for New Zealand estuaries. Potential nutrient concentrations are the concentrations that would occur in the absence of nutrient uptake or losses through biogeochemical processes, and so represent the pressure on a system due to nutrient loading. The dilution modelling approach gives a single time- and space-averaged concentration as a function of flow and nutrient input, with the capability to include seasonal nutrient and flow differences. This tool is intended to be used to identify estuaries likely to be highly sensitive to current nutrient loads based on their physical attributes, or to quickly compare the effects of different land-use scenarios on estuaries. The dilution modelling approach is applied both to a case study of a single New Zealand estuary, and used in a New Zealand-wide assessment of 415 estuaries. For the NZ-wide assessment, annual nutrient loads to each estuary were obtained from a GIS-based land-use model. Comparison with measured data shows that the predicted potential nitrate concentrations are significantly correlated with, but higher than, measured nitrate values from water quality sampling time series. This is consistent with expectations given that the measured concentrations include the effects of nitrogen uptake and loss. The estuary dilution modelling approach is currently incorporated into the GIS-land use model, and is also available as a web-app for assessing eutrophication susceptibility of New Zealand estuaries.     

Floodgate Operations and Fish Communities in Tidal Creeks of the Lower Fraser River (British Columbia,Canada)

Tidal creeks in large coastal deltas can be important habitat for fish but are often highly modified by human activities. Connectivity between tributary creeks and mainstem channels is often constrained by structures such as dikes and floodgates, designed to protect urban and agricultural areas from flooding. While they play important roles in flood mitigation, floodgates can diminish habitat quality and block fish from accessing tidal creeks. It is likely that floodgates differ in their operations and may consequently open for different amounts of time; however, floodgate operations and their effects are not well quantified. We asked the question: how does the mechanical functioning of these floodgates affect fish communities in tidal creeks? We used time-lapse cameras and quantified the timing of gate openings for 22 tributaries of the Lower Fraser River in British Columbia, Canada, and related these operational data to differences in fish communities above and below floodgates. Floodgate operations varied substantially, with some floodgates opening daily while others opened less than 20% of the day, on average. Sites with floodgates that seldom opened were associated with greater differences in fish communities and with reduced upstream native species richness by about one species on average. Where floodgates opened infrequently, we also found lower upstream dissolved oxygen concentrations than at sites where floodgates opened for longer periods of time. Thus, floodgate operations can influence fish communities as well as water quality. These data indicate a large scope for improving floodgate operations for connectivity.     

Formation Mechanism of Mud Bank Along the Southwest Coast of India

Mud bank formation during the southwest monsoon along the southwest coast of India remains an enigma to the researchers and coastal community in spite of several earlier studies. The present study attempts to unravel the mystery through a high-frequency, season-long time-series observation at Alappuzha, located at the southern part of the west coast of India, a region of frequent occurrence of mud bank. Using 7-month-long weekly time-series observation, we identified strong winds and high waves associated with onset of the southwest monsoon and subsequent three episodic atmospheric low-pressure events (LPEs).With the help of in situ time-series data, we show that the strong winds and high waves associated with southwest monsoon pre-conditions the near shore bottom sediment to bring it into suspension. The high amplitude waves associated with the southwest monsoon, while propagating from the deep water to shallow water region, interact with the bottom initiating bottom-sediment movement and its suspension due to wave refraction and shoaling. The sporadic occurrence of the atmospheric LPEs enhances the process of suspension of bottom sediment in the near shore region leading to the formation of fluid mud. Simulations with a cohesive sediment transport model yielded realistic estimates of sediment transport, in the presence of an onshore current, a pre-requisite for transporting the fluid mud toward the coast. The prevailing onshore upwelling current during the southwest monsoon provides the favorable pre-requisite conditions for transporting the fluid mud through depression channel network towards the coast. Once sufficient quantity and thickness of fluid mud is accumulated in the near shore region, it acts as a wave damper for subsequent high monsoon waves, as indicated by the time-series wave data, leading to the formation of tranquil mud bank region. Depression channel networks extending from the shelf to the coast off Alappuzha, Kochi, Ponnani, Beypore, and Ullal were found in the bathymetric charts, thus explaining why mud banks occur only at few locations in spite of the prevalence of similar monsoon conditions.     

Potential Salinity and Temperature Futures for the Chesapeake Bay Using a Statistical Downscaling Spatial Disaggregation Framework

Estuaries are productive and ecologically important ecosystems, incorporating environmental drivers from watersheds, rivers, and the coastal ocean. Climate change has potential to modify the physical properties of estuaries, with impacts on resident organisms. However, projections from general circulation models (GCMs) are generally too coarse to resolve important estuarine processes. Here, we statistically downscaled near-surface air temperature and precipitation projections to the scale of the Chesapeake Bay watershed and estuary. These variables were linked to Susquehanna River streamflow using a water balance model and finally to spatially resolved Chesapeake Bay surface temperature and salinity using statistical model trees. The low computational cost of this approach allowed rapid assessment of projected changes from four GCMs spanning a range of potential futures under a high CO₂ emission scenario, for four different downscaling methods. Choice of GCM contributed strongly to the spread in projections, but choice of downscaling method was also influential in the warmest models. Models projected a ~2–5.5 °C increase in surface water temperatures in the Chesapeake Bay by the end of the century. Projections of salinity were more uncertain and spatially complex. Models showing increases in winter-spring streamflow generated freshening in the Upper Bay and tributaries, while models with decreased streamflow produced salinity increases. Changes to the Chesapeake Bay environment have implications for fish and invertebrate habitats, as well as migration, spawning phenology, recruitment, and occurrence of pathogens. Our results underline a potentially expanded role of statistical downscaling to complement dynamical approaches in assessing climate change impacts in dynamically challenging estuaries.     

Upstream—Downstream Shifts in Peak Recruitment of the Native Olympia Oyster in San Francisco Bay During Wet and Dry Years

Understanding the conditions that drive variation in recruitment of key estuarine species can be important for effective conservation and management of their populations. The Olympia oyster (Ostrea lurida) is native to the Pacific coast of North America and has been a target of conservation efforts, though relatively little information on larval recruitment exists across much of its range. This study examined the recruitment of Olympia oysters at biweekly to monthly intervals at four sites in northern San Francisco Bay from 2010 to 2015 (except 2013). Mean monthly temperatures warmed at all sites during the study, while winter (January–April) mean monthly salinity decreased significantly during a wet year (2011), but otherwise remained high as a result of a drought. A recurring peak in oyster recruitment was identified in mid-estuary, in conditions corresponding to a salinity range of 25–30 and >16 °C at the time of settlement (April–November). Higher average salinities and temperatures were positively correlated with greater peak recruitment. Interannual variation in the timing of favorable conditions for recruitment at each site appears to explain geographic and temporal variation in recruitment onset. Higher winter/spring salinities and warmer temperatures at the time of recruitment corresponded with earlier recruitment onset within individual sites. Across all sites, higher winter/spring salinities were also correlated with earlier onset and earlier peak recruitment. Lower winter salinities during 2011 also resulted in a downstream shift in the location of peak recruitment.     

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.