Search

Strontium and barium to calcium ratios are often used as proxies for tracking animal movement across salinity gradients. As sea level rise continues, many estuarine rivers in the world face saltwater intrusion, which may cause changes in mobility and distribution of these metals upstream. Despite intensive research on metal adsorption and desorption in marine systems, knowledge of the spatiotemporal distribution of these elements along estuarine rivers is still limited. In this study, we conducted an intensive monitoring of Sr and Ba dynamics along an 88 km long estuary, the Calcasieu River in South Louisiana, USA, which has been strongly affected by saltwater intrusion. Over the period from May 2013 to August 2015, we collected monthly water samples and performed in-situ water quality measurements at six sites from the upstream to the river mouth, with a salinity range from 0.02 to 29.50 ppt. Water samples were analyzed for Sr, Ba, and Ca concentrations. In-situ measurements were made on salinity, pH, water temperature, dissolved oxygen concentration, and specific conductance. We found that the Sr and Ca concentrations and the Sr =Ca ratio all increased significantly with increasing salinity. The average Sr concentration at the site closest to the Gulf of Mexico (site 6) was 46.21 μmol L-1, which was about 130 times higher than that of the site furthest upstream (site 1, 0.35 ìmol L-1). The average Ca concentration at site 6 was 8.19 mmol L-1, which was about 60 times higher than that of site 1 (0.13 mmol L-1). The average Sr =Ca ratio at site 6 (8.41mmolmol-1) was about 3 times the average Sr = Ca ratio at site 1 (2.89 mmol mol-1). However, the spatial variation in Ba concentration was marginal, varying from 0.36 μmol L-1 at site 6 to 0.47 at site 5. The average Ba = Ca ratio at site 1 (4.82 mmol mol-1) was about 54 times the average Ba = Ca ratio at site 6 (0.09 mmol mol-1), showing a clear negative relation between the Ba = Ca ratio and increasing salinity. All the elemental concentrations and ratios had considerable seasonal variations, with significant differences among sampling months for the Sr, Ba concentrations and the Ba = Ca ratio (p < 0.01). The results from this study suggest that concentrations of Sr and Ca in the world's estuaries will very likely increase in the future as sea level rise continues. For low-gradient estuarine rivers such as the Calcasieu River in South Louisiana, USA, water chemistry upstream would experience substantial Sr and Ca enrichment, which could affect aquatic environments and biological communities. [ABSTRACT FROM AUTHOR]

We present evidence of temporal variation in nitrous oxide (N2O) concentrations in the bottom waters of the northern Gulf of Mexico (nGOM) hypoxic zone. The analysis is based on a conceptual model simulating N2O biogeochemical processes in conjunction with water-column O2 levels, derived from summer Texas-Louisiana shelf-wide hydrographic data for twenty Julys between 1985 and 2007. The mean modeled nGOM N2O concentration was 7.7±6.7 nmol L-1, and was significantly correlated with the areal extent of hypoxia. Our modeling analysis indicates that the nGOM is a persistent summer source of N2O, and nitrification is a primary factor leading to its production in this region. Based on the ongoing increase in the areal extent of hypoxia in the nGOM, we conclude that N2O emission from this environmentally stressed region will continue to increase into the future contributing to the global increase in greenhouse gases. [ABSTRACT FROM AUTHOR]

With an extension of > 40 km2 the recently discovered Campeche cold-water coral province located at the northeastern rim of the Campeche Bank in the southern Gulf of Mexico belongs to the largest coherent cold-water coral areas discovered so far. The Campeche province consists of numerous 20 to 40m high coral ridges that are developed in intermediate water depths of 500 to 600 m. The ridges are colonized by a vivid cold-water coral ecosystem that covers the upper flanks and summits. The rich coral community is dominated by the framework-building scleractinia Enallopsammia profunda and Lophelia pertusa while the associated benthic megafauna shows a rather scarce occurrence. The recent environmental setting is characterized by a high surface water production caused by a local upwelling center and a dynamic bottom water regime comprising vigorous bottom currents, internal waves and strong density contrasts, which all together provide optimal conditions for the growth of cold-water corals. The strong hydrodynamics -- potentially supported by the diel vertical migration of zooplankton in the Campeche area -- drive the delivering of food particles to the corals. The Campeche cold-water coral province is, thus, an excellent example highlighting the importance of the hydrographic setting in securing the food supply for the development of large and vivid cold-water coral ecosystems. [ABSTRACT FROM AUTHOR]

Deep-sea bivalves found at hydrothermal vents, cold seeps and organic falls are sustained by chemosynthetic bacteria which ensure part or all of their carbon nutrition. These symbioses are of prime importance for the functioning of the ecosystems. Similar symbioses occur in other bivalve species living in shallow and coastal reduced habitats worldwide. In recent years, several deep-sea species have been investigated from continental margins around Europe, West Africa, East America, the Gulf of Mexico, and from hydrothermal vents on the Mid-Atlantic Ridge. In parallel, numerous more easily accessible shallow marine species were studied. We here provide a summary of the current knowledge available on chemosymbiotic bivalves in the area ranging west-to-east from the Gulf of Mexico to Marmara Sea, and north-to-south from the Arctic to the Gulf of Guinea. Characteristics of symbioses in 51 species from the area are summarized for each of the five bivalve families documented to harbor chemosynthetic symbionts (Mytilidae, Vesicomyidae, Solemyidae, Thyasiridae and Lucinidae), and compared among families with special emphasis on ecology, life cycle, and connectivity. Chemosynthetic symbioses are a major adaptation to ecosystems and habitats exposed to reducing conditions, yet relatively little is known regarding their diversity and functioning apart from a few "model species" on which effort has focused over the last 30 yr. In the context of increasing concern about biodiversity and ecosystems, and increasing anthropogenic pressure on Oceans, we advocate for a better assessment of bivalve symbioses diversity in order to evaluate the capacities of these remarkable ecological and evolutionary units to withstand environmental change [ABSTRACT FROM AUTHOR]

The Texas-Louisiana shelf in the Northern Gulf of Mexico receives large inputs of nutrients and freshwater from the Mississippi/Atchafalaya River system. The nutrients stimulate high rates of primary production in the river plume, which contributes to the development of a large and recurring hypoxic area in summer. The mechanistic links between hypoxia and river discharge of freshwater and nutrients are complex as the accumulation and vertical export of organic matter, the establishment and maintenance of vertical stratification, and the microbial degradation of organic matter are controlled by a non-linear interplay of factors. We present results from a realistic, 3-dimensional, physical-biological model that includes the processes thought to be of first order importance to hypoxia formation and demonstrate that the model realistically reproduces many features of observed nitrate and phytoplankton dynamics including observed property distributions and rates. We then contrast the environmental factors and phytoplankton source and sink terms characteristic of three model subregions that represent an ecological gradient from eutrophic to oligotrophic conditions. We analyze specifically the reasons behind the counterintuitive observation that primary production in the light-limited plume region near the Mississippi River delta is positively correlated with river nutrient input. We find that, while primary production and phytoplankton biomass are positively correlated with nutrient load, phytoplankton growth rate is not. This suggests that accumulation of biomass in this region is not primarily controlled bottom up by nutrient-stimulation, but top down by systematic differences in the loss processes. We hypothesize that increased retention of river water in high discharge years explains this phenomenon. [ABSTRACT FROM AUTHOR]

Presents a response by S. B. Joye et al to the comments on their paper "Geohysical and Geochemical Signatures of Gulf of Mexico Seafloor Brines," which appeared in the 2005 issue of "Biogeosciences Discussions." Questions raised about the thermal stratification and chemical distribution data for the brine pools; Purpose of brine trapper positioning; Implication of the temperature profiles shown for the brine pool site.

This paper presents the first comprehensive comparison of several different dynamical-systems-based measures of stirring and Lagrangian coherence, computed from real ocean drifters. Seven commonly used methods (finite-time Lyapunov exponent (FTLE), trajectory path length, trajectory correlation dimension, trajectory encounter volume, Lagrangian-averaged vorticity deviation, dilation, and spectral clustering) were applied to 144 surface drifters in the Gulf of Mexico in order to map out the dominant Lagrangian coherent structures. Among the detected structures were regions of hyperbolic nature resembling stable manifolds from classical examples, divergent and convergent zones, and groups of drifters that moved more coherently and stayed closer together than the rest of the drifters. Many methods highlighted the same structures, but there were differences too. Overall, five out of seven methods provided useful information about the geometry of transport within the domain spanned by the drifters, whereas the path length and correlation dimension methods were less useful than others. [ABSTRACT FROM AUTHOR]

The ability to constrain the mechanisms that transport organic carbon into the deep ocean is complicated by the multiple physical, chemical, and ecological processes that intersect to create, transform, and transport particles in the ocean. In this paper we develop and parameterize a data-assimilative model of the multiple pathways of the biological carbon pump (NEMURO BCP). The mechanistic model is designed to represent sinking particle flux, active transport by vertically migrating zooplankton, and passive transport by subduction and vertical mixing, while also explicitly representing multiple biological and chemical properties measured directly in the field (including nutrients, phytoplankton and zooplankton taxa, carbon dioxide and oxygen, nitrogen isotopes, and 234 Thorium). Using 30 different data types (including standing stock and rate measurements related to nutrients, phytoplankton, zooplankton, and non-living organic matter) from Lagrangian experiments conducted on 11 cruises from four ocean regions, we conduct an objective statistical parameterization of the model and generate 1 million different potential parameter sets that are used for ensemble model simulations. The model simulates in situ parameters that were assimilated (net primary production and gravitational particle flux) and parameters that were withheld (234 Thorium and nitrogen isotopes) with reasonable accuracy. Model results show that gravitational flux of sinking particles and vertical mixing of organic matter from the euphotic zone are more important biological pump pathways than active transport by vertically migrating zooplankton. However, these processes are regionally variable, with sinking particles most important in oligotrophic areas of the Gulf of Mexico and California Current, sinking particles and vertical mixing roughly equivalent in productive coastal upwelling regions and the subtropical front in the Southern Ocean, and active transport an important contributor in the eastern tropical Pacific. We further find that mortality at depth is an important component of active transport when mesozooplankton biomass is high, but it is negligible in regions with low mesozooplankton biomass. Our results also highlight the high degree of uncertainty, particularly amongst mesozooplankton functional groups, that is derived from uncertainty in model parameters. Indeed, variability in BCP pathways between simulations for a specific location using different parameter sets (all with approximately equal misfit relative to observations) is comparable to variability in BCP pathways between regions. We discuss the implications of these results for other data-assimilation approaches and for studies that rely on non-ensemble model outputs. [ABSTRACT FROM AUTHOR]

