Original Research ARTICLE
Scenarios of Deoxygenation of the Eastern Tropical North Pacific During the Past Millennium as a Window Into the Future of Oxygen Minimum Zones
- 1Department of Earth and Planetary Sciences, University of California, Riverside, Riverside, CA, United States
- 2Instituto de Geofisica, Universidad Nacional Autonóma de México, Mexico City, Mexico
- 3Department of Environmental Sciences, University of California, Riverside, Riverside, CA, United States
Diverse studies predict global expansion of Oxygen Minimum Zones (OMZs) as a consequence of anthropogenic global warming. While the observed dissolved oxygen concentrations in many coastal regions are slowly decreasing, sediment core paleorecords often show contradictory trends. This is the case for numerous high-resolution reconstructions of oxygenation in the Eastern Tropical North Pacific (ETNP). While major shifts in redox conditions of the ETNP are dominated by glacial-interglacial cycling, important fluctuations also occur in response to minor climatic and oceanographic perturbations. It is important to understand these scenarios of past redox variation, as they are the closest analog for near future climate and oceanographic change. We present recently collected sediment core proxy records from the Gulf of California in which we reproduce the variability of productivity and oxygenation of the ETNP OMZ during the past millennium. We emphasize paleoproductivity (Corg, Ni, Baexcess) and paleoredox indicators (Mo, Cd, V, Uauth) in sediment cores collected in Alfonso and La Paz basins and compare these OMZ records with other archives of the Eastern Pacific. Our findings indicate that the OMZ expanded in response to increased upwelling and productivity during cold intervals of the early 1400s, early 1500s, late 1600s, and early 1800s AD (evidenced by higher Ni, V, Cd, Mo, and Uauth). The most hypoxic times corresponded to the beginning of the Little Ice Age (expressed in elevated Mo). Significant OMZ contractions occurred around late 1300s, early 1700s, and late 1900s AD after reoxygenation events that were instigated by low productivity (lower Ni, V, Cd, Mo, and Uauth). The mechanisms that control decadal-to-centennial oxygen variability in the ETNP remain unidentified but are likely influenced by solar forcing not only driving migrations of the Intertropical Convergence Zone (ITCZ) but more importantly changes in the intensity of the Pacific Walker Circulation (PWC). During the Little Ice Age solar irradiance was at its lowest for the past millennium, which strengthened the PWC. This pattern contributed to more frequent La Niña-like conditions, which enhanced upwelling of nutrient-rich waters in the west coast of North America, driving productivity and reducing bottom oxygen levels, as seen in our ETNP records.
Numerous experimental, observational, and numerical model studies predict global expansion of Oxygen Minimum Zones (OMZ) linked to human induced climate change (Keeling et al., 2010; Gilly et al., 2013; Howes et al., 2015; Praetorius et al., 2015; Levin, 2018). The Eastern Tropical North Pacific (ETNP) hosts the largest OMZ in the world, and its economic importance for the coastal regions of North and South America is indisputable (Levin, 2018). Numerous in-situ measurements during recurring oceanographic campaigns (e.g., CalCOFI, IMECOCAL; Lynn and Simpson, 1987; Bograd et al., 2003; Durazo-Arvizu and Gaxiola-Castro, 2010) revealed declining O2 trends during the last few decades (Bograd et al., 2003, 2008; Lavín et al., 2013; Breitburg et al., 2018; Levin, 2018). It is important to note that some of the decreasing patterns could be biased by improvements in O2 measurement technologies that have decreased minimum detection thresholds. Despite this caveat, high-resolution reconstructions of changes in dissolved oxygen in the ETNP OMZ over the past century yield contrasting observations, with warming and increasing extent of anoxia (<0.1 mL/L O2) not always correlated (Choumiline, 2011; Deutsch et al., 2014; Tems et al., 2016; Fu et al., 2018; Ontiveros-Cuadras et al., 2019). Ongoing discussions about the causes of OMZ intensification over the past decades typically invoke decreased O2 solubility due to warming of the upper ocean, increased stratification slowing ventilation (Oschlies et al., 2008; Shaffer et al., 2009; Moffitt et al., 2015; Praetorius et al., 2015; Fu et al., 2018), and the strong effect of enhanced marine productivity on the formation of settling organic matter that undergoes decay and consumes more O2 during remineralization (Wright et al., 2012; Levin, 2018). The latter relationship is the main mechanism by which some OMZ basins in the Eastern Pacific become anoxic (Moffitt et al., 2015; Breitburg et al., 2018; Levin, 2018). However, what causes increases in primary productivity is still debated, with initial results pointing to a complex interaction of factors involving increased upwelling rates and terrigenous discharge—both of which can provide the necessary nutrients for proliferation of primary producers (Howes et al., 2015).
Because the first analytical measurements of dissolved oxygen in seawater started in the late 1800s (Carpenter, 1965), data are not available for preceding portions of the nineteenth century. Nevertheless, multiple applications of well-established paleoceanographic proxies allow us to infer changes in earlier productivity and oxygenation by reconstructing past conditions preserved in the elemental composition of marine sedimentary records emphasizing organic carbon (Corg), Ba, Cd, Mo, Mn, Ni, U, and V (Calvert and Pedersen, 1993; Helz et al., 1996; Morford and Emerson, 1999; Zheng et al., 2000; Lyons and Severmann, 2006; Tribovillard et al., 2006; Algeo and Rowe, 2012; Schoepfer et al., 2015). Multiple lines of proxy evidence show that major oxygenation shifts in the ETNP are coupled to glacial-interglacial cycling, often leading to a decrease in dissolved O2 concentrations in North Pacific Intermediate and North Pacific Deep waters during glacial stages (e.g., Last Glacial Maximum) as well as cold Dansgaard-Oeschger and/or Heinrich-like events (e.g., the Younger Dryas) (Ganeshram and Pedersen, 1998; Hendy and Kennett, 2000; Nameroff et al., 2004; Hendy and Pedersen, 2006; Tetard et al., 2017). Paleoceanographic changes are not only observed in sediment core records at OMZ-depths (400–1,000 m) but also in sediments collected at greater depths of the Pacific and shallower waters of the glacial Southern Ocean, which confirms the likelihood of widespread impacts on oxygenation far beyond intra-OMZ variability (Lu et al., 2016; Hoogakker et al., 2018). Important fluctuations of the upper OMZ also occurred. Rather than the traditional view of these fluctuations as lesser responses to minor climate perturbations, they may have dominated certain regions such as the ETNP. The Pacific Decadal Oscillation (PDO)—detected as decadal cyclic changes in climatic regime from warm to cold and vice versa—is one of these atmospheric effects with strong oceanographic implications (Mantua and Hare, 2002). The North Pacific Gyre Oscillation (NPGO) manifests in decadal fluctuations of nutrients, chlorophyll, and consequentially primary productivity that very often go hand to hand with atmospheric PDO changes in the Pacific (Di Lorenzo et al., 2008). Despite the similarities between the PDO and NPGO, they do not always co-vary.
High-resolution variability in paleoceanographic records, previously categorized by the scientific community as regional background noise, is now elucidating decadal- and centennial-scale OMZ fluctuations (Dean et al., 2004; Barron and Bukry, 2007; Choumiline, 2011; Deutsch et al., 2014; Tems et al., 2015, 2016). In the ETNP, cold PDO and/or NPGO phases cause intense upwelling off the western mainland of North America and lead to enhanced productivity, denitrification, and extreme water column anoxia (Di Lorenzo et al., 2008). However, numerous reports of marine sedimentary records in the literature typically lack the highly resolved age models needed to study decadal phenomena, which coupled with the absence of published multi-proxy records for denitrification (δ15N) and water column redox (Cd, Mn, V, Mo, U) make ETNP-wide redox reconstructions challenging (Moffitt et al., 2015; Borreggine et al., 2017). The importance of reconstructing OMZ reoxygenation and deoxygenation events and mechanistically explaining their climatic causes is essential in light of current and potential future responses of the ocean to modern global warming. These efforts encounter further challenges as redox-sensitive data can be compromised when sediment cores from anoxic settings are not carefully sampled or stored with adequate care in core repositories (Zheng et al., 2000; McManus et al., 2006; Scholz et al., 2011; Costa et al., 2018).
Here we reconstruct changes in the ETNP OMZ based on sediment cores collected at various depths in the Gulf of California. To achieve an ETNP-wide reconstruction, we integrate our multiproxy record of paleoproductivity and paleoredox with published sediment cores from the region. Further, to avoid misinterpreting local variability as global patterns, we apply geostatistical and numerical approaches to existing datasets to unmix local from larger-scale controls.
Strengths and Weaknesses of Geochemical Proxies as Applied to Oxygen Minimum Zone Reconstructions
In this section we provide background details for productivity and redox proxies as archived in marine sediments and as applied to OMZ reconstructions.
Productivity Proxies (Corg, Baexcess, Ni)
Organic carbon (Corg) has long been used as a proxy for exported organic matter because of its analytical simplicity and wide range of applications (Schoepfer et al., 2015). While Corg could directly represent superficial primary productivity, this possibility is likely confounded by multiple controls on remineralization and burial that are tightly related to reactivity of organic matter, bioturbation, and oxygen availability (Calvert et al., 1996; Burdige, 2007; Bianchi et al., 2016).
Barite, Ba/Al, and excess barium (Baexcess) have been used widely as productivity proxies (Dymond et al., 1992; Eagle et al., 2003; Anderson and Winckler, 2005; Schoepfer et al., 2015), although barite and Baexcess relationships with productivity may be complicated by post-burial reprecipitation (McManus et al., 1998; Riedinger et al., 2006; Schoepfer et al., 2015). In some well-oxygenated marine settings, sedimentary Baexcess corresponds well to changes in productivity established by other proxies (e.g., Ni, P, and biogenic silica) (Schoepfer et al., 2015). However, post-depositional diagenetic remobilization of Ba has been found both in regions with transient anoxia and where bottom waters are well-oxygenated (McManus et al., 1998; Riedinger et al., 2006; Schoepfer et al., 2015). Sedimentary records of Baexcess can be further compromised over glacial-interglacial cycles when significant reoxygenation can cause redissolution and reprecipitation of barite along reaction fronts (McManus et al., 1998; Riedinger et al., 2006; Schoepfer et al., 2015).
Nickel is actively involved in biological cycling related to diatom export (Dupont et al., 2010; Twining et al., 2012). Dissolved Ni is removed from the water column by organic matter and supplied to the bottom sediments through particle settling (Böning et al., 2015). In sediments, Ni is the best-known analog for fresh organic matter because it mirrors the input of chlorins—one of the degradation products of chlorophyll (Böning et al., 2015). An advantage of Ni over other biologically active trace elements is its detachment from Mn and S cycling, compared to Cd and Zn, which more readily form insoluble sulfides (Morse and Luther, 1999; Böning et al., 2015). Potential complications for Ni as a productivity proxy include its elevated concentration in the detrital fraction (Wedepohl, 1995) and the possibility of anthropogenic contribution of Ni to marine sediments (Reck et al., 2008).
Redox Proxies (Mo, Cd, V, U)
Molybdenum enrichment processes are linked with incorporation into pyrite and associations with organic matter, which often require abundant hydrogen sulfide (Calvert and Pedersen, 1993; Tribovillard et al., 2008; Chappaz et al., 2014; Scholz et al., 2017). An important pathway for Mo enrichment is the transformation of molybdate to thiomolybdates , where x = 0 ~ 4), which are highly particle reactive (Helz et al., 1996; Erickson and Helz, 2000; Scholz et al., 2017). Additionally, Mo can be scavenged by Fe-Mo-S clusters (Helz et al., 2011) or adsorbed onto sulfide-rich organic matter (Helz et al., 1996). Molybdenum enrichments can occur diagenetically (within the sediments), in the water column, or in some cases near the sediment-water interface, which was reported for OMZs (Zheng et al., 2000). While widely used as an indicator of past anoxic and euxinic (anoxic and sulfidic) conditions, several complications have been reported for this proxy (Morford and Emerson, 1999; Chappaz et al., 2014), including localized Mo enriched in turbidites despite fully oxygenated water column conditions (McKay and Pedersen, 2014) and uptake exclusively within the sediments when significant amounts of sulfide are present solely in pore waters (Zheng et al., 2000; Scholz et al., 2017; Hardisty et al., 2018)—although these enrichments tend to fall below those typical of euxinic settings (e.g., Scott and Lyons, 2012).
Cadmium is a biologically active trace element supplied to marine sediments in association with organic matter (Morford and Emerson, 1999). It has been used as a proxy for strongly reducing waters (Calvert and Pedersen, 1993; Tribovillard et al., 2006). Cadmium can be precipitated in the presence of H2S, which makes it a particularly good indicator of euxinia (Morse and Luther, 1999).
Vanadium is present in well-oxygenated seawater in its V(+5) form, usually as H2, which is reduced to V(+4) as VO2+ under oxygen-deficient conditions (Calvert and Pedersen, 1993; Morford and Emerson, 1999). This V(+4) form is more effectively scavenged and can be adsorbed to Fe-(oxyhydr)oxides, which is one of its removal mechanisms from seawater (Bauer et al., 2017; Ho et al., 2018). Over sub-millennial and glacial-interglacial timescales, V is potentially one of the most reliable redox proxies, tightly following changes in oxygenation (Nameroff et al., 2004; Costa et al., 2018). The contribution of non-lithogenic V through particulate settling has not been fully constrained; however, slight enrichments in suspended and settling matter have been found (Rodríguez-Castañeda, 2008; Hakspiel-Segura et al., 2016; Bauer et al., 2017; Ho et al., 2018; Conte et al., 2019).
Uranium is widely used as paleoredox proxy; however, the pathways of its sedimentary enrichment are still not fully understood. Uranium is present in seawater in the dissolved U(+6) form as uranyl carbonate complexes (UO2)2+ that are reduced under anoxic conditions to the insoluble U(+4) form by inorganic and microbially mediated processes (Ku et al., 1977; Klinkhammer and Palmer, 1991; Dunk et al., 2002) resulting in the formation of solid particles. Uranium can also be adsorbed to organic matter and to Fe- and Mn-(oxyhydr)oxides (Anderson et al., 1989). While increases in particulate U are thought to mostly occur diagenetically in sediment pore waters underlying strongly anoxic and sulfidic (euxinic) waters (Klinkhammer and Palmer, 1991), some particulate non-lithogenic U can form very slowly in the water column (Zheng et al., 2002; Holmden et al., 2015). The extent to which non-lithogenic U can be enriched in marine particles is not well-established and may be more common in highly productive OMZ-type settings than previously thought. We propose an updated model for the behavior of U in OMZs, which will be discussed further in the paper.
The Gulf of California (GC) is located within the ETNP OMZ and is well-connected to the Eastern Pacific through its southern entrance (Lavín et al., 2013). Circulation is not known to be limited in most GC basins, including Guaymas, Pescadero, and La Paz. Alfonso Basin is the only partially restricted setting in the GC, with a bathymetric sill extending today to a depth of 250 m (Nava-Sánchez et al., 2001). The sill restricts the incursion of North Pacific Intermediate Water (NPIW) and North Pacific Deep Water (NPDW) (Figures 1A,B). Primary productivity is higher on the eastern side of the GC because it is nourished by nutrients upwelled predominantly during dry and cold seasons with strong northwesterly winds (Lavín et al., 2013; Silverberg et al., 2014). Mesoscale eddies together with wind-driven coastal upwelling (primarily along the mainland margin during winter) promote primary productivity in the region (Badan-Dangon et al., 1991; Trasviña-Castro et al., 2003; Deutsch et al., 2014). A strong east-to-west gradient persists from autumn through spring (Lavín et al., 2013), which results in organic material produced in the upwelling system on the eastern margin being deposited as biogenic remains on the west side of the basin.
Figure 1. (A) Location of the studied marine setting in the eastern margin of the Gulf of California including the anoxic Alfonso and La Paz basins. The line between SW and NE represents a cross section shown in (B). (B) Cross section from the Baja California peninsula (SW) to the Gulf of California (NE) representing the basin morphology and dissolved oxygen (mL/L) concentration spanning Alfonso and La Paz basins.
Particularly in Alfonso Basin, productivity is regulated by cyclonic flow that produces water divergence—transporting nutrients from the subsurface waters (>60 m) into the photic zone. This cyclonic gyre is semi-permanent and has been linked by past researchers to the increase in productivity (Monreal-Gómez et al., 2001; Coria-Monter et al., 2014). Satellite and in-situ measurements of Net Primary Productivity (NPP) in Alfonso Basin range from 0.5 to 3.5 gC m−2d−1 over the last decade and show an inverse correlation with surface temperature. Mesoscale transport processes explain why the Gulf of California, like most of the Eastern Pacific, is influenced by ENSO, PDO, and/or NPGO phenomena (Mantua and Hare, 2002; Di Lorenzo et al., 2008). During cold La Niña years, the Walker Circulation triggers upwelling and productivity in the Eastern Pacific—weakening the typical upper ocean stratification that dominates during normal and El Niño years (Thunell, 1998; McPhaden, 2004; Bustos-Serrano and Castro-Valdez, 2006). Atmospheric circulation is also an important mechanism that seems to control upwelling in the Gulf of California and is dictated by the North American Monsoon (NAM) (Metcalfe et al., 2015). When the North Pacific High (NPH) and the Intertropical Convergence Zone (ITCZ) migrate southward, intense upwelling is triggered by northwesterly winds (Schneider et al., 2014; Metcalfe et al., 2015; Staines-Urías et al., 2015). Longer-scale variability of similar nature is attributed to cold PDO and sometimes NPGO phases, but the atmospheric and oceanographic links are still unclear.
The water columns of Alfonso and La Paz basins show a pronounced decrease in oxygen with depth, reaching values below 0.1 ml/L O2 (Figure 1B). To avoid the imprecise term “suboxic” (Canfield and Thamdrup, 2009), we classify the water column using the threshold O2 values from Hofmann et al. (2011) and Moffitt et al. (2015), which yielded the following categories: oxic ([bottom O2] > 1.5 mL/L), moderately hypoxic (0.5 mL/L < [bottom O2] <1.5 mL/L), severely hypoxic (0.1 mL/L < [bottom O2] <0.5 mL/L), and anoxic ([bottom O2] <0.1 mL/L). According to this classification, both Alfonso and La Paz basins fall within the “anoxic” category (Figure 1B).
Due to their locations near the open Pacific Ocean, the southern basins are a sensitive recorder of regional variation in the GC and the larger-scale climate circulation of the ETNP, capturing changes in productivity and redox conditions useful for tracking fluctuations of the OMZ (Dean et al., 2004; Barron and Bukry, 2007; Dean, 2007; Pérez-Cruz, 2013; Deutsch et al., 2014; Tems et al., 2015, 2016; Ontiveros-Cuadras et al., 2019).
Laminated sediments from marine basins of the Gulf of California with varying degrees of deoxygenation were collected using a Reineck box corer during the PALEOMAR-I campaign in November 2014 aboard R/V “El Puma” owned by UNAM, Mexico (Figures 1A,B). The 50 cm sediment core PALEOMAR-I C37A was collected at 408 m water depth in the silled Alfonso Basin. Core PALEOMAR-I C52A represents a sedimentary sequence of the upper 24 cm of La Paz Basin at 785 m depth. Both cores show a clear, continuously laminated fabric. Additionally, the core top of PALEOMAR-I C42A was sampled from the shallower part of Alfonso Basin (265 m). This core did not show clear laminations. Poor recovery during retrieval limited its utility for paleoceanographic reconstructions. From this point forward, we will refer to these cores as C37A, C42A, and C52A.
Each core was extruded and subsampled immediately after collection, and pore waters were extracted from each section in a glove bag under N2 atmosphere using CSS 5 cm Rhizon samplers attached to plastic syringes (Lyons and Berner, 1992; Dickens et al., 2007; Raiswell and Canfield, 2012; Riedinger et al., 2017). Each core segment yielded approximately 8 mL of pore waters, which were filtered at 0.2 μm, acidified with trace grade pure nitric acid (10 μL of acid per each 1 mL of sample), and stored at 4°C (Riedinger et al., 2014, 2017). Prior to analysis, solid-phase samples were dried, weighed, and homogenized.
We use the present-day sedimentation rates previously reported in the literature for Alfonso Basin to generate the age model for the studied core C37A. Due to the extended timespan of the sedimentary section, the composite age model is based on 210Pb and 14C dating methods (Figure 2; Table 1). For this purpose, we selected box core DIPAL-II C44 (Choumiline, 2011) and gravity core DIPAL-III T43 (Pérez-Cruz, 2013). The dated materials were collected at the same location within the basin and span a matching sedimentary package (Supplementary Figure 1). High-resolution sections of the uppermost 10 cm (subsampled at <3 mm) of core DIPAL-II C44 were analyzed for 210Pb at Flett Research Ltd. Laboratory (Winnipeg, Canada) (Supplementary Figure 2). The mean sedimentation rate for the upper 10 cm is 0.61 mm/year in Alfonso Basin (Choumiline et al., 2010; Choumiline, 2011). The sedimentary section of the gravity core DIPAL-III T43 corresponding to the lower part of C37 (between 15 and 50 cm) was dated using the 14C method (Pérez-Cruz, 2013). Measurements were performed in benthic foraminifera Bolivina subadvena by accelerator mass spectrometry at Beta Analytic (Pérez-Cruz, 2013) and the radiocarbon values were calibrated using the Marine13 curve (Reimer et al., 2013). This part of the core yielded a mean sedimentation rate of 0.8 mm/year (Pérez-Cruz, 2013). In our age model for C37A we extrapolate the sedimentation rate of 0.61 mm/year to the uppermost core section (0–15 cm) and use 0.8 mm/year for the lowermost segment (15–50 cm) as shown in Table 1 and Supplementary Figure 1. To demonstrate that the core sections are equivalent, we compared the total carbon (Ctot), calcium carbonate (CaCO3) and Ca changes through time for cores C37A (this study), DIPAL-II C44 (Choumiline et al., 2010, 2017; Choumiline, 2011) and DIPAL-III T43 (Pérez-Cruz, 2013). The curves covary, confirming the consistency between age models for different core records (Supplementary Figure 1). Additionally, selected dry splits of sediment core C37A (Alfonso Basin) were dated using the 210Pb and 137Cs methods via gamma-ray spectroscopy in the Gray Lab (University of California, Riverside). The subsampling resolution of this core was too coarse to produce a robust age model using these methods. Despite known difficulties with radiocesium detection in the region (Ruiz-Fernández et al., 2009) it was still possible to find slightly elevated 137Cs activity values around the 2–4 cm core horizon, very likely corresponding to the mid-1960s peak of nuclear testing in the atmosphere (Foster and Walling, 1994). The resulting linear sedimentation rate using the radiocesium method for core C37A was roughly 0.6 mm/year. The sedimentation rates determined for Alfonso Basin from different age models are consistent. Further, the age model for Alfonso Basin shows that the deepest sample of core C37A at 500 mm extends to 1300 AD (also called the Common Era or CE).
Table 1. Age model for core C37A, constructed by a combination of extrapolated average sedimentation rates produced using the 210Pb method on core DIPAL-II C44 (Choumiline et al., 2010; Choumiline, 2011) and acceleration mass spectrometry 14C on core DIPAL-III T43 (Pérez-Cruz, 2013).
Total carbon (Ctot) and total inorganic carbon (Cinorg) in dry sediment subsamples were measured on an ELTRA CS-500. The CO2 produced during subsample combustion at 1350°C (Ctot) or acidification with 20% HCl (Cinorg) was quantified with an infrared analyzer. The organic carbon (Corg) fraction was calculated by numerical difference: Corg = Ctot -Cinorg.
Prior to sample digestion, 100 mg splits were ashed in ceramic vessels for 8 h at 550°C in a high temperature oven. A two-step digestion procedure involved a mixture of HNO3 and HF, followed by aqua regia (HCl and HNO3) and acidification/drying cycles on temperature-controlled hot plates. A very low HF:HNO3 ratio of 1:8 was employed during the first step to prevent insoluble fluoride formation. Careful evaluation of each sample revealed that no undissolved particles were detected, suggesting complete digestion (Durand et al., 2016). The final digested extracts were diluted in 2% HNO3 and measured for major and trace elements (Al, Ba, Cd, Fe, Ni, Mo, V, U) on an Agilent 7900 Inductively Coupled Mass Spectrometer (ICP-MS).
Pore waters were analyzed for elemental composition on an Agilent 7900 ICP-MS, emphasizing dissolved Fe, Mn, and selected trace elements. The quality of our data was assessed by analyzing duplicates and certified reference materials USGS SGR-1B, USGS SDO-1, NIST 2702, AR4007, and AR4024. Analytical precision was always below 3%. Recoveries of measured sample relative to certified standard value ranged between 96 and 105%.
Numerical and Statistical Approaches
The productivity and redox proxy approaches used here rely heavily on authigenic enrichment of trace elements relative to detrital inputs through uptake in marine sediments during diagenesis and via scavenging from seawater during primary production and during particle settling. To assess enrichments relative to detrital contributions, each element was normalized to Al as a reliable crustal proxy (Wedepohl, 1995). Further, to account for dilution by CaCO3, which can cause apparent elemental variability, we recalculated the data on a carbonate free basis (CFB) by dividing each concentration by the fraction of non-carbonate material. Elemental ratios were also calculated for selected elements. The authigenic uranium (Uauth) fraction was calculated as Uauth = Utotal – (Th/3) (Wignall and Myers, 1988).
Due to limited sample sizes, it was not possible to perform the barite sequential leach sometimes used as a productivity proxy in well-oxygenated marine settings (Dymond et al., 1992; Eagle et al., 2003; Schoepfer et al., 2015). Instead, we used our measurements for total Ba and Al to determine excess barium (Baexcess) as the closest calculated alternative (Eagle et al., 2003; Anderson and Winckler, 2005): Baexcess = Batotal – (Baterrigenous/Alterrigenous) × Altotal. We assumed a Baterrigenous/Alterrigenous ratio of 0.0045 as suggested by Dymond et al. (1992) and Eagle et al. (2003).
In order to define proxy associations, we applied a multivariate statistical approach for the proxy datasets from Alfonso and La Paz basins. Specifically, we performed R-mode Factor Analysis with Varimax rotation using the software STATISTICA by StatSoft. The analysis included the following variables: Corg, Baexcess, Ni, U, Uauth, Cd, Mo, and V.
Geostatistical analyses and data extraction techniques were performed using ArcGIS 10.6 software by ESRI. We combined our data with previously published geochemical records obtained in part from the EarthChem SedDB database via the Interdisciplinary Earth Data Alliance (IEDA) repository. Estimated oceanographic data (dissolved O2 and chlorophyll) were obtained for each core location in GIS and used in our statistical analysis. Global oceanic oxygen data are based on the World Ocean Atlas (2005) and Garcia et al. (2006); surficial chlorophyll estimates from the SeaWIFS satellite sensor were published by OBPG/NASA.
Paleoceanographic Proxies Used in This Study
To interpret past OMZ conditions preserved in marine sediments we used a combination of proxies assuming Corg as an indicator for organic carbon burial (Schoepfer et al., 2015); Baexcess (despite the caveats described above) and Ni for export production (Böning et al., 2015; Schoepfer et al., 2015); Cd, Mo, and V for bottom water anoxia (Calvert and Pedersen, 1993; McManus et al., 2006; Tribovillard et al., 2006; Lyons et al., 2009). Due to a strong contribution of Uauth via settling particulates (Choumiline et al., 2010; Choumiline, 2011) potentially after its reduction or adsorption to organic matter, this proxy could trace both redox and productivity—as will be explained further in the paper.
While each proxy has distinct biogeochemistry and enrichment pathways, they are all affected by a circular relationship between productivity and redox (Tribovillard et al., 2006; Schoepfer et al., 2015; Bianchi et al., 2016). In transiently anoxic settings high productivity often leads to the development of water column anoxia, which favors organic preservation and reduces the mobility of productivity proxies (Schoepfer et al., 2015). Due to the latter, it is often common to expect a covariance between productivity and redox proxies (Tribovillard et al., 2006; Bianchi et al., 2016). The differences in proxy relationships are often helpful to unravel more specific biogeochemical processes tied to preservation biases, such as organic matter reactivity and post-depositional diagenetic transformations.
Results and Discussion
Geochemistry of Modern Sediments From the Southwestern Gulf of California at Various OMZ Depths
The chemical compositions of core tops (upper 2 cm) represent the most recent sediments from Alfonso and La Paz basins (Table 2). A water depth transect from the 265 to 785 m through the study region revealed that dissolved O2 values decreased from 0.3 mL/L (uppermost OMZ) to <0.1 mL/L (OMZ core) (Figure 1B). Most redox-sensitive elements increase in concentration in the deeper sites C37A and C52A relative to the shallower C42A. For all three sites the elemental concentrations range as follow: Corg (2.96–6.41%), Baexcess (175.50–682.77 mg/kg), Ni (26.84–49.06 mg/kg), and Mo (0.75–6.75 mg/kg). In contrast, Cd and Uauth are enriched at the shallow, more oxygenated site C42A—reaching 13.17 mg/kg for Cd and 5.34 mg/kg for Uauth, compared to 7.24 mg/kg for Uauth in C52A and 2.48 mg/kg for Cd in C37A (Table 2). Although circulation in Alfonso Basin is slightly restricted due to the presence of a bathymetric sill at 250 m depth, restriction was apparently not severe enough to limit trace element enrichments in La Paz Basin relative to the more open ocean. We calculated Uauth values for the CFB fraction as 6.43 mg/kg in C42A, 7.26 mg/kg in C37A, and 7.24 mg/kg in C52A. Carbonate-free values were also calculated for the other proxies but did not produce major differences relative to either total or Al-normalized data (see Supplementary Material).
Table 2. Composition of modern sedimentary material represented by core tops from different OMZ depths in Alfonso (C42A and C37A) and La Paz basin (C52A).
Bioturbation and Post-sampling Effects as Complicating Overprints on Primary Proxy Records
Visual observation confirmed that the sediments from Alfonso (C37A) and La Paz (C52A) basins are clearly laminated, with no evidence for bioturbation, which could obscure primary sedimentary features. Pore water profiles were used to establish the chemical zonation of electron acceptors at the time of sample collection (Canfield and Thamdrup, 2009). It is essential to evaluate diagenetic conditions before applying paleoreconstruction techniques because many proxies can be overprinted during diagenesis, obscuring records of primary water column conditions (Calvert and Pedersen, 1993; Tribovillard et al., 2006). Further, pore water data shed light on bioturbational overprints via burrow irrigation and sediment homogenization and disturbance during sediment collection and processing. At both Alfonso Basin (C37A) and La Paz Basin (C52A), concentrations of dissolved Fe (Fediss) and dissolved Mn (Mndiss) show downcore decreases (Figure 3). The observed trend for Fediss is similar to those from analogous OMZ settings of the Gulf of California (Guaymas Basin), offshore California (San Pedro and Catalina basins), and the Peruvian margin—that is, pore water Fediss reaches its maximum a few centimeters below the sediment-water interface and gradually declines down core (Brumsack and Gieskes, 1983; McManus et al., 1997; Canfield and Thamdrup, 2009; Scholz et al., 2011; Scholz, 2018). Our recovery of expected trends suggests that proper techniques for sediment collection and preservation were employed in this study.
Figure 3. Pore water chemical profiles in cores PALEOMAR-I C37A from Alfonso Basin and PALEOMAR-I C52A from La Paz Basin (Data: Choumiline et al., 2017).
Downcore Records of Changing Productivity and Redox
Here we present downcore variations of productivity and redox proxies recorded in our sedimentary material from the ETNP (Figure 1). We found almost no difference between the trends for absolute concentrations and corrected data, so we present only absolute values. However, Al-normalized and CFB data are presented in Supplementary Materials for the interested reader.
Geochemistry of Sediments From the Upper OMZ in Alfonso Basin
Elemental profiles for C37A from Alfonso Basin are shown in Table 3. Data for Corg are in good agreement with Baexcess, Ni, Uauth, and V as proxies for past oxygenation and/or primary productivity. Values range from 5.60 to 6.76% for Corg, 372.22–491.54 mg/kg for Baexcess, 40.25–51.61 mg/kg for Ni, 9.23–14.15 mg/kg for U, 7.71–12.48 mg/kg for Uauth, and 115.81–172.40 mg/kg for V. Trends for Cd and Mo differ slightly from those for the other proxies and range from 2.48 to 7.04 mg/kg and 6.34–15.16 mg/kg, respectively (Figure 4).
Table 3. Elemental composition of the solid phase from cores C37A, C42A, and C52A collected during PALEOMAR-I expedition to the Gulf of California.
Figure 4. Organic carbon, productivity, and redox proxies in core PALEOMAR-I C37A from Alfonso Basin. Data for temperature anomaly reconstructions of the last millennium taken from IPCC (2014).
Geochemistry of Sediment From the Mid-Depth OMZ of La Paz Basin
Downcore concentrations for the productivity and redox proxies in core C52A from La Paz Basin are presented in Table 3. Despite differences in core length and depositional setting, the records for Corg, Ni, Uauth, V, and Cd in C52A are similar to those from core C37A (Supplementary Figure 3). These elements range as follow: 6.51–7.12% for Corg, 47.34–53.25 mg/kg for Ni, 5.97–9.32 mg/kg for Uauth, and 74.31–102.80 mg/kg for V. In C52A, Baexcess does not covary with Ni and Uauth, pointing to its potential remobilization during diagenesis. Excess-Ba concentrations range from 419.76 to 763.24 mg/kg but are mostly uniform over the length of the core—averaging 670.89 mg/kg. The trends for Cd and Mo decrease down core in the upper 6 cm of laminated sediments (Supplementary Figure 3; Table 3) and range from 2.38 to 2.97 and 4.90 to 12.55 mg/kg, respectively.
The lowest concentrations for most proxy data are found in the upper 5 cm, at ~14 cm, and below 23 cm—perhaps corresponding to well-ventilated past conditions. The highest values, indicating enhanced anoxia, were registered around 20 cm for all proxies, while the second anoxic horizon close to 10 cm was only captured by Corg, Ni, Uauth, V, and more subtly by Cd and Mo.
Variability of the Upper OMZ in Alfonso Basin Over the Last 700 Years
We present down-core elemental concentrations in the context of the age model described earlier in the paper. Due to multiple challenges during our efforts to date cores C42A and C52A, such as poor recovery, only core C37A from Alfonso Basin produced a robust age model. Core C37 will be emphasized for paleoceanographic reconstructions for the remainder of the report.
The geochemical profiles for C37A from Alfonso Basin record changes in the upper OMZ boundary dating back to 1300 AD (Figure 4; Table 3). Data for Corg are in good agreement with Baexcess, Ni, Uauth, and V as proxies for past oxygenation and/or primary productivity. Unlike our other proxies, Mo and Cd remain high from 1300 to 1700 AD and then slowly transition toward lower values after 1700 AD. This slight decrease toward the present corresponds to the uppermost 20 cm of the record in the Alfonso Basin. The lowest concentrations for our proxies are generally concurrent with the late 1900s, early 1700s, and late 1300s AD. These times correspond to warm periods of the past millennium as reconstructed by multiple records (IPCC, 2014), including the anthropogenically caused warming trend that started around the beginning of the nineteenth century and extends to present day (Figure 4). Lower proxy values indicate decreased primary productivity and correspond to slightly more oxygenated times. In contrast, high proxy values occurred in the early 1800s, late 1600s, early 1500s, and early 1400s, with maxima between roughly 1700 and 1400 AD are likely associated with the prolonged cold period of the Little Ice Age. Higher values usually correspond to increased primary productivity in surface waters and enhanced anoxia.
Multivariate Statistical Associations of Paleoproxies in Alfonso and La Paz Basins
Factor analysis results for the combined geochemical datasets from Alfonso Basin (core C37A) and La Paz Basin (C52A) distinguished three elemental associations between Factor 1 (48.23% variance) and Factor 2 (27.85% variance) (Table 4; Figure 5). The “productivity” association is based on Corg, Ni and Baexcess–elements related to organic matter production at shallow depths. The “redox with weak particulate contribution” defined by Cd and Mo, which require strongly reducing conditions in bottom waters for strong authigenic enrichments. The “redox with strong particulate contribution” links U, Uauth, and V—elements that can be supplied through particulate fluxes in oxygen-deficient water column depths and/or enriched diagenetically either under anoxic bottom waters or at the sediment-water interface.
Table 4. Factor loadings of the R-mode Factor Analysis with Varimax rotation performed on a combined dataset of cores C37A (Alfonso Basin) and C52A (La Paz Basin).
Figure 5. Productivity and redox proxy associations defined by R-mode factor analysis of a combined geochemical dataset from Alfonso Basin (core C37A) and La Paz Basin (C52A).
Comparison of Productivity and Redox Proxies in Core Tops Across Varying OMZ Depths and Productivity Regimes in the ETNP
Here we present a compilation of data from core tops from this and previous studies (Table 5; Figure 6 and Supplementary Figure 4). Comparing core tops, most of the cores collected at upper- and mid-OMZ depths have significantly stronger Corg, Baexcess, Ni, Mo, and V enrichments compared to shallower cores (e.g., core C42A)—with the latter perhaps representing the fluctuating upper OMZ boundary often exposed to high O2 levels (Table 5). These enrichments correspond to regions where high chlorophyll is reported in the surface waters associated with eastern North Pacific upwelling (Table 5; Figure 6 and Supplementary Figure 4) (Badan-Dangon et al., 1991; Trasviña-Castro et al., 2003). The elevated Corg, Baexcess, and Ni in these high-chlorophyll regions are likely to indicate enhanced primary productivity in the photic zone but could also reflect strong preservation of refractory organics delivered from land (Calvert et al., 1996; Burdige, 2007; Bianchi et al., 2016) (Table 5). These productive coastal sectors correspond to severe hypoxia/anoxia in which trace element enrichments reach 75.05–149.08 mg/kg for V and 6.75–7.35 mg/kg for Mo in the deepest cores from La Paz and Pescadero basins, respectively (Table 5; Figure 6). Detrital contributions in such settings are relatively minor (Rodríguez-Castañeda, 2008; Choumiline, 2011). These enrichments are favored by oxygen-deficient conditions in the water column and, in the case of Mo, sulfide accumulation in seawater (Calvert and Pedersen, 1993; Zheng et al., 2000; Tribovillard et al., 2006). The studied core tops from the Eastern Pacific OMZ show Mo values sufficiently high to indicate anoxia but are too low to suggest euxinia (presence of H2S) (Table 5; Figure 6A). For comparison, only certain places off the Peruvian margin, where water conditions are extremely reducing even for an OMZ, show Mo values as high as 70 mg/kg (McManus et al., 2006). In contrast, at 265 m depth in the upper portion of the OMZ, core C42A displays low V and Mo values of 26.05 and 0.75 mg/kg, respectively, indicating the presence of some oxygen in the overlying waters and that V and Mo are mostly supplied through detrital inputs.
Table 5. Comparison of productivity and/or redox proxy values for selected sedimentary material from different OMZ depths in the Eastern Pacific.
Figure 6. Molybdenum and organic carbon in core tops from various locations across the Gulf of California and adjacent Pacific, (A) bottom dissolved oxygen, and (B) superficial chlorophyll estimate from the studied marine OMZ settings of ETNP (Alfonso Basin and La Paz Basin) (Oxygen based on World Ocean Atlas, 2005; Garcia et al., 2006; superficial chlorophyll estimates by SeaWIFS satellite sensor published by OBPG/NASA).
Uranium is the exception. The Uauth concentrations are high in the shallower and more-oxygenated core (C42A), reaching 13.17 mg/kg—with values similar to those of the deeper, anoxic core tops (C37A and C52A). We hypothesize that the added contribution of uranium is through its uptake in settling particles and associated anoxic microenvironments. Under oxic and weakly oxic conditions, dissolved uranium is present in its +6 form but is reduced to the insoluble +4 form once oxygen declines further (moderate hypoxia) or disappears (anoxia). This authigenic U is preserved as the organic-rich marine aggregates settle through the water column, in the absence of Uauth reoxidation even in ventilated waters. This exported authigenic uranium is delivered to the sediments—often showing a stronger correlation with Corg compared to the relationships between the organic matter and V or Mo (Choumiline, 2011). Evidence for authigenic U enrichments in settling particles was reported by Zheng et al. (2002) in marine basins of the California Margin and also from Alfonso Basin, with particulate U concentrations reaching 40 mg/kg (Choumiline et al., 2010; Choumiline, 2011)—more than an order of magnitude higher than the crustal average of 2.5 mg/kg (Wedepohl, 1995).
Like Uauth, V is enriched in sediments from regions of the Pacific where bottom water anoxia is present. Additionally, its concentration in suspended particles also peaks seasonally at depths of 150–250 m in response to intra-annual changes in water column O2 (Hakspiel-Segura et al., 2016). The contribution of V to the bottom via settling particles is usually lower than the average sedimentary concentration under low O2 conditions but still within the same order of magnitude (Rodríguez-Castañeda, 2008).
To assess the contribution of authigenic U through settling, we contrast U and Mo sedimentary records published in the literature with ours from the ETNP (Ivanochko and Pedersen, 2004; McManus et al., 2006). Typically, Mo is used as an indicator for water euxinia (Calvert and Pedersen, 1993; Algeo and Lyons, 2006; Algeo and Tribovillard, 2009; Lyons et al., 2009; Algeo and Rowe, 2012; Scott and Lyons, 2012), but most of the enrichments in modern OMZs are likely tied to the presence of sulfide only in pore waters (Scholz et al., 2017; Hardisty et al., 2018). We are not using V because, like U, there could also be important contribution of the redox sensitive element through settling particles, associated to reducing microenvironments (Table 5; Choumiline, 2011). We expect good correlation between U and Mo records if both elements are enriched during burial (i.e., diagenetically or at the sediment-water interface; Klinkhammer and Palmer, 1991; Calvert and Pedersen, 1993; Zheng et al., 2002; Algeo and Tribovillard, 2009). A weak correlation is likely to be found for samples where significant amounts of non-lithogenic U are scavenged into particles and settle through the water column. In contrast, Mo uptake is typically limited to the sediments when oxygen is present in the water column. For the previously published data, we estimated bottom O2 values for each sediment core in a GIS using an interpolated global oxygen database (World Ocean Atlas, 2005; Garcia et al., 2006). This approach allowed us to classify data in two main groups depending on the dominant redox conditions: oxic to moderately hypoxic (0.5 mL/L < [bottom O2] <1.5 mL/L) (Figure 7A) and severely hypoxic (0.1 mL/L < [bottom O2] <0.5 mL/L) to anoxic ([bottom O2] <0.1 mL/L) (Figure 7B). We found that sediments underlying oxic–moderately hypoxic waters do not show correlation between U and Mo—with an R2 value close to 0.00 (Figure 7A). Some samples do, however, show significant enrichment above crustal values of 2.5 mg/kg that could be due to water column U authigenesis (Zheng et al., 2002; Choumiline, 2011). In contrast, sediments deposited under severely hypoxic to anoxic waters show strong correlation between Mo and U, with an R2 value of 0.80. While most of the U concentrations for the oxygen-deficient sites correlate with Mo, some reach even higher values of 17 mg/kg U with no significant Mo increase. In many cases, the lack of Mo enrichments could be attributed to bottom water redox potentials that were not sufficiently reducing, especially for sites in the upper OMZ. These observations of high U and low Mo could also indicate a water column supply of authigenic U through settling particles. In this context of U-Mo relationships, the studied sites from Alfonso and La Paz basins fall above the trend line, showing potential surpluses in Uauth enrichments perhaps tied to high levels of productivity and associated particle scavenging (Figure 7).
Figure 7. Correlation between U and Mo in sediments collected under (A) oxic or moderately hypoxic waters and (B) severely hypoxic or anoxic waters from different basins of the Eastern Tropical North Pacific (Data from: this study; Ivanochko and Pedersen, 2004; McManus et al., 2006; Choumiline et al., 2017; EarthChem database with n = 183).
Chronology of Reoxygenation and Deoxygenation Events Over the Past Millennium in a Predominantly Anoxic Marine Setting
In the following discussing, we use the term reoxygenation to imply a relative increase in dissolved O2 at OMZ depths, usually corresponding to deepening of the upper OMZ boundary. Reoxygenation need not indicate disappearance of the OMZ, as no proxy record indicates fully oxygenated water conditions (Tems et al., 2016), including the data produced in our study—specifically, Mo, V, and U values are several times higher than their crustal average values (Wedepohl, 1995). The oxygenation chronology of the last millennium starts with the predominantly warm climatic trends of the Medieval Warm Period (or Medieval Climate Optimum), which ended as cooling began with gradual change in the insolation regime after the Wolf Minimum (Figure 8). This cooling transitioned the planet into the colder Little Ice Age, which started around 1300 AD (Bard et al., 2000; IPCC, 2014; Griffiths et al., 2016) and was the most anoxic period of the last millennium as indicated by our V, Cd, Mo, and Uauth data from Alfonso. The Little Ice Age coincides with the highest values of the SOI index, indicating more La Niña-like years and enhanced upwelling (Yan et al., 2011; Griffiths et al., 2016). Most marine basins of the ETNP experienced short-term reoxygenation events and corresponding decreases in water column denitrification during warm times around late 1300s, early 1700s, and the mid to late 1900s (Figure 8). These warm reoxygenation periods lasted from 100 to 150 years but reversed several times toward anoxic conditions in a matter of decades due to cooling and strengthening of the OMZ (Tems et al., 2016). During the Little Ice Age, denitrification, and bottom anoxia were the strongest during the Spörer Minimum, leading to OMZ expansion (this study; Tems et al., 2016). Subsequent periods of lower solar irradiance (Maunder and Dalton minimums) also corresponded with expanded OMZs but to a lesser degree than during the Spörer (Figure 8). The most recent stage was marked by a gradual reoxygenation during the twentieth century toward the mid-1990s, followed by a slight possible OMZ intensification over the last decade as shown by some of our proxy V and Uauth data from Alfonso and La Paz basins in accordance with the findings of Tems et al. (2016), Ontiveros-Cuadras et al. (2019) and a high-resolution U record of the last century by Choumiline et al. (2010) and Choumiline (2011).
Figure 8. Reconstruction of OMZ variability of the Eastern Pacific during the last 800 years based on proxies for water column denitrification, redox and productivity compared to major climatic forcings. Note that the solar irradiance plot is inverted for better visualization.
Climatic Drivers of Historical OMZ Fluctuations in the ETNP Over the Last Millennium
As suggested by our data from marine basins of the southwestern Gulf of California, patterns of productivity tied to upwelling are one of the key drivers of oxygenation in coastal OMZs. Importantly, most of the variability in the ETNP is indirectly controlled by short and long-term variations of solar irradiance (Bard et al., 2000). La Niña years are known to positively impact productivity in coastal regions of the Eastern Pacific by increasing upwelling rates (Thunell, 1998; Chang et al., 2013). Atmospheric dynamics, also triggered by insolation changes, play a substantial role in controlling coastal wind-driven upwelling, caused by migrations of NPH and ITCZ (Schneider et al., 2014; Metcalfe et al., 2015; Staines-Urías et al., 2015). Southward position of the NPH and ITCZ would trigger more upwelling in the Eastern Pacific (Metcalfe et al., 2015). Our research suggests—based on the chronology of Ni, V, and Uauth variability—that one of the most important controls on upwelling in the ETNP leading to OMZ fluctuations is the position of the Walker Circulation and its effect on El Niño or La Niña years on longer time scales (Figure 8). Evidence for long-term ENSO trends are manifested in PDO and/or NPGO variability in the Eastern Tropical North Pacific (Mantua and Hare, 2002; Moy et al., 2002; Di Lorenzo et al., 2010; Choumiline, 2011; Tems et al., 2016). These changes control delivery of nutrients to the marine basins of the Eastern Pacific, subsequently enhancing productivity and favoring oxygen loss (deoxygenation) of these oceanic settings (and OMZs). The mechanistic links between changes in PWC and PDO and their impacts on upwelling variability are still not clear, mainly due to difficulties in the paleoreconstruction of these phenomena.
Our record of productivity and redox over the last millennium can be compared to changes in COCOC (Composite Coccolith Proxy) presented by Beaufort and Grelaud (2017) in an attempt to reconstruct long-term variation in the PDO (Figure 8). The Southern Oscillation Index (SOI) presented by Yan et al. (2011) and Griffiths et al. (2016) is proposed to represent changes in Walker Circulation through time (Figure 8). On decadal and centennial timescales, the driving force for the major and minor fluctuations in the PWC is solar irradiance (Bard et al., 2000). Colder periods corresponding to the solar irradiance minimums—known as Dalton, Maunder, Spörer, and Wolf minimums—are also times of higher productivity and lower oxygen concentrations (seen in the Alfonso, La Paz, and Pescadero records). Importantly, denitrification was also higher during times of low insolation, as evidenced by higher sedimentary δ15N values on the Pescadero Slope (eastern side of the Gulf of California) published by Tems et al. (2016). Periods of high productivity are captured in records for Baexcess, Ni, and Uauth in Alfonso and La Paz basins during the Little Ice Age—with maximum peaks of productivity during Spörer, Maunder, and Dalton solar minimums (Bard et al., 2000; Barron and Bukry, 2007). Independent records of biogenic silica from Guaymas Basin, central Gulf of California, confirm elevated export production during colder solar episodes (Barron and Bukry, 2007). These productive times led to enhanced bottom water anoxia manifested in high V and Uauth concentrations in marine sediments of Alfonso and La Paz basins (Figures 8, 9A). The drivers of climatic change over the last millennium remain unclear, but the leading theories involve a combination of factors from orbital change to periods of strong volcanism that triggered a cooling phase during the Little Ice Age (Atwood et al., 2016). That cooling phase affected the strength of the PWC and led to lesser El Niño years.
Figure 9. Two dominating oceanographic regimes: (A) high productivity and enhanced anoxia during colder periods, (B) low productivity and reoxygenation during warm periods.
Drivers of Redox in the Eastern Pacific and Apparent Contradictions With the Glacial-Interglacial Model for OMZ Expansions
The mechanisms that control decadal to centennial oxygen variability in the ETNP remain unresolved but are likely influenced by solar forcing. That forcing not only drives migration of the ITCZ (Schneider et al., 2014) but also, more importantly, variation in the intensity of the PWC. During the Little Ice Age, solar irradiance was at its lowest over the past millennium, which strengthened Walker Circulation (Bard et al., 2000; Yan et al., 2011; Griffiths et al., 2016). The result was La Niña-like conditions, which enhanced upwelling of colder nutrient-rich waters and favored elevated productivity and reduced bottom oxygen levels (Figure 9A), as seen in Alfonso and La Paz basins and on the Pescadero Slope (Figure 8). A smaller-scale form of PWC variability is expressed in changes in the PDO and/or NPGO indices, which control rates of denitrification in OMZ waters of the Eastern Pacific (Di Lorenzo et al., 2010; Tems et al., 2016; Beaufort and Grelaud, 2017). As indicated by Tems et al. (2016) and confirmed by our findings, despite a decline in measured O2 values in the Eastern Pacific over the last few decades (Bograd et al., 2003, 2008; Levin, 2018), sedimentary records of the ETNP do not indicate increasing anoxia. Our geochemical proxy data from Alfonso from this study, Pescadero (Tems et al., 2016), and Mazatlan Margin (Ontiveros-Cuadras et al., 2019) also confirm that despite the recent warming trend and related anthropogenically triggered ocean-scale O2 decline (Bograd et al., 2003, 2008; Breitburg et al., 2018; Levin, 2018), the sedimentary records of benthic foraminifera, δ15N, Cd, Mo, Uauth, and V do not suggest a significant increase in anoxia after the mid-1800s (Figures 4, 8). Our evidence for gradually increasing dissolved oxygen levels are also supported by a decline in export production toward the beginning of the twenty-first century, evidenced by our Corg and Ni data from Alfonso and La Paz, as well as biogenic silica from Guaymas Basin (Barron and Bukry, 2007).
Despite a longstanding decrease of productivity and redox during the last century (Figure 8), a subtle regime change toward increased anoxia could have begun around the 1990s but is represented only by a few data points for the redox proxies (this study; Dean et al., 2004; IPCC, 2014; Tems et al., 2016). Remarkably, this relationship is only evidenced by δ15N, V, and Mo and not by productivity proxies, pointing at the global decline of O2 (Tems et al., 2016; Breitburg et al., 2018; Levin, 2018). Rather than reflecting global warming, reoxygenation and deoxygenation scenarios of the last millennium in the ETNP are primarily controlled by coastal productivity and dominated by long-term persistence of La Niña conditions more often occurring during stronger PWC. Deutsch et al. (2014) proposed another mechanism for the apparent OMZ contraction despite declining O2 trends in the North Pacific. The authors argued that global warming is resulting in a weakening of easterly trade winds and will continue to do so, leading to contraction of the ETNP OMZ rather than the previously predicted expansion (Deutsch et al., 2014; Tems et al., 2016). Regardless of the dominant mechanisms, our results (Figure 8) also support the interpretation that OMZ expansion occurred during cooling periods rather than being a consequence of warming—at least over the last millennium (Figure 9).
A clear distinction should be made between the dominant mechanisms that control oxygenation in the ETNP over decadal-centennial timescales (Deutsch et al., 2014; Tems et al., 2015, 2016) and those that govern long-term glacial-interglacial variability (Ganeshram and Pedersen, 1998; Nameroff et al., 2004; Hendy and Pedersen, 2006; Ivanochko et al., 2008; Pichevin et al., 2012, 2014; Cartapanis et al., 2014; Umling and Thunell, 2017). High-amplitude variability (such as the one described in our results) is superimposed over the long-term millennial-multimillenial O2 trends (Nameroff et al., 2004; Cartapanis et al., 2014; Chang et al., 2014; Moffitt et al., 2015; Praetorius et al., 2015; Tetard et al., 2017). In contrast, the long-term tendencies (e.g., transition from the Last Glacial Maximum toward the Holocene) are driven by global changes in temperature and circulation that control intermediate and deep-water oxygenation (Ganeshram and Pedersen, 1998; Hendy and Pedersen, 2006; Cartapanis et al., 2014; Praetorius et al., 2015; Lu et al., 2016; Choumiline et al., 2017; Breitburg et al., 2018; Levin, 2018). In other words, short-term decreases in anoxia during warm times in response to the weakening of the Walker Circulation and increase abundance of ENSO years can counteract long-term warming of the upper ocean, which causes stratification and reduced ventilation and O2 concentration on the global scale ocean.
Lagged Response of the Ocean to Climate Change and Future Scenarios
The effect of global warming on the world ocean is not expected to be spatially uniform. For example, latitudinal differences in the change of O2 volume have been observed, with faster declines near the equator (0–20°N) compared to higher latitudes (20–90°N) (Schmidtko et al., 2017). Some vulnerable regions have experienced a rapid response (e.g., ETSP), while others (e.g., ETNP) show a lag of several decades (Chavez et al., 2011; Howes et al., 2015; Long et al., 2016). It has been suggested that the response time of the OMZs to ocean deoxygenation caused by a decrease in O2 solubility and enhanced stratification linked to global warming (Praetorius et al., 2015) could be quicker than 25 years (Tems et al., 2016). However, the effects of climate change on OMZ oxygenation seem to only begin toward the end of the 20th century as preserved in multiple sediment records (Dean et al., 2004; Tems et al., 2016)—including ours—wherein the gradual oxygen increase in the beginning of the twentieth century shifts to a decline in O2 toward the twenty-first century. The observed lagged response in O2 trends to warming is consistent with current estimates of oceanic oxygen volume (Chavez et al., 2011; Long et al., 2016; Schmidtko et al., 2017).
Overall, there is still no consensus in the community about predicted OMZ behavior for the near future due to differences in the volume estimates among the models (Cocco et al., 2013; Long et al., 2016; Battaglia and Joos, 2018). A more recent study concludes that after the gradual increase of anoxia, the trend will be reversed globally, and the overall volume of OMZs will start declining (Fu et al., 2018). The mechanism proposed by Fu et al. (2018) involves a reduction of biological export linked to surface warming.
Lessons Learned About the Utility of Productivity and Redox Proxies From Our ETNP OMZ Records
Our findings indicate that the efficacy of productivity and redox proxies in modern OMZs depends somewhat on the supply of non-detrital trace elements via settling due to biological uptake or scavenging, as well as their preservation potential during diagenetic overprinting. We conclude that Ni, U, and V are the most reliable, while Ba and Mo should be used with caution in OMZ-type settings. Barium is significantly enriched in sediment trap material, which could elevate its proxy potential in the region (Table 5; Rodríguez-Castañeda, 2008; Choumiline, 2011). However, its long-term preservation in sedimentary records is an important issue. For La Paz Basin, we can confirm that Baexcess did not behave as a reliable productivity proxy—showing trends that differ significantly from those of the other indicators (Table 3). However, our Baexcess record in Alfonso Basin does mirror other productivity proxies (Corg and Ni) and is statistically associated with them (Figure 5). This agreement could be because our Alfonso Basin core was not exposed to severe diagenetic remobilization of sedimentary barite as much as the La Paz basin core (McManus et al., 1998; Riedinger et al., 2006). Despite the good agreement between proxy records and Baexcess in Alfonso Basin, long-term preservation of the Ba fraction remains a concern (Schoepfer et al., 2015).
In our studied cores, Ni shows good temporal agreement with both the productivity (Corg and Baexcess) and water column redox proxies (V and Uauth). In Alfonso Basin, this element served as the best indicator for productivity—mirroring an independent biogenic SiO2 record from Guaymas Basin generated by Barron and Bukry (2007) (Figure 8). Nickel does not follow Mo and Cd, which are bottom water and/or diagenetic redox proxies (Figure 5), which could be explained by a significant contribution of particulate Ni (Table 5), enriched in settling matter by biological uptake (Böning et al., 2015) and potentially related to diatom export (Twining et al., 2012).
In our sediment cores from Alfonso and La Paz basins, Mo covaries with Cd (Figure 5) and reaches high values. However, due to the low and transient concentrations of water column H2S in most modern OMZs, it is unlikely that Mo enrichments are related to bottom water euxinia in such settings (Scott and Lyons, 2012; McKay and Pedersen, 2014; Costa et al., 2018). Nevertheless, Mo is commonly elevated in ETNP sediment cores, reaching maxima several times higher than the average crustal value of 1.4 mg/kg (Wedepohl, 1995) (Figure 7; Table 5). These values are likely linked to the presence of sulfide only in the pore waters, leading to diagenetic enrichments (Zheng et al., 2000; Scholz et al., 2017; Hardisty et al., 2018). Because of this diagenetic enrichment, Mo is a less sensitive indicator of very low and transient sulfide in the water column.
We confirm the utility of V as a water column redox proxy as expressed through strong correlation with water column denitrification records reconstructed from sedimentary δ15N in the ETNP OMZ (Tems et al., 2016). Moreover, V is sensitive to O2 changes occurring not only at upper OMZ depths but also those closer to the OMZ-core, which have less pronounced redox variations. It is worth noting that, as for Uauth, there is evidence for an important contribution of V through settling matter (Rodríguez-Castañeda, 2008; Hakspiel-Segura et al., 2016; Bauer et al., 2017; Ho et al., 2018). The source of such enrichments is yet unknown, but we hypothesize that it could be related to redox processes causing authigenic V enrichment in organic matter aggregates (Ho et al., 2018).
An important methodological distinction needs to be made for proper numerical treatment of sediments with highly variable CaCO3 and biogenic silica contents. As previously suggested (Pichevin et al., 2012, 2014), proxy records can, in some cases, be significantly impacted by dilution with major components. In this study we have CaCO3 values for the studied sediment cores, so CFB fraction calculations were possible. By comparing elemental concentrations against CFB-corrected values, no substantial differences were found. Most importantly, the paleoceanographic and paleoclimatic trends described in this paper also prevailed for the carbonate-corrected values. We have no data to corroborate the dilution effect with biogenic opal, as it was not measured for the cores studied here. The consequences of the dilution effects with biogenic detritus need to be further investigated for Alfonso Basin and similar marine settings.
New Insights on Uranium Biogeochemistry in OMZs
In our study, uranium shows a consistently good correlation with Corg and the micronutrient Ni in both the deeper (mid OMZ) and shallow cores (upper OMZ) (Figure 4 and Supplementary Figure 3). We attribute this observation to seasonal variation in dissolved O2 and to the particulate authigenic U contribution during settling (Rodríguez-Castañeda, 2008; Choumiline, 2011). As particles settle deeper, there is more reaction time for dissolved U (+6) to be reduced to its particulate phase U (+4) and/or for adsorption to occur. This relationship is a consistent feature of vertical sediment trap moorings—that is, higher authigenic U enrichment in the particles collected in deeper waters (Zheng et al., 2002; Rodríguez-Castañeda, 2008; Huang and Conte, 2009; Choumiline, 2011). The extent to which Uauth enrichments are microbially mediated remains unknown, but evidence suggests that these biological interactions could be among the principal controls (Lovley et al., 1991; Cumberland et al., 2016). There is growing evidence that anaerobic metabolism can occur within aggregates of settling matter if conditions are energetically favorable (Wright et al., 2012; Lehto et al., 2014; Bianchi et al., 2018). Authigenic U formation is favored after the ferruginous (Fe reduction) and sulfidic ( reduction) chemical zones (Zheng et al., 2002; Canfield and Thamdrup, 2009). This range of redox conditions is common in marine aggregates of modern OMZs (Bianchi et al., 2018), which could explain the strong authigenic U enrichments in settling matter of the Gulf of California and California Margin (Zheng et al., 2002; Choumiline, 2011). Other marine settings, such as Saanich Inlet, also reveal the possibility of water column authigenic U enrichments, where up to one third of the estimated authigenic fraction is captured in particles in the water column (Holmden et al., 2015). It is worth mentioning that while ideal microenvironments for particulate-Uauth formation exist in strongly reducing and restricted settings in the Black and Baltic seas, as well as the Cariaco Basin (Dellwig et al., 2010; Calvert et al., 2015; Bauer et al., 2017), significant Uauth enrichments have not been detected (Anderson et al., 1989; Calvert et al., 2015). Perhaps the necessary conditions for authigenic U enrichments are favored in coastal upwelling sites, where high productivity and thus abundant particulate organic matter formation are common.
We propose an updated model for Uauth enrichments (Figure 10). This model points to a significant contribution of Uauth, beyond that linked to diagenesis, via settling matter in regions of high productivity and upwelling—such as the coastal regions of the ETNP.
In summary, the current understanding of U geochemistry in OMZs supports a significant contribution via settling particulates during times of high productivity. The actual enrichment mechanism could either be linked to localized reduction within organic-rich microenvironments or simply adsorption to organic matter. Both processes could explain the correlation with Corg and Ni and excessive U enrichments during well-ventilated times. Improved understanding of the details of U uptake will enhance its utility as a paleoproxy.
Suggestions for Future Work
This research highlights the importance of comparing high-resolution OMZ records of both modern and past oceanic conditions with global climatic trends. We have demonstrated it is possible to reconstruct high-resolution OMZ variability using a multiproxy approach. However, there are multiple caveats when interpreting deoxygenation and productivity trends because of the complex mechanisms by which trace elements are transferred to sedimentary records. We encourage particular caution for elements that demand an understanding of the relative contributions from diagenetic vs. water column processes and how specifically those inputs translate to regional patterns of redox and productivity. The importance of particulate supply of scavenged non-detrital trace elements in OMZ-type settings is clear, and further research, when possible, should include isotopic data as a source of additional mechanistic perspective.
In the modern Eastern Tropical North Pacific, marine sediments are mostly enriched in Cd, Mo, V at sites underlying severely hypoxic and anoxic waters—in contrast to shallower and oxygenated sites where these elements are present in low, mostly crustal, abundances. The exception is U, which is found to be enriched not only in sediments underlying anoxic waters but also in those of well-ventilated regions, suggesting a strong contribution of authigenic uranium through settling of marine snow aggregates.
Most of the sites located at present Oxygen Minimum Zone depths remained anoxic over the last millennium. The Eastern Tropical North Pacific Oxygen Minimum Zone expanded throughout the Little Ice Age, with the most anoxic events during the Spörer, Maunder, and Dalton insolation minima. Reoxygenation in the Eastern Tropical North Pacific occurred during warm periods around late 1300s, early 1700s, and late 1900s AD.
There is no clear effect of the temperature increase of twentieth century on the intensification of the Eastern Tropical North Pacific Oxygen Minimum Zone. The prevailing trend shows an increase in dissolved oxygen from the 1900s AD toward present time, with a subtle decrease in the beginning of the twenty-first century.
Redox and productivity changes reconstructed for the last millennium from marine basins of the Eastern Tropical North Pacific indicate a strong relationship with insolation fluctuations and the Pacific Walker Circulation that drive upwelling rates, productivity, and oceanic anoxia at Oxygen Minimum Zone depths.
The most reliable redox proxies for water column oxygenation in Oxygen Minimum Zones are vanadium and uranium. Uranium also shows strong potential to serve as a productivity proxy because it can be exported from the water column to sediments through settling particulate matter during periods of high primary productivity. Molybdenum is the least useful proxy for Oxygen Minimum Zone oxygenation because it is not sensitive to small redox variations and can be enriched diagenetically when pore water sulfide is present.
KC wrote the paper including comments by LP-C, AG, SB, and TL. Project planning logistics were performed by KC, TL, and LP-C. Samples used in this research were collected, preserved, and handled by KC and LP-C. Analytical measurements were performed by KC, SB, and AG.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors acknowledge financial support granted by UC MEXUS - CONACYT Collaborative Research Grant “High resolution geochemical reconstructions of recent climate and oxygenation history in La Paz Bay, Gulf of California” that covered travel and analytical costs. The authors express gratitude to Coordinación de Plataformas Oceanográficas (UNAM, Mexico) for funding the PALEOMAR-I campaign and the crew members of R/V “El Puma” for exceptional support during the field stage. KC was thankful for the UC MEXUS – CONACYT Doctoral Fellowship for a fully funded scholarship during his Ph.D. studies and for the 2016 UC MEXUS Dissertation Grant award for analytical expenses. Partial support was received by TL, KC, and SB from the NASA Astrobiology Institute, and AG from the USDA NIFA Hatch program (project # CA-R-ENS-5120-H).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2019.00237/full#supplementary-material
Algeo, T. J., and Lyons, T. W. (2006). Mo-total organic carbon covariation in modern anoxic marine environments: implications for analysis of paleoredox and paleohydrographic conditions. Paleoceanogr. Paleoclimatol. 21:PA1016 doi: 10.1029/2004PA001112
Anderson, R. F., Fleisher, M. Q., and Lehuray, A. P. (1989). Concentration, oxidation state, and particulate flux of uranium in the Black Sea. Geochim. Cosmochim. Acta 53, 2215–2224. doi: 10.1016/0016-7037(89)90345-1
Atwood, A. R., Wu, E., Frierson, D. M. W., Battisti, D. S., and Sachs, J. P. (2016). Quantifying climate forcings and feedbacks over the Last Millennium in the CMIP5-PMIP3 Models. J. Clim. 29, 1161–1178. doi: 10.1175/JCLI-D-15-0063.1
Barron, J. A., and Bukry, D. (2007). Solar forcing of Gulf of California climate during the past 2000 yr suggested by diatoms and silicoflagellates. Mar. Micropaleontol. 62, 115–139. doi: 10.1016/j.marmicro.2006.08.003
Battaglia, G., and Joos, F. (2018). Hazards of decreasing marine oxygen: the near-term and millennial-scale benefits of meeting the Paris climate targets. Earth Syst. Dyn. 9, 797–816. doi: 10.5194/esd-9-797-2018
Bauer, S., Blomqvist, S., and Ingri, J. (2017). Distribution of dissolved and suspended particulate molybdenum, vanadium, and tungsten in the Baltic Sea. Mar. Chem. 196, 135–147. doi: 10.1016/j.marchem.2017.08.010
Beaufort, L., and Grelaud, M. (2017). A 2700-year record of ENSO and PDO variability from the Californian margin based on coccolithophore assemblages and calcification. Prog. Earth Planet. Sci. 4:5. doi: 10.1186/s40645-017-0123-z
Bianchi, D., Weber, T. S., Kiko, R., and Deutsch, C. (2018). Global niche of marine anaerobic metabolisms expanded by particle microenvironments. Nat. Geosci. 11, 263–268. doi: 10.1038/s41561-018-0081-0
Bianchi, T. S., Schreiner, K. M., Smith, R. W., Burdige, D. J., Woodard, S., and Conley, D. J. (2016). Redox effects on organic matter storage in coastal sediments during the holocene: a biomarker/proxy perspective. Annu. Rev. Earth Planet. Sci. 44, 295–319. doi: 10.1146/annurev-earth-060614-105417
Bograd, S. J., Castro, C. G., Di Lorenzo, E., Palacios, D. M., Bailey, H., Gilly, W., et al. (2008). Oxygen declines and the shoaling of the hypoxic boundary in the California Current. Geophys. Res. Lett. 35:L12607. doi: 10.1029/2008GL034185
Bograd, S. J., Checkley, D. A., and Wooster, W. S. (2003). CalCOFI: a half century of physical, chemical, and biological research in the California Current System. Deep Sea Res. Part II Top. Stud. Oceanogr. 50, 2349–2353. doi: 10.1016/S0967-0645(03)00122-X
Borreggine, M., Myhre, S. E., Mislan, K. A. S., Deutsch, C., and Davis, C. V. (2017). A database of paleoceanographic sediment cores from the North Pacific, 1951-2016. Earth Syst. Sci. Data 9, 739–749. doi: 10.5194/essd-9-739-2017
Breitburg, D., Levin, L. A., Oschlies, A., Grégoire, M., Chavez, F. P., Conley, D. J., et al. (2018). Declining oxygen in the global ocean and coastal waters. Science 359:eaam7240. doi: 10.1126/science.aam7240
Calvert, S. E., Bustin, R. M., and Ingall, E. D. (1996). Influence of water column anoxia and sediment supply on the burial and preservation of organic carbon in marine shales. Geochimica Cosmochimica Acta 60, 1577–1593.
Canfield, D. E., and Thamdrup, B. (2009). Towards a consistent classification scheme for geochemical environments, or, why we wish the term 'suboxic' would go away. Geobiology 7, 385–392. doi: 10.1111/j.1472-4669.2009.00214.x
Cartapanis, O., Tachikawa, K., Romero, O. E., and Bard, E. (2014). Persistent millennial-scale link between Greenland climate and northern Pacific Oxygen Minimum Zone under interglacial conditions. Clim. Past 10, 405–418. doi: 10.5194/cp-10-405-2014
Chang, A. S., Bertram, M. A., Ivanochko, T., Calvert, S. E., Dallimore, A., and Thomson, R. E. (2013). Annual record of particle fluxes, geochemistry and diatoms in Effingham Inlet, British Columbia, Canada, and the impact of the 1999 La Niña event. Mar. Geol. 337, 20–34. doi: 10.1016/j.margeo.2013.01.003
Chang, A. S., Pedersen, T. F., and Hendy, I. L. (2014). Effects of productivity, glaciation, and ventilation on late Quaternary sedimentary redox and trace element accumulation on the Vancouver Island margin, western Canada. Paleoceanography 29, 730–746. doi: 10.1002/2013PA002581
Chappaz, A., Lyons, T. W., Gregory, D. D., Reinhard, C. T., Gill, B. C., Li, C., et al. (2014). Does pyrite act as an important host for molybdenum in modern and ancient euxinic sediments? Geochim. Cosmochim. Acta 126, 112–122. doi: 10.1016/j.gca.2013.10.028
Chavez, F. P., Messié, M., and Pennington, J. T. (2011). Marine primary production in relation to climate variability and change. Ann. Rev. Mar. Sci. 3, 227–260. doi: 10.1146/annurev.marine.010908.163917
Choumiline, K. (2011). Geoquímica de la materia particulada en hundimiento y de los sedimentos recientes de Cuenca Alfonso, Bahía de La Paz. (Master's thesis). Centro Interdisciplinario de Ciencias Marinas – Instituto Politécnico Nacional, La Paz, Mexico.
Choumiline, K., Lyons, T. W., Pérez-Cruz, L., Carriquiry, J. D., Raiswell, R., and Beaufort, L. (2017). “Sensitivity of redox proxies to rapid variability in Oxygen Minimum Zones. Goldschmidt 2017,” in Abstract Retrieved from Goldschmidt Conference Archive (Paris).
Choumiline, K., Rodríguez-Castañeda, A. P., Silverberg, N., Shumilin, E., Aguirre-Bahena, F., Sapozhnikov, D., et al. (2010). “Arsenic and uranium in the settling particulate matter and sediments of Alfonso Basin, La Paz Bay,” in Proceedings of the 13th International Conference on Water-Rock Interaction. Reviewed Extended Abstract (Guanajuato).
Cocco, V., Joos, F., Steinacher, M., Frölicher, T. L., Bopp, L., Dunne, J., et al. (2013). Oxygen and indicators of stress for marine life in multi-model global warming projections. Biogeosciences 10, 1849–1868. doi: 10.5194/bg-10-1849-2013
Conte, M. H., Carter, A. M., Koweek, D. A., Huang, S., and Weber, J. C. (2019). The elemental composition of the deep particle flux in the Sargasso Sea. Chem. Geol. 511, 279–313. doi: 10.1016/j.chemgeo.2018.11.001
Coria-Monter, E., Monreal-Gómez, M. A., Salas-de-León, D. A., Aldeco-Ramírez, J., and Merino-Ibarra, M. (2014). Differential distribution of diatoms and dinoflagellates in a cyclonic eddy confined in the Bay of La Paz, Gulf of California. J. Geophys. Res. 119, 6258–6268. doi: 10.1002/2014JC009916
Costa, K. M., Anderson, R. F., McManus, J. F., Winckler, G., Middleton, J. L., and Langmuir, C. H. (2018). Trace element (Mn, Zn, Ni, V) and authigenic uranium (aU) geochemistry reveal sedimentary redox history on the Juan de Fuca Ridge, North Pacific Ocean. Geochim. Cosmochim. Acta 236, 79–98. doi: 10.1016/j.gca.2018.02.016
Cumberland, S. A., Douglas, G., Grice, K., and Moreau, J. W. (2016). Uranium mobility in organic matter-rich sediments: a review of geological and geochemical processes. Earth Sci. Rev. 159, 160–185. doi: 10.1016/j.earscirev.2016.05.010
Dean, W., Pride, C., and Thunell, R. (2004). Geochemical cycles in sediments deposited on the slopes of the Guaymas and Carmen Basins of the Gulf of California over the last 180 years. Quat. Sci. Rev. 23, 1817–1833. doi: 10.1016/j.quascirev.2004.03.010
Dean, W. E. (2007). Sediment geochemical records of productivity and oxygen depletion along the margin of western North America during the past 60,000 years: teleconnections with Greenland Ice and the Cariaco Basin. Quat. Sci. Rev. 26, 98–114. doi: 10.1016/j.quascirev.2006.08.006
Dellwig, O., Leipe, T., März, C., Glockzin, M., Pollehne, F., Schnetger, B., et al. (2010). A new particulate Mn–Fe–P-shuttle at the redoxcline of anoxic basins. Geochim. Cosmochim. Acta 74, 7100–7115. doi: 10.1016/j.gca.2010.09.017
Deutsch, C., Berelson, W., Thunell, R., Weber, T., Tems, C., McManus, J., et al. (2014). Oceanography. Centennial changes in North Pacific anoxia linked to tropical trade winds. Science 345, 665–668. doi: 10.1126/science.1252332
Di Lorenzo, E., Cobb, K. M., Furtado, J. C., Schneider, N., Anderson, B. T., Bracco, A., et al. (2010). Central Pacific El Niño and decadal climate change in the North Pacific Ocean. Nat. Geosci. 3, 762–765. doi: 10.1038/ngeo984
Di Lorenzo, E., Schneider, N., Cobb, K. M., Franks, P. J. S., Chhak, K., Miller, A. J., et al. (2008). North Pacific Gyre Oscillation links ocean climate and ecosystem change. Geophys. Res. Lett. 35:L08607. doi: 10.1029/2007GL032838
Dickens, G. R., Koelling, M., Smith, D. C., and Schnieders, L. (2007). Rhizon sampling of pore waters on scientific drilling expeditions: an example from the IODP Expedition 302, Arctic Coring Expedition (ACEX). Sci. Drilling 4, 22–25. doi: 10.5194/sd-4-22-2007
Dupont, C. L., Buck, K. N., Palenik, B., and Barbeau, K. (2010). Nickel utilization in phytoplankton assemblages from contrasting oceanic regimes. Deep Sea Res. Part I Oceanogr. Res. Pap. 57, 553–566. doi: 10.1016/j.dsr.2009.12.014
Durand, A., Chase, Z., Townsend, A. T., Noble, T., Panietz, E., and Goemann, K. (2016). Improved methodology for the microwave digestion of carbonate-rich environmental samples. Int. J. Environ. Anal. Chem. 96, 119–136. doi: 10.1080/03067319.2015.1137904
Durazo-Arvizu, R., and Gaxiola-Castro, G. (2010). Dinámica del ecosistema pelágico frente a Baja California, 1997–2007. Diez años de investigaciones mexicanas de la Corriente de California. Mexico, Distrito Federal: INE/CICESE/UABC/SEMARNAT.
Eagle, M., Paytan, A., Arrigo, K. R., van Dijken, G., and Murray, R. W. (2003). A comparison between excess barium and barite as indicators of carbon export. Paleoceanogr. Paleoclimatol. 18:1021. doi: 10.1029/2002PA000793
Erickson, B. E., and Helz, G. R. (2000). Molybdenum(VI) speciation in sulfidic waters:: Stability and lability of thiomolybdates. Geochim. Cosmochim. Acta 64, 1149–1158. doi: 10.1016/S0016-7037(99)00423-8
Foster, I. D. L., and Walling, D. E. (1994). Using reservoir deposits to reconstruct changing sediment yields and sources in the catchment of the Old Mill Reservoir, South Devon, UK, over the past 50 years. Hydrol. Sci. J. 39, 347–368. doi: 10.1080/02626669409492755
Fu, W., Primeau, F., Moore, J. K., Lindsay, K., and Randerson, J. T. (2018). Reversal of increasing tropical ocean hypoxia trends with sustained climate warming. Global Biogeochem. Cycles 32, 551–564. doi: 10.1002/2017GB005788
Ganeshram, R. S., and Pedersen, T. F. (1998). Glacial-interglacial variability in upwelling and bioproductivity off NW Mexico: implications for quaternary paleoclimate. Paleoceanography Paleoclimatol. 13, 634–645. doi: 10.1029/98PA02508
Garcia, H. E., Locarnini, R. A., Boyer, T. P., and Antonov, J. I. (2006). World Ocean Atlas 2005, Volume 3: Dissolved Oxygen, Apparent Oxygen Utilization, and Oxygen Saturation. S. Levitus, Ed. NOAA Atlas NESDIS 63, U.S. Government Printing Office, Washington, D.C., 342 pp.
Gilly, W. F., Beman, J. M., Litvin, S. Y., and Robison, B. H. (2013). Oceanographic and biological effects of shoaling of the oxygen minimum zone. Ann. Rev. Mar. Sci. 5, 393–420. doi: 10.1146/annurev-marine-120710-100849
Griffiths, M. L., Kimbrough, A. K., Gagan, M. K., Drysdale, R. N., Cole, J. E., Johnson, K. R., et al. (2016). Western Pacific hydroclimate linked to global climate variability over the past two millennia. Nat. Commun. 7:11719. doi: 10.1038/ncomms11719
Hakspiel-Segura, C., Martínez-López, A., Pinedo-González, P., Verdugo-Díaz, G., and Acevedo-Acosta, J. D. (2016). Composition of metals in suspended particulate matter of Alfonso basin, southern Gulf of California. Reg. Stud. Mar. Sci. 3, 144–153. doi: 10.1016/j.rsma.2015.07.001
Hardisty, D. S., Lyons, T. W., Riedinger, N., Isson, T. T., Owens, J. D., Aller, R. C., et al. (2018). An evaluation of sedimentary molybdenum and iron as proxies for pore fluid paleoredox conditions. Am. J. Sci. 318, 527–556. doi: 10.2475/05.2018.04
Helz, G. R., Miller, C. V., Charnock, J. M., Mosselmans, F. W., Pattrick, R., Garner, D., et al. (1996). Mechanism of molybdenum removal from the sea and its concentration in black shales: EXAFS evidence. Geochim. Cosmochim. Acta 60, 3631–3642. doi: 10.1016/0016-7037(96)00195-0
Hendy, I. L., and Kennett, J. P. (2000). Dansgaard-Oeschger cycles and the California current system: Planktonic foraminiferal response to rapid climate change in Santa Barbara Basin, Ocean Drilling Program Hole 893A. Paleoceanogr. Paleoclimatol. 15, 30–42. doi: 10.1029/1999PA000413
Ho, P., Lee, J.-M., Heller, M. I., Lam, P. J., and Shiller, A. M. (2018). The distribution of dissolved and particulate Mo and V along the U.S. GEOTRACES East Pacific Zonal Transect (GP16): the roles of oxides and biogenic particles in their distributions in the oxygen deficient zone and the hydrothermal plume. Mar. Chem. 201, 242–255. doi: 10.1016/j.marchem.2017.12.003
Hofmann, A. F., Peltzer, E. T., Walz, P. M., and Brewer, P. G. (2011). Hypoxia by degrees: establishing definitions for a changing ocean. Deep Sea Res. Part I Oceanogr. Res. Pap. 58, 1212–1226. doi: 10.1016/j.dsr.2011.09.004
Holmden, C., Amini, M., and Francois, R. (2015). Uranium isotope fractionation in Saanich Inlet: a modern analog study of a paleoredox tracer. Geochim. Cosmochim. Acta 153, 202–215. doi: 10.1016/j.gca.2014.11.012
Hoogakker, B. A. A., Lu, Z., Umling, N., Jones, L., Zhou, X., Rickaby, R. E. M., et al. (2018). Glacial expansion of oxygen-depleted seawater in the eastern tropical Pacific. Nature 562, 410–413. doi: 10.1038/s41586-018-0589-x
Howes, E. L., Joos, F., Eakin, C. M., and Gattuso, J.-P. (2015). An updated synthesis of the observed and projected impacts of climate change on the chemical, physical and biological processes in the oceans. Front. Mar. Sci. 2:36. doi: 10.3389/fmars.2015.00036
Huang, S., and Conte, M. H. (2009). Source/process apportionment of major and trace elements in sinking particles in the Sargasso sea. Geochim. Cosmochim. Acta 73, 65–90. doi: 10.1016/j.gca.2008.08.023
IPCC (2014). Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Core Writing Team, R.K. Pachauri and L.A. Meyer (eds.)]. IPCC, Geneva, Switzerland, 151pp.
Ivanochko, T. S., Calvert, S. E., Southon, J. R., Enkin, R. J., Baker, J., Dallimore, A., et al. (2008). Determining the post-glacial evolution of a northeast Pacific coastal fjord using a multiproxy geochemical approach. Can. J. Earth Sci. 45, 1331–1344. doi: 10.1139/E08-030
Ivanochko, T. S., and Pedersen, T. F. (2004). Determining the influences of Late Quaternary ventilation and productivity variations on Santa Barbara Basin sedimentary oxygenation: a multi-proxy approach. Quat. Sci. Rev. 23, 467–480. doi: 10.1016/j.quascirev.2003.06.006
Lavín, M. F., Castro, R., Beier, E., and Godínez, V. M. (2013). Mesoscale eddies in the southern Gulf of California during summer: characteristics and interaction with the wind stress. J. Geophys. Res. Oceans 118, 1367–1381. doi: 10.1002/jgrc.20132
Lehto, N., Glud, R. N.Á, Norð*i, G., Zhang, H., and Davison, W. (2014). Anoxic microniches in marine sediments induced by aggregate settlement: biogeochemical dynamics and implications. Biogeochemistry 119, 307–327. doi: 10.1007/s10533-014-9967-0
Lu, Z., Hoogakker, B. A. A., Hillenbrand, C.-D., Zhou, X., Thomas, E., et al. (2016). Oxygen depletion recorded in upper waters of the glacial Southern Ocean. Nat. Commun. 7:11146. doi: 10.1038/ncomms11146
Lyons, T. W., and Severmann, S. (2006). A critical look at iron paleoredox proxies: new insights from modern euxinic marine basins. Geochim. Cosmochim. Acta 70, 5698–5722. doi: 10.1016/j.gca.2006.08.021
McManus, J., Berelson, W. M., Coale, K. H., Johnson, K. S., and Kilgore, T. E. (1997). Phosphorus regeneration in continental margin sediments. Geochim. Cosmochim. Acta 61, 2891–2907. doi: 10.1016/S0016-7037(97)00138-5
McManus, J., Berelson, W. M., Klinkhammer, G. P., Johnson, K. S., Coale, K. H., Anderson, R. F., et al. (1998). Geochemistry of barium in marine sediments: implications for its use as a paleoproxy. Geochim. Cosmochim. Acta 62, 3453–3473. doi: 10.1016/S0016-7037(98)00248-8
McManus, J., Berelson, W. M., Severmann, S., Poulson, R. L., Hammond, D. E., Klinkhammer, G. P., et al. (2006). Molybdenum and uranium geochemistry in continental margin sediments: paleoproxy potential. Geochim. Cosmochim. Acta 70, 4643–4662. doi: 10.1016/j.gca.2006.06.1564
Metcalfe, S. E., Barron, J. A., and Davies, S. J. (2015). The Holocene history of the North American Monsoon: ‘known knowns' and ‘known unknowns' in understanding its spatial and temporal complexity. Quat. Sci. Rev. 120, 1–27. doi: 10.1016/j.quascirev.2015.04.004
Moffitt, S. E., Moffitt, R. A., Sauthoff, W., Davis, C. W., Hewett, K., and Hill, T. M. (2015). Paleoceanographic insights on recent oxygen minimum zone expansion: lessons for modern oceanography. PLoS ONE 10:e0115246. doi: 10.1371/journal.pone.0115246
Monreal-Gómez, M. A., Molina-Cruz, A., and Salas-de-León, D. A. (2001). Water masses and cyclonic circulation in Bay of la Paz, Gulf of California, during June 1998. J. Mar. Syst. 30, 305–315. doi: 10.1016/S0924-7963(01)00064-1
Moy, C. M., Seltzer, G. O., Rodbell, D. T., and Anderson, D. M. (2002). Variability of El Niño/Southern Oscillation activity at millennial timescales during the Holocene epoch. Nature 420, 162–165. doi: 10.1038/nature01194
Nameroff, T. J., Calvert, S. E., and Murray, J. W. (2004). Glacial-interglacial variability in the eastern tropical North Pacific oxygen minimum zone recorded by redox-sensitive trace metals. Paleoceanogr. Paleoclimatol. 19:PA1010. doi: 10.1029/2003PA000912
Nava-Sánchez, E. H., Gorsline, D. S., and Molina-Cruz, A. (2001). The Baja California peninsula borderland: structural and sedimentological characteristics. Sediment. Geol. 144, 63–82. doi: 10.1016/S0037-0738(01)00135-X
Ontiveros-Cuadras, J. F., Thunell, R., Ruiz-Fernández, A. C., Benitez-Nelson, C., Machain-Castillo, M. L., Tappa, E., et al. (2019). Centennial OMZ changes in the NW Mexican Margin from geochemical and foraminiferal sedimentary records. Cont. Shelf Res. 176, 64–75. doi: 10.1016/j.csr.2019.02.009
Oschlies, A., Schulz, K. G., Riebessel, U., and Schmittner, A. (2008). Simulated 21st century's increase in oceanic suboxia by CO2-enhanced biotic carbon export. Global Biogeochem. Cycles 22:GB4008. doi: 10.1029/2007GB003147
Pérez-Cruz, L. (2013). Hydrological changes and paleoproductivity in the Gulf of California during middle and late Holocene and their relationship with ITCZ and North American Monsoon variability. Quat. Res. 79, 138–151. doi: 10.1016/j.yqres.2012.11.007
Pichevin, L., Ganeshram, R. S., Reynolds, B. C., Prahl, F., Pedersen, T. F., Thunell, R., et al. (2012). Silicic acid biogeochemistry in the Gulf of California: Insights from sedimentary Si isotopes. Paleoceanogr. Paleoclimatol. 27:PA2201. doi: 10.1029/2011PA002237
Pichevin, L. E., Ganeshram, R. S., Geibert, W., Thunell, R., and Hinton, R. (2014). Silica burial enhanced by iron limitation in oceanic upwelling margins. Nat. Geosci. 7, 541–546. doi: 10.1038/ngeo2181
Praetorius, S. K., Mix, A. C., Walczak, M. H., Wolhowe, M. D., Addison, J. A., and Prahl, F. G. (2015). North Pacific deglacial hypoxic events linked to abrupt ocean warming. Nature 527, 362–366. doi: 10.1038/nature15753
Reck, B. K., Müller, D. B., Rostkowski, K., and Graedel, T. E. (2008). Anthropogenic nickel cycle: insights into use, trade, and recycling. Environ. Sci. Technol. 42, 3394–3400. doi: 10.1021/es072108l
Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Ramsey, C. B., et al. (2013). IntCal13 and Marine13 Radiocarbon Age Calibration Curves 0–50,000 Years cal BP. Radiocarbon 55, 1869–1887. doi: 10.2458/azu_js_rc.55.16947
Riedinger, N., Brunner, B., Krastel, S., Arnold, G. L., Wehrmann, L. M., Formolo, M. J., et al. (2017). Sulfur cycling in an iron oxide-dominated, dynamic marine depositional system: the Argentine continental margin. Front. Earth Sci. 5:33. doi: 10.3389/feart.2017.00033
Riedinger, N., Formolo, M. J., Lyons, T. W., Henkel, S., Beck, A., and Kasten, S. (2014). An inorganic geochemical argument for coupled anaerobic oxidation of methane and iron reduction in marine sediments. Geobiology 12, 172–181. doi: 10.1111/gbi.12077
Riedinger, N., Kasten, S., Gröger, J., Franke, C., and Pfeifer, K. (2006). Active and buried authigenic barite fronts in sediments from the Eastern Cape Basin. Earth Planet. Sci. Lett. 241, 876–887. doi: 10.1016/j.epsl.2005.10.032
Rodríguez-Castañeda, A. P. (2008). Variación de flujos de los elementos particulados en Cuenca Alfonso, Bahía de La Paz, en el periodo 2002–2005. [Ph.D. dissertation]. [La Paz, Mexico]: Centro Interdisciplinario de Ciencias Marinas – Instituto Politécnico Nacional
Ruiz-Fernández, A. C., Hillaire-Marcel, C., De Vernal, A., Machain-Castillo, M. L., Vásquez, L., Ghaleb, B., et al. (2009). Changes of coastal sedimentation in the Gulf of Tehuantepec, South Pacific Mexico, over the last 100 years from short-lived radionuclide measurements. Estuar. Coast. Shelf Sci. 82, 525–536. doi: 10.1016/j.ecss.2009.02.019
Schoepfer, S. D., Shen, J., Wei, H., Tyson, R. V., Ingall, E., and Algeo, T. J. (2015). Total organic carbon, organic phosphorus, and biogenic barium fluxes as proxies for paleomarine productivity. Earth Sci. Rev. 149, 23–52. doi: 10.1016/j.earscirev.2014.08.017
Scholz, F., Hensen, C., Noffke, A., Rohde, A., Liebetrau, V., and Wallmann, K. (2011). Early diagenesis of redox-sensitive trace metals in the Peru upwelling area – response to ENSO-related oxygen fluctuations in the water column. Geochim. Cosmochim. Acta 75, 7257–7276. doi: 10.1016/j.gca.2011.08.007
Scholz, F., Siebert, C., Dale, A. W., and Frank, M. (2017). Intense molybdenum accumulation in sediments underneath a nitrogenous water column and implications for the reconstruction of paleo-redox conditions based on molybdenum isotopes. Geochim. Cosmochim. Acta 213, 400–417. doi: 10.1016/j.gca.2017.06.048
Scott, C., and Lyons, T. W. (2012). Contrasting molybdenum cycling and isotopic properties in euxinic versus non-euxinic sediments and sedimentary rocks: refining the paleoproxies. Chem. Geol. 324–325, 19–27. doi: 10.1016/j.chemgeo.2012.05.012
Silverberg, N., Aguirre-Bahena, F., and Mucci, A. (2014). Time-series measurements of settling particulate matter in Alfonso Basin, La Paz Bay, southwestern Gulf of California. Cont. Shelf Res. 84, 169–187. doi: 10.1016/j.csr.2014.05.005
Staines-Urías, F., González-Yajimovich, O., and Beaufort, L. (2015). Reconstruction of past climate variability and ENSO-like fluctuations in the southern Gulf of California (Alfonso Basin) since the last glacial maximum. Q. Res. 83, 488–501. doi: 10.1016/j.yqres.2015.03.007
Tems, C. E., Berelson, W. M., and Prokopenko, M. G. (2015). Particulate δ15N in laminated marine sediments as a proxy for mixing between the California Undercurrent and the California Current: a proof of concept. Geophys. Res. Lett. 42, 419–427. doi: 10.1002/2014GL061993
Tems, C. E., Berelson, W. M., Thunell, R., Tappa, E., Xu, X., Khider, D., et al. (2016). Decadal to centennial fluctuations in the intensity of the eastern tropical North Pacific oxygen minimum zone during the last 1200 years. Paleoceanography Paleoclimatol. 31, 1138–1151. doi: 10.1002/2015PA002904
Thunell, R. (1998). Seasonal and annual variability in particle fluxes in the Gulf of California. A response to climate forcing. Deep Sea Res. Part I Oceanogr. Res. Pap. 45, 2059–2083. doi: 10.1016/S0967-0637(98)00053-3
Trasviña-Castro, A., Gutierrez De Velasco, G., Valle-Levinson, A., González-Armas, R., Muhlia, A., and Cosio, M. A. (2003). Hydrographic observations of the flow in the vicinity of a shallow seamount top in the Gulf of California. Estuar. Coast. Shelf Sci. 57, 149–162. doi: 10.1016/S0272-7714(02)00338-4
Tribovillard, N., Bout-Roumazeilles, V., Algeo, T., Lyons, T. W., Sionneau, T., Montero-Serrano, J. C., et al. (2008). Paleodepositional conditions in the Orca Basin as inferred from organic matter and trace metal contents. Mar. Geol. 254, 62–72. doi: 10.1016/j.margeo.2008.04.016
Yan, H., Sun, L., Wang, Y., Huang, W., Qiu, S., and Yang, C. (2011). A record of the Southern Oscillation Index for the past 2,000 years from precipitation proxies. Nat. Geosci. 4, 611–614. doi: 10.1038/ngeo1231
Zheng, Y., Anderson, R. B., Van Geen, A., and Fleisher, M. Q. (2002). Preservation of particulate non-lithogenic uranium in marine sediments. Geochim. Cosmochim. Acta 66, 3085–3092. doi: 10.1016/S0016-7037(01)00632-9
Keywords: Oxygen Minimum Zones (OMZs), Eastern Tropical North Pacific (ETNP), primary productivity, deoxygenation, hypoxia, paleoredox, redox-sensitive trace elements, authigenic uranium
Citation: Choumiline K, Pérez-Cruz L, Gray AB, Bates SM and Lyons TW (2019) Scenarios of Deoxygenation of the Eastern Tropical North Pacific During the Past Millennium as a Window Into the Future of Oxygen Minimum Zones. Front. Earth Sci. 7:237. doi: 10.3389/feart.2019.00237
Received: 20 September 2018; Accepted: 28 August 2019;
Published: 26 September 2019.
Edited by:Babette Hoogakker, Heriot-Watt University, United Kingdom
Reviewed by:Nicolas Tribovillard, Lille University of Science and Technology, France
Laetitia Pichevin, University of Edinburgh, United Kingdom
Copyright © 2019 Choumiline, Pérez-Cruz, Gray, Bates and Lyons. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Konstantin Choumiline, email@example.com