ACS Publications. Most Trusted. Most Cited. Most Read
My Activity
Recently Viewed
You have not visited any articles yet, Please visit some articles to see contents here.
CONTENT TYPES

Figure 1Loading Img

Retrospective Analysis of Midsummer Hypoxic Area and Volume in the Northern Gulf of Mexico, 1985–2011

View Author Information
School of Natural Resources & Environment, University of Michigan, Ann Arbor, Michigan 48109-1041, United States
Department of Civil & Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109-2125, United States ,
§ Louisiana Universities Marine Consortium, 8124 Highway 56, Chauvin, Louisiana 70344, United States
Department of Oceanography and Coastal Sciences, Louisiana State University, Baton Rouge, Louisiana 70803, United States
Department of Global Ecology, Carnegie Institution for Science, Stanford, California 94305, United States
*Phone: 734-647-2464; fax: 734-647-6730; e-mail: [email protected]
Cite this: Environ. Sci. Technol. 2013, 47, 17, 9808–9815
Publication Date (Web):July 29, 2013
https://doi.org/10.1021/es400983g
Copyright © 2013 American Chemical Society
ACS AuthorChoiceACS AuthorChoice
Article Views
1687
Altmetric
-
Citations
LEARN ABOUT THESE METRICS
PDF (2 MB)
Supporting Info (1)»

Abstract

Robust estimates of hypoxic extent (both area and volume) are important for assessing the impacts of low dissolved oxygen on aquatic ecosystems at large spatial scales. Such estimates are also important for calibrating models linking hypoxia to causal factors, such as nutrient loading and stratification, and for informing management decisions. In this study, we develop a rigorous geostatistical modeling framework to estimate the hypoxic extent in the northern Gulf of Mexico from data collected during midsummer, quasi-synoptic monitoring cruises (1985–2011). Instead of a traditional interpolation-based approach, we use a simulation-based approach that yields more robust extent estimates and quantified uncertainty. The modeling framework also makes use of covariate information (i.e., trend variables such as depth and spatial position), to reduce estimation uncertainty. Furthermore, adjustments are made to account for observational bias resulting from the use of different sampling instruments in different years. Our results suggest an increasing trend in hypoxic layer thickness (p = 0.05) from 1985 to 2011, but less than significant increases in volume (p = 0.12) and area (p = 0.42). The uncertainties in the extent estimates vary with sampling network coverage and instrument type, and generally decrease over the study period.

1 Introduction

ARTICLE SECTIONS
Jump To

A large hypoxic zone has formed nearly every year in the northern Gulf of Mexico over at least the last few decades. (1-3) It is the largest human-caused coastal hypoxic zone in the western hemisphere and one of the largest worldwide. (1, 4) Because of its negative ecological impacts and potential fisheries impacts, (5) the hypoxic zone has received attention from various stakeholders and policy makers who have developed plans to reduce its average size. (6, 7)
The hypoxic zone is generally largest in summer due to the seasonal cycles of phytoplankton production and water column stratification, which are the two primary factors leading to dissolved oxygen (DO) depletion. (8-10) Both factors are related to the outflow from the Mississippi River Basin, which typically peaks in March–May. Nutrient loads from the Basin stimulate phytoplankton production in offshore waters, and much of the resulting organic matter eventually settles along the Louisiana–Texas shelf. Aerobic bacteria decompose the organic matter, consuming large quantities of DO. At the same time, freshwater outflows from the Basin promote stratification through a fresher surface layer, which warms in spring and summer, and overlays a colder, saltier water layer, inhibiting reoxygenation of bottom waters. (1)
The hypoxic extent, operationally defined as the region (area or volume) where DO concentrations are below 2 mg L–1, is often assessed using data from “shelfwide” sampling cruises. The cruises are considered “quasi-synoptic” because changes in weather conditions (e.g., tropical storms) preceding or during the cruises can disrupt typical seasonal patterns in hypoxia. Cruises to document the occurrence of hypoxia, as well as related physical and biological parameters, have been performed by the Louisiana Universities Marine Consortium (LUMCON) beginning in 1985. (1, 2) The estimated bottom-water hypoxic area (2) has been determined by interpolating between sampling locations and hand-contouring parallel to isobaths over a calibrated (planimeter) grid. These estimates have been used in multiple modeling studies linking hypoxia to nutrient loads from the Mississippi River basin and other environmental factors (11-14) and for setting policy goals to reduce the severity of hypoxia. (6, 7) However, these estimates are generally conservative in that the most inshore, most offshore, and most western extents of hypoxia are not always captured because of logistical constraints. (15, 16) Also, these estimates do not quantify the uncertainty inherent in making estimates from limited observations. (17, 18) Quantified uncertainties are useful for assessing the adequacy of existing sampling programs and for improving models that link hypoxic extent to environmental causes and effects.
Hypoxic area estimates only partially characterize the hypoxic extent. The thickness and volume of hypoxia are also important, as they relate to how hypoxia affects pelagic organisms. (19, 20) In addition, hypoxic volume, rather than area, should be more closely related to the total oxygen deficit of the system, potentially making it a more useful metric for biogeochemical modeling studies. (4) While hypoxic volume measurements are available for other systems, such as Chesapeake Bay, (21, 22) estimates for the Gulf have only been available for a subset of years, using an unpublished methodology. (4)
The primary goal of this study is to improve our knowledge of the Gulf’s hypoxic extent over the 27-year study period (1985–2011) by systematically estimating midsummer hypoxic area and volume. These are the first Gulf hypoxic extent estimates to include “instrument bias adjustments” that account for the use of different sampling instruments, capable of being lowered to within different proximities of the sea floor, in different years. (23) Our approach also uses trends between DO, hypoxic fraction, depth, and spatial position to account for consistent, large-scale spatial patterns in DO, thereby improving the explanatory power of the model. (18) Our extent estimates are developed using a Monte Carlo-type simulation approach, rather than traditional interpolation, allowing for uncertainty quantification (i.e., confidence intervals). (17, 18) These confidence intervals reflect the spatial stochasticity of the system, uncertainties in trends among variables, and uncertainties in instrument bias adjustments. Hypoxic volume is estimated using a novel two-step approach where DO is simulated first, and hypoxic fraction (i.e., the fraction of the water column that is hypoxic) is simulated second using DO as a trend variable.

2 Methods

ARTICLE SECTIONS
Jump To

2.1 Data and Study Boundaries

We use DO data from LUMCON midsummer sampling cruises conducted between 1985 and 2011. Data for 1998–2008 are retrieved from the National Ocean Data Center, while data for other years are obtained directly from LUMCON. (24) Sampling locations are geo-referenced using the Universal Transverse Mercator (UTM) Zone 15 projection, and bathymetry is determined from a 3-arc-second digital elevation model obtained from the National Oceanic and Atmospheric Administration (NOAA). (25)
During cruises, DO is sampled using one or two types of instruments (Figure 1): a rosette-mounted DO probe and a hand-held DO probe, with the latter capable of being lowered closer to the sea floor. (9, 23) (Hereafter, we refer to these instruments as the “rosette sampler” and “handheld sampler”.) While instrument technology has changed over time, all instruments were calibrated against Winkler titrations, and postcruise corrections were made as necessary. For sampling events where both instruments were used, data are combined into synthesized profiles. (9) (Here, a “sampling event” refers to the depth profile of data collected at a specific latitude/longitude and time.)

Figure 1

Figure 1. Number of locations sampled during the annual midsummer shelfwide cruises using hand-held and rosette instruments.

We extract bottom water dissolved oxygen (BWDO) and minimum dissolved oxygen (MinDO) concentrations from the DO profiles. BWDO and MinDO concentrations greater than 7 mg L–1 (less than 1% of total samples), are treated as 7 mg L–1, because higher concentrations are outliers representing supersaturation conditions that are not of interest in this study. For sampling events where hypoxia was present, we also extract the fraction of the water column that is within the hypoxic bottom layer (hereafter, bottom water hypoxic fraction, or BWHF) and, where other hypoxic layers are present, the total hypoxic fraction (THF) of all hypoxic layers. Bottom and upper layers of hypoxia are observed at 49% and 12% of sampling events, respectively. Here, we focus on the models for BWDO and BWHF because bottom layer results are most comparable to previous studies. (2, 16) The models for MinDO and THF are formulated in the same way and yield comparable results, but for brevity are discussed only in Supporting Information (SI) Section S7.
In this study, adjustments are made to address biases that arise from the use of different sampling instruments. For cruises when only the rosette sampler was used (Figure 1), one would expect the hypoxic extent to be underestimated because the rosette sampler does not reach as close to the sea floor as the hand-held sampler. We quantify this bias using data from sampling events where both instruments were used together. For these cases, observed BWDO and BWHF are calculated for both the synthesized profile and rosette-only profile. Probabilistic relationships are then developed between the synthesized results and the rosette-only results (Section S1 of the SI). When performing simulations (described below), we adjust the rosette-only observations by sampling from these relationships. A larger bias adjustment is required for the first 38 sampling events in 1991 because the ship’s fathometer was not properly calibrated, causing the rosette sampler to be lowered 1.5 m less than it would have been otherwise (N.N. Rabalais, cruise records). For these observations, the probabilistic relationships between synthesized and rosette-only results are again determined based on sampling events where both instruments were functioning properly (as described above) but with the bottom 1.5 m of the rosette-only profiles removed.
Geostatistical modeling is performed on a 5 × 5 km2 grid of estimation points (Figure 2), covering 342.5–837.5 km UTM easting and 3122.5–3292.5 km UTM northing, reflecting the general extent of sampling. The estimation grid is limited to depths between 3 and 80 m, which are typical of the Louisiana–Texas shelf region where hypoxia occurs.

Figure 2

Figure 2. Study area bathymetry, sampling, and estimation locations.

2.2 Model Formulation