During the last interglacial (LIG) the volume of additional water in the world's oceans was large enough to raise global sea levels about 6–9 m higher than present levels. However, LIG sea levels vary regionally and those regional differences hold clues about the past distribution of ice sheets and local rates of subsidence and tectonic uplift. In this study, I used a standardized database template to review and summarize the existing constraints on LIG sea levels across the northern Gulf of Mexico and Caribbean shoreline of the Yucatán Peninsula. In total, I extracted 32 sea-level indicators including the insertion of 16 U-series ages on corals, 1 electron spin resonance age, 2 amino acid racemization ages, and 26 luminescence ages. Most dated sea-level indicators for the northern Gulf of Mexico are based on optically stimulated luminescence (OSL) ages of beach deposits of a mappable LIG shoreline. This shoreline extends from the Florida Panhandle through south Texas but is buried or removed by the Mississippi River across most of Louisiana. A similar feature is observed in satellite images south of the Rio Grande within the Mexican portions of the Gulf of Mexico but has yet to be dated. Elevations measured on portions of this feature close to the modern coast point to sea levels less than 1 m to ∼5 m higher than present for much of the northern Gulf of Mexico. However, a few, albeit undated, portions of the same shoreline located at more inland locations point to sea levels up to +7.2 m, attesting to up to 7 m of differential subsidence between the inland and coastal sites. Across the Yucatán Peninsula, U-series dating of corals has provided the main index points for LIG sea levels. Other carbonate coastal features such as beach ridges and eolianites have also been described but rely on corals for their dating. The maximum elevation of the LIG coral-based relative sea-level (RSL) estimates decrease from around +6 m across the Caribbean shoreline of the Yucatán Peninsula near Cancún, Mexico, to as low as -6 m to the south beneath the southern atolls of Belize, although discussion continues as to the validity of the ages for these southern corals. If these lower-elevation corals are LIG in age, their below-present elevations may be a result of vertical motion along faults dipping into the Cayman Trough. South of Belize only one purported LIG coral has been dated on the Isla de Roatán off the coast of Honduras at a likely tectonically uplifted elevation of 37.2 m. Thus the elevation of LIG sea levels within the inland siliciclastic shorelines of Guatemala and Honduras as well as the southwestern Gulf of Mexico remains poorly constrained and a potential venue for future research. The database described in this paper is available open access in spreadsheet format as Simms (2020), at this link: 10.5281/zenodo.4556163. [ABSTRACT FROM AUTHOR]

The design of coastal protection measures and the quantification of coastal risks at locations affected by tropical cyclones (TCs) are often based solely on the analysis of historical cyclone tracks. Due to data scarcity and the random nature of TCs, the assumption that a hypothetical TC could hit a neighboring area with equal likelihood to past events can potentially lead to over- and/or underestimations of extremes and associated risks. The simulation of numerous synthetic TC tracks based on (historical) data can overcome this limitation. In this paper, a new method for the generation of synthetic TC tracks is proposed. The method has been implemented in the highly flexible open-source Tropical Cyclone Wind Statistical Estimation Tool (TCWiSE). TCWiSE uses an empirical track model based on Markov chains and can simulate thousands of synthetic TC tracks and wind fields in any oceanic basin based on any (historical) data source. Moreover, the tool can be used to determine the wind extremes, and the output can be used for the reliable assessment of coastal hazards. Validation results for the Gulf of Mexico show that TC patterns and extreme wind speeds are well reproduced by TCWiSE. [ABSTRACT FROM AUTHOR]

Comments on the paper "Geophysical and Geochemical Signatures of Gulf of Mexico Seafloor Brines" by S. B. Joye and colleagues. Comparison of the GC233 with the GB425 site; Background on the brine discharge process; Role of microbial activity in the geochemical budget.

A three-dimensional coupled hydrodynamic–biogeochemical model with multiple nutrient and plankton functional groups was developed and adapted to the Gulf of Mexico to investigate the role of nutrients and the complexity of plankton community in dissolved oxygen (DO) dynamics. A 15-year hindcast was achieved covering the period of 2006–2020. Extensive model validation against in situ data demonstrates that the model was capable of reproducing vertical distributions of DO and spatial distributions of bottom DO concentration, as well as their interannual variations. The study demonstrates that bottom DO dynamics and hypoxia evolution are significantly influenced by both physical processes and local biochemistry, with sedimentary oxygen consumption and vertical diffusion identified as key contributors. Summer hydrodynamics play a critical role in nutrient distribution and limitation: a notable expansion of Si limitation was simulated when coastal currents shifted eastward or northward. This effect, especially pronounced on the western part of the Louisiana–Texas shelf, underscores the importance of nutrient limitation in shaping DO dynamics. The model identifies a bi-peak primary production pattern in spring and early summer, aligned with satellite chlorophyll a variations, attributed to the complexity of the plankton community and interactions among different plankton groups. Our findings emphasize the necessity of integrating sophisticated plankton community dynamics into biogeochemical models to understand primary production variability and its impact on bottom hypoxia. [ABSTRACT FROM AUTHOR]

Presents comments on the study "Geophysical and Geochemical Signatures of Gulf of Mexico Seafloor Brines" by S. B. Joyce et al. Quality of the study; Critique of the authors' discussion of their observations; Limitations of the study.

Fires have great ecological, social, and economic impact. However, fire prediction and management remain challenges due to a limited understanding of their roles in the Earth system. Fires over southern Mexico and Central America (SMCA) are a good example of this, greatly impacting local air quality and regional climate. Here we report that the spring peak (April–May) of fire activities in this region has a distinct quasi-biennial signal based on multiple satellite datasets measuring different fire characteristics. The variability is initially driven by quasi-biennial variations in precipitation. Composite analysis indicates that strong fire years correspond to suppressed ascending motion and weakened precipitation over the SMCA. The anomalous precipitation over the SMCA is further found to be mostly related to the East Pacific–North Pacific (EP-NP) pattern 2 months prior to the fire season. The positive phase of the EP-NP leads to enhanced precipitation over the eastern US but suppressed precipitation over the SMCA, similar to the spatial pattern of precipitation differences between strong and weak fire years. Meanwhile, the quasi-biennial signals in precipitation and fires appear to be amplified by their interactions through a positive feedback loop at short timescales. Model simulations show that in strong fire years, more aerosol particles are released and transported downstream over the Gulf of Mexico and the eastern US, where suspended light-absorbing aerosols warm the atmosphere and cause the ascending motion of the air aloft. Subsequently, a compensating downward motion is formed over the region of the fire source and ultimately suppresses precipitation and intensifies fires. Statistical analysis shows the different durations of the two-way interaction, where the fire suppression effect of precipitation lasts for more than 20 d, while fire leads to a decrease in precipitation at shorter timescales (3–5 d). This study demonstrates the importance of fire–climate interactions in shaping the fire activities on an interannual scale and highlights how precipitation–fire interactions at short timescales contribute to the interannual variability in both fire and precipitation. [ABSTRACT FROM AUTHOR]

​​​​​​​The ability to detect and quantify methane emissions from offshore platforms is of considerable interest in providing actionable feedback to industrial operators. While satellites offer a distinctive advantage for remote sensing of offshore platforms which may otherwise be difficult to reach, offshore measurements of methane from satellite instruments in the shortwave infrared are challenging due to the low levels of diffuse sunlight reflected from water surfaces. Here, we use the GHGSat satellite constellation in a sun glint configuration to detect and quantify methane emissions from offshore targets around the world. We present a variety of examples of offshore methane plumes, including the largest single emission at (84 000 ± 24 000) kg h -1 observed by GHGSat from the Nord Stream 2 pipeline leak in 2022 and the smallest offshore emission measured from space at (180 ± 130) kg h -1 in the Gulf of Mexico. In addition, we provide an overview of the constellation's offshore measurement capabilities. We measure a median column precision of 2.1 % of the background methane column density and estimate a detection limit, from analytical modelling and orbital simulations, that varies between 160 and 600 kg h -1 depending on the latitude and season. [ABSTRACT FROM AUTHOR]

