skip to main content


Title: Reconstructing primary production in a changing estuary: A mass balance modeling approach
Abstract

Estuarine primary production (PP) is a critical rate process for understanding ecosystem function and response to environmental change. PP is fundamentally linked to estuarine eutrophication, and as such should respond to ongoing efforts to reduce nutrient inputs to estuaries globally. However, concurrent changes including warming, altered hydrology, reduced input of sediments, and emergence of harmful algal blooms (HABs) could interact with nutrient management to produce unexpected changes in PP. Despite its fundamental importance, estuarine PP is rarely measured. We reconstructed PP in the York River Estuary with a novel mass balance model based on dissolved inorganic nitrogen (DIN) for the period 1994–2018. Modeled PP compared well to previous estimates and demonstrated a long‐term increase and down‐estuary shift over the study period. This increase occurred despite reductions in discharge, flushing time, DIN loading, and DIN standing stock over the same period. Increased PP corresponded to increased water temperature, decreased turbidity and light attenuation, and increased photic depth and assimilation ratio, suggesting that phytoplankton in the York River Estuary have become more efficient at converting nutrients into biomass primarily due to a release from light limitation. The increase in PP also coincided with the increasing occurrence of late summer HABs in the lower York River Estuary, including the emergence of a second bloom‐forming dinoflagellate in 2007. Results demonstrate how changes concurrent with nutrient management could alter expected system responses and illustrate the utility of the mass balance approach for estimating critical rate processes like PP in the absence of observations.

 
more » « less
Award ID(s):
1737258 1655221
NSF-PAR ID:
10443378
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Limnology and Oceanography
Volume:
66
Issue:
6
ISSN:
0024-3590
Page Range / eLocation ID:
p. 2535-2546
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Modifications for navigation since the late 1800s have increased channel depth (H) in the lower Hudson River estuary by 10–30%, and at the mouth the depth has more than doubled. Observations along the lower estuary show that both salinity and stratification have increased over the past century. Model results comparing predredging bathymetry from the 1860s with modern conditions indicate an increase in the salinity intrusion of about 30%, which is roughly consistent with theH5/3scaling expected from theory for salt flux dominated by steady exchange. While modifications including a recent deepening project have been concentrated near the mouth, the changes increase salinity and threaten drinking water supplies more than 100 km landward. The deepening has not changed the responses to river discharge (Qr) of the salinity intrusion (~Qr−1/3) or mean stratification (Qr2/3). Surprisingly, the increase in salinity intrusion with channel deepening results in almost no change in the estuarine circulation. This contrasts sharply with local scaling based on local dynamics of anH2dependence, but it is consistent with a steady state salt balance that allows scaling of the estuarine circulation based on external forcing factors and is independent of depth. In contrast, the observed and modeled increases in stratification are opposite of expectations from the steady state balance, which could be due to reduction in mixing with loss of shallow subtidal regions. Overall, the mean shift in estuarine parameter space due to channel deepening has been modest compared with the monthly‐to‐seasonal variability due to tides and river discharge.

     
    more » « less
  2. Abstract

    Subterranean estuaries (STEs) form at the land‐sea boundary where groundwater and seawater mix. These biogeochemically reactive zones influence groundwater‐borne nutrient concentrations and speciation prior to export via submarine groundwater discharge (SGD). We examined a STE located along the York River Estuary (YRE) to determine if SGD delivers dissolved inorganic nitrogen (DIN) and phosphorus (DIP) to the overlying water. We assessed variations in STE geochemical profiles with depth across locations, times, and tidal stages, estimated N removal along the STE flow path, measured hydraulic gradients to estimate SGD, and calculated potential nutrient fluxes. Salinity, dissolved oxygen (DO), DIN, and DIP varied significantly with depth and season (p < 0.05), but not location or tidal stage. Ammonium dominated the DIN pool deep in the STE. Moving toward the sediment surface, ammonium concentrations decreased as nitrate and DO concentrations increased, suggesting nitrification. Potential sediment N removal rates mediated by denitrification were <8 mmoles N m−2 d−1. The total groundwater discharge rate was 38 ± 11 L m−2 d−1; discharge followed tidal and seasonal patterns. Net SGD nutrient fluxes were 0.065–3.2 and 0.019–0.093 mmoles m−2 d−1for DIN and DIP, respectively. However, microbial N removal in the STE may attenuate 0.58% to >100% of groundwater DIN. SGD fluxes were on the same order of magnitude as diffusive benthic fluxes but accounted for <10% of the nutrients delivered by fluvial advection in the YRE. Our results indicate the importance of STE biogeochemical transformations to SGD flux estimations and their role in coastal eutrophication and nutrient dynamics.

     
    more » « less
  3. Abstract

    Dissolved inorganic nutrient concentrations in the surface waters (0 to 5 m) of the Northern Gulf of Mexico (NGoM) were analyzed from 1985 to 2019 (> 10,000 observations) to determine spatiotemporal trends and their connection to nutrients supplied from the Mississippi/Atchafalaya River (MAR). In the NGoM, annual mean dissolved inorganic P (DIP) concentrations increased significantly over time, while dissolved inorganic N (DIN) concentrations showed no temporal trend. With greater salinity, mean DIN:DIP decreased from above the Redfield ratio of 16 to below it, reflecting DIN losses and the more conservative behavior of DIP with salinity. Over the same time period, annual mean P (total dissolved P, DIP, dissolved organic P) loading from the MAR to the NGoM significantly increased, annual mean DIN and total dissolved N loading showed no temporal trend, and dissolved organic N loading significantly decreased. Though DIP increased in the MAR, MAR DIP alone was insufficient to explain the surface distribution of DIP with salinity. Therefore, increases in surface DIP in the NGoM are not simply a reflection of increasing MAR DIP, pointing to temporal changes in other DIP sources. The increase in NGoM DIP suggests greater N limitation for phytoplankton, with implications for N fixation and nutrient management.

     
    more » « less
  4. Abstract

    Water temperature is a critical ecological indicator; however, few studies have statistically modeled century‐scale trends in riverine or estuarine water temperature, or their cause. Here, we recover, digitize, and analyze archival temperature measurements from the 1850s onward to investigate how and why water temperatures in the lower Columbia River are changing. To infill data gaps and explore changes, we develop regression models of daily historical Columbia River water temperature using time‐lagged river flow and air temperature as the independent variables. Models were developed for three time periods (mid‐19th, mid‐20th, and early 21st century), using archival and modern measurements (1854–1876; 1938–present). Daily and monthly averaged root‐mean‐square errors overall are 0.89°C and 0.77°C, respectively for the 1938–2018 period. Results suggest that annual averaged water temperature increased by 2.2°C ± 0.2°C since the 1850s, a rate of 1.3°C ± 0.1°C/century. Increased water temperatures are seasonally dependent. An increase of approximately 2.0°C ± 0.2°C/century occurs in the July–Dec time‐frame, while springtime trends are statistically insignificant. Rising temperatures change the probability of exceeding ecologically important thresholds; since the 1850s, the number of days with water temperatures over 20°C increased from ~5 to 60 per year, while the number below 2°C decreased from ~10 to 0 days/per year. Overall, the modern system is warmer, but exhibits less temperature variability. The reservoir system reduces sensitivity to short‐term atmospheric forcing. Statistical experiments within our modeling framework suggest that increased water temperature is driven by warming air temperatures (~29%), altered river flow (~14%), and water resources management (~57%).

     
    more » « less
  5. Abstract River deltas all over the world are sinking beneath sea-level rise, causing significant threats to natural and social systems. This is due to the combined effects of anthropogenic changes to sediment supply and river flow, subsidence, and sea-level rise, posing an immediate threat to the 500–1,000 million residents, many in megacities that live on deltaic coasts. The Mississippi River Deltaic Plain (MRDP) provides examples for many of the functions and feedbacks, regarding how human river management has impacted source-sink processes in coastal deltaic basins, resulting in human settlements more at risk to coastal storms. The survival of human settlement on the MRDP is arguably coupled to a shifting mass balance between a deltaic landscape occupied by either land built by the Mississippi River or water occupied by the Gulf of Mexico. We developed an approach to compare 50 % L:W isopleths (L:W is ratio of land to water) across the Atchafalaya and Terrebonne Basins to test landscape behavior over the last six decades to measure delta instability in coastal deltaic basins as a function of reduced sediment supply from river flooding. The Atchafalaya Basin, with continued sediment delivery, compared to Terrebonne Basin, with reduced river inputs, allow us to test assumptions of how coastal deltaic basins respond to river management over the last 75 years by analyzing landward migration rate of 50 % L:W isopleths between 1932 and 2010. The average landward migration for Terrebonne Basin was nearly 17,000 m (17 km) compared to only 22 m in Atchafalaya Basin over the last 78 years (p\0.001), resulting in migration rates of 218 m/year (0.22 km/year) and\0.5 m/year, respectively. In addition, freshwater vegetation expanded in Atchafalaya Basin since 1949 compared to migration of intermediate and brackish marshes landward in the Terrebonne Basin. Changes in salt marsh vegetation patterns were very distinct in these two basins with gain of 25 % in the Terrebonne Basin compared to 90 % decrease in the Atchafalaya Basin since 1949. These shifts in vegetation types as L:W ratio decreases with reduced sediment input and increase in salinity also coincide with an increase in wind fetch in Terrebonne Bay. In the upper Terrebonne Bay, where the largest landward migration of the 50 % L:W ratio isopleth occurred, we estimate that the wave power has increased by 50–100 % from 1932 to 2010, as the bathymetric and topographic conditions changed, and increase in maximum storm-surge height also increased owing to the landward migration of the L:W ratio isopleth. We argue that this balance of land relative to water in this delta provides a much clearer understanding of increased flood risk from tropical cyclones rather than just estimates of areal land loss. We describe how coastal deltaic basins of the MRDP can be used as experimental landscapes to provide insights into how varying degrees of sediment delivery to coastal deltaic floodplains change flooding risks of a sinking delta using landward migrations of 50 % L:W isopleths. The nonlinear response of migrating L:W isopleths as wind fetch increases is a critical feedback effect that should influence human river-management decisions in deltaic coast. Changes in land area alone do not capture how corresponding landscape degradation and increased water area can lead to exponential increase in flood risk to human populations in low-lying coastal regions. Reduced land formation in coastal deltaic basins (measured by changes in the land:water ratio) can contribute significantly to increasing flood risks by removing the negative feedback of wetlands on wave and storm-surge that occur during extreme weather events. Increased flood risks will promote population migration as human risks associated with living in a deltaic landscape increase, as land is submerged and coastal inundation threats rise. These system linkages in dynamic deltaic coasts define a balance of river management and human settlement dependent on a certain level of land area within coastal deltaic basins (L). 
    more » « less