Geostatistical methods provide an effective means to model data that exhibit spatial correlation. (17, 26) The efficacy of these methods has been demonstrated in previous environmental analyses of rainfall, snow depth, soil phosphorus, and water quality indicators. (18, 21, 27-29) Many of these studies have focused on mapping through spatial interpolation, and comparison studies have demonstrated that geostatistical interpolation outperforms simpler interpolation methods, such as inverse distance weighting. (21, 28)
In a geostatistical model, the response variable, z, is represented as a combination of deterministic and stochastic components: (26)(1)
This formulation is similar to linear regression, but includes an additional term, η, representing spatially correlated stochasticity, in addition to the more commonly modeled uncorrelated stochasticity, ε. As in linear regression, is the portion of z that can be expressed as a deterministic function of categorical and/or trend variables (X) and their corresponding regression coefficients (β). If the X term includes only cruise-specific categorical variables (vectors of zeros and ones that bin the data by cruise), then the model can essentially be used to perform ordinary kriging (OK). If X also includes trend variables, then it can be used to perform universal kriging (UK). (17)
In this study, the stochastic components (η and ε) are found to be well represented by the commonly used exponential covariance function with a nugget effect, (26) defined as follows:(2)where Q is the covariance between two observations (zi and zj) separated by distance (hi,j); σε2 and ση2 are parameters representing the variances of the stochasticity that is spatially correlated and is not spatially correlated, respectively; and r is a range parameter (3r is approximately the distance at which observations are no longer spatially correlated). We allow for anisotropy in the covariance model by scaling hi,j using parameter α, which represents the ratio of east–west to north–south correlation ranges. Covariance function parameters are estimated using restricted maximum likelihood (30, 31) and deterministic component parameters (β) are estimated using generalized least-squares, as outlined in a previous study by Obenour et al. (9) that used geostatistical methods to explore biophysical drivers of DO depletion.
Geostatistical models are developed for both BWDO and BWHF to determine bottom layer hypoxic area and volume, respectively. The model for BWDO uses all observations, while the model for BWHF uses only observations from locations where BWDO is below the hypoxic threshold. Both models use the UK formulation to take advantage of potential relationships between the response variables and available trend variables. (18)
The deterministic components of the models include two types of trends. First, “constant trends” represent relationships between response and trend variables using regression coefficients that are the same for all cruises. For the BWDO model, potential constant trends include linear and quadratic trends with depth (Depth and Depth2), easting (Easting and Easting2), and northing (Northing and Northing2). On the basis of an examination of model residuals, it was noted that the overall trend between depth and BWDO does not continue for depths greater than 40 m, and so depths of greater than 40 m are treated as 40 m in the final BWDO model development. This approach is justified because model residuals are evenly distributed across all depths when depths greater than 40 m are modified as described here (see Figure S3 of the SI). For BWHF, potential constant trends include linear and quadratic trends with easting, northing, and BWDO.
Second, “cruise-specific trends” (with cruise-specific regression coefficients for Easting) represent relationships that are specific to individual cruises. Allowing these trends was motivated by previous studies (16, 23) that indicate the east–west distribution of hypoxia is influenced by alongshore current velocity, which can vary interannually in response to prevailing winds. When selected, these cruise-specific trends modify the constant trend.
To prevent overparameterization of the model, we use only trend variables selected through a geostatistical adaptation of the Bayesian information criterion (BIC), (9, 32, 33) wherein models with different subsets of trend variables are compared based on their BIC score, and the model with the lowest BIC score is optimal in terms of its explanatory power relative to its complexity. Because of the large number of variables considered in this study, an initial optimal model is first selected among those formed from constant trend variables only. Cruise-specific trend variables are then individually tested relative to that model, and only cruise-specific trends that improve the BIC score are included in the final model.
The deterministic portion of the model also includes categorical variables that bin data by cruise, reflecting the fact that the mean BWDO varies from cruise to cruise. In UK, the categorical variables essentially allow for different “intercepts” (as in linear regression), that shift the trend up or down to best fit the observations from each cruise. In OK, these categorical variables simply allow for a different mean for each cruise.
The covariance model, selected deterministic variables, and categorical variables are used to determine a set of geostatistical weights, Λ, that are applied to observations when performing geostatistical interpolation and simulation. For each cruise, y, a unique set of weights, Λy, are determined by solving a system of linear equations:(3)where Qoo is an n × n covariance matrix for the n observation locations (all cruises), with elements determined from eq 2. Because we did not assume correlation among stochasticity from different cruises, intercruise covariances are assigned a value of zero. Similarly, Qoe,y is an n × m covariance matrix of n observation locations and m estimation locations, and the rows of Qoe,y that correspond to observations from cruises other than cruise y are assigned a value of zero. The matrix Xo is n × p and includes the p deterministic variables (trend and categorical) for the observation locations, and the matrix Xe (m × p) includes the same variables for the estimation locations. The trend variables are normalized to a mean of zero and variance of one. Note that terms with a “y” subscript are cruise-specific.
In eq 3, Λy is an n × m matrix of cruise-specific weights that reflect both the spatial correlation structure and the deterministic trends. Also, Gy is a p × m matrix of Lagrange multipliers that can be used with Λy to determine location-specific estimation uncertainties. Using the weights, along with the observations zo (an n × 1 vector), one can develop estimates of the response variable across the estimation grid:(4)where ze,y0 is an m × 1 vector of interpolated BWDO or BWHF values for cruise y.
Although we use a UK model formulation, hypoxic area and volume are not determined by kriging (i.e., spatial interpolation), but instead by developing “spatially consistent Monte Carlo simulations” (17) which are often referred to as conditional realizations (CRs). While both kriging and CRs provide equivalent information for individual estimation locations, CRs are necessary to estimate spatially aggregated quantities (e.g., area and volume) probabilistically. A CR is performed by first creating an “unconditional realization” and then “conditioning” it to the observed data and deterministic trend. (17) Unconditional realizations (eq 5) include simulated values at the estimation locations (ze,yu) as well as at the observation locations (zo,yu):(5)Note that the vector zo,yu includes simulated values corresponding to the observations from all cruises, but only observations from cruise y have their stochasticity correlated with that of the estimation locations. Here, Qee is the m × m covariance matrix between estimation locations, and u is an (m + n) × 1 vector of random independent samples from the standard normal distribution. The operator C( ) returns the triangular matrix resulting from Cholesky decomposition of the subject matrix.
The unconditional realizations are then conditioned to the observed data and deterministic trends through:(6)Here, zo is the n × 1 vector of the observed values, and ze,yc is the resulting cruise-specific CR, an m × 1 vector of values corresponding to the estimation locations.
The CRs are performed in two steps. First, BWDO concentration is simulated across the entire estimation grid. Then, BWHF is simulated over those locations where the simulated BWDO is below the hypoxic threshold. At hypoxic estimation locations, the simulated BWDO values are used as a trend variable in the BWHF model. Simulated values are limited to within realistic ranges (determined to be 0–7 mg L–1 for BWDO and 0.001–0.8 for BWHF, based on a review of the observed data). Without these constraints, the estimated hypoxic volume would be 6% smaller, on average, as a result of (unrealistic) negative hypoxic thicknesses.
The two-step CR process is repeated 1000 times for each cruise, resulting in 1000 realizations of both BWDO and BWHF. For each realization of BWDO, the hypoxic area is calculated as the number of estimation locations simulated to be hypoxic multiplied by the grid cell area (25 km2); and for each realization of BWHF, the hypoxic volume is calculated as the vector of simulated BWHF values multiplied by the corresponding vector of water column depths, all multiplied by the grid cell area. From this ensemble of results, we determine the mean and 95% confidence intervals for the hypoxic area and volume for each cruise.

3 Results

ARTICLE SECTIONS
Jump To

The BWDO and BWHF models include several parameters that characterize the deterministic and stochastic model components (eq 1). Regression coefficients for the BIC-selected trend variables, explaining a portion of the spatial variability in observed BWDO and BWHF, are provided in Table 1. The standard errors of these coefficients are low (i.e., p < 0.05), suggesting these trends are statistically significant. The coefficients for the cruise-specific categorical variables, accounting for year-to-year variability in the responses, are included in Section S3 of the SI. Overall, the deterministic model components explain 28% and 32% of the total (spatial plus interannual) variability in BWDO and BWHF, respectively, while the stochastic components of the models explain the remainder of the spatial variability. If trend variables are omitted (i.e., the OK formulation), then the deterministic components only account for interannual variability, and only 12% and 11% of the total variability in BWDO and BWHF is explained, respectively. A check of the linearity assumption implicit in these deterministic relationships is included in Section S2 of the SI.
Table 1. Regression Coefficients (β̂) with Standard Errors (σβ̂) for Normalized, BIC-Selected Trend Variables in BWDO and BWHF Modelsa,

b

variableBWDO (mg L–1)BWHF
β̂σβ̂β̂σβ̂
Easting–0.620.090.0180.007
Easting20.250.07–0.0200.006
Northing–0.360.09n.s.
Depth–2.310.18n.a.
Depth22.450.17n.a.
BWDOn.a.–0.0650.005
c.s.E 1998–1.350.45n.s.
a

Parameters optimized by generalized least squares.

b

c.s.E = cruise specific trend for Easting, n.s.=not selected, n.a.=not available.

The BWDO model includes trends with easting, northing, and depth. The trend between BWDO and easting is quadratic with a minimum at 794 km easting (UTM coordinate, Figure 2), between the Mississippi and Atchafalaya river outfalls, which is reasonable given that these rivers provide the freshwater flows and nutrients that are important for hypoxia formation. In 1998, the only year for which a cruise-specific east–west trend with BWDO is selected, the minimum is shifted to 1111 km easting, outside of the study area, indicating that BWDO concentrations decrease monotonically with easting within the study area. The unique spatial distribution in 1998 has been noted previously, and is generally attributed to unusually persistent eastward currents. (16) A similar BWDO pattern has been noted for 2009, (23) but it was not sufficiently strong to result in a cruise-specific trend in this analysis. The trend between BWDO and northing is linear, and suggests that BWDO concentrations are higher to the south. The trend between BWDO and bathymetry (depth) is quadratic with a minimum at 22 m. The overall deterministic trends in BWDO (Figure S7 of the SI) compare well with spatial patterns in hypoxic frequency, as mapped by Walker and Rabalais. (34)
Because BWHF is modeled as a function of BWDO, it effectively inherits the trends from the BWDO model. The trend with BWDO is linear and suggests that BWHF is larger where BWDO concentrations are lower. The model for BWHF also includes a quadratic trend with easting, suggesting a maximum at 701 km, about 30 km east of the Atchafalaya outfall location (although the overall trend in BWHF is also affected by the trends in BWDO). No cruise-specific east–west trends were selected for the BWHF model.
The models also include covariance parameters that characterize their stochastic components. For BWDO, σε2 and ση2, commonly referred to as the nugget and partial sill, were 0.39 and 2.33 (mg2 L–2), respectively. Since σε2 is smaller than ση2, the majority of the stochasticity in the model is spatially correlated. The approximate range of spatial correlation (3r, per eq 2) is 94 km in the east–west direction and 53 km in the north–south direction. The greater correlation distance in the east–west direction was expected due to the dominant east–west current pattern. (35, 36) For BWHF, the spatial correlation of the stochasticity is somewhat weaker with σε2 and ση2 similar in magnitude, having values of 0.009 and 0.011 (unit-less), respectively. Also, the correlation range was 63 km in all directions (anisotropy was negligible). Overall, spatially correlated stochasticity accounts for the greatest portion of the variability in both BWDO and BWHF. As discussed in Obenour et al., (9) the spatial correlation of the stochasticity is consistent with the effects of varying coastal current patterns, influencing the distribution of hypoxia over the spatial scales described above. In general, the correlation ranges are considerably longer than the typical distances between sampling locations (Figure 2), especially on the eastern shelf, suggesting that the sampling network is adequate for resolving the spatially correlated stochasticity of the system.
Using CR, we determined the mean and 95% confidence intervals for hypoxic area and volume of each cruise (Figure 3). The largest estimated hypoxic area was for 1996, but it was not significantly different from 1991, 1993, 1995, 1997, 1999, 2001, 2002, 2007, or 2008 (p > 0.05), given the uncertainty in the area estimates. The largest hypoxic volume estimate was for 2008, but it was not significantly different from 1993, 1996, 1999, 2004, or 2007 (p > 0.05). Both hypoxic area and volume were lowest in 1988, a drought year. (16)