Atmospheric methane (CH 4) concentrations have more than doubled since the beginning of the industrial age, making CH 4 the second most important anthropogenic greenhouse gas after carbon dioxide (CO 2). The oil and gas sector represents one of the major anthropogenic CH 4 emitters as it is estimated to account for 22 % of global anthropogenic CH 4 emissions. An airborne field campaign was conducted in April–May 2019 to study CH 4 emissions from offshore gas facilities in the southern North Sea with the aim of deriving emission estimates using a top-down (measurement-led) approach. We present CH 4 fluxes for six UK and five Dutch offshore platforms or platform complexes using the well-established mass balance flux method. We identify specific gas production emissions and emission processes (venting and fugitive or flaring and combustion) using observations of co-emitted ethane (C 2 H 6) and CO 2. We compare our top-down estimated fluxes with a ship-based top-down study in the Dutch sector and with bottom-up estimates from a globally gridded annual inventory, UK national annual point-source inventories, and operator-based reporting for individual Dutch facilities. In this study, we find that all the inventories, except for the operator-based facility-level reporting, underestimate measured emissions, with the largest discrepancy observed with the globally gridded inventory. Individual facility reporting, as available for Dutch sites for the specific survey date, shows better agreement with our measurement-based estimates. For all the sampled Dutch installations together, we find that our estimated flux of (122.9 ± 36.8) kg h -1 deviates by a factor of 0.64 (0.33–12) from reported values (192.8 kg h -1). Comparisons with aircraft observations in two other offshore regions (the Norwegian Sea and the Gulf of Mexico) show that measured, absolute facility-level emission rates agree with the general distribution found in other offshore basins despite different production types (oil, gas) and gas production rates, which vary by 2 orders of magnitude. Therefore, mitigation is warranted equally across geographies. [ABSTRACT FROM AUTHOR]

The design of coastal protection measures and the quantification of coastal risks at locations affected by tropical cyclone (TC) are often based solely on the analysis of historical cyclone tracks. Due to data scarcity and the random nature of TCs, the assumption that future TCs could hit a neighboring area with equal likelihood than past events can potentially lead to over- and/or underestimations of extremes and associated risks. The simulation of numerous synthetic TC tracks based on (historical) data can overcome this limitation. In this paper, a new method for the generation of synthetic TC tracks is proposed. The method has been implemented in the highly flexible open-source Tropical Cyclone Wind Statistical Estimation Tool (TCWiSE). TCWiSE uses an Empirical Track Model based on Markov-chains and can simulate thousands of synthetic TC tracks and wind fields in any oceanic basin based on any (historical) data source. Moreover, the tool can be used to determine the wind extremes and the output can be used for the reliable assessment of coastal hazards. Validation results for the Gulf of Mexico show that TC patterns and extreme wind speeds are well reproduced by TCWiSE. [ABSTRACT FROM AUTHOR]

Intense tropical cyclones (TCs) can cause catastrophic damage to coastal regions after landfall. Recent studies have linked the devastation associated with TCs to climate change, which induces favorable conditions, such as increasing sea-surface temperature, to supercharge storms. Meanwhile, environmental factors, such as atmospheric aerosols, also impact the development and intensity of TCs, but their effects remain poorly understood, particularly coupled with ocean dynamics. Here, we quantitatively assess the aerosol microphysical effects and aerosol-modified ocean feedbacks during Hurricane Katrina using a cloud-resolving atmosphere–ocean coupled model: Weather Research and Forecasting (WRF) in conjunction with the Regional Ocean Model System (ROMS). Our model simulations reveal that an enhanced storm destructive power, as reflected by larger integrated kinetic energy, heavier precipitation, and higher sea-level rise, is linked to the combined effects of aerosols and ocean feedbacks. These effects further result in an expansion of the storm circulation with a reduced intensity because of a decreasing moist static energy supply and enhancing vorticity Rossby wave outward propagation. Both accumulated precipitation and storm surge are enhanced during the mature stage of the TC with elevated aerosol concentrations, implying exacerbated flood damage over the polluted coastal region. The ocean feedback following the aerosol microphysical effects tends to mitigate the vertical mixing cooling in the ocean mixing layer and offsets the aerosol-induced storm weakening by enhancing cloud and precipitation near the eyewall region. Our results highlight the importance of accounting for the effects of aerosol microphysics and ocean-coupling feedbacks to improve the forecast of TC destructiveness, particularly near the heavily polluted coastal regions along the Gulf of Mexico. [ABSTRACT FROM AUTHOR]

The chemical composition of foraminiferal shells is a well-known tool in paleoceanography to reconstruct past environments and climate. Their application is based on the relation between environmental variables and the concentration of elements incorporated or stable isotope fractionation during calcification. The vast majority of these so-called proxy relationships are based on the foraminiferal order of the Rotaliida, which, for example, encompasses all living planktonic species. However, there are more orders of foraminifera with calcifying members, some of which have fundamentally different biomineralization pathways, such as the Nodosariida, the Polymorphinida and the Vaginulinida. All these belong to the class of the Nodosariata and produce calcite shells, which may serve as carriers of paleoenvironmental and climate signals. The microstructures of these shells and overall morphology of these foraminifera strongly deviate from the Rotaliida, suggesting that their elemental and stable isotopic composition do not necessarily respond similarly to environmental parameters. A potential advantage of the Nodosariata is that they appear considerably earlier in the fossil record (Carboniferous) than the Rotaliida (Jurassic), thereby possibly extending the range of foraminifer-based paleoceanographic reconstructions considerably. To test the potential application of Nodosariata foraminifera as paleoproxies, we investigated incorporation of 5 elements in 11 species as a function of environmental parameters from a transect sampled in the Gulf of Mexico. Their element composition (B / Ca, Na / Ca, Mg / Ca, Sr / Ca and Ba / Ca) shows a distinct geochemical signature for these foraminifera, different to that of members of other foraminiferal orders. Results also show an increase in Mg / Ca values with increasing temperature, similar to that known for the Rotaliida, which suggest that Nodosariata shells might be useful for paleotemperature reconstructions. The difference in Mg / Ca–temperature calibration in Nodosariata compared to Rotaliida, with the large differences in their morphology, shell microstructures and overall geochemical composition, suggests that the Mg / Ca-to-temperature relationship is partly independent of the exact calcification mechanism. We compare Mg / Ca–temperature sensitivities across foraminiferal orders and describe a relationship between the average Mg / Ca and the sensitivity of the Mg / Ca–temperature calibration. For other elements, the variability across orders is smaller compared to that in Mg / Ca, which results in more similar El / Ca–environmental calibrations. [ABSTRACT FROM AUTHOR]

The Paleocene–Eocene Thermal Maximum (PETM) represents the most pronounced hyperthermal of the Cenozoic era and is hypothesized to have resulted in an intensification of the paleohydrologic cycle, including enhanced seasonality and increased sediment discharge to the coastal ocean. Although the PETM has been widely documented, there are few records from deposits that form the distal, deepwater components of large sediment-routing systems. This study presents new constraints on the stratigraphic placement of the PETM in the deepwater Gulf of Mexico basin through analysis of geochemical, carbon isotopic, and biostratigraphic data within a ∼124 m cored interval of the Wilcox Group. Biostratigraphic and carbon isotopic data indicate that the PETM extends over ∼13 m based on acmes in the dinoflagellate Apectodinium homomorphum and calcareous nannoplankton Rhomboaster cuspis as well as a ∼-2 ‰ shift in bulk organic δ13 C values. A decrease in bioturbation and benthic foraminifera suggests that a reduction in oxygen of Gulf of Mexico bottom waters and/or an increase in sedimentation rates were coincident with the onset of the PETM. A ∼2 m lag in the depositional record separates the onset of the PETM negative carbon isotope excursion (CIE) and deposition of a 5.7 m thick interval of organic-lean claystone and marlstone that reflects a shut-off of the supply of sand, silt, and terrestrial palynomorphs to the basin. We interpret deposits of the PETM in the deepwater Gulf of Mexico to reflect the combined effects of increased erosional denudation and rising sea level that resulted in sequestration of sand and silt near the coastline but that allowed delivery of terrigenous mud to the deep sea. The similarity of oceanographic changes observed in the Gulf of Mexico and Atlantic Ocean during the PETM supports the inference that these water masses were connected during latest Paleocene–earliest Eocene times. Although deposition of typical Wilcox Group facies resumed during and after the PETM recovery, an increased influx of terrestrial detritus (i.e., pollen, spores, terrestrial organic debris) relative to marine dinoflagellates is suggestive of long-lasting effects of the PETM. This study illustrates the profound and prolonged effects of climatic warming on even the most distal reaches of large (≥1×106 km 2) sediment-routing systems. [ABSTRACT FROM AUTHOR]

Both marine dinoflagellate cysts and terrestrially derived pollen and spores are abundant in coastal sediments close to river mouths, making sediment records from such settings ideal to simultaneously study land–ocean climate interactions, marine productivity patterns and freshwater input over time. However, few studies consider the combined calibration of these palynological proxies in modern coastal sediments offshore from rivers, which is needed to strengthen the interpretation of paleoreconstructions. Here, we analyze the palynological content of marine surface sediments along land–sea transects off the Mississippi and Atchafalaya river mouths in the northern Gulf of Mexico (GoM) and test three palynological indices which are often employed in the paleo-domain: (i) the abundance of cysts of heterotrophic and autotrophic dinoflagellates (dinocysts) as a tracer for primary productivity (H/A ratio) and (ii) the ratio between non-bisaccate pollen and bisaccates (P/B) as well as (iii) the ratio between pollen (excluding bisaccates) and dinocysts (P/D), which are both tracers for river input and distance to the coast. Our results show that dinoflagellate cysts are most abundant on the shelf, where heterotrophic dinocyst species dominate coastal assemblages in reach of the river plume, while autotrophic taxa are more present in the oligotrophic open ocean. This is clearly reflected in decreasing H/A values further offshore. Individual dinocyst taxa also seem to inhabit specific niches along an onshore–offshore transect, linked to nutrient availability and proximity to the turbid river plume. The highest pollen concentrations are found close to the Mississippi river mouth and mostly represent a mixture of local coastal and upstream vegetation, whereas bisaccate pollen was most abundant further offshore of the Mississippi river. Multivariate redundancy analysis (RDA) performed on both pollen and dinocyst assemblages, a set of environmental parameters, and the three palynological ratios showed that net primary productivity was the most important variable influencing the dinocyst assemblages, likely as the result of nutrient input. Additionally, the RDA confirmed that the H/A ratio indeed seems to track primary productivity, while the P/B ratio results in a robust indicator for distance to the coast, and the P/D ratio better reflects river input. Together, our data confirm and further specify the suitability of these three palynological ratios in river-dominated coastal margins as proxies for (past) marine productivity and distance to the coast and river. [ABSTRACT FROM AUTHOR]

This study describes and implements an integrated, multimedia, process-based system-level approach to estimating nitrogen (N) fate and transport in large river basins. The modeling system includes the following components: 1) Community Multi-Scale Air Quality (CMAQ); 2) Water Research and Forecasting (WRF); 3) Environmental Policy Integrated Climate (EPIC); and 4) Soil and Water Assessment Tool (SWAT). The previously developed Fertilizer Emission Scenario Tool for the Community Multiscale Air Quality (FEST-C) system integrated EPIC with the WRF model and CMAQ. FEST-C, driven by process-based WRF weather simulations, includes atmospheric N additions to agricultural cropland, and agricultural cropland contributions to ammonia emissions. Watershed hydrology and water quality models need to be integrated with the system (FEST-C), however, so it can be used in large river basins to address impacts of fertilization, meteorology, and atmospheric N deposition on water quality. Objectives of this paper are to describe how to expand the previous effort by integrating a watershed model with the FEST-C (CMAQ/WRF/EPIC) modeling system, as well as demonstrate application of the Integrated Modeling System (IMS) to the Mississippi River Basin (MRB) to simulate streamflow and dissolved N loadings to the Gulf of Mexico (GOM). IMS simulation results generally agree with USGS observations/estimations; the annual simulated streamflow is 218.9 mm and USGS observation is 211.1 mm and the annual simulated dissolved N is 2.1 kg/ha. and the USGS estimation is 2.8 kg/ha. Integrating SWAT with the CMAQ/WRF/EPIC modeling system allows for its use within large river basins without losing EPIC’s more detailed biogeochemistry processes, which will strengthen assessment of impacts of future climate scenarios, regulatory and voluntary programs for nitrogen oxide air emissions, and land use and land management on N transport and transformation in large river basins. [ABSTRACT FROM AUTHOR]

Rivers play a key role in the global carbon cycle by transporting terrestrial organic matter (TerrOM) from land to the ocean. Upon burial in marine sediments, this TerrOM may be a significant long-term carbon sink, depending on its composition and properties. However, much remains unknown about the dispersal of different types of TerrOM in the marine realm upon fluvial discharge since the commonly used bulk organic matter (OM) parameters do not reach the required level of source- and process-specific information. Here, we analyzed bulk OM properties, lipid biomarkers (long-chain n-alkanes , sterols, long-chain diols, alkenones, branched and isoprenoid glycerol dialkyl glycerol tetraethers (brGDGTs and isoGDGTs)), pollen, and dinoflagellate cysts in marine surface sediments along two transects offshore the Mississippi–Atchafalaya River (MAR) system, as well as one along the 20 m isobath in the direction of the river plume. We use these biomarkers and palynological proxies to identify the dispersal patterns of soil–microbial organic matter (SMOM), fluvial, higher plant, and marine-produced OM in the coastal sediments of the northern Gulf of Mexico (GoM). The Branched and Isoprenoid Tetraether (BIT) index and the relative abundance of C32 1,15-diols indicative for freshwater production show high contributions of SMOM and fluvial OM near the Mississippi River (MR) mouth (BIT = 0.6, FC321,15 > 50 %), which rapidly decrease further away from the river mouth (BIT < 0.1, FC321,15 < 20 %). In contrast, concentrations of long-chain n-alkanes and pollen grains do not show this stark decrease along the path of transport, and especially n-alkanes are also found in sediments in deeper waters. Proxy indicators show that marine productivity is highest close to shore and reveal that marine producers (diatoms, dinoflagellates, coccolithophores) have different spatial distributions, indicating their preferred niches. Close to the coast, where food supply is high and waters are turbid, cysts of heterotrophic dinoflagellates dominate the assemblages. The dominance of heterotrophic taxa in shelf waters in combination with the rapid decrease in the relative contribution of TerrOM towards the deeper ocean suggest that TerrOM input may trigger a priming effect that results in its rapid decomposition upon discharge. In the open ocean far away from the river plume, autotrophic dinoflagellates dominate the assemblages, indicating more oligotrophic conditions. Our combined lipid biomarker and palynology approach reveals that different types of TerrOM have distinct dispersal patterns, suggesting that the initial composition of this particulate OM influences the burial efficiency of TerrOM on the continental margin. [ABSTRACT FROM AUTHOR]

Low bottom water dissolved oxygen conditions (hypoxia) occur almost every summer in the northern Gulf of Mexico due to a combination of nutrient loadings and water column stratification. Several statistical and mechanistic models have been used to forecast the midsummer hypoxic area, based on spring nitrogen loading from major rivers. However, sub-seasonal forecasts are needed to fully characterize the dynamics of hypoxia over the summer season, which is important for informing fisheries and ecosystem management. Here, we present an approach to forecasting hypoxic conditions at a daily resolution through Bayesian mechanistic modeling that allows for rigorous uncertainty quantification. Within this framework, we develop and test different representations and projections of hydrometeorological model inputs. We find that May precipitation over the Mississippi River basin is a key predictor of summer discharge and loading that substantially improves forecast performance. Accounting for spring wind conditions also improves forecast performance, though to a lesser extent. The proposed approach generates forecasts for two different sections of the Louisiana–Texas shelf (east and west), and it explains about 50 % of the variability in the total hypoxic area when tested against historical observations (1985–2016). Results also show how forecast uncertainties build over the summer season, with longer lead times from the nominal forecast release date of 1 June, due to increasing stochasticity in riverine and meteorological inputs. Consequently, the portion of overall forecast variance associated with uncertainties in data inputs increases from 26 % to 41 % from June–July to August–September, respectively. Overall, the study demonstrates a unique approach to assessing and reducing uncertainties in temporally resolved hypoxia forecasting. [ABSTRACT FROM AUTHOR]

A long short-term memory (LSTM) neural network is proposed to predict hurricane-forced significant wave heights (SWHs) in the Caribbean Sea (CS) based on a dataset of 20 CS, Gulf of Mexico, and western Atlantic hurricane events collected from 10 buoys from 2010–2020. SWH nowcasting and forecasting are initiated using LSTM on 0, 3, 6, 9, and 12 h horizons. Through examining study cases Hurricanes Dorian (2019), Sandy (2012), and Igor (2010), results illustrate that the model is well suited to forecast hurricane-forced wave heights much more rapidly at a significantly cheaper computational cost compared to numerical wave models, with much less required expertise. Forecasts are highly accurate with regards to observations. For example, Hurricane Dorian nowcasts had correlation (R), root mean square error (RMSE), and mean absolute percentage error (MAPE) values of 0.99, 0.16 m, and 2.6 %, respectively. Similarly, on the 3, 6, 9, and 12 h forecasts, results produced R (RMSE; MAPE) values of 0.95 (0.51 m; 7.99 %), 0.92 (0.74 m; 10.83 %), 0.85 (1 m; 13.13 %), and 0.84 (1.24 m; 14.82 %), respectively. In general, the model can provide accurate predictions within 12 h (R≥0.8) and errors can be maintained at under 1 m within 6 h of forecast lead time. However, the model also consistently overpredicted the maximum observed SWHs. From a comparison of LSTM with a third-generation wave model, Simulating Waves Nearshore (SWAN), it was determined that when using Hurricane Dorian as a case example, nowcasts were far more accurate with regards to the observations. This demonstrates that LSTM can be used to supplement, but perhaps not replace, computationally expensive numerical wave models for forecasting extreme wave heights. As such, addressing the fundamental problem of phase shifting and other errors in LSTM or other data-driven forecasting should receive greater scrutiny from Small Island Developing States. To improve models results, additional research should be geared towards improving single-point LSTM neural network training datasets by considering hurricane track and identifying the hurricane quadrant in which buoy observations are made. [ABSTRACT FROM AUTHOR]

Flooding is of particular concern in low-lying coastal zones that are prone to flooding impacts from multiple drivers, such as oceanographic (storm surge and wave), fluvial (excessive river discharge), and/or pluvial (surface runoff). In this study, we analyse, for the first time, the compound flooding potential along the contiguous United States (CONUS) coastline from all flooding drivers, using observations and reanalysis data sets. We assess the overall dependence from observations by using Kendall's rank correlation coefficient (τ) and tail (extremal) dependence (χ). Geographically, we find the highest dependence between different drivers at locations in the Gulf of Mexico, southeastern, and southwestern coasts. Regarding different driver combinations, the highest dependence exists between surge–waves, followed by surge–precipitation, surge–discharge, waves–precipitation, and waves–discharge. We also perform a seasonal dependence analysis (tropical vs. extra-tropical season), where we find higher dependence between drivers during the tropical season along the Gulf and parts of the East Coast and stronger dependence during the extra-tropical season on the West Coast. Finally, we compare the dependence structure of different combinations of flooding drivers, using observations and reanalysis data, and use the Kullback–Leibler (KL) divergence to assess significance in the differences of the tail dependence structure. We find, for example, that models underestimate the tail dependence between surge–discharge on the East and West coasts and overestimate tail dependence between surge–precipitation on the East Coast, while they underestimate it on the West Coast. The comprehensive analysis presented here provides new insights on where the compound flooding potential is relatively higher, which variable combinations are most likely to lead to compounding effects, during which time of the year (tropical versus extra-tropical season) compound flooding is more likely to occur, and how well reanalysis data capture the dependence structure between the different flooding drivers. [ABSTRACT FROM AUTHOR]

Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability, such as the Younger Dryas and Dansgaard–Oeschger events, but this relationship is not straightforward. Large-scale glacial runoff events, such as Meltwater Pulse 1a (MWP1a), are not always temporally proximal to subsequent large-scale cooling. Moreover, the typical design of hosing experiments that support this relationship tends to artificially amplify the climate response. This study explores the impact that limitations in the representation of runoff in conventional "hosing" simulations has on our understanding of this relationship by examining where coastally released freshwater is transported when it reaches the ocean. We particularly focus on the impact of (1) the injection of freshwater directly over sites of deep-water formation (DWF) rather than at runoff locations (i.e. hosing), (2) excessive freshwater injection volumes (often by a factor of 5), and (3) the use of present-day (rather than palaeo) ocean gateways. We track the routing of glaciologically constrained freshwater volumes from four different inferred injection locations in a suite of eddy-permitting glacial ocean simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm) under both open and closed Bering Strait conditions. Restricting freshwater forcing values to realistic ranges results in less spreading of freshwater across the North Atlantic and indicates that the freshwater anomalies over DWF sites depend strongly on the geographical location of meltwater input. In particular, freshwater released into the Gulf of Mexico generates a very weak freshwater signal over DWF regions as a result of entrainment by the turbulent Gulf Stream. In contrast, freshwater released into the Arctic with an open Bering Strait or from the Eurasian ice sheet is found to generate the largest salinity anomalies over DWF regions in the North Atlantic and GIN (Greenland–Iceland–Norwegian) seas region respectively. Experiments show that when the Bering Strait is open, the Mackenzie River source exhibits more than twice as much freshening of the North Atlantic deep-water formation regions as when the Bering Strait is closed. Our results illustrate that applying freshwater hosing directly into the North Atlantic with even "realistic" freshwater amounts still overestimates the amount of terrestrial runoff reaching DWF regions. Given the simulated salinity anomaly distributions and the lack of reconstructed impact on deep-water formation during the Bølling–Allerød, our results support that the majority of the North American contribution to MWP1a was not routed through the Mackenzie River. [ABSTRACT FROM AUTHOR]

The coastal ocean is affected by an array of co-occurring biogeochemical and anthropogenic processes, resulting in substantial heterogeneity in water chemistry, including carbonate chemistry parameters such as pH and partial pressure of CO 2 (p CO 2). To better understand coastal and estuarine acidification and air-sea CO 2 fluxes, it is important to study baseline variability and driving factors of carbonate chemistry. Using both discrete bottle sample collection (2014–2020) and hourly sensor measurements (2016–2017), we explored temporal variability, from diel to interannual scales, in the carbonate system (specifically pH and p CO 2) at the Aransas Ship Channel located in the northwestern Gulf of Mexico. Using other co-located environmental sensors, we also explored the driving factors of that variability. Both sampling methods demonstrated significant seasonal variability at the location, with highest pH (lowest p CO 2) in the winter and lowest pH (highest p CO 2) in the summer. Significant diel variability was also evident from sensor data, but the time of day with elevated p CO 2 and depressed pH was not consistent across the entire monitoring period, sometimes reversing from what would be expected from a biological signal. Though seasonal and diel fluctuations were smaller than many other areas previously studied, carbonate chemistry parameters were among the most important environmental parameters for distinguishing between time of day and between seasons. It is evident that temperature, biological activity, freshwater inflow, and tide level (despite the small tidal range) are all important controls on the system, with different controls dominating at different timescales. The results suggest that the controlling factors of the carbonate system may not be exerted equally on both pH and p CO 2 on diel timescales, causing separation of their diel or tidal relationships during certain seasons. Despite known temporal variability on shorter timescales, discrete sampling was generally representative of the average carbonate system and average air-sea CO 2 flux on a seasonal and annual basis when compared with sensor data. [ABSTRACT FROM AUTHOR]

Given current threats to ocean ecosystem health, there is a growing demand for accurate biogeochemical hindcasts, nowcasts, and predictions. Provision of such products requires data assimilation, i.e., a comprehensive strategy for incorporating observations into biogeochemical models, but current data streams of biogeochemical observations are generally considered insufficient for the operational provision of such products. This study investigates to what degree the assimilation of satellite observations in combination with a priori model calibration by sparse BGC-Argo profiles can improve subsurface biogeochemical properties. The multivariate deterministic ensemble Kalman filter (DEnKF) has been implemented to assimilate physical and biological observations into a three-dimensional coupled physical–biogeochemical model, the biogeochemical component of which has been calibrated by BGC-Argo float data for the Gulf of Mexico. Specifically, observations of sea surface height, sea surface temperature, and surface chlorophyll were assimilated, and profiles of both physical and biological variables were updated based on the surface information. We assessed whether this leads to improved subsurface distributions, especially of biological properties, using observations from five BGC-Argo floats that were not assimilated. An alternative light parameterization that was tuned a priori using BGC-Argo observations was also applied to test the sensitivity of data assimilation impact on subsurface biological properties. Results show that assimilation of the satellite data improves model representation of major circulation features, which translate into improved three-dimensional distributions of temperature and salinity. The multivariate assimilation also improves the agreement of subsurface nitrate through its tight correlation with temperature, but the improvements in subsurface chlorophyll were modest initially due to suboptimal choices of the model's optical module. Repeating the assimilation run by using the alternative light parameterization greatly improved the subsurface distribution of chlorophyll. Therefore, even sparse BGC-Argo observations can provide substantial benefits for biogeochemical prediction by enabling a priori model tuning. Given that, so far, the abundance of BGC-Argo profiles in the Gulf of Mexico and elsewhere has been insufficient for sequential assimilation, updating 3D biological properties in a model that has been well calibrated is an intermediate step toward full assimilation of the new data types. [ABSTRACT FROM AUTHOR]

Seismic oceanography is a new cross-discipline between geophysics and oceanography that uses seismic reflection data to image and study the oceanic water column. Previous work on seismic oceanography was largely limited to two-dimensional (2D) seismic data and methods. Here we explore and quantify temporal and spatial variations in three-dimensional (3D) seismic oceanography to address whether 3D seismic imaging is meaningful in all directions and how one can take advantage of the variations. From a 3D multichannel seismic survey acquired for oil and gas exploration in the Gulf of Mexico over a 6-month period, a 3D oceanic seismic volume was derived. The 3D seismic images exhibit both temporal and spatial variations of the ocean, and theoretical and data analyses were used to quantify their contribution. Our results suggest that temporal variation is more prominent in the crossline direction than in the inline direction, causing discontinuities in crossline images. However, a series of 3D inline images can be seen as snapshots of the water column at different times, capturing temporal variation of thermohaline structures induced by ocean dynamics. Our findings suggest the potential uses of marine 3D seismic data in studying time-evolving mesoscale ocean dynamics. [ABSTRACT FROM AUTHOR]

A method for objectively extracting the displacement signals associated with coherent eddies from Lagrangian trajectories is presented, refined, and applied to a large dataset of 3770 surface drifters from the Gulf of Mexico. The method, wavelet ridge analysis, is a general method for the analysis of modulated oscillations, here modified to be more suitable to the eddy-detection problem. A means for formally assessing statistical significance is introduced, addressing the issue of false positives arising by chance from an unstructured turbulent background and opening the door to confident application of the method to very large datasets. Significance is measured through a frequency-dependent comparison with a stochastic dataset having statistical and spectral properties that match the original, but lacking organized oscillations due to eddies or waves. The application to the Gulf of Mexico reveals major asymmetries between cyclones and anticyclones, with anticyclones dominating at radii larger than about 50 km, but an unexpectedly rich population of highly nonlinear cyclones dominating at smaller radii. Both the method and the Gulf of Mexico eddy dataset are made freely available to the community for noncommercial use in future research. [ABSTRACT FROM AUTHOR]

A large set of historical surface drifter data from the Gulf of Mexico – 3770 trajectories spanning 28 years and more than a dozen data sources – are collected, uniformly processed and quality controlled, and assimilated into a spatially and temporally gridded dataset called GulfFlow. This dataset is available in two versions, with 1/4 ∘ or 1/12 ∘ spatial resolution respectively, both of which have overlapping monthly temporal bins with semimonthly spacing and which extend from the years 1992 through 2020. Together these form a significant resource for studying the circulation and variability in this important region. The uniformly processed historical drifter data from all publicly available sources, interpolated to hourly resolution, are also distributed in a separate product called GulfDriftersOpen. Forming a mean surface current map by directly bin-averaging the hourly drifter data is found to lead to severe artifacts, a consequence of the extremely inhomogeneous temporal distribution of the drifters. Averaging instead the already monthly-averaged data in GulfFlow avoids these problems, resulting in the highest-resolution map of the mean Gulf of Mexico surface currents yet produced. The consolidated drifter dataset is freely available at 10.5281/zenodo.3985916 , while the gridded products are available for noncommercial use only (for reasons discussed herein) at 10.5281/zenodo.3978793. [ABSTRACT FROM AUTHOR]

Agricultural watersheds are significant contributors to downstream nutrient excess issues. The timing and magnitude of nutrient mobilization in these watersheds are driven by a combination of anthropogenic, hydrologic, and biogeochemical factors that operate across a range of spatial and temporal scales. However, how, when, and where these complex factors drive nutrient mobilization has previously been difficult to capture with low-frequency or spatially limited data sets. To address this knowledge gap, we analyzed daily nitrate concentration (c) and discharge (Q) data for a 4-year period (2016–2019) from five nested, agricultural watersheds in the midwestern United States that contribute nutrient loads to the Gulf of Mexico. These records allow us to investigate nutrient mobilization patterns at a temporal and spatial resolution not previously possible. The watersheds span two distinct landforms shaped by differences in glacial history, resulting in natural soil properties that necessitated different drainage infrastructure across the study area. To investigate nutrient export patterns under different hydrologic conditions, we partitioned the hydrograph into stormflow and baseflow periods and examined those periods separately through the analysis of their concentration–discharge (c – Q) relationships on annual, seasonal, and event timescales. Stormflow showed consistent chemostatic patterns across all seasons, while baseflow showed seasonally dynamic c – Q patterns. Baseflow exhibited chemodynamic conditions in the summer and fall and more chemostatic conditions in the winter and spring, suggesting that water source contributions during baseflow were nonstationary. Baseflow chemodynamic behavior was driven by low-flow, low-NO 3- conditions during which in-stream and near-stream biological processing likely moderated in-stream NO 3- concentrations. Additionally, inputs from deeper groundwater with longer residence times and lower-NO 3- concentration likely contributed to low-NO 3- conditions in stream, particularly in the larger watersheds. Stormflow c – Q behavior was consistent across watersheds, but baseflow c – Q behavior was linked to the intensity of agriculture and the density of built drainage infrastructure, with more drainage infrastructure associated with higher loads and more chemostatic export patterns across the watersheds. This suggests that the way humans replumb the subsurface in response to geologic conditions has implications for hydrologic connectivity, homogenization of source areas, and, subsequently, nutrient export during both baseflow and stormflow. Our analysis also showed that anomalous flow periods greatly influenced overall c – Q patterns, suggesting that the analysis of high-resolution records at multiple scales is critical when interpreting seasonal or annual patterns. [ABSTRACT FROM AUTHOR]

The hypoxic zone in the northern Gulf of Mexico varies spatially (area, location) and temporally (onset, duration) on multiple scales. Exposure of fish to hypoxic dissolved oxygen (DO) concentrations (< 2 mg L -1) is often lethal and avoided, while exposure to 2 to 4 mg L -1 occurs readily and often causes the sublethal effects of decreased growth and fecundity for individuals of many species. We simulated the movement of individual fish within a high-resolution 3-D coupled hydrodynamic water quality model (FVCOM-WASP) configured for the northern Gulf of Mexico to examine how spatial variability in DO concentrations would affect fish exposure to hypoxic and sublethal DO concentrations. Eight static snapshots (spatial maps) of DO were selected from a 10 d FVCOM-WASP simulation that showed a range of spatial variation (degree of clumpiness) in sublethal DO for when total sublethal area was moderate (four maps) and for when total sublethal area was high (four maps). An additional case of allowing DO to vary in time (dynamic DO) was also included. All simulations were for 10 d and were performed for 2-D (bottom layer only) and 3-D (allows for vertical movement of fish) sets of maps. Fish movement was simulated every 15 min with each individual switching among three algorithms: tactical avoidance when exposure to hypoxic DO was imminent, strategic avoidance when exposure had occurred in the recent past, and default (independent of DO) when avoidance was not invoked. Cumulative exposure of individuals to hypoxia was higher under the high sublethal area snapshots compared to the moderate sublethal area snapshots but spatial variability in sublethal concentrations had little effect on hypoxia exposure. In contrast, relatively high exposures to sublethal DO concentrations occurred in all simulations. Spatial variability in sublethal DO had opposite effects on sublethal exposure between moderate and high sublethal area maps: the percentage of fish exposed to 2–3 mg L -1 decreased with increasing variability for high sublethal area but increased for moderate sublethal area. There was also a wide range of exposures among individuals within each simulation. These results suggest that averaging DO concentrations over spatial cells and time steps can result in underestimation of sublethal effects. Our methods and results can inform how movement is simulated in larger models that are critical for assessing how management actions to reduce nutrient loadings will affect fish populations. [ABSTRACT FROM AUTHOR]

We introduced a sediment-induced light attenuation algorithm into a biogeochemical model of the Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) modeling system. A fully coupled ocean–atmospheric–sediment–biogeochemical simulation was carried out to assess the impact of sediment-induced light attenuation on primary production in the northern Gulf of Mexico during the passage of Hurricane Gustav in 2008. When compared with model results without sediment-induced light attenuation, our new model showed a better agreement with satellite data on both the magnitude of nearshore chlorophyll concentration and the spatial distribution of offshore bloom. When Hurricane Gustav approached, resuspended sediment shifted the inner shelf ecosystem from a nutrient-limited one to a light-limited one. Only 1 week after Hurricane Gustav's landfall, accumulated nutrients and a favorable optical environment induced a posthurricane algal bloom in the top 20 m of the water column, while the productivity in the lower water column was still light-limited due to slow-settling sediment. Corresponding with the elevated offshore NO3 flux (38.71 mmol N m -1 s -1) and decreased chlorophyll flux (43.10 mg m -1 s -1), the outer shelf posthurricane bloom should have resulted from the cross-shelf nutrient supply instead of the lateral dispersed chlorophyll. Sensitivity tests indicated that sediment light attenuation efficiency affected primary production when sediment concentration was moderately high. Model uncertainties due to colored dissolved organic matter and parameterization of sediment-induced light attenuation are also discussed. [ABSTRACT FROM AUTHOR]

Oceanic primary production forms the basis of the marine food web and provides a pathway for carbon sequestration. Despite its importance, spatial and temporal variations of primary production are poorly observed, in large part because the traditional measurement techniques are laborious and require the presence of a ship. More efficient methods are emerging that take advantage of miniaturized sensors integrated into autonomous platforms such as gliders and profiling floats. One such method relies on determining the diurnal cycle of dissolved oxygen in the mixed layer and has been applied successfully to measurements from gliders and mixed-layer floats. This study is the first documented attempt to estimate primary production from diurnal oxygen changes measured by Argo-type profiling floats, thus accounting for the whole euphotic zone. We first present a novel method for correcting measurement errors that result from the relatively slow response time of the oxygen optode sensor. This correction relies on an in situ determination of the sensor's effective response time. The method is conceptually straightforward and requires only two minor adjustments in current Argo data transmission protocols: (1) transmission of measurement time stamps and (2) occasional transmission of downcasts in addition to upcasts. Next, we present oxygen profiles collected by 10 profiling floats in the northern Gulf of Mexico, evaluate whether community production and respiration can be detected, and show evidence of internal oscillations influencing the diurnal oxygen signal. Our results show that profiling floats are capable of measuring diurnal oxygen variations although the confounding influence of physical processes does not permit a reliable estimation of biological rates in our dataset. We offer suggestions for recognizing and removing the confounding signals. [ABSTRACT FROM AUTHOR]

Biogeochemical ocean models are useful tools but subject to uncertainties arising from simplifications, inaccurate parameterization of processes, and poorly known model parameters. Parameter optimization is a standard method for addressing the latter but typically cannot constrain all biogeochemical parameters because of insufficient observations. Here we assess the trade-offs between satellite observations of ocean color and biogeochemical (BGC) Argo profiles and the benefits of combining both observation types for optimizing biogeochemical parameters in a model of the Gulf of Mexico. A suite of optimization experiments is carried out using different combinations of satellite chlorophyll and profile measurements of chlorophyll, phytoplankton biomass, and particulate organic carbon (POC) from autonomous floats. As parameter optimization in 3D models is computationally expensive, we optimize the parameters in a 1D model version and then perform 3D simulations using these parameters. We show first that the use of optimal 1D parameters, with a few modifications, improves the skill of the 3D model. Parameters that are only optimized with respect to surface chlorophyll cannot reproduce subsurface distributions of biological fields. Adding profiles of chlorophyll in the parameter optimization yields significant improvements for surface and subsurface chlorophyll but does not accurately capture subsurface phytoplankton and POC distributions because the parameter for the maximum ratio of chlorophyll to phytoplankton carbon is not well constrained in that case. Using all available observations leads to significant improvements of both observed (chlorophyll, phytoplankton, and POC) and unobserved (e.g., primary production) variables. Our results highlight the significant benefits of BGC-Argo measurements for biogeochemical parameter optimization and model calibration. [ABSTRACT FROM AUTHOR]

Uncertainties in carbon chemistry variability still remain large in the Gulf of Mexico (GoM), as data gaps limit our ability to infer basin-wide patterns. Here we configure and validate a regional high-resolution ocean biogeochemical model for the GoM to describe seasonal patterns in surface pressure of CO2 (pCO2), aragonite saturation state (ΩAr), and sea–air CO2 flux. Model results indicate that seasonal changes in surface pCO2 are strongly controlled by temperature across most of the GoM basin, except in the vicinity of the Mississippi–Atchafalaya river system delta, where runoff largely controls dissolved inorganic carbon (DIC) and total alkalinity (TA) changes. Our model results also show that seasonal patterns of surface ΩAr are driven by seasonal changes in DIC and TA, and reinforced by the seasonal changes in temperature. Simulated sea–air CO2 fluxes are consistent with previous observation-based estimates that show CO2 uptake during winter–spring, and CO2 outgassing during summer–fall. Annually, our model indicates a basin-wide mean CO2 uptake of 0.35 molm-2yr-1 , and a northern GoM shelf (< 200 m) uptake of 0.93 molm-2yr-1. The observation and model-derived patterns of surface pCO2 and CO2 fluxes show good correspondence; thus this study contributes to improved constraints of the carbon budget in the region. [ABSTRACT FROM AUTHOR]

The coastal Gulf of Mexico (GOM) and coastal sea off the Korean Peninsula (CSK) both suffer from human-induced eutrophication. We used a nitrogen (N) mass balance model in two different regions with different nitrogen input sources to estimate organic carbon fluxes and predict future carbon fluxes under different model scenarios. The coastal GOM receives nitrogen predominantly from the Mississippi and Atchafalaya rivers and atmospheric nitrogen deposition is only a minor component in this region. In the CSK, groundwater and atmospheric nitrogen deposition are more important controlling factors. Our model includes the fluxes of nitrogen to the ocean from the atmosphere, groundwater and rivers, based on observational and literature data, and identifies three zones (brown, green and blue waters) in the coastal GOM and CSK with different productivity and carbon fluxes. Based on our model results, the potential primary production rate in the inner (brown water) zone are over 2 gC m -2 d -1 (GOM) and 1.5 gC m -2 d -1 (CSK). In the middle (green water) zone, potential production is from 0.1 to 2 (GOM) and 0.3 to 1.5 gC m -2 d -1 (CSK). In the offshore (blue water) zone, productivity is less than 0.1 (GOM) and 0.3 (CSK) gC m -2 d -1. Through our model scenario results, overall oxygen demand in the GOM will increase approximately 21 % if we fail to reduce riverine N input, likely increasing considerably the area affected by hypoxia. Comparing the results from the USA with those from the Korean Peninsula shows the importance of considering both riverine and atmospheric inputs of nitrogen. This has direct implications for investigating how changes in energy technologies can lead to changes in the production of various atmospheric contaminants that affect air quality, climate and the health of local populations. [ABSTRACT FROM AUTHOR]

Assessments of ocean data assimilation (DA) systems and observing system design experiments typically rely on identical or nonidentical twin experiments. The identical twin approach has been recognized as yielding biased impact assessments in atmospheric predictions, but these shortcomings are not sufficiently appreciated for oceanic DA applications. Here we present the first direct comparison of the nonidentical and identical twin approaches in an ocean DA application. We assess the assimilation impact for both approaches in a DA system for the Gulf of Mexico that uses the ensemble Kalman filter. Our comparisons show that, despite a reasonable error growth rate in both approaches, the identical twin produces a biased skill assessment, overestimating the improvement from assimilating sea surface height and sea surface temperature observations while underestimating the value of assimilating temperature and salinity profiles. Such biases can lead to an undervaluation of some observing assets (in this case profilers) and thus a misguided distribution of observing system investments. [ABSTRACT FROM AUTHOR]

Very-near-surface ocean currents are dominated by wind and wave forcing and have large impacts on the transport of buoyant materials in the ocean. Surface currents, however, are under-resolved in most operational ocean models due to the difficultly of measuring ocean currents close to, or directly at, the air–sea interface with many modern instrumentations. Here, observations of ocean currents at two depths within the first meter of the surface are made utilizing trajectory data from both drogued and undrogued Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) drifters, which have draft depths of 60 and 5 cm, respectively. Trajectory data of dense, colocated drogued and undrogued drifters were collected during the Lagrangian Submesoscale Experiment (LASER) that took place from January to March of 2016 in the northern Gulf of Mexico. Examination of the drifter data reveals that the drifter velocities become strongly wind- and wave-driven during periods of high wind, with the pre-existing regional circulation having a smaller, but non-negligible, influence on the total drifter velocities. During these high wind events, we deconstruct the total drifter velocities of each drifter type into their wind- and wave-driven components after subtracting an estimate for the regional circulation, which pre-exists each wind event. In order to capture the regional circulation in the absence of strong wind and wave forcing, a Lagrangian variational method is used to create hourly velocity field estimates for both drifter types separately, during the hours preceding each high wind event. Synoptic wind and wave output data from the Unified Wave INterface-Coupled Model (UWIN-CM), a fully coupled atmosphere, wave and ocean circulation model, are used for analysis. The wind-driven component of the drifter velocities exhibits a rotation to the right with depth between the velocities measured by undrogued and drogued drifters. We find that the average wind-driven velocity of undrogued drifters (drogued drifters) is ∼3.4 %–6.0 % (∼2.3 %–4.1 %) of the wind speed and is deflected ∼5 –55 ∘ (∼30 –85 ∘) to the right of the wind, reaching higher deflection angles at higher wind speeds. Results provide new insight on the vertical shear present in wind-driven surface currents under high winds, which have vital implications for any surface transport problem. [ABSTRACT FROM AUTHOR]

Net community production (NCP) in the surface water of the northern Gulf of Mexico (nGOM) and its coupling with the CO2 system were examined during the productive spring season. NCP was estimated using multiple approaches: (1) underway O2 and Ar ratio, (2) oxygen changes during light/dark bottle oxygen incubations, and (3) non-conservative changes in dissolved inorganic carbon or nutrients. These methods all showed high spatial variability of NCP and displayed similar patterns along the river–ocean mixing gradient, showing high production rates in plume regions. NCP O2Ar estimated from high-resolution O2 and Ar underway measurement indicated heterotrophic conditions at the high-nutrient and high-turbidity Mississippi River end (-51.3±11.9 mmol C m -2 d -1 when salinity < 2) resulting from the influence of terrestrial carbon input and light limitation on photosynthesis. High NCP O2Ar rates (105.0±59.2 mmol C m -2 d -1 , up to 235.4 mmol C m -2 d -1) were observed in the Mississippi and Atchafalaya plumes at intermediate salinities between 15 and 30 where light and nutrients were both favorable for phytoplankton production. NCP O2Ar rates observed in the high-salinity, oligotrophic offshore waters (salinity > 35.5) were close to zero due to nutrient limitation. Air–sea CO2 fluxes generally showed corresponding changes, from being a strong CO2 source in the river channel (55.5±7.6 mmol C m -2 d -1), to a CO2 sink in the plume (-13.4±5.5 mmol C m -2 d -1), and to being nearly in equilibrium with the atmosphere in offshore waters. Overall, the surface water of the nGOM was net autotrophic during spring 2017, with an area-weighted mean NCP O2Ar of 21.2 mmol C m -2 d -1 , and was a CO2 sink of - 6.7 mmol C m -2 d -1. A temporal mismatch between in situ biological production and gas exchange of O2 and CO2 was shown through a box model to result in decoupling between NCP O2Ar and CO2 flux (e.g., autotrophic water as a CO2 source outside the Mississippi River mouth and heterotopic water as a CO2 sink in the Atchafalaya coastal water). This decoupling was a result of in situ biological production superimposed on the lingering background pCO2 from the source water because of the slow air–sea CO2 exchange rate and the buffering effect of the carbonate system. [ABSTRACT FROM AUTHOR]

Tehuantepecers or Tehuanos are extreme winds produced in the Isthmus of Tehuantepec, blowing south through Chivela Pass, the mountain gap across the isthmus, from the Gulf of Mexico into the Pacific Ocean. They are the result of the complex interaction between large-scale meteorological conditions and local orographic forcings around Chivela Pass, and occur mainly in winter months due to cold air damming in the wake of cold fronts that reach as far south as the Bay of Campeche. Even though the name refers mostly to the intense mountain gap outflow, Tehuantepecer episodes can also generate other localized extreme wind conditions across the region, such as downslope windstorms and hydraulic jumps, which are strong turbulent flows that have a direct effect on the Pacific side of the isthmus and the Gulf of Tehuantepec. This study focuses on investigating these phenomena using high horizontal and vertical resolution WRF (Weather Research and Forecasting) model simulations. In particular, we employ a 4-nested grid configuration with up to 444 m horizontal spacing in the innermost domain and 70 hybrid-sigma vertical levels, 8 of which lie within the first 200 m above ground. We select one 36 h period in December 2013, when favorable conditions for a strong gap wind situation were observed. The high-resolution WRF experiment reveals a significant fine-scale structure in the strong Tehuano wind flow, beyond the well known gap jet. Depending on the Froude number upstream of the topographic barrier, different downslope windstorm conditions and hydraulic jumps with rotor circulations develop simultaneously at different locations east of Chivela Pass with varied crest height. A comparison with observations suggests that the model accurately represents the spatially heterogeneous intense downslope windstorm and the formation of mountain wave clouds for several hours, with low errors in wind speed, wind direction and temperature. [ABSTRACT FROM AUTHOR]

A 3-D unstructured-grid hydrodynamic model for the northern Gulf of Mexico was developed, with a hybrid s – z vertical grid and high-resolution horizontal grid for the main estuarine systems along the Texas–Louisiana coast. This model, based on the Semi-implicit Cross-scale Hydroscience Integrated System Model (SCHISM), is driven by the observed river discharge, reanalysis atmospheric forcing, and open boundary conditions from global HYCOM output. The model reproduces the temporal and spatial variation of observed water level, salinity, temperature, and current velocity in Galveston Bay and on the shelf. The validated model was applied to examine the remote influence of neighboring large rivers, specifically the Mississippi–Atchafalaya River (MAR) system, on salinity, stratification, vertical mixing, and longshore transport along the Texas coast. Numerical experiments reveal that the MAR discharge could significantly decrease the salinity and change the stratification and vertical mixing on the inner Texas shelf. It would take about 25 and 50 d for the MAR discharge to reach the mouth of Galveston Bay and Port Aransas, respectively. The influence of the MAR discharge is sensitive to the wind field. Winter wind constrains the MAR freshwater to form a narrow lower-salinity band against the shore from the Mississippi Delta all the way to the southwestern Texas coast, while summer wind reduces the downcoast longshore transport significantly, weakening the influence of the MAR discharge on surface salinity along Texas coast. However, summer wind causes a much stronger stratification on the Texas shelf, leading to a weaker vertical mixing. The decrease in salinity of up to 10 psu at the mouth of Galveston Bay due to the MAR discharge results in a decrease in horizontal density gradient, a decrease in the salt flux, and a weakened estuarine circulation and estuarine–ocean exchange. We highlight the flexibility of the model and its capability to simulate not only estuarine dynamics and shelf-wide transport, but also the interactions between them. [ABSTRACT FROM AUTHOR]

Atmospheric aerosol particles that can nucleate ice are referred to as ice-nucleating particles (INPs). Recent studies have confirmed that aerosol particles emitted by the oceans can act as INPs. This very relevant information can be included in climate and weather models to predict the formation of ice in clouds, given that most of them do not consider oceans as a source of INPs. Very few studies that sample INPs have been carried out in tropical latitudes, and there is a need to evaluate their availability to understand the potential role that marine aerosol may play in the hydrological cycle of tropical regions. This study presents results from the first measurements obtained during a field campaign conducted in the tropical village of Sisal, located on the coast of the Gulf of Mexico of the Yucatan Peninsula in Mexico in January–February 2017, and one of the few data sets currently available at such latitudes (i.e., 21 ∘ N). Aerosol particles sampled in Sisal are shown to be very efficient INPs in the immersion freezing mode, with onset freezing temperatures in some cases as high as -3 ∘ C, similarly to the onset temperature from Pseudomonas syringae. The results show that the INP concentration in Sisal was higher than at other locations sampled with the same type of INP counter. Air masses arriving in Sisal after the passage of cold fronts have surprisingly higher INP concentrations than the campaign average, despite their lower total aerosol concentration. The high concentrations of INPs at warmer ice nucleation temperatures (T>-15 ∘ C) and the supermicron size of the INPs suggest that biological particles may have been a significant contributor to the INP population in Sisal during this study. However, our observations also suggest that at temperatures ranging between -20 and -30 ∘ C mineral dust particles are the likely source of the measured INPs. [ABSTRACT FROM AUTHOR]

This work examines the origin of atmospheric water vapor arriving to the western North American monsoon (WNAM) region over a 34-year period (1981–2014) using a Lagrangian approach. This methodology computes budgets of evaporation minus precipitation (E-P) by calculating changes in the specific humidity of thousands of air particles advected into the study area by the observed winds. The length of the period analyzed (34 years) allows the method to identify oceanic and terrestrial sources of moisture to the WNAM region from a climatological perspective. During the wet season, the WNAM region itself is on average the main evaporative source, followed by the Gulf of California. However, water vapor originating from the Caribbean Sea, the Gulf of Mexico, and terrestrial eastern Mexico is found to influence regional-scale rainfall generation. Enhanced (reduced) moisture transport from the Caribbean Sea and the Gulf of Mexico from 4 to 6 days before precipitation events seems to be responsible for increased (decreased) rainfall intensity on regional scales during the monsoon peak. Westward propagating mid- to upper-level inverted troughs (IVs) seem to favor these water vapor fluxes from the east. In particular, a 200 % increase in the moisture flux from the Caribbean Sea to the WNAM region is found to be followed by the occurrence of heavy precipitation in the WNAM area a few days later. Low-level troughs off the coast of northwestern Mexico and upper-level IVs over the Gulf of Mexico are also related to these extreme rainfall events. [ABSTRACT FROM AUTHOR]

Low-level jets (LLJs) can be defined as wind corridors of anomalously high wind speed values located within the first kilometre of the troposphere. These structures are one of the major meteorological systems in the meridional transport of moisture on a global scale. In this work, we focus on the southerly Great Plains low-level jet, which plays an important role in the moisture transport balance over the central United States. The Gulf of Mexico is the main moisture source for the Great Plains low-level jet (GPLLJ), which has been identified as a key factor for rainfall modulation over the eastern and central US. The relationship between moisture transport from the Gulf of Mexico to the Great Plains and precipitation has been well documented in previous studies. Nevertheless, a large uncertainty still remains in the quantification of the moisture amount actually carried by the GPLLJ. The main goal of this work is to address this question. For this purpose, a relatively new tool, the regional atmospheric Weather Research and Forecasting Model with 3-D water vapour tracers (WRF-WVT; Insua-Costa and Miguez-Macho, 2018) is used together with the Lagrangian model FLEXPART to estimate the load of precipitable water advected within the GPLLJ. Both models were fed with data from ERA Interim. From a climatology of jet intensity over a 37-year period, which follows a Gaussian distribution, we select five cases for study, representing the mean and 1 and 2 standard deviations above and below it. Results show that the jet is responsible for roughly 70 %–80 % of the moisture transport occurring in the southern Great Plains when a jet event occurs. Furthermore, moisture transport by the GPLLJ extends to the north-east US, accounting for 50 % of the total in areas near the Great Lakes. Vertical distributions show the maximum of moisture advected by the GPLLJ at surface levels and maximum values of moisture flux about 500 m above, in coincidence with the wind speed profile. [ABSTRACT FROM AUTHOR]

The beach-ridge sequence of the Usumacinta–Grijalva delta borders a 300 km long section of the southern Gulf of Mexico coast. With around 500 beach ridges formed in the last 6500 years, the sequence is unsurpassed in the world in terms of numbers of individual ridges preserved, continuity of the record, and temporal resolution. We mapped and dated the most extensively accreted part of the sequence, linking six phases of accretion to river mouth reconfigurations and constraining their ages with 14C and OSL dating. The geomorphological and sedimentological reconstruction relied on lidar data, coring transects, GPR measurements, grain-size analyses, and chemical fingerprinting of volcanic glass and pumice encountered within the beach and dune deposits. We demonstrate that the beach-ridge complex was formed under ample long-term fluvial sediment supply and shorter-term wave- and aeolian-modulated sediment reworking. The abundance of fluvially supplied sand is explained by the presence of easily weatherable Los Chocoyos ignimbrites from the ca. 84 ka eruption of the Atitlán volcano (Guatemala) in the catchment of the Usumacinta River. Autocyclic processes seem responsible for the formation of ridge–swale couplets. Fluctuations in their periodicity (ranging from 6–19 years) are governed by progradation rate, and are therefore not indicative of sea level fluctuations or variability in storm activity. The fine sandy beach ridges are mainly swash built. Ridge elevation, however, is strongly influenced by aeolian accretion during the time the ridge is located next to the beach. Beach-ridge elevation is negatively correlated with progradation rate, which we relate to the variability in sediment supply to the coastal zone, reflecting decadal-scale precipitation changes within the river catchment. In the southern Mexican delta plain, the coastal beach ridges therefore appear to be excellent recorders of hinterland precipitation. [ABSTRACT FROM AUTHOR]

This study investigates current ground deformation derived from the GPS geodesy infrastructure in the Gulf of Mexico region. The positions and velocity vectors of 161 continuous GPS (CGPS) stations are presented with respect to a newly established local reference frame, the Stable Gulf of Mexico Reference Frame (SGOMRF). Thirteen long-term (>5 years) CGPS are used to realize the local reference frame. The root mean square (RMS) of the velocities of the 13 SGOMRF reference stations achieves 0.2mmyr-1 in the horizontal and 0.3mmyr-1 in the vertical directions. GPS observations presented in this study indicate significant land subsidence in the coastal area of southeastern Louisiana, the greater Houston metropolitan area, and two cities in Mexico (Aguascalientes and Mexico City). The most rapid subsidence is recorded at the Mexico City International airport, which is up to 26.6 cm yr-1 (2008-2014). Significant spatial variation of subsidence rates is observed in both Mexico City and the Houston area. The overall subsidence rate in the Houston area is decreasing. The subsidence rate in southeastern Louisiana is relatively smaller (4.0-6.0mmyr-1) but tends to be steady over time. This poses a potential threat to the safety of coastal infrastructure in the long-term. [ABSTRACT FROM AUTHOR]

The Paleocene-Eocene Thermal Maximum (PETM, ~56 Ma) was a ~200 kyr episode of global warming, associated with massive injections of 13C-depleted carbon into the ocean–atmosphere system. Although climate change during the PETM is relatively well constrained, effects on marine oxygen concentrations and nutrient cycling remain largely unclear. We identify the PETM in a sediment core from the US margin of the Gulf of Mexico. Biomarker-based paleotemperature proxies (methylation of branched tetraether-cyclization of branched tetraether (MBT-CBT) and TEX86) indicate that continental air and sea surface temperatures warmed from 27-29 to ~35 °C, although variations in the relative abundances of terrestrial and marine biomarkers may have influenced these estimates. Vegetation changes, as recorded from pollen assemblages, support this warming. The PETM is bracketed by two unconformities. It overlies Paleocene silt- and mudstones and is rich in angular (thus in situ produced; autochthonous) glauconite grains, which indicate sedimentary condensation. A drop in the relative abundance of terrestrial organic matter and changes in the dinoflagellate cyst assemblages suggest that rising sea level shifted the deposition of terrigenous material landward. This is consistent with previous findings of eustatic sea level rise during the PETM. Regionally, the attribution of the glauconite-rich unit to the PETM implicates the dating of a primate fossil, argued to represent the oldest North American specimen on record. The biomarker isorenieratene within the PETM indicates that euxinic photic zone conditions developed, likely seasonally, along the Gulf Coastal Plain. A global data compilation indicates that O2 concentrations dropped in all ocean basins in response to warming, hydrological change, and carbon cycle feedbacks. This culminated in (seasonal) anoxia along many continental margins, analogous to modern trends. Seafloor deoxygenation and widespread (seasonal) anoxia likely caused phosphorus regeneration from suboxic and anoxic sediments. We argue that this fueled shelf eutrophication, as widely recorded from microfossil studies, increasing organic carbon burial along many continental margins as a negative feedback to carbon input and global warming. If properly quantified with future work, the PETM offers the opportunity to assess the biogeochemical effects of enhanced phosphorus regeneration, as well as the timescales on which this feedback operates in view of modern and future ocean deoxygenation. [ABSTRACT FROM AUTHOR]