Figure 3

Figure 3. Bottom layer hypoxic extent estimates with 95% confidence intervals by year; estimates prior to making adjustments for instrument bias as triangles; previous LUMCON area estimates as open squares; revised LUMCON area estimates as solid squares.

Two sets of hypoxic area estimates, determined by LUMCON, are also included in Figure 3 for comparison. The first set are the original LUMCON estimates, determined by hand contouring, as described above. (2, 19) The second set are revised estimates developed as part of this study; also by hand contouring, but using the updated BWDO values for this study (which are superior because they reflect postcruise DO calibrations). In most cases, the original LUMCON estimates (Figure 3, open squares) overlap the revised estimates (solid squares), and the majority of the LUMCON estimates fall within the 95% confidence intervals determined by this study. However, for the first third of the study period (1985–1993), the new estimates are consistently higher than the LUMCON estimates.
LUMCON did not make an estimate for 1989, because only a portion of the shelf was sampled that year. Our method does allow for extent estimates for 1989, and the wide confidence intervals (Figure 3) reflect the relatively large uncertainties. Even though data for 1989 were limited, the overall spatial trends (based on the data from all years) help to constrain the variability in the CRs to within realistic ranges, across the study area. We note that the 1989 estimate presented here is consistent with previous estimates developed from nutrient loading models. (37, 38)
There was considerable interannual variability in the spatial distribution of hypoxia and in the thickness of the hypoxic bottom layer (Figure 4). Years with similar hypoxic areas may have very different average hypoxic thicknesses and thus volumes. For example, 2002 and 2008 have similar hypoxic areas, but they are very different in terms of hypoxic thickness, such that 2008 has approximately twice the hypoxic volume of 2002. The results indicate that the average bottom layer hypoxic thickness for the 27-year study period was 3.9 m, with the thickest average layers of approximately 6.2 and 6.3 m occurring in 2008 and 2009, respectively. Over the 27-year study period, hypoxic volume is correlated with area (r2 = 0.77). Interestingly, hypoxic area and thickness are also somewhat positively correlated over most years (r2 = 0.37, not including 1998 and 2009, which were subject to unusually strong eastward currents (16, 23)), such that volume increases exponentially relative to area. Tabulated results and additional maps showing estimated BWDO, estimated hypoxic thickness, probability of hypoxia, and example CRs for the entire study period are included in Sections S4 and S5 of the SI.

Figure 4

Figure 4. Example maps of estimated bottom layer hypoxic thickness (median values from CRs), 2001–2008; observation locations shown as white dots.

4 Discussion

ARTICLE SECTIONS
Jump To

The geostatistical modeling results can be used to assess temporal trends in hypoxic zone size for the 27-year study period. Hypoxic volume increased by an average of 2.3% per year as a linear trend with time (1.4 km3 yr–1, percent increase determined by dividing by the 27-year mean hypoxic extent), but this trend was not significant (p = 0.12). Hypoxic area increased at a lesser rate of 0.9% per year (140 km2 yr–1) and was also not significant (p = 0.42). However, the relatively large increase in volume relative to area reflects a significant increasing trend in hypoxic layer thickness (1.8% per year, 0.069 m yr–1, p = 0.05). Note that trend significance is affected by the uncertainty in the geostatistical estimates; trend significance is determined using a Monte Carlo approach where temporal trend coefficients are developed for each of the 1000 sets of CRs, and an overall probability distribution is developed by sampling from the uncertainty in these trend coefficients (note p-values are based on a two-sided test). Without accounting for uncertainty, the volume, area, and thickness trends would have p-values of 0.09, 0.38, and 0.01, respectively.
The new hypoxic area estimates can be compared to the previous hypoxic area estimates developed by LUMCON. (2, 16) For the 27-year study period, the previous area estimates increased at a highly significant rate of 2.6% per year (360 km2 yr–1, p = 0.01), substantially greater than the 0.9% rate for the new estimates. This dissimilarity is due primarily to differences in estimates for the earlier years of the study period (Figure 3). For 1985–1987, the new estimates are consistently higher than the previous estimates primarily because the geostatistical methodology accounts for the possibility of hypoxia occurring outside the envelopes of the cruises, which were relatively small in these years. For 1990, 1991, and 1993, the new estimates are also higher because of the instrument bias adjustments developed in this study. For 1985–1993 (except 1989, for which there is no previous estimate), the mean geostatistically determined hypoxic area is 39% (3650 km2) greater than that of the previous estimates.
For 1994–2011, the new area estimates are in general agreement with the previous area estimates. Over this time period, the two sets of estimates are highly correlated (r= 0.88), the means of the two data sets are of negligible difference, and only three of the previous estimates (1996, 2003, and 2010) fall outside the 95% confidence intervals of the new estimates (Figure 3). This suggests that when the biasing issues noted above are avoided, there is approximate agreement between the geostatistical and hand-contouring estimates.
The new extent estimates have implications for our understanding of how hypoxia is changing over long temporal scales. Multiple studies using nutrient loads to predict the previous hypoxic area estimates have suggested the Gulf is becoming increasingly susceptible to hypoxia, based on increasing hypoxic area relative to nutrient loading. (11, 12, 39, 40) (Nutrient loading increased greatly in the 1970s but remained relatively stable throughout the study period.) (37, 41) Increases in hypoxic susceptibility have also been suggested in other studies of the Gulf and other coastal systems. (42-44) The new hypoxic area estimates, however, exhibit relatively little increase over time, especially when compared to the previous estimates, potentially suggesting less system change during the study period than previously thought. But, the new hypoxic thickness and volume estimates do increase to a greater degree, potentially suggesting a more vertically oriented increase in hypoxic extent. Future studies could focus on recalibrating existing nutrient-loading models to these new hypoxic extent estimates to develop a refined understanding of how Gulf hypoxia may be changing over time.
Uncertainties in the new hypoxic extent estimates, represented by the 95% confidence intervals (Figure 3), reflect the limited spatial scope and resolution of shelfwide cruise sampling. Uncertainties are generally greatest in the earlier years when cruises were smaller and did not always use instruments that reached the sea floor. From 1985 to 1993, the mean relative standard error for hypoxic area was 23%, but it decreased to 11% for 1994–2011. Uncertainties also appear to be larger in years with relatively severe hypoxia (i.e., lower average BWDO and larger hypoxic area), likely because more of the estimation grid is subject to the possibility of hypoxia (in years with higher average BWDO, most of the estimation grid is determined to be well above the hypoxic threshold, such that simulated values rarely fall below the hypoxic threshold). In the future, the modeling framework presented here could be used to evaluate different sampling designs based on how well they constrain estimate uncertainty. The modeling approach could also be used to compare hypoxic extent results and associated uncertainties using hypoxic threshold choices other than 2 mg L–1.
Regardless of the precision of the estimates, it is important to remember that they represent conditions at only certain points in time (i.e., during the shelfwide cruises), and that hypoxic extent can vary substantially throughout the summer due to changes in organic matter production, wind-driven mixing events, and fluctuating current patterns. (10, 45, 46) As a result, these estimates do not necessarily reflect hypoxic conditions over the entire summer (although back-to-back cruises have demonstrated fairly consistent hypoxic areas under stable weather conditions). (15) More detailed monitoring and biophysical modeling are needed to better understand the short- and long-term dynamics of hypoxia formation, and to mechanistically interpret the temporal variability in the extent estimates presented here.
The primary feature of the geostatistical approach, when compared to more traditional, interpolation-based approaches, is the use of simulations (i.e., CRs). (18) The most obvious benefit is the ability to quantify the uncertainty in extent estimates. A second benefit is that CRs provide more realistic extent estimates than can be derived from kriging (or largely equivalent methods such as Gauss-Markov smoothing (47)) alone. In this study, hypoxic area estimates derived from kriged maps (UK interpolation, Figure S14 of the SI) were substantially lower than both our CR estimates and the LUMCON estimates. This is because the kriged maps (Figure S9 of the SI) tend to characterize large portions of the estimation grid as slightly above the hypoxic threshold. However, due to the stochasticity of the system, these locations still have some probability of being hypoxic. The CR approach, which samples from the uncertainty in the system, accounts for this possibility, and thus results in larger hypoxic area estimates. This is consistent with Chiles and Delfiner, (17) who argue that CR, rather than kriging, is the more appropriate approach for determining spatially aggregated quantities. A third advantage of the CR approach is that it provides a framework for performing instrument bias adjustments probabilistically, such that adjustment uncertainty is propagated to the hypoxic extent estimates. Finally, CR was fundamental to probabilistically determining hypoxic volume (in addition to area). The two-step CR approach, developed here, can be compared to other volume estimation methods, such as multilayer kriging (21) and three-dimensional CR, (22) which have been applied in Chesapeake Bay. A benefit of our approach is that it allows uncertainty quantification of hypoxic volume within a relatively simple (two-dimensional) geostatistical framework.
This study also demonstrates the benefits of including trend variables within the geostatistical model (i.e., the UK formulation) because the deterministic trends help reduce model uncertainty and result in more realistic extent estimates (Section S6 of the SI). The advantages of UK have been demonstrated previously for Lake Erie and Chesapeake Bay hypoxia, (18, 21) but are perhaps even more salient for an open system such as the Gulf shelf. When trends in DO are not modeled, it has been necessary to limit the estimation grid around the bounds of the sampling cruise, (15, 48) such that the size of the cruise can potentially bias the inferred hypoxic area. By including trend variables that explain the large-scale spatial patterns in BWDO and BWHF, it is possible to develop realistic CRs of DO across the entire study area, so the same estimation grid can be used for all cruises.
Finally, the results confirm that the hypoxic area on the Louisiana-Texas shelf greatly exceeds the Hypoxia Task Force goal of 5000 km2 as a five-year running average. (6, 7) The most recent five-year period from our study, 2007–2011, has a mean hypoxic area of 16 600 km2 with a 95% confidence range of 15 100–18 000 km2. Clearly, additional management measures are required if the hypoxic extent is to be reduced to comply with the Task Force goal.

Supporting Information

ARTICLE SECTIONS
Jump To

A description of the instrument bias adjustments, a graphical test of linearity for the deterministic components of the models, the intercept parameters for the categorical variables, a set of results maps, tabulated bottom layer hypoxic extent results, and a summary of the models for MinDO and THF. This material is available free of charge via the Internet at http://pubs.acs.org.

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

ARTICLE SECTIONS
Jump To

  • Corresponding Author
    • Daniel R. Obenour - School of Natural Resources & Environment, University of Michigan, Ann Arbor, Michigan 48109-1041, United StatesDepartment of Civil & Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109-2125, United States , Email: [email protected]
  • Authors
    • Donald Scavia - School of Natural Resources & Environment, University of Michigan, Ann Arbor, Michigan 48109-1041, United StatesDepartment of Civil & Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109-2125, United States ,
    • Nancy N. Rabalais - Louisiana Universities Marine Consortium, 8124 Highway 56, Chauvin, Louisiana 70344, United States
    • R. Eugene Turner - Department of Oceanography and Coastal Sciences, Louisiana State University, Baton Rouge, Louisiana 70803, United States
    • Anna M. Michalak - Department of Global Ecology, Carnegie Institution for Science, Stanford, California 94305, United States
  • Notes
    The authors declare no competing financial interest.

Acknowledgment

ARTICLE SECTIONS
Jump To

We are grateful to the LUMCON personnel who collected and provided quality control for the data used in this study. The work was supported in part by the US EPA STAR Fellowship program, NOAA’s Center for Sponsored Coastal Ocean Research grants NA09NOS4780204, NA06NPS4780197, NA09NOS4780204, and NA16OP2670, and the University of Michigan’s Graham Sustainability Institute. This is NGOMEX Contribution 176.

References

ARTICLE SECTIONS
Jump To

This article references 48 other publications.

  1. 1
    Rabalais, N. N.; Turner, R. E.; Wiseman, W. J. Gulf of Mexico hypoxia, aka “The dead zone” Annu. Rev. Ecol. Syst. 2002, 33, 235 263
  2. 2
    LUMCON hypoxia in the northern Gulf of Mexico: Research; http://www.gulfhypoxia.net/Research/. (Accessed October, 2012.)
  3. 3
    Rabalais, N. N.; Turner, R. E.; Scavia, D. Beyond science into policy: Gulf of Mexico hypoxia and the Mississippi River Bioscience 2002, 52 (2) 129 142
  4. 4
    Rabalais, N. N.; Diaz, R. J.; Levin, L. A.; Turner, R. E.; Gilbert, D.; Zhang, J. Dynamics and distribution of natural and human-caused hypoxia Biogeosciences 2010, 7 (2) 585 619
  5. 5
    Rabalais, N. N.; Turner, R. E.Coastal Hypoxia: Consequences for Living Resources and Ecosystems. Coastal and Estuarine Studies 58; American Geophysical Union: Washington, D.C., 2001.
  6. 6
    Mississippi River/Gulf of Mexico Watershed Nutrient Task Force Action Plan for Reducing, Mitigating, And Controlling Hypoxia in the Northern Gulf of Mexico; Office of Wetlands, Oceans, and Watersheds, U.S. Environmental Protection Agency: Washington, DC, 2001, http://water.epa.gov/type/watersheds/named/msbasin/history.cfm.
  7. 7
    Mississippi River/Gulf of Mexico Watershed Nutrient Task Force Gulf Hypoxia Action Plan 2008 for Reducing Mitigating, And Controlling Hypoxia in the Northern Gulf of Mexico and Improving Water Quality in the Mississippi River Basin; U.S. Environmental Protection Agency, Office of Wetlands, Oceans, and Watersheds: Washington, D.C., 2008, http://water.epa.gov/type/watersheds/named/msbasin/actionplan.cfm.
  8. 8
    Hypoxia in the Northern Gulf of Mexico, An Update by the EPA Science Advisory Board; EPA Science Advisory Board: Washington, DC, 2007, EPA-SAB-08–003.
  9. 9
    Obenour, D. R.; Michalak, A. M.; Zhou, Y.; Scavia, D. Quantifying the impacts of stratification and nutrient loading on hypoxia in the northern Gulf of Mexico Environ. Sci. Technol. 2012, 46 (10) 5489 5496
  10. 10
    Wiseman, W. J.; Rabalais, N. N.; Turner, R. E.; Dinnel, S. P.; MacNaughton, A. Seasonal and interannual variability within the Louisiana coastal current: stratification and hypoxia J. Mar. Syst. 1997, 12 (1–4) 237 248
  11. 11
    Forrest, D. R.; Hetland, R. D.; DiMarco, S. F. Multivariable statistical regression models of the areal extent of hypoxia over the Texas–Louisiana continental shelf Environ. Res. Lett. 2011, 6 (4) 045002
  12. 12
    Greene, R. M.; Lehrter, J. C.; Hagy, J. D. Multiple regression models for hindcasting and forecasting midsummer hypoxia in the Gulf of Mexico Ecol. Appl. 2009, 19 (5) 1161 1175
  13. 13
    Scavia, D.; Rabalais, N. N.; Turner, R. E.; Justic, D.; Wiseman, W. J. Predicting the response of Gulf of Mexico hypoxia to variations in Mississippi River nitrogen load Limnol. Oceanogr. 2003, 48 (3) 951 956
  14. 14
    Turner, R. E.; Rabalais, N. N.; Justic, D. Predicting summer hypoxia in the northern Gulf of Mexico: Riverine N, P, and Si loading Mar. Pollut. Bull. 2006, 52 (2) 139 148
  15. 15
    Rabalais, N. N.; Turner, R.E.; Justic, D.; Dortch, Q.; Wiseman, W.Characterization of Hypoxia: Coastal Ocean Program Decision Analysis Series No. 15; National Oceanic and Atmospheric Administration (NOAA): Silver Springs, MD, 1999.
  16. 16
    Rabalais, N. N.; Turner, R. E.; Sen Gupta, B. K.; Boesch, D. F.; Chapman, P.; Murrell, M. C. Hypoxia in the northern Gulf of Mexico: Does the science support the plan to reduce, mitigate, and control hypoxia? Estuaries Coasts 2007, 30 (5) 753 772
  17. 17
    Chiles, J.-P.; Delfiner, P. Geostatistics: Modeling Spatial Uncertainty; Wiley: New York, 1999.
  18. 18
    Zhou, Y.; Obenour, D. R.; Scavia, D.; Johengen, T. H.; Michalak, A. M. Spatial and temporal trends in Lake Erie hypoxia, 1987–2007 Environ. Sci. Technol. 2013, 47 (2) 899 905
  19. 19
    Zhang, H. Y.; Ludsin, S. A.; Mason, D. M.; Adamack, A. T.; Brandt, S. B.; Zhang, X. S.; Kimmel, D. G.; Roman, M. R.; Boicourt, W. C. Hypoxia-driven changes in the behavior and spatial distribution of pelagic fish and mesozooplankton in the northern Gulf of Mexico J. Exp. Mar. Biol. Ecol. 2009, 381, S80 S91
  20. 20
    Kimmel, D. G.; Boicourt, W. C.; Pierson, J. J.; Roman, M. R.; Zhang, X. S. A comparison of the mesozooplankton response to hypoxia in Chesapeake Bay and the northern Gulf of Mexico using the biomass size spectrum J. Exp. Mar. Biol. Ecol. 2009, 381, S65 S73
  21. 21
    Murphy, R. R.; Curriero, F. C.; Ball, W. P. Comparison of spatial interpolation methods for water quality evaluation in the Chesapeake Bay J. Environ. Eng.–ASCE 2010, 136 (2) 160 171
  22. 22
    Zhou, Y.; Scavia, D.; Michalak, A. M. Nutrient loading and meteorological conditions explain interannual variability of hypoxia in the Chesapeake Bay, (submitted manuscript).
  23. 23
    Turner, R. E.; Rabalais, N. N.; Justic, D. Predicting summer hypoxia in the northern Gulf of Mexico: Redux Mar. Pollut. Bull. 2012, 64 (2) 319 324
  24. 24
    Rabalais, N. N. Louisiana Hypoxia Surveys, available from NOAA National Oceanographic Data Center (NODC); http://www.nodc.noaa.gov/General/getdata.html. (Accessed January, 2011.)
  25. 25
    NOAA National Geophysical Data Center (NGDC) Coastal ReliefModel; http://www.ngdc.noaa.gov/mgg/coastal/startcrm.htm. (Accessed June, 2010.)
  26. 26
    Zimmerman, D. L.; Stein, M. Classical Geostatistical Methods. In Handbook of Spatial Statistics; Gelfand, A. E.; Diggle, P. J.; Fuentes, M.; Guttorp, P., Eds.; CRC Press: Boca Raton, FL, 2010.
  27. 27
    Erickson, T. A.; Williams, M. W.; Winstral, A. Persistence of topographic controls on the spatial distribution of snow in rugged mountain terrain, Colorado, United States. Water Resour. Res. 2005, 41 (4).
  28. 28
    Goovaerts, P. Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall J. Hydrol. 2000, 228 (1–2) 113 129
  29. 29
    Qian, S. S. Estimating the area affected by phosphorus runoff in an Everglades wetland: A comparison of universal kriging and Bayesian kriging Environ. Ecol. Stat. 1997, 4 (1) 1 26
  30. 30
    Patterson, H. D.; Thompson, R. Recovery of inter-block information when block sizes are unequal Biometrika 1971, 58 (3) 545 554
  31. 31
    Zimmerman, D. L. Likelihood-Based Methods. In Handbook of Spatial Statistics; Gelfand, A. E.; Diggle, P. J.; Fuentes, M.; Guttorp, P., Eds.; CRC Press: Boca Raton, FL, 2010.
  32. 32
    Mueller, K. L.; Yadav, V.; Curtis, P. S.; Vogel, C.; Michalak, A. M. Attributing the variability of eddy-covariance CO2 flux measurements across temporal scales using geostatistical regression for a mixed northern hardwood forest. Global Biogeochem. Cycles 2010, 24.
  33. 33
    Schwarz, G. Estimating the dimension of a model Ann. Stat. 1978, 6 (2) 461 464
  34. 34
    Walker, N. D.; Rabalais, N. N. Relationships among satellite chlorophyll a, river inputs, and hypoxia on the Louisiana continental shelf, gulf of Mexico Estuaries Coasts 2006, 29 (6B) 1081 1093
  35. 35
    Wiseman, W. J.; Kelly, F. J. Salinity variability within the Louisiana coastal current during the 1982 flood season Estuaries 1994, 17 (4) 732 739
  36. 36
    Zavala-Hidalgo, J.; Morey, S. L.; O’Brien, J. J. Seasonal circulation on the western shelf of the Gulf of Mexico using a high-resolution numerical model J. Geophys. Res., [Oceans] 2003, 108 (C12) C123389
  37. 37
    Scavia, D.; Donnelly, K. A. Reassessing hypoxia forecasts for the Gulf of Mexico Environ. Sci. Technol. 2007, 41 (23) 8111 8117
  38. 38
    Turner, R. E.; Rabalais, N. N.; Swenson, E. M.; Kasprzak, M.; Romaire, T. Summer hypoxia in the northern Gulf of Mexico and its prediction from 1978 to 1995 Mar. Environ. Res. 2005, 59 (1) 65 77
  39. 39
    Liu, Y.; Evans, M. A.; Scavia, D. Gulf of Mexico hypoxia: Exploring increasing sensitivity to nitrogen loads Environ. Sci. Technol. 2010, 44 (15) 5836 5841
  40. 40
    Turner, R. E.; Rabalais, N. N.; Justic, D. Gulf of Mexico hypoxia: Alternate states and a legacy Environ. Sci. Technol. 2008, 42 (7) 2323 2327
  41. 41
    Goolsby, D. A.; Battaglin, W. A. Long-term changes in concentrations and flux of nitrogen in the Mississippi River Basin, USA Hydrol. Processes 2001, 15 (7) 1209 1226
  42. 42
    Stow, C. A.; Qian, S. S.; Craig, J. K. Declining threshold for hypoxia in the Gulf of Mexico Environ. Sci. Technol. 2005, 39 (3) 716 723
  43. 43
    Kemp, W. M.; Testa, J. M.; Conley, D. J.; Gilbert, D.; Hagy, J. D. Temporal responses of coastal hypoxia to nutrient loading and physical controls Biogeosciences 2009, 6 (12) 2985 3008
  44. 44
    Conley, D. J.; Carstensen, J.; Vaquer-Sunyer, R.; Duarte, C. M. Ecosystem thresholds with hypoxia Hydrobiologia 2009, 629 (1) 21 29
  45. 45
    Hetland, R. D.; DiMarco, S. F. How does the character of oxygen demand control the structure of hypoxia on the Texas–Louisiana continental shelf? J. Mar. Syst. 2008, 70 (1–2) 49 62
  46. 46
    Wang, L. X.; Justic, D. A modeling study of the physical processes affecting the development of seasonal hypoxia over the inner Louisiana–Texas shelf: Circulation and stratification Cont. Shelf Res. 2009, 29 (11–12) 1464 1476
  47. 47
    Emery, W. J.; Thomson, R. E. Data Analysis Methods in Physical Oceanography, 2nd and rev. ed.; Elsevier: Amsterdam/New York, 2001.
  48. 48
    NOAA Gulf of Mexico Hypoxia Watch: Hypoxia Contour Process; http://www.ncddc.noaa.gov/hypoxia/. (Accessed October, 2012.)

Cited By


This article is cited by 74 publications.

  1. Venkata Rohith Reddy Matli, Arnaud Laurent, Katja Fennel, Kevin Craig, Jacob Krause, Daniel R. Obenour. Fusion-Based Hypoxia Estimates: Combining Geostatistical and Mechanistic Models of Dissolved Oxygen Variability. Environmental Science & Technology 2020, 54 (20) , 13016-13025. https://doi.org/10.1021/acs.est.0c03655
  2. Arnaud Laurent, Katja Fennel. Time-Evolving, Spatially Explicit Forecasts of the Northern Gulf of Mexico Hypoxic Zone. Environmental Science & Technology 2019, 53 (24) , 14449-14458. https://doi.org/10.1021/acs.est.9b05790
  3. V. Rohith Reddy Matli, Shiqi Fang, Joseph Guinness, Nancy. N. Rabalais, J. Kevin Craig, Daniel R. Obenour. Space-Time Geostatistical Assessment of Hypoxia in the Northern Gulf of Mexico. Environmental Science & Technology 2018, 52 (21) , 12484-12493. https://doi.org/10.1021/acs.est.8b03474
  4. Donald Scavia, Mary Anne Evans, and Daniel R. Obenour . A Scenario and Forecast Model for Gulf of Mexico Hypoxic Area and Volume. Environmental Science & Technology 2013, 47 (18) , 10423-10428. https://doi.org/10.1021/es4025035
  5. Richard A Snyder, Joseph A Moss, Luciana Santoferrara, Marie Head, Wade H Jeffrey, . Ciliate microzooplankton from the Northeastern Gulf of Mexico. ICES Journal of Marine Science 2021, 78 (9) , 3356-3371. https://doi.org/10.1093/icesjms/fsab002
  6. Donald Scavia, Isabella Bertani, Jeremy M. Testa, Aaron J. Bever, Joel D. Blomquist, Marjorie A. M. Friedrichs, Lewis C. Linker, Bruce D. Michael, Rebecca R. Murphy, Gary W. Shenk. Advancing estuarine ecological forecasts: seasonal hypoxia in Chesapeake Bay. Ecological Applications 2021, 31 (6) https://doi.org/10.1002/eap.2384
  7. Elizabeth D. LaBone, Kenneth A. Rose, Dubravko Justic, Haosheng Huang, Lixia Wang. Effects of spatial variability on the exposure of fish to hypoxia: a modeling analysis for the Gulf of Mexico. Biogeosciences 2021, 18 (2) , 487-507. https://doi.org/10.5194/bg-18-487-2021
  8. Veronica Ruiz Xomchuk, Robert D. Hetland, Lixin Qu. Small‐Scale Variability of Bottom Oxygen in the Northern Gulf of Mexico. Journal of Geophysical Research: Oceans 2021, 126 (1) https://doi.org/10.1029/2020JC016279
  9. Ryan T. Munnelly, David B. Reeves, Edward J. Chesney, Donald M. Baltz. Spatial and Temporal Influences of Nearshore Hydrography on Fish Assemblages Associated with Energy Platforms in the Northern Gulf of Mexico. Estuaries and Coasts 2021, 44 (1) , 269-285. https://doi.org/10.1007/s12237-020-00772-7
  10. Jiatang Hu, Zhongren Zhang, Bin Wang, Jia Huang. Long-term spatiotemporal variations in and expansion of low-oxygen conditions in the Pearl River estuary: a study synthesizing observations during 1976–2017. Biogeosciences 2021, 18 (18) , 5247-5264. https://doi.org/10.5194/bg-18-5247-2021
  11. Serghei A. Bocaniov, Kevin G. Lamb, Wentao Liu, Yerubandi R. Rao, Ralph E. H. Smith. High Sensitivity of Lake Hypoxia to Air Temperatures, Winds, and Nutrient Loading: Insights From a 3‐D Lake Model. Water Resources Research 2020, 56 (12) https://doi.org/10.1029/2019WR027040
  12. Wenfei Ni, Ming Li, Jeremy M. Testa. Discerning effects of warming, sea level rise and nutrient management on long-term hypoxia trends in Chesapeake Bay. Science of The Total Environment 2020, 737 , 139717. https://doi.org/10.1016/j.scitotenv.2020.139717
  13. Hongjie Wang, John Lehrter, Kanchan Maiti, Katja Fennel, Arnaud Laurent, Nancy Rabalais, Najid Hussain, Qian Li, Baoshan Chen, K. Michael Scaboo, Wei‐Jun Cai. Benthic Respiration in Hypoxic Waters Enhances Bottom Water Acidification in the Northern Gulf of Mexico. Journal of Geophysical Research: Oceans 2020, 125 (10) https://doi.org/10.1029/2020JC016152
  14. Elena Vasechkina. Object-Based Modeling of Marine Phytoplankton and Seaweeds. Journal of Marine Science and Engineering 2020, 8 (9) , 685. https://doi.org/10.3390/jmse8090685
  15. Wenjing Li, Huaiyang Fang, Guangxiong Qin, Xiuqin Tan, Zhiwei Huang, Fantang Zeng, Hongwei Du, Shuping Li. Concentration estimation of dissolved oxygen in Pearl River Basin using input variable selection and machine learning techniques. Science of The Total Environment 2020, 731 , 139099. https://doi.org/10.1016/j.scitotenv.2020.139099
  16. Brandon M. Jarvis, John C. Lehrter, Lisa L. Lowe, James D. Hagy, Yongshan Wan, Michael C. Murrell, Dong S. Ko, Bradley Penta, Richard W. Gould. Modeling Spatiotemporal Patterns of Ecosystem Metabolism and Organic Carbon Dynamics Affecting Hypoxia on the Louisiana Continental Shelf. Journal of Geophysical Research: Oceans 2020, 125 (4) https://doi.org/10.1029/2019JC015630
  17. Dario Del Giudice, V. R. R. Matli, Daniel R. Obenour. Bayesian mechanistic modeling characterizes Gulf of Mexico hypoxia: 1968–2016 and future scenarios. Ecological Applications 2020, 30 (2) https://doi.org/10.1002/eap.2032
  18. Todd V. Royer. Stoichiometry of nitrogen, phosphorus, and silica loads in the Mississippi‐Atchafalaya River basin reveals spatial and temporal patterns in risk for cyanobacterial blooms. Limnology and Oceanography 2020, 65 (2) , 325-335. https://doi.org/10.1002/lno.11300
  19. Xiuqin Li, Chuqian Lu, Yafeng Zhang, Huade Zhao, Juying Wang, Hongbin Liu, Kedong Yin. Low dissolved oxygen in the Pearl River estuary in summer: Long-term spatio-temporal patterns, trends, and regulating factors. Marine Pollution Bulletin 2020, 151 , 110814. https://doi.org/10.1016/j.marpolbul.2019.110814
  20. Taavi Liblik, Yijing Wu, Daidu Fan, Dinghui Shang. Wind-driven stratification patterns and dissolved oxygen depletion off the Changjiang (Yangtze) Estuary. Biogeosciences 2020, 17 (10) , 2875-2895. https://doi.org/10.5194/bg-17-2875-2020
  21. Fabian A. Gomez, Sang-Ki Lee, Frank J. Hernandez, Luciano M. Chiaverano, Frank E. Muller-Karger, Yanyun Liu, John T. Lamkin. ENSO-induced co-variability of Salinity, Plankton Biomass and Coastal Currents in the Northern Gulf of Mexico. Scientific Reports 2019, 9 (1) https://doi.org/10.1038/s41598-018-36655-y
  22. Shiqi Fang, Dario Del Giudice, Donald Scavia, Caren E. Binding, Thomas B. Bridgeman, Justin D. Chaffin, Mary Anne Evans, Joseph Guinness, Thomas H. Johengen, Daniel R. Obenour. A space-time geostatistical model for probabilistic estimation of harmful algal bloom biomass and areal extent. Science of The Total Environment 2019, 695 , 133776. https://doi.org/10.1016/j.scitotenv.2019.133776
  23. Jingbo Gao, Yongli Lu, Zhujun Chen, Lei Wang, Jianbin Zhou. Land‐use change from cropland to orchard leads to high nitrate accumulation in the soils of a small catchment. Land Degradation & Development 2019, 30 (17) , 2150-2161. https://doi.org/10.1002/ldr.3412
  24. Nancy N. Rabalais, R. Eugene Turner. Gulf of Mexico Hypoxia: Past, Present, and Future. Limnology and Oceanography Bulletin 2019, 28 (4) , 117-124. https://doi.org/10.1002/lob.10351
  25. Elizabeth LaBone, Dubravko Justic, Kenneth Rose, Lixia Wang, Haosheng Huang. Modeling Fish Movement in 3-D in the Gulf of Mexico Hypoxic Zone. Estuaries and Coasts 2019, 42 (6) , 1662-1685. https://doi.org/10.1007/s12237-019-00601-6
  26. Fabian Große, Katja Fennel, Arnaud Laurent. Quantifying the Relative Importance of Riverine and Open‐Ocean Nitrogen Sources for Hypoxia Formation in the Northern Gulf of Mexico. Journal of Geophysical Research: Oceans 2019, 124 (8) , 5451-5467. https://doi.org/10.1029/2019JC015230
  27. Ryan T. Munnelly, David B. Reeves, Edward J. Chesney, Donald M. Baltz. Summertime hydrography of the nearshore Louisiana Continental Shelf: Effects of riverine outflow, shelf morphology, and the presence of sand shoals on water quality. Continental Shelf Research 2019, 179 , 18-36. https://doi.org/10.1016/j.csr.2019.04.002
  28. David B. Reeves, Edward J. Chesney, Ryan T. Munnelly, Donald M. Baltz, Kanchan Maiti. Trophic Ecology of Sheepshead and Stone Crabs at Oil and Gas Platforms in the Northern Gulf of Mexico's Hypoxic Zone. Transactions of the American Fisheries Society 2019, 148 (2) , 324-338. https://doi.org/10.1002/tafs.10135
  29. Donald Scavia, Dubravko Justić, Daniel R Obenour, J Kevin Craig, Lixia Wang. Hypoxic volume is more responsive than hypoxic area to nutrient load reductions in the northern Gulf of Mexico—and it matters to fish and fisheries. Environmental Research Letters 2019, 14 (2) , 024012. https://doi.org/10.1088/1748-9326/aaf938
  30. RT Munnelly, DB Reeves, EJ Chesney, DM Baltz, BD Marx. Habitat suitability for oil and gas platform-associated fishes in Louisiana’s nearshore waters. Marine Ecology Progress Series 2019, 608 , 199-219. https://doi.org/10.3354/meps12772
  31. Zong-Pei Jiang, Wei-Jun Cai, John Lehrter, Baoshan Chen, Zhangxian Ouyang, Chengfeng Le, Brian J. Roberts, Najid Hussain, Michael K. Scaboo, Junxiao Zhang, Yuanyuan Xu. Spring net community production and its coupling with the CO2 dynamics in the surface water of the northern Gulf of Mexico. Biogeosciences 2019, 16 (18) , 3507-3525. https://doi.org/10.5194/bg-16-3507-2019
  32. David B. Reeves, Edward J. Chesney, Ryan T. Munnelly, Donald M. Baltz, Brian D. Marx. Abundance and Distribution of Reef-Associated Fishes Around Small Oil and Gas Platforms in the Northern Gulf of Mexico’s Hypoxic Zone. Estuaries and Coasts 2018, 41 (7) , 1835-1847. https://doi.org/10.1007/s12237-017-0349-4
  33. Jason K. Jolliff, Ewa Jarosz, Sherwin Ladner, Travis Smith, Stephanie Anderson, James Dykes. The Optical Signature of a Bottom Boundary Layer Ventilation Event in the Northern Gulf of Mexico's Hypoxic Zone. Geophysical Research Letters 2018, 45 (16) , 8390-8398. https://doi.org/10.1029/2018GL078228
  34. Il-Nam Kim. Estimating “Mean-State” July (1985–2007) N2O Fluxes in the Northern Gulf of Mexico Hypoxic Region: Variation, Distribution, and Implication. Frontiers in Marine Science 2018, 5 https://doi.org/10.3389/fmars.2018.00249
  35. Brian Fry. Respiration and Metabolic Age as Controls of Bottom Water Hypoxia on the Louisiana Continental Shelf; 18Δ as the Ghost of Respiration Past. Estuaries and Coasts 2018, 41 (5) , 1297-1313. https://doi.org/10.1007/s12237-018-0365-z
  36. K. J. Van Meter, P. Van Cappellen, N. B. Basu. Legacy nitrogen may prevent achievement of water quality goals in the Gulf of Mexico. Science 2018, 360 (6387) , 427-430. https://doi.org/10.1126/science.aar4462
  37. Kenneth A. Rose, Sean Creekmore, Peter Thomas, J. Kevin Craig, Md Saydur Rahman, Rachael Miller Neilan. Modeling the Population Effects of Hypoxia on Atlantic Croaker (Micropogonias undulatus) in the Northwestern Gulf of Mexico: Part 1—Model Description and Idealized Hypoxia. Estuaries and Coasts 2018, 41 (1) , 233-254. https://doi.org/10.1007/s12237-017-0266-6
  38. Kenneth A. Rose, Sean Creekmore, Dubravko Justić, Peter Thomas, J. Kevin Craig, Rachael Miller Neilan, Lixia Wang, Md Saydur Rahman, David Kidwell. Modeling the Population Effects of Hypoxia on Atlantic Croaker (Micropogonias undulatus) in the Northwestern Gulf of Mexico: Part 2—Realistic Hypoxia and Eutrophication. Estuaries and Coasts 2018, 41 (1) , 255-279. https://doi.org/10.1007/s12237-017-0267-5
  39. Nancy N. Rabalais, Leslie M. Smith, R. Eugene Turner. The Deepwater Horizon oil spill and Gulf of Mexico shelf hypoxia. Continental Shelf Research 2018, 152 , 98-107. https://doi.org/10.1016/j.csr.2017.11.007
  40. S. Kent Hoekman, Amber Broch, Xiaowei (Vivian) Liu. Environmental implications of higher ethanol production and use in the U.S.: A literature review. Part I – Impacts on water, soil, and air quality. Renewable and Sustainable Energy Reviews 2018, 81 , 3140-3158. https://doi.org/10.1016/j.rser.2017.05.050
  41. Katja Fennel, Arnaud Laurent. N and P as ultimate and proximate limiting nutrients in the northern Gulf of Mexico: implications for hypoxia reduction strategies. Biogeosciences 2018, 15 (10) , 3121-3131. https://doi.org/10.5194/bg-15-3121-2018
  42. Kevin M. Purcell, J. Kevin Craig, James M. Nance, Martin D. Smith, Lori S. Bennear, . Fleet behavior is responsive to a large-scale environmental disturbance: Hypoxia effects on the spatial dynamics of the northern Gulf of Mexico shrimp fishery. PLOS ONE 2017, 12 (8) , e0183032. https://doi.org/10.1371/journal.pone.0183032
  43. Donald Scavia, Isabella Bertani, Daniel R. Obenour, R. Eugene Turner, David R. Forrest, Alexey Katin. Ensemble modeling informs hypoxia management in the northern Gulf of Mexico. Proceedings of the National Academy of Sciences 2017, 114 (33) , 8823-8828. https://doi.org/10.1073/pnas.1705293114
  44. Xinping Hu, Qian Li, Wei-Jen Huang, Baoshan Chen, Wei-Jun Cai, Nancy N. Rabalais, R. Eugene Turner. Effects of eutrophication and benthic respiration on water column carbonate chemistry in a traditional hypoxic zone in the Northern Gulf of Mexico. Marine Chemistry 2017, 194 , 33-42. https://doi.org/10.1016/j.marchem.2017.04.004
  45. E. Sinha, A. M. Michalak, V. Balaji. Eutrophication will increase during the 21st century as a result of precipitation changes. Science 2017, 357 (6349) , 405-408. https://doi.org/10.1126/science.aan2409
  46. David B. Reeves, Ryan T. Munnelly, Edward J. Chesney, Donald M. Baltz, Brian D. Marx. Stone Crab Menippe spp. Populations on Louisiana’s Nearshore Oil and Gas Platforms: Higher Density and Size at Maturity on a Sand Shoal. Transactions of the American Fisheries Society 2017, 146 (3) , 371-383. https://doi.org/10.1080/00028487.2017.1281164
  47. Martin D. Smith, Atle Oglend, A. Justin Kirkpatrick, Frank Asche, Lori S. Bennear, J. Kevin Craig, James M. Nance. Seafood prices reveal impacts of a major ecological disturbance. Proceedings of the National Academy of Sciences 2017, 114 (7) , 1512-1517. https://doi.org/10.1073/pnas.1617948114
  48. Arnaud Laurent, Katja Fennel, Wei‐Jun Cai, Wei‐Jen Huang, Leticia Barbero, Rik Wanninkhof. Eutrophication‐induced acidification of coastal waters in the northern Gulf of Mexico: Insights into origin and processes from a coupled physical‐biogeochemical model. Geophysical Research Letters 2017, 44 (2) , 946-956. https://doi.org/10.1002/2016GL071881
  49. Elizabeth LaBone, Dubravko Justic, Kenneth A. Rose, Lixia Wang, Haosheng Huang. Comparing Default Movement Algorithms for Individual Fish Avoidance of Hypoxia in the Gulf of Mexico. 2017,,, 239-278. https://doi.org/10.1007/978-3-319-54571-4_10
  50. Kenneth A. Rose, Sean Creekmore, Shaye Sable. Simulation of the Population-Level Responses of Fish to Hypoxia: Should We Expect Sampling to Detect Responses?. 2017,,, 359-376. https://doi.org/10.1007/978-3-319-54571-4_13
  51. Kenneth A. Rose, Dubravko Justic, Katja Fennel, Robert D. Hetland. Numerical Modeling of Hypoxia and Its Effects: Synthesis and Going Forward. 2017,,, 401-421. https://doi.org/10.1007/978-3-319-54571-4_15
  52. Arnaud Laurent, Katja Fennel. Modeling River-Induced Phosphorus Limitation in the Context of Coastal Hypoxia. 2017,,, 149-171. https://doi.org/10.1007/978-3-319-54571-4_7
  53. John C. Lehrter, Dong S. Ko, Lisa L. Lowe, Bradley Penta. Predicted Effects of Climate Change on Northern Gulf of Mexico Hypoxia. 2017,,, 173-214. https://doi.org/10.1007/978-3-319-54571-4_8
  54. Brian J. Langseth, Amy M. Schueller, Kyle W. Shertzer, J.Kevin Craig, Joseph W. Smith. Management implications of temporally and spatially varying catchability for the Gulf of Mexico menhaden fishery. Fisheries Research 2016, 181 , 186-197. https://doi.org/10.1016/j.fishres.2016.04.013
  55. Katja Fennel, Arnaud Laurent, Robert Hetland, Dubravko Justić, Dong S. Ko, John Lehrter, Michael Murrell, Lixia Wang, Liuqian Yu, Wenxia Zhang. Effects of model physics on hypoxia simulations for the northern Gulf of Mexico: A model intercomparison. Journal of Geophysical Research: Oceans 2016, 121 (8) , 5731-5750. https://doi.org/10.1002/2015JC011577
  56. Serghei A. Bocaniov, Donald Scavia. Temporal and spatial dynamics of large lake hypoxia: Integrating statistical and three-dimensional dynamic models to enhance lake management criteria. Water Resources Research 2016, 52 (6) , 4247-4263. https://doi.org/10.1002/2015WR018170
  57. Dong Ko, Richard Gould, Bradley Penta, John Lehrter. Impact of Satellite Remote Sensing Data on Simulations of Coastal Circulation and Hypoxia on the Louisiana Continental Shelf. Remote Sensing 2016, 8 (5) , 435. https://doi.org/10.3390/rs8050435
  58. Chengfeng Le, John C. Lehrter, Chuanmin Hu, Daniel R. Obenour. Satellite‐based empirical models linking river plume dynamics with hypoxic area and volume. Geophysical Research Letters 2016, 43 (6) , 2693-2699. https://doi.org/10.1002/2015GL067521
  59. Jochen Kämpf, Piers Chapman. Seasonal Wind-Driven Coastal Upwelling Systems. 2016,,, 315-361. https://doi.org/10.1007/978-3-319-42524-5_8
  60. Jordon S. Beckler, Nicole Kiriazis, Christophe Rabouille, Frank J. Stewart, Martial Taillefert. Importance of microbial iron reduction in deep sediments of river-dominated continental-margins. Marine Chemistry 2016, 178 , 22-34. https://doi.org/10.1016/j.marchem.2015.12.003
  61. Jianhong Xue, Wei‐Jun Cai, Xinping Hu, Wei‐Jen Huang, Steven E. Lohrenz, Kjell Gundersen. Temporal variation and stoichiometric ratios of organic matter remineralization in bottom waters of the northern G ulf of M exico during late spring and summer. Journal of Geophysical Research: Oceans 2015, 120 (12) , 8304-8326. https://doi.org/10.1002/2015JC011453
  62. Mark D. Rowe, Daniel R. Obenour, Thomas F. Nalepa, Henry A. Vanderploeg, Foad Yousef, W. Charles Kerfoot. Mapping the spatial distribution of the biomass and filter-feeding effect of invasive dreissenid mussels on the winter-spring phytoplankton bloom in Lake Michigan. Freshwater Biology 2015, 60 (11) , 2270-2285. https://doi.org/10.1111/fwb.12653
  63. Songjie He, Y. Jun Xu. Three Decadal Inputs of Nitrogen and Phosphorus from Four Major Coastal Rivers to the Summer Hypoxic Zone of the Northern Gulf of Mexico. Water, Air, & Soil Pollution 2015, 226 (9) https://doi.org/10.1007/s11270-015-2580-6
  64. Liuqian Yu, Katja Fennel, Arnaud Laurent. A modeling study of physical controls on hypoxia generation in the northern Gulf of Mexico. Journal of Geophysical Research: Oceans 2015, 120 (7) , 5019-5039. https://doi.org/10.1002/2014JC010634
  65. Daniel R. Obenour, Anna M. Michalak, Donald Scavia. Assessing biophysical controls on Gulf of Mexico hypoxia through probabilistic modeling. Ecological Applications 2015, 25 (2) , 492-505. https://doi.org/10.1890/13-2257.1
  66. Gerald Whittaker, Bradley L. Barnhart, Raghavan Srinivasan, Jeffrey G. Arnold. Cost of areal reduction of gulf hypoxia through agricultural practice. Science of The Total Environment 2015, 505 , 149-153. https://doi.org/10.1016/j.scitotenv.2014.09.101
  67. MH Monk, JE Powers, EN Brooks. Spatial patterns in species assemblages associated with the northwestern Gulf of Mexico shrimp trawl fishery. Marine Ecology Progress Series 2015, 519 , 1-12. https://doi.org/10.3354/meps11150
  68. Songjie He, Y.Jun Xu. Three decadal inputs of total organic carbon from four major coastal river basins to the summer hypoxic zone of the Northern Gulf of Mexico. Marine Pollution Bulletin 2015, 90 (1-2) , 121-128. https://doi.org/10.1016/j.marpolbul.2014.11.010
  69. L. Yu, K. Fennel, A. Laurent, M. C. Murrell, J. C. Lehrter. Numerical analysis of the primary processes controlling oxygen dynamics on the Louisiana shelf. Biogeosciences 2015, 12 (7) , 2063-2076. https://doi.org/10.5194/bg-12-2063-2015
  70. Eurico J. D’Sa. Assessment of chlorophyll variability along the Louisiana coast using multi-satellite data. GIScience & Remote Sensing 2014, 51 (2) , 139-157. https://doi.org/10.1080/15481603.2014.895578
  71. Yang Feng, Katja Fennel, George A. Jackson, Steven F. DiMarco, Robert D. Hetland. A model study of the response of hypoxia to upwelling-favorable wind on the northern Gulf of Mexico shelf. Journal of Marine Systems 2014, 131 , 63-73. https://doi.org/10.1016/j.jmarsys.2013.11.009
  72. Arnaud Laurent, Katja Fennel, , . Simulated reduction of hypoxia in the northern Gulf of Mexico due to phosphorus limitation. Elementa: Science of the Anthropocene 2014, 2 https://doi.org/10.12952/journal.elementa.000022
  73. Dubravko Justić, Lixia Wang. Assessing temporal and spatial variability of hypoxia over the inner Louisiana–upper Texas shelf: Application of an unstructured-grid three-dimensional coupled hydrodynamic-water quality model. Continental Shelf Research 2014, 72 , 163-179. https://doi.org/10.1016/j.csr.2013.08.006
  74. Brian J. Langseth, Kevin M. Purcell, J. Kevin Craig, Amy M. Schueller, Joseph W. Smith, Kyle W. Shertzer, Sean Creekmore, Kenneth A. Rose, Katja Fennel. Effect of Changes in Dissolved Oxygen Concentrations on the Spatial Dynamics of the Gulf Menhaden Fishery in the Northern Gulf of Mexico. Marine and Coastal Fisheries 2014, 6 (1) , 223-234. https://doi.org/10.1080/19425120.2014.949017
  • Abstract

    Figure 1

    Figure 1. Number of locations sampled during the annual midsummer shelfwide cruises using hand-held and rosette instruments.

    Figure 2

    Figure 2. Study area bathymetry, sampling, and estimation locations.

    Figure 3

    Figure 3. Bottom layer hypoxic extent estimates with 95% confidence intervals by year; estimates prior to making adjustments for instrument bias as triangles; previous LUMCON area estimates as open squares; revised LUMCON area estimates as solid squares.

    Figure 4

    Figure 4. Example maps of estimated bottom layer hypoxic thickness (median values from CRs), 2001–2008; observation locations shown as white dots.

  • References

    ARTICLE SECTIONS
    Jump To

    This article references 48 other publications.

    1. 1
      Rabalais, N. N.; Turner, R. E.; Wiseman, W. J. Gulf of Mexico hypoxia, aka “The dead zone” Annu. Rev. Ecol. Syst. 2002, 33, 235 263
    2. 2
      LUMCON hypoxia in the northern Gulf of Mexico: Research; http://www.gulfhypoxia.net/Research/. (Accessed October, 2012.)
    3. 3
      Rabalais, N. N.; Turner, R. E.; Scavia, D. Beyond science into policy: Gulf of Mexico hypoxia and the Mississippi River Bioscience 2002, 52 (2) 129 142
    4. 4
      Rabalais, N. N.; Diaz, R. J.; Levin, L. A.; Turner, R. E.; Gilbert, D.; Zhang, J. Dynamics and distribution of natural and human-caused hypoxia Biogeosciences 2010, 7 (2) 585 619
    5. 5
      Rabalais, N. N.; Turner, R. E.Coastal Hypoxia: Consequences for Living Resources and Ecosystems. Coastal and Estuarine Studies 58; American Geophysical Union: Washington, D.C., 2001.
    6. 6
      Mississippi River/Gulf of Mexico Watershed Nutrient Task Force Action Plan for Reducing, Mitigating, And Controlling Hypoxia in the Northern Gulf of Mexico; Office of Wetlands, Oceans, and Watersheds, U.S. Environmental Protection Agency: Washington, DC, 2001, http://water.epa.gov/type/watersheds/named/msbasin/history.cfm.
    7. 7
      Mississippi River/Gulf of Mexico Watershed Nutrient Task Force Gulf Hypoxia Action Plan 2008 for Reducing Mitigating, And Controlling Hypoxia in the Northern Gulf of Mexico and Improving Water Quality in the Mississippi River Basin; U.S. Environmental Protection Agency, Office of Wetlands, Oceans, and Watersheds: Washington, D.C., 2008, http://water.epa.gov/type/watersheds/named/msbasin/actionplan.cfm.
    8. 8
      Hypoxia in the Northern Gulf of Mexico, An Update by the EPA Science Advisory Board; EPA Science Advisory Board: Washington, DC, 2007, EPA-SAB-08–003.
    9. 9
      Obenour, D. R.; Michalak, A. M.; Zhou, Y.; Scavia, D. Quantifying the impacts of stratification and nutrient loading on hypoxia in the northern Gulf of Mexico Environ. Sci. Technol. 2012, 46 (10) 5489 5496
    10. 10
      Wiseman, W. J.; Rabalais, N. N.; Turner, R. E.; Dinnel, S. P.; MacNaughton, A. Seasonal and interannual variability within the Louisiana coastal current: stratification and hypoxia J. Mar. Syst. 1997, 12 (1–4) 237 248
    11. 11
      Forrest, D. R.; Hetland, R. D.; DiMarco, S. F. Multivariable statistical regression models of the areal extent of hypoxia over the Texas–Louisiana continental shelf Environ. Res. Lett. 2011, 6 (4) 045002
    12. 12
      Greene, R. M.; Lehrter, J. C.; Hagy, J. D. Multiple regression models for hindcasting and forecasting midsummer hypoxia in the Gulf of Mexico Ecol. Appl. 2009, 19 (5) 1161 1175
    13. 13
      Scavia, D.; Rabalais, N. N.; Turner, R. E.; Justic, D.; Wiseman, W. J. Predicting the response of Gulf of Mexico hypoxia to variations in Mississippi River nitrogen load Limnol. Oceanogr. 2003, 48 (3) 951 956
    14. 14
      Turner, R. E.; Rabalais, N. N.; Justic, D. Predicting summer hypoxia in the northern Gulf of Mexico: Riverine N, P, and Si loading Mar. Pollut. Bull. 2006, 52 (2) 139 148
    15. 15
      Rabalais, N. N.; Turner, R.E.; Justic, D.; Dortch, Q.; Wiseman, W.Characterization of Hypoxia: Coastal Ocean Program Decision Analysis Series No. 15; National Oceanic and Atmospheric Administration (NOAA): Silver Springs, MD, 1999.
    16. 16
      Rabalais, N. N.; Turner, R. E.; Sen Gupta, B. K.; Boesch, D. F.; Chapman, P.; Murrell, M. C. Hypoxia in the northern Gulf of Mexico: Does the science support the plan to reduce, mitigate, and control hypoxia? Estuaries Coasts 2007, 30 (5) 753 772
    17. 17
      Chiles, J.-P.; Delfiner, P. Geostatistics: Modeling Spatial Uncertainty; Wiley: New York, 1999.
    18. 18
      Zhou, Y.; Obenour, D. R.; Scavia, D.; Johengen, T. H.; Michalak, A. M. Spatial and temporal trends in Lake Erie hypoxia, 1987–2007 Environ. Sci. Technol. 2013, 47 (2) 899 905
    19. 19
      Zhang, H. Y.; Ludsin, S. A.; Mason, D. M.; Adamack, A. T.; Brandt, S. B.; Zhang, X. S.; Kimmel, D. G.; Roman, M. R.; Boicourt, W. C. Hypoxia-driven changes in the behavior and spatial distribution of pelagic fish and mesozooplankton in the northern Gulf of Mexico J. Exp. Mar. Biol. Ecol. 2009, 381, S80 S91
    20. 20
      Kimmel, D. G.; Boicourt, W. C.; Pierson, J. J.; Roman, M. R.; Zhang, X. S. A comparison of the mesozooplankton response to hypoxia in Chesapeake Bay and the northern Gulf of Mexico using the biomass size spectrum J. Exp. Mar. Biol. Ecol. 2009, 381, S65 S73
    21. 21
      Murphy, R. R.; Curriero, F. C.; Ball, W. P. Comparison of spatial interpolation methods for water quality evaluation in the Chesapeake Bay J. Environ. Eng.–ASCE 2010, 136 (2) 160 171
    22. 22
      Zhou, Y.; Scavia, D.; Michalak, A. M. Nutrient loading and meteorological conditions explain interannual variability of hypoxia in the Chesapeake Bay, (submitted manuscript).
    23. 23
      Turner, R. E.; Rabalais, N. N.; Justic, D. Predicting summer hypoxia in the northern Gulf of Mexico: Redux Mar. Pollut. Bull. 2012, 64 (2) 319 324
    24. 24
      Rabalais, N. N. Louisiana Hypoxia Surveys, available from NOAA National Oceanographic Data Center (NODC); http://www.nodc.noaa.gov/General/getdata.html. (Accessed January, 2011.)
    25. 25
      NOAA National Geophysical Data Center (NGDC) Coastal ReliefModel; http://www.ngdc.noaa.gov/mgg/coastal/startcrm.htm. (Accessed June, 2010.)
    26. 26
      Zimmerman, D. L.; Stein, M. Classical Geostatistical Methods. In Handbook of Spatial Statistics; Gelfand, A. E.; Diggle, P. J.; Fuentes, M.; Guttorp, P., Eds.; CRC Press: Boca Raton, FL, 2010.
    27. 27
      Erickson, T. A.; Williams, M. W.; Winstral, A. Persistence of topographic controls on the spatial distribution of snow in rugged mountain terrain, Colorado, United States. Water Resour. Res. 2005, 41 (4).
    28. 28
      Goovaerts, P. Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall J. Hydrol. 2000, 228 (1–2) 113 129
    29. 29
      Qian, S. S. Estimating the area affected by phosphorus runoff in an Everglades wetland: A comparison of universal kriging and Bayesian kriging Environ. Ecol. Stat. 1997, 4 (1) 1 26
    30. 30
      Patterson, H. D.; Thompson, R. Recovery of inter-block information when block sizes are unequal Biometrika 1971, 58 (3) 545 554
    31. 31
      Zimmerman, D. L. Likelihood-Based Methods. In Handbook of Spatial Statistics; Gelfand, A. E.; Diggle, P. J.; Fuentes, M.; Guttorp, P., Eds.; CRC Press: Boca Raton, FL, 2010.
    32. 32
      Mueller, K. L.; Yadav, V.; Curtis, P. S.; Vogel, C.; Michalak, A. M. Attributing the variability of eddy-covariance CO2 flux measurements across temporal scales using geostatistical regression for a mixed northern hardwood forest. Global Biogeochem. Cycles 2010, 24.
    33. 33
      Schwarz, G. Estimating the dimension of a model Ann. Stat. 1978, 6 (2) 461 464
    34. 34
      Walker, N. D.; Rabalais, N. N. Relationships among satellite chlorophyll a, river inputs, and hypoxia on the Louisiana continental shelf, gulf of Mexico Estuaries Coasts 2006, 29 (6B) 1081 1093
    35. 35
      Wiseman, W. J.; Kelly, F. J. Salinity variability within the Louisiana coastal current during the 1982 flood season Estuaries 1994, 17 (4) 732 739
    36. 36
      Zavala-Hidalgo, J.; Morey, S. L.; O’Brien, J. J. Seasonal circulation on the western shelf of the Gulf of Mexico using a high-resolution numerical model J. Geophys. Res., [Oceans] 2003, 108 (C12) C123389
    37. 37
      Scavia, D.; Donnelly, K. A. Reassessing hypoxia forecasts for the Gulf of Mexico Environ. Sci. Technol. 2007, 41 (23) 8111 8117
    38. 38
      Turner, R. E.; Rabalais, N. N.; Swenson, E. M.; Kasprzak, M.; Romaire, T. Summer hypoxia in the northern Gulf of Mexico and its prediction from 1978 to 1995 Mar. Environ. Res. 2005, 59 (1) 65 77
    39. 39
      Liu, Y.; Evans, M. A.; Scavia, D. Gulf of Mexico hypoxia: Exploring increasing sensitivity to nitrogen loads Environ. Sci. Technol. 2010, 44 (15) 5836 5841
    40. 40
      Turner, R. E.; Rabalais, N. N.; Justic, D. Gulf of Mexico hypoxia: Alternate states and a legacy Environ. Sci. Technol. 2008, 42 (7) 2323 2327
    41. 41
      Goolsby, D. A.; Battaglin, W. A. Long-term changes in concentrations and flux of nitrogen in the Mississippi River Basin, USA Hydrol. Processes 2001, 15 (7) 1209 1226
    42. 42
      Stow, C. A.; Qian, S. S.; Craig, J. K. Declining threshold for hypoxia in the Gulf of Mexico Environ. Sci. Technol. 2005, 39 (3) 716 723
    43. 43
      Kemp, W. M.; Testa, J. M.; Conley, D. J.; Gilbert, D.; Hagy, J. D. Temporal responses of coastal hypoxia to nutrient loading and physical controls Biogeosciences 2009, 6 (12) 2985 3008
    44. 44
      Conley, D. J.; Carstensen, J.; Vaquer-Sunyer, R.; Duarte, C. M. Ecosystem thresholds with hypoxia Hydrobiologia 2009, 629 (1) 21 29
    45. 45
      Hetland, R. D.; DiMarco, S. F. How does the character of oxygen demand control the structure of hypoxia on the Texas–Louisiana continental shelf? J. Mar. Syst. 2008, 70 (1–2) 49 62
    46. 46
      Wang, L. X.; Justic, D. A modeling study of the physical processes affecting the development of seasonal hypoxia over the inner Louisiana–Texas shelf: Circulation and stratification Cont. Shelf Res. 2009, 29 (11–12) 1464 1476
    47. 47
      Emery, W. J.; Thomson, R. E. Data Analysis Methods in Physical Oceanography, 2nd and rev. ed.; Elsevier: Amsterdam/New York, 2001.
    48. 48
      NOAA Gulf of Mexico Hypoxia Watch: Hypoxia Contour Process; http://www.ncddc.noaa.gov/hypoxia/. (Accessed October, 2012.)
  • Supporting Information

    Supporting Information

    ARTICLE SECTIONS
    Jump To

    A description of the instrument bias adjustments, a graphical test of linearity for the deterministic components of the models, the intercept parameters for the categorical variables, a set of results maps, tabulated bottom layer hypoxic extent results, and a summary of the models for MinDO and THF. This material is available free of charge via the Internet at http://pubs.acs.org.


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

You’ve supercharged your research process with ACS and Mendeley!

STEP 1:
Click to create an ACS ID

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

MENDELEY PAIRING EXPIRED
Your Mendeley pairing has expired. Please reconnect

This website uses cookies to improve your user experience. By continuing to use the site, you are accepting our use of cookies. Read the ACS privacy policy.

CONTINUE