ALBERT

Description: Vegetation biomass is a globally important climate-relevant terrestrial carbon pool. Landsat, Sentinel-2 and Sentinel-1 satellite missions provide a landscape-level opportunity to upscale tundra vegetation communities and biomass in high latitude terrestrial environments. We assessed the applicability of landscape-level remote sensing for the low Arctic Lena Delta region in Northern Yakutia, Siberia, Russia. The Lena Delta is the largest delta in the Arctic and is located North of the treeline and the 10 °C July isotherm at 72° Northern Latitude in the Laptev Sea region. During the LENA2018 expedition, we set up plots for plant projective cover and Above Ground Biomass (AGB) and sampled shrubs for shrub-ring analyses. AGB is providing the magnitude of the carbon flux, whereas stand age is irreplaceable to provide the cycle rate. AGB data and shrub age data clearly show a separation between i) low disturbance landscape types with dominant AGB moss contribution, but always low vascular plant AGB (〈0.5 kg m-2) characterised by old shrubs of several decades of stand age versus ii) a much higher vascular plant AGB contribution (〉 0.5 kg m-2) with only young shrubs in high disturbance regimes. The low disturbance regimes are represented on the Holocene and Pleistocene delta terraces in form of azonal polygonal tundra complexes and softly dissected valleys with zonal tussock tundra. In contrast, the high disturbance regimes are sites of thermo-erosion such as along thermo-erosional valleys and on floodplains. We upscaled AGB and above ground carbon pool ages using a Sentinel-2 satellite acquisition from early August 2018. We classified via classification training using Elementary Sampling Units that are the 30 m x 30 m vegetation field plots. We then used the land cover classes and grouped them according to their settings either in high disturbance or low disturbance regimes with each associated AGB value ranges and shrub age regimes. We also evaluated circum-Arctic harmonized ESA GlobPermafrost land cover and vegetation height remote sensing products covering subarctic to Arctic land cover types for the central Lena Delta. The products are freely available and published in the PANGAEA data repository under https://doi.org/10.1594/PANGAEA.897916 and https://doi.org/10.1594/PANGAEA.897045. ESA GlobPermafrost land cover and vegetation height remote sensing products and our Sentinel-2 derived AGB product for the central Lena Delta shows realistic spatial patterns of landcover classes and biomass distribution at landscape level. However, in all products, the high biomass patches of high shrubs in the tundra landscape could not spatially be resolved as they are confined to patchy and linear distribution, not representing large enough areas suitable for upscaling. We found that high disturbance regimes with linked high and rapid AGB fluxes are distributed mainly on the floodplains and as patches along thermoerosioal features, e.g. valleys. Whereas the low disturbance landscapes on Yedoma upland tundra and Holocene terraces occur with larger area coverage representing decades slower and in magnitude smaller AGB fluxes.

Description: As air temperatures rise and sea ice cover declines in the Arctic, permafrost coastal cliffs thaw more rapidly and wave energy rises. Thus, as the open water season continues to lengthen, climate change triggers a large part of the Arctic shoreline to become increasingly vulnerable to erosion. Arctic erosion supplies nutrient-laden and carbon-rich sediment into nearshore ecosystems. A retreating coastline also has consequences for residential, cultural, and industrial infrastructure. Despite its importance, erosion is currently neglected in global climate models, and existing physics-based numerical models of Arctic shoreline erosion are too complex and regionally-focused to be applied on a pan-Arctic scale. Here, we apply our simplified numerical erosion model, ArcticBeach v1.0, to the entire Arctic coastline. ArcticBeach v1.0 has previously been shown to simulate retreat rates at two sites that differ substantially in their main mechanisms of retreat (sub-aerial erosion/thaw slumping versus notch/block erosion). The model uses heat and sediment volume balances in order to predict horizontal cliff retreat and vertical erosion of a fronting beach. It contains an erosion module that uses empirical equations to estimate cross-shore sediment transport, coupled to a storm surge module forced by wind. We present Arctic maps of regional variation in trends in 2-meter air temperature, sea ice concentration, and wind speed.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 22(3), (2021): e2020GC009472, https://doi.org/10.1029/2020GC009472.

Description: Carbonatite volcanism remains poorly understood compared to silicic volcanism due to the scarcity of carbonatite volcanoes worldwide and because volcanic H2O and CO2—major components in carbonatite volcanic systems—are not well preserved in the rock record. To further our understanding of carbonatite genesis, we utilize the non-traditional thallium (Tl) isotope system in Khanneshin carbonatites in Afghanistan. These carbonatites contain 250–30,000 ng/g Tl and have ε205Tl values (−4.6 to +4.6) that span much of the terrestrial igneous range. We observe that δ18OVSMOW (+8.6‰ to +23.5‰) correlates positively with δ13CVPDB (−4.6‰ to +3.5‰) and ε205Tl up to δ18O = 15‰. Rayleigh fractionation of calcite from an immiscible CO2-H2O fluid with a mantle-like starting composition can explain the δ18O and δ13C—but not ε205Tl—trends. Biotite fractionates Tl isotopes in other magmatic settings, so we hypothesize that a Tl-rich hydrous brine caused potassic metasomatism (i.e., biotite fenitization) of wall rock that increased the ε205Tl of the residual magma-fluid reservoir. Our results imply that, in carbonatitic volcanic systems, simultaneous igneous differentiation and potassic metasomatism increase ε205Tl, δ18O, δ13C, and light rare earth element concentrations in residual fluids. Our fractionation models suggest that the Tl isotopic compositions of the primary magmas were among the isotopically lightest (less than or equal to ε205Tl = −4.6) material derived from the mantle for which Tl isotopic constraints exist. If so, the ultimate source of Tl in Khanneshin lavas—and perhaps carbonatites elsewhere—may be recycled ocean crust.

Description: This project was supported by funding from Woods Hole Oceanographic Institution Independent Research & Development funds and the National Science Foundation (Award #1911699).

Description: 2021-07-27

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Lattaud, J., Broder, L., Haghipour, N., Rickli, J., Giosan, L., & Eglinton, T., I. Influence of hydraulic connectivity on carbon burial efficiency in Mackenzie Delta lake sediments. Journal of Geophysical Research: Biogeosciences, 126(3), (2021): e2020JG006054, https://doi.org/10.1029/2020JG006054.

Description: The Arctic is undergoing accelerated changes in response to ongoing modifications to the climate system, and there is a need for local to regional scale records of past climate variability in order to put these changes into context. The Mackenzie Delta region in northern Canada is populated by numerous small shallow lakes. They are classified as no-, low-, and high-closure (NC, LC, and HC, respectively) lakes, reflecting varying degrees of connection to the river main stem, and have different sedimentation characteristics. This study examines sedimentological (mineral surface area, grain size), carbon isotopic (bulk and molecular-level) and inorganic isotopic (neodymium) characteristics of sediment cores from three lakes representing each class. We find that HC lake sediments exhibit strikingly different properties from the other lake sediments. Specifically, they are characterized by higher organic carbon loadings per unit mineral surface area and with relatively minor influence from allochthonous, petrogenic (rock-derived) organic carbon. In contrast, LC and NC lakes have the potential to record basin-scale climatic changes at a high resolution by virtue of enhanced detrital sedimentation. Overall the delta lakes have the capacity to bury about 2 MtC year−1, with little changes in the last 200 years. However, in the (near) future, an increased number of high closure lakes might change the carbon burial efficiency of the Mackenzie Delta as they seem to retain less carbon than NC and LC lakes.

Description: J. Lattaud was funded by a Rubicon grant (019.183EN.002) from NWO, Netherlands Organization for scientific research.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Jackson, R. L., Gabric, A. J., Matrai, P. A., Woodhouse, M. T., Cropp, R., Jones, G. B., Deschaseaux, E. S. M., Omori, Y., McParland, E. L., Swan, H. B., & Tanimoto, H. Parameterizing the impact of seawater temperature and irradiance on dimethylsulfide (DMS) in the Great Barrier Reef and the contribution of coral reefs to the global sulfur cycle. Journal of Geophysical Research:Oceans, 126(3), (2021): e2020JC016783, https://doi.org/10.1029/2020JC016783.

Description: Biogenic emissions of dimethylsulfide (DMS) are an important source of sulfur to the atmosphere, with implications for aerosol formation and cloud albedo over the ocean. Natural aerosol sources constitute the largest uncertainty in estimates of aerosol radiative forcing and climate and thus, an improved understanding of DMS sources is needed. Coral reefs are strong point sources of DMS; however, this coral source of biogenic sulfur is not explicitly included in climatologies or in model simulations. Consequently, the role of coral reefs in local and regional climate remains uncertain. We aim to improve the representation of tropical coral reefs in DMS databases by calculating a climatology of seawater DMS concentration (DMSw) and sea-air flux in the Great Barrier Reef (GBR), Australia. DMSw is calculated from remotely sensed observations of sea surface temperature and photosynthetically active radiation using a multiple linear regression model derived from field observations of DMSw in the GBR. We estimate that coral reefs and lagoon waters in the GBR (∼347,000 km2) release 0.03–0.05 Tg yr−1 of DMS (0.02 Tg yr−1 of sulfur). Based on this estimate, global tropical coral reefs (∼600,000 km2) could emit 0.08 Tg yr−1 of DMS (0.04 Tg yr−1 of sulfur), with the potential to influence the local radiative balance.

Description: Australian Research Council. Grant Number: DP150101649 National Science Foundation (NSF). Grant Number: 1543450 Ministry of Education, Culture, Sports, Science and Technology Grants-in-Aid for Scientific Research. Grant Number: 23310016,16H02967,24241010,15H01732 Ministry of Education, Culture, Sports, Science and Technology Grant-in-Aid for Young Scientists. Grant Number: 17K12812

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(3), (2021): e2020JB021136, https://doi.org/10.1029/2020JB021136.

Description: The Geomagnetic Polarity Time Scale (GPTS) provides a basis for the geological timescale, quantifies geomagnetic field behavior, and gives a time framework for geologic studies. We build a revised Middle to Late Jurassic GPTS by using a new multiscale magnetic profile, combining sea surface, midwater, and autonomous underwater vehicle near-bottom magnetic anomaly data from the Hawaiian lineation set in the Pacific Jurassic Quiet Zone (JQZ). We correlate the new profile with a previously published contemporaneous magnetic sequence from the Japanese lineation set. We then establish geomagnetic polarity block models as a basis for our interpretation of the origin and nature of JQZ magnetic anomalies and a GPTS. A significant level of coherency between short-wavelength anomalies for both the Japanese and Hawaiian lineation magnetic anomaly sequences suggests the existence of a regionally coherent field during this period of rapid geomagnetic reversals. Our study implies the rapid onset of the Mesozoic Dipole Low from M42 through M39 and then a subsequent gradual recovery in field strength into the Cenozoic. The new GPTS, together with the Japanese sequence, extends the magnetic reversal history from M29 back in time to M44. We identify a zone of varying, difficult-to-correlate anomalies termed the Hawaiian Disturbed Zone, which is similar to the zone of low amplitude, difficult-to-correlate anomalies in the Japanese sequence termed the Low Amplitude Zone (LAZ). We suggest that the LAZ, bounded by M39–M41 isochrons, may in fact represent the core of what is more commonly known as the JQZ crust.

Description: This study is funded by National Science Foundation grants OCE-1029965 (Tominaga, Tivey, and Lizarralde) and OCE-1233000 (Tominaga and Tivey) and OCE-1029573 (Sager).

Description: 2021-07-21

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Catunda, M. C. A., Bahr, A., Kaboth-Bahr, S., Zhang, X., Foukal, N. P., & Friedrich, O. Subsurface heat channel drove sea surface warming in the high-latitude North Atlantic during the Mid-Pleistocene Transition. Geophysical Research Letters, 48(11), (2021): e2020GL091899, https://doi.org/10.1029/2020GL091899.

Description: The Mid-Pleistocene Transition (MPT, 1,200–600 ka) marks the rapid expansion of Northern Hemisphere (NH) continental ice sheets and stronger precession pacing of glacial/interglacial cyclicity. Here, we investigate the relationship between thermocline depth in the central North Atlantic, subsurface northward heat transport and the initiation of the 100-kyr cyclicity during the MPT. To reconstruct deep-thermocline temperatures, we generated a Mg/Ca-based temperature record of deep-dwelling (∼800 m) planktonic foraminifera from mid-latitude North Atlantic at Site U1313. This record shows phases of pronounced heat accumulation at subsurface levels during the mid-MPT glacial driven by increased outflow of the Mediterranean Sea. Concurrent warming of the subtropical thermocline and subpolar surface waters indicates enhanced (subsurface) inter-gyre transport of warm water to the subpolar North Atlantic, which provided moisture for ice-sheet growth. Precession-modulated variability in the northward transport of subtropical waters imprinted this orbital cyclicity into NH ice-sheets after Marine Isotope Stage 24.

Description: Catunda and A. Bahr were funded by DFG project BA 3809/8, O.F. by DFG project FR 2544/11. S. Kaboth-Bahr acknowledges an Open-Topic Post-Doc Grant from the University of Potsdam. X.Z. was funded via the Lanzhou University (project 225000–830006) and National Science Foundation of China (Grant 42075047). N.F. was funded by the NSF Grant 1756361. Open access funding enabled and organized by Projekt DEAL.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(5), (2021): e2020GL091461, https://doi.org/10.1029/2020GL091461.

Description: We investigate how the near-surface chlorophyll (CHL)-a evolves in Gulf Stream (GS) warm-core rings (WCRs) and cold-core rings (CCRs) using multi-platform satellite observations. Averaged CHL anomaly (CHLA) within the rings exhibits both positive and negative linear trends during the evolution of the WCRs while negative trends dominate in CCRs. This difference is associated with a variety of physical processes occurring during the evolution process. Meanwhile, eddy-centric analysis reveals four spatial patterns of CHLA long-term trends, some of which highlights the importance of rings in shaping surface CHL. Short-term fluctuations of CHLA in WCRs and CCRs are closely correlated with mixed layer depth and sea surface temperature anomaly and highlight the complex interplay between multiple mechanisms. In addition, we find higher concentration CHL in some WCRs than that in CCRs during the same season, providing an alternative view of the characteristics of the surface ecosystem in Gulf Stream rings.

Description: This work was supported by the National Science Foundation Ocean Science Division under grant OCE-1558960. JN was supported by the Fundamental Research Funds for the Central Universities (Hohai University) under grant B200203005 and the China Scholarship Council.

Description: 2021-08-16

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(4), (2021): e2020JC016789, https://doi.org/10.1029/2020JC016789.

Description: Argo profiling floats and L-band passive microwave remote sensing have significantly improved the global sampling of sea surface salinity (SSS) in the past 15 years, allowing the study of the range of SSS seasonal variability using concurrent satellite and in situ platforms. Here, harmonic analysis was applied to four 0.25° satellite products and two 1° in situ products between 2016 and 2018 to determine seasonal harmonic patterns. The 0.25° World Ocean Atlas (WOA) version 2018 was referenced to help assess the harmonic patterns from a long-term perspective based on the 3-year period. The results show that annual harmonic is the most characteristic signal of the seasonal cycle, and semiannual harmonic is important in regions influenced by monsoon and major rivers. The percentage of the observed variance that can be explained by harmonic modes varies with products, with values ranging between 50% and 72% for annual harmonic and between 15% and 19% for semiannual harmonic. The large spread in the explained variance by the annual harmonic reflects the large disparity in nonseasonal variance (or noise) in the different products. Satellite products are capable of capturing sharp SSS features on meso- and frontal scales and the patterns agree well with the WOA 2018. These products are, however, subject to the impacts of radiometric noises and are algorithm dependent. The coarser-resolution in situ products may underrepresent the full range of high-frequency small scale SSS variability when data record is short, which may have enlarged the explained SSS variance by the annual harmonic.

Description: L. Yu was funded by NASA Ocean Salinity Science Team (OSST) activities through Grant 80NSSC18K1335. FMB was funded by the NASA OSST through Grant 80NSSC18K1322. E. P. Dinnat was funded by NASA through Grant 80NSSC18K1443. This research is carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA.

Description: 2021-09-17

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Schultz, C., Doney, S. C., Hauck, J., Kavanaugh, M. T., & Schofield, O. Modeling phytoplankton blooms and inorganic carbon responses to sea-ice variability in the West Antarctic Peninsula. Journal of Geophysical Research: Biogeosciences, 126(4), (2021): e2020JG006227, https://doi.org/10.1029/2020JG006227.

Description: The ocean coastal-shelf-slope ecosystem west of the Antarctic Peninsula (WAP) is a biologically productive region that could potentially act as a large sink of atmospheric carbon dioxide. The duration of the sea-ice season in the WAP shows large interannual variability. However, quantifying the mechanisms by which sea ice impacts biological productivity and surface dissolved inorganic carbon (DIC) remains a challenge due to the lack of data early in the phytoplankton growth season. In this study, we implemented a circulation, sea-ice, and biogeochemistry model (MITgcm-REcoM2) to study the effect of sea ice on phytoplankton blooms and surface DIC. Results were compared with satellite sea-ice and ocean color, and research ship surveys from the Palmer Long-Term Ecological Research (LTER) program. The simulations suggest that the annual sea-ice cycle has an important role in the seasonal DIC drawdown. In years of early sea-ice retreat, there is a longer growth season leading to larger seasonally integrated net primary production (NPP). Part of the biological uptake of DIC by phytoplankton, however, is counteracted by increased oceanic uptake of atmospheric CO2. Despite lower seasonal NPP, years of late sea-ice retreat show larger DIC drawdown, attributed to lower air-sea CO2 fluxes and increased dilution by sea-ice melt. The role of dissolved iron and iron limitation on WAP phytoplankton also remains a challenge due to the lack of data. The model results suggest sediments and glacial meltwater are the main sources in the coastal and shelf regions, with sediments being more influential in the northern coast.

Description: C. Schultz, S. C. Doney, M. T. Kavanaugh, and O. Schofield acknowledge support by the US National Science Foundation (Grant no. PLR-1440435), and C. Schultz and S. C. Doney acknowledge support from the University of Virginia. This research has also received funding from the Helmholtz Young Investigator Group Marine Carbon and Ecosystem Feedbacks in the Earth System (MarESys), Grant number VH-NG-1301.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Grooms, I., Loose, N., Abernathey, R., Steinberg, J. M., Bachman, S. D., Marques, G., Guillaumin, A. P., & Yankovsky, E. Diffusion-Based smoothers for spatial filtering of gridded geophysical data. Journal of Advances in Modeling Earth Systems, 13(9), (2021): e2021MS002552, https://doi.org/10.1029/2021MS002552.

Description: We describe a new way to apply a spatial filter to gridded data from models or observations, focusing on low-pass filters. The new method is analogous to smoothing via diffusion, and its implementation requires only a discrete Laplacian operator appropriate to the data. The new method can approximate arbitrary filter shapes, including Gaussian filters, and can be extended to spatially varying and anisotropic filters. The new diffusion-based smoother's properties are illustrated with examples from ocean model data and ocean observational products. An open-source Python package implementing this algorithm, called gcm-filters, is currently under development.

Description: I.G. and N.L. are supported by NSF OCE 1912332. R.A. is supported by NSF OCE 1912325. J.S. is supported by NSF OCE 1912302. S.B. and G.M. are supported by NSF OCE 1912420. A.G. and E.Y. are supported by NSF GEO 1912357 and NOAA CVP NA19OAR4310364.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(10),(2021): e2021JB022228, https://doi.org/10.1029/2021JB022228.

Description: Seafloor massive sulfide deposits form in remote environments, and the assessment of deposit size and composition through drilling is technically challenging and expensive. To aid the evaluation of the resource potential of seafloor massive sulfide deposits, three-dimensional inverse modeling of geophysical potential field data (magnetic and gravity) collected near the seafloor can be carried out to further enhance geologic models interpolated from sparse drilling. Here, we present inverse modeling results of magnetic and gravity data collected from the active mound at the Trans-Atlantic Geotraverse hydrothermal vent field, located at 26°08′N on the Mid-Atlantic Ridge, using autonomous underwater vehicle and submersible surveying. Both minimum-structure and surface geometry inverse modeling methods were utilized. Through deposit-scale magnetic modeling, the outer extent of a chloritized alteration zone within the basalt host rock below the mound was resolved, providing an indication of the angle of the rising hydrothermal fluid and the depth and volume of seawater/hydrothermal mixing zone. The thickness of the massive sulfide mound was determined by modeling the gravity data, enabling the tonnage of the mound to be estimated at 2.17 ± 0.44 Mt through this geophysics-based, noninvasive approach.

Description: The authors would like to thank the captain, crew, and scientific team from the 2016 R/V Meteor M127 and 1994 R/V Yokosuka MODE'94 cruises for all their work collecting the data modeled in this study. C. Galley is funded through an NSERC Discovery Grant and Memorial University's School of Graduate Studies Grant.

Description: 2022-03-29

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Beasley, C., Kender, S., Giosan, L., Bolton, C. T., Anand, P., Leng, M. J., Nilsson-Kerr, K., Ullmann, C. V., Hesselbo, S. P., & Littler, K. Evidence of a South Asian proto-monsoon during the Oligocene-Miocene transition. Paleoceanography and Paleoclimatology, 36(9), (2021): e2021PA004278, https://doi.org/10.1029/2021PA004278.

Description: The geological history of the South Asian monsoon (SAM) before the Pleistocene is not well-constrained, primarily due to a lack of available continuous sediment archives. Previous studies have noted an intensification of SAM precipitation and atmospheric circulation during the middle Miocene (∼14 Ma), but no records are available to test how the monsoon changed prior to this. In order to improve our understanding of monsoonal evolution, geochemical and sedimentological data were generated for the Oligocene-early Miocene (30–20 Ma) from Indian National Gas Hydrate Expedition 01 Site NGHP-01-01A in the eastern Arabian Sea, at 2,674 m water depth. We find the initial glaciation phase (23.7–23.0 Ma) of the Oligocene-Miocene transition (OMT) to be associated with an increase in water column ventilation and water mass mixing, suggesting an increase in winter monsoon type atmospheric circulation, possibly driven by a relative southward shift of the intertropical convergence zone. During the latter part of the OMT, or “deglaciation” phase (23.0–22.7 Ma), a long-term decrease in Mn (suggestive of deoxygenation), increase in Ti/Ca and dissolution of the biogenic carbonate fraction suggest an intensification of a proto-summer SAM system, characterized by the formation of an oxygen minimum zone in the eastern Arabian Sea and a relative increase of terrigenous material delivered by runoff to the site. With no evidence at this site for an active SAM prior to the OMT we suggest that changes in orbital parameters, as well as possibly changing Tethyan/Himalayan tectonics, caused this step change in the proto-monsoon system at this intermediate-depth site.

Description: This research forms part of a PhD study funded by the Natural Environment Research Council (NERC) Centre for Doctoral Training in Oil & Gas (grant number NE/M00578X/1) awarded to C. Beasley, and was also supported by a NERC National Environmental Isotope Facility Steering Committee grant (IP-1865-1118) awarded to S. Kender. L. Giosan acknowledges funding from USSP and WHOI and thanks colleagues from the NGHP-01 expedition. C. Ullmann acknowledges funding via NERC grant NE/N018508/1.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 22(5), (2021): e2020GC009608, https://doi.org/10.1029/2020GC009608.

Description: Thallium (Tl) isotope ratios are an emerging tool that can be used to trace crustal recycling processes in arc lavas and ocean island basalts (OIBs). Thallium is a highly volatile metal that is enriched in volcanic fumaroles, but it is unknown whether degassing of Tl from subaerial lavas has a significant effect on their residual Tl isotope compositions. Here, we present Tl isotope and concentration data from degassing experiments that are best explained by Rayleigh kinetic isotope fractionation during Tl loss. Our data closely follow predicted isotope fractionation models in which TlCl is the primary degassed species and where Tl loss is controlled by diffusion and natural convection, consistent with the slow gas advection velocity utilized during our experiments. We calculate that degassing into air should be associated with a net Tl isotope fractionation factor of αnet = 0.99969 for diffusion and natural gas convection (low gas velocities) and αnet = 0.99955 for diffusion and forced gas convection (high gas velocities). We also show that lavas from three volcanoes in the Kamchatka arc exhibit Tl isotope and concentration patterns that plot in between the two different gas convection regimes, implying that degassing played an important role in controlling the observed Tl isotope compositions in these three volcanoes. Literature inspection of Tl isotope data for subaerial lavas reveals that the majority of these appear only minorly affected by degassing, although a few samples from both OIBs and arc volcanoes can be identified that likely experienced some Tl degassing.

Description: National Science Foundation (NSF). Grant Numbers: EAR 1829546

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(15), (2021): e2021GL092779, https://doi.org/10.1029/2021GL092779.

Description: Double diffusion refers to a variety of turbulent processes in which potential energy is released into kinetic energy, made possible in the ocean by the difference in molecular diffusivities between salinity and temperature. Here, we present a new method for estimating the kinetic energy dissipation rates forced by double-diffusive convection using temperature and salinity data alone. The method estimates the up-gradient diapycnal buoyancy flux associated with double diffusion, which is hypothesized to balance the dissipation rate. To calculate the temperature and salinity gradients on small scales we apply a canonical scaling for compensated thermohaline variance (or ‘spice’) on sub-measurement scales with a fixed buoyancy gradient. Our predicted dissipation rates compare favorably with microstructure measurements collected in the Chukchi Sea. Fine et al. (2018), https://doi.org/10.1175/jpo-d-18-0028.1, showed that dissipation rates provide good estimates for heat fluxes in this region. Finally, we show the method maintains predictive skill when applied to a sub-sampling of the Conductivity, Temperature, Depth (CTD) data.

Description: This work was supported by the Natural Environment Research Council (grant number NE/L002507/1).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(5), (2021): e2020JC017042, https://doi.org/10.1029/2020JC017042.

Description: In frontal zones, water masses that are tens of kilometers in extent with origins in the mixed layer can be identified in the pycnocline for days to months. Here, we explore the pathways and mechanisms of subduction, the process by which water from the surface mixed layer makes its way into the pycnocline, using a submesoscale-resolving numerical model of a mesoscale front. By identifying Lagrangian trajectories of water parcels that exit the mixed layer, we study the evolution of dynamical properties from a statistical standpoint. Velocity- and buoyancy-gradients increase as water parcels experience both mesoscale (geostrophic) and submesoscale (ageostrophic) frontogenesis and subduct beneath the mixed layer into the stratified pycnocline along isopycnals that outcrop in the mixed layer. Subduction is transient and occurs in coherent regions along the front, the spatial and temporal scales of which influence the scales of the subducted water masses in the pycnocline. An examination of specific subduction events reveals a range of submesoscale features that support subduction. Contrary to the forced submesoscale processes that sequester low potential vorticity (PV) anomalies in the interior, we find that PV can be elevated in subducting water masses. The rate of subduction is of similar magnitude to previous studies (∼100 m/year), but the Lagrangian evolution of properties on water parcels and pathways that are unraveled in this study emphasize the role of submesoscale dynamics coupled with mesoscale frontogenesis.

Description: This research was funded by the ONR CALYPSO DRI grant N00014-16-1-3130. MAF was partially funded by a Martin Fellowship from MIT.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Long, M. H. Aquatic biogeochemical eddy covariance fluxes in the presence of waves. Journal of Geophysical Research: Oceans, 126(2), (2021): e2020JC016637, https://doi.org/10.1029/2020JC016637.

Description: The eddy covariance (EC) technique is a powerful tool for measuring atmospheric exchange rates that was recently adapted by biogeochemists to measure aquatic oxygen fluxes. A review of aquatic biogeochemical EC literature revealed that the majority of studies were conducted in shallow waters where waves were likely present, and that waves biased sensor and turbulence measurements. This review identified that larger measurement heights shifted turbulence to lower frequencies, producing a spectral gap between turbulence and wave frequencies. However, some studies sampled too close to the boundary to allow for a spectral turbulence‐wave gap, and a change in how EC measurements are conducted and analyzed is needed to remove wave‐bias. EC fluxes have only been derived from the time‐averaged product of vertical velocity and oxygen, often resulting in wave‐bias. Presented is a new analysis framework for removing wave‐bias by accumulation of cross‐power spectral densities below wave frequencies. This analysis framework also includes new measurement guidelines based on wave period, currents, and measurement heights. This framework is applied to sand, seagrass, and reef environments where traditional EC analysis resulted in wave‐bias of 7.0% ± 9.2% error in biogeochemical (oxygen and H+) fluxes, while more variable and higher error was evident in momentum fluxes (10.5% ± 21.0% error). It is anticipated that this framework will lead to significant changes in how EC measurements are conducted and evaluated, and help overcome the major limitations caused by wave‐sensitive and slow‐response sensors, potentially expanding new chemical tracer applications and more widespread use of the EC technique.

Description: This work was supported by the Independent Research & Development Program at WHOI grant 25307and NSF OCE grants 1657727 and 1633951.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(2), (2021): e2020JB019962, https://doi.org/10.1029/2020JB019962.

Description: The largest slip in great megathrust earthquakes often occurs in the 10–30 km depth range, yet seismic imaging of the material properties in this region has proven difficult. We utilize a dense onshore‐offshore passive seismic dataset from the southernmost Cascadia subduction zone where seismicity in the mantle of the subducted Gorda Plate produces S‐to‐P and P‐to‐S conversions generated within a few km of the plate interface. These conversions typically occur in the 10–20 km depth range at either the top or bottom of a ∼5 km thick layer with a high Vp/Vs that we infer to be primarily the subducted crust. We use their arrival times and amplitudes to infer the location of the top and bottom of the subducted crust as well as the velocity contrasts across these discontinuities. Comparing with both the Slab1.0 and the updated Slab2 interface models, the Slab2 model is generally consistent with the converted phases, while the Slab1.0 model is 1–2 km deeper in the 2–20 km depth range and ∼6–8 km too deep in the 10–20 km depth range between 40.25°N and 40.4°N. Comparing the amplitudes of the converted phases to synthetics for simplified velocity structures, the amplitude of the converted phases requires models containing a ∼5 km thick zone with at least a ∼10%–20% reduction in S wave velocity. Thus, the plate boundary is likely contained within or at the top of this low velocity zone, which potentially indicates a significant porosity and fluid content within the seismogenic zone.

Description: This work is funded by National Science Foundation Award Numbers EAR‐1520690.

Description: 2021-07-25

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(11), (2020): e2019JC015851, doi:10.1029/2019JC015851.

Description: Influences of the ocean mixed layer (OML) dynamics on intensity, pathway, and landfall of October 2012 Hurricane Sandy were examined through an experiment using the Weather Research and Forecasting (WRF) model. The WRF model was run for two cases with or without coupling to the OML. The OML in the WRF was calculated by an oceanic mixed layer submodel. The initial conditions of the depth and mean water temperature of the OML were specified using Global‐FVCOM and Global‐HYCOM fields. The comparison results between these two cases clearly show that including the OML dynamics enhanced the contribution of vertical mixing to the air‐sea heat flux. When the hurricane moved toward the coast, the local OML rapidly deepened with an increase of storm wind. Intense vertical mixing brought cold water in the deep ocean toward the surface to produce a cold wake underneath the storm, with the lowest sea temperature at the maximum wind zone. This process led to a significant latent heat loss from the ocean within the storm and hence rapid drops of the air temperature and vapor mixing ratio above the sea surface. As a result, the storm was intensified as the central sea level pressure dropped. Improving air pressure simulation with OML tended to reduce the storm size and strengthened the storm intensity and hence provided a better simulation of hurricane pathway and landfall.

Description: This work was supported by the MIT Sea Grant College Program through grant 2017‐R/RCM‐49C and 2012‐R/RC‐127, the NSF grants OCE1459096, OCE1332207, and OCE1332666, the NOAA‐funded IOOS NERACOOS program for NECOFS with subcontract numbers NA16NOS0120023 and NERACOOS A002 and A007, and the NOAA‐CINAR Hurricane Sandy fund. The development of the Global‐FVCOM system has been supported by NSF grants OCE1603000. S. Li was supported partially by the oversea Ph.D. fellowship from the China Scholarship Council (No. 1409010025).

Description: 2021-04-07

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Evans, D. G., Frajka-Williams, E., Garabato, A. C. N., Polzin, K. L., & Forryan, A. Mesoscale eddy dissipation by a "zoo" of submesoscale processes at a western boundary. Journal of Geophysical Research: Oceans, 125(11), (2020): e2020JC016246, doi:10.1029/2020JC016246.

Description: Mesoscale eddies are ubiquitous dynamical features that tend to propagate westward and disappear along ocean western boundaries. Using a multiscale observational study, we assess the extent to which eddies dissipate via a direct cascade of energy at a western boundary. We analyze data from a ship‐based microstructure and velocity survey, and an 18‐month mooring deployment, to document the dissipation of energy in anticyclonic and cyclonic eddies impinging on the topographic slope east of the Bahamas, in the North Atlantic Ocean. These observations reveal high levels of turbulence where the steep and rough topographic slope modified the intensified northward flow associated with, in particular, anticyclonic eddies. Elevated dissipation was observed both near‐bottom and at mid depths (200–800 m). Near‐bottom turbulence occurred in the lee of a protruding escarpment, where elevated Froude numbers suggest hydraulic control. Energy was also radiated in the form of upward‐propagating internal waves. Elevated dissipation at mid depths occurred in regions of strong vertical shear, where the topographic slope modified the vertical structure of the northward eddy flow. Here, low Richardson numbers and a local change in the isopycnal gradient of potential vorticity (PV) suggest that the elevated dissipation was associated with horizontal shear instability. Elevated mid‐depth dissipation was also induced by topographic steering of the flow. This led to large anticyclonic vorticity and negative PV adjacent to the topographic slope, suggesting that centrifugal instability underpinned the local enhancement in dissipation. Our results provide a mechanistic benchmark for the realistic representation of eddy dissipation in ocean models.

Description: The MeRMEED project, DGE, EFW, ACNG and AF were funded under Natural Environment Research Council standard grant NE/N001745/2. ACNG further acknowledges the support of the Royal Society and the Wolfson Foundation.

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 21(11), (2020): e2020GC009074, doi:10.1029/2020GC009074.

Description: Marine ferromanganese deposits, often called the scavengers of the sea, adsorb and coprecipitate with a wide range of metals of great interest for paleo‐environmental reconstructions and economic geology. The long (up to ∼75 Ma), near‐continuous record of seawater chemistry afforded by ferromanganese deposits offers much historical information about the global ocean and surface earth including crustal processes, mantle processes, ocean circulation, and biogeochemical cycles. The extent to which the ferromanganese deposits hosting these geochemical proxies undergo diagenesis on the seafloor, however, remains an important and challenging factor in assessing the fidelity of such records. In this study, we employ multiple X‐ray techniques including micro–X‐ray fluorescence, bulk and micro–X‐ray absorption spectroscopy, and X‐ray powder diffraction to probe the structural, compositional, redox, and mineral changes within a single ferromanganese crust. These techniques illuminate a complex two‐dimensional structure characterized by crust growth controlled by the availability of manganese (Mn), a dynamic range in Mn oxidation state from +3.4 to +4.0, changes in Mn mineralogy over time, and recrystallization in the lower phosphatized portions of the crust. Iron (Fe) similarly demonstrates spatial complexity with respect to concentration and mineralogy, but lacks the dynamic range of oxidation state seen for Mn. Micrometer‐scale measurements of metal abundances reveal complex element associations between trace elements and the two major oxide phases, which are not typically resolvable via bulk analytical methods. These findings provide evidence of post‐depositional processes altering chemistry and mineralogy, and provide important geochemical context for the interpretation of element and isotopic records in ferromanganese crusts.

Description: This research is supported by NASA Exobiology NNX15AM046 to Scott D. Wankel and Colleen M. Hansel, NASA NESSF NNX15AR62H to Kevin M. Sutherland, and WHOI Ocean Exploration Institute to Colleen M. Hansel. The Stanford Synchrotron Radiation Lightsource was utilized in this study. Use of the Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE‐AC02‐76SF00515.

Description: 2021-04-26

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Wallace, E. J., Donnelly, J. P., van Hengstum, P. J., Winkler, T. S., McKeon, K., MacDonald, D., d'Entremont, N. E., Sullivan, R. M., Woodruff, J. D., Hawkes, A. D., & Maio, C. 1,050 years of hurricane strikes on long island in the Bahamas. Paleoceanography and Paleoclimatology, 36(3), (2021): e2020PA004156, https://doi.org/10.1029/2020PA004156.

Description: Sedimentary records of past hurricane activity indicate centennial-scale periods over the past millennium with elevated hurricane activity. The search for the underlying mechanism behind these active hurricane periods is confounded by regional variations in their timing. Here, we present a new high resolution paleohurricane record from The Bahamas with a synthesis of published North Atlantic records over the past millennium. We reconstruct hurricane strikes over the past 1,050 years in sediment cores from a blue hole on Long Island in The Bahamas. Coarse-grained deposits in these cores date to the close passage of seven hurricanes over the historical interval. We find that the intensity and angle of approach of these historical storms plays an important role in inducing storm surge near the site. Our new record indicates four active hurricane periods on Long Island that conflict with published records on neighboring islands (Andros and Abaco Island). We demonstrate these three islands do not sample the same storms despite their proximity, and we compile these reconstructions together to create the first regional compilation of annually resolved paleohurricane records in The Bahamas. Integrating our Bahamian compilation with compiled records from the U.S. coastline indicates basin-wide increased storminess during the Medieval Warm Period. Afterward, the hurricane patterns in our Bahamian compilation match those reconstructed along the U.S. East Coast but not in the northeastern Gulf of Mexico. This disconnect may result from shifts in local environmental conditions in the North Atlantic or shifts in hurricane populations from straight-moving to recurving storms over the past millennium.

Description: This work was funded by the National Science Foundation Graduate Research Fellowship (to E. J. W.), the Dalio Explore Foundation, and National Science Foundation grant OCE-1356708 (to J. P. D. and P. J. vH.).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(3), (2021): e2021JB021709, https://doi.org/10.1029/2021JB021709.

Description: Serpentinites are increasingly recognized as playing an important role in the global geochemical cycle. However, discriminating the contributions of serpentinites to arc magmas from those of other subduction components is challenging. The Okinawa Trough is a back-arc basin developed behind the Ryukyu subduction zone, where magmas are extensively affected by sediment subduction. In this study, we reported the F-Cl concentrations and Sr-Nd-Pb-B isotopes of basaltic andesites from the Yaeyama Graben, Yonaguni Graben, and Irabu Knoll in the southern Okinawa Trough. The Irabu Knoll lavas show the most enrichment of fluid-mobile elements and F ± Cl, and have the heaviest B isotopes (δ11B: +6.6 ± 1.5‰). They also have decoupled Sr-Nd isotopes: higher 87Sr/86Sr (∼0.7049) but have no obvious decrease of 143Nd/144Nd (∼0.5128). Results from slab dehydration modeling and mixing calculations suggest that the heavy δ11B in the Irabu Knoll lavas is not consistent with fluids derived from altered oceanic crust (AOC), sediments, or wedge serpentinites (formed in the mantle wedge), but rather from slab serpentinites (formed within the subducting plate); sediments control the subduction input of Nd, whereas the decoupled Sr-Nd isotopes are most likely due to the excess radiogenic Sr carried by AOC fluids. Our results imply that recycling of serpentinite fluids and AOC fluids are usually coupled in subduction zones, as the arc lavas influenced by subducted serpentinite generally show Sr-Nd isotopes decoupling. The large variation of Sr-Nd-B isotopes observed in a relatively localized area is consistent with a focused migration through the mantle wedge of components from multiple sources.

Description: This study was funded by the National Natural Science Foundation of China (91958213), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB42020402), the China Postdoctoral Science Foundation (2019M662454), the Shandong Provincial Natural Science Foundation, China (ZR2020QD068 and ZR2020MD068), the International Partnership Program of the Chinese Academy of Sciences (133137KYSB20170003), the Special Fund for the Taishan Scholar Program of Shandong Province (ts201511061), and the China Scholarship Council (201709410550).

Description: 2021-09-12

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(15), (2021): e2021GL094469, https://doi.org/10.1029/2021GL094469.

Description: The magnitude of natural oceanic dissolved oxygen (DO) variability remains poorly understood due to the short duration of the observational record. Here we present a high-resolution (4–9 years) reconstruction of the Southern California oxygen minimum zone (OMZ) through the Common Era using redox-sensitive metals. Rapid OMZ intensification on multidecadal timescales reveals greater DO variability than observed in instrumental records. An anomalous interval of intensified OMZ between 1600–1750 CE contradicts the expectation of better-ventilated mid-depth North Pacific during cool climates. Although the influence of low-DO Equatorial Pacific Intermediate Water on the Southern California Margin was likely weaker during this interval, we attribute the observed rapid deoxygenation to reduced North Pacific Intermediate Water (NPIW) ventilation. NPIW ventilation thus appears very sensitive to atmospheric circulation reorganization (e.g., a weakened Siberian High and Aleutian Low). In addition to temperature-induced gas solubility, atmospheric forcing under future anthropogenic influences could amplify OMZ variability.

Description: The authors are grateful for financial supports from NSF (OCE-1851242), SMAST, and UMass Dartmouth. GG was supported by NSF under grants OCE-1657853 and OCE-1558521.

Description: 2022-01-16

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 47(22), (2020): e2020GL088692, doi:10.1029/2020GL088692.

Description: Rapid increases in upper 700‐m Indian Ocean heat content (IOHC) since the 2000s have focused attention on its role during the recent global surface warming hiatus. Here, we use ocean model simulations to assess distinct multidecadal IOHC variations since the 1960s and explore the relative contributions from wind stress and buoyancy forcing regionally and with depth. Multidecadal wind forcing counteracted IOHC increases due to buoyancy forcing from the 1960s to the 1990s. Wind and buoyancy forcing contribute positively since the mid‐2000s, accounting for the drastic IOHC change. Distinct timing and structure of upper ocean temperature changes in the eastern and western Indian Ocean are linked to the pathway how multidecadal wind forcing associated with the Interdecadal Pacific Oscillation is transmitted and affects IOHC through local and remote winds. Progressive shoaling of the equatorial thermocline—of importance for low‐frequency variations in Indian Ocean Dipole occurrence—appears to be dominated by multidecadal variations in wind forcing.

Description: This work was supported by the Alexander von Humboldt Foundation (CCU and SR), The Investment in Science Fund given primarily by WHOI Trustee and Corporation Members (CCU), James E. and Barbara V. Moltz Fellowship for climate‐related research (CCU), the ARC Centre of Excellence for Climate Extremes (CE170100023; CCU and MHE), ARC DP150101331 (CCU and MHE), and PW was supported through grant IndoArchipel from the Deutsche Forschungsgemeinschaft (DFG) as part of the Special Priority Program (SPP)‐1889”Regional Sea Level Change and Society” (SeaLevel).

Description: 2021-04-26

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(6), (2021): e2020JB021341, https://doi.org/10.1029/2020JB021341.

Description: Constraining how the physical properties and seismic responses of recently deposited sands change with time is important for understanding earthquake site response, subsurface fluid flow, and early stages of lithification. Currently, however, there is no detailed (cm-scale) assessment of how sand's physical properties and associated seismic velocities evolve over the first two centuries after deposition. Here, we integrate sedimentation rates with seismic velocity and sediment physical properties data to assess how the vadose zone sands at Port Royal Beach, Jamaica, change within 180 years after deposition. We show that compressional and shear wave velocities increase with sediment age, whereas porosity, grain size, sorting, mineralogy, and cementation fraction remain relatively unchanged during the same period. Rock physics models (constrained by the measured physical properties) predict constant seismic velocities at all sites regardless of sediment age, though misfits between modeled and observed velocities increase with sediment age. We explain these misfits by proposing that shallow sands undergo microstructural grain reorganization that leads to a more uniform distribution of grain contact forces with time. Our results imply that beach sands undergo a previously undocumented lithification process that occurs before compaction.

Description: The Society of Exploration Geophysicists Geoscientists without Borders Grant and the Institute for Earth, Science, and Man at Southern Methodist University partially supported this work.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Johnson, J. E., Phillips, S. C., Clyde, W. C., Giosan, L., & Torres, M. E. Isolating detrital and diagenetic signals in magnetic susceptibility records from methane-bearing marine sediments. Geochemistry Geophysics Geosystems, 22(9), (2021): e2021GC009867, https://doi.org/10.1029/2021GC009867.

Description: Volume-dependent magnetic susceptibility (κ) is commonly used for paleoenvironmental reconstructions in both terrestrial and marine sedimentary environments where it reflects a mixed signal between primary deposition and secondary diagenesis. In the marine environment, κ is strongly influenced by the abundance of ferrimagnetic minerals regulated by sediment transport processes. Post-depositional alteration by H2S, however, can dissolve titanomagnetite, releasing reactive Fe that promotes pyritization and subsequently decreases κ. Here, we provide a new approach for isolating the detrital signal in κ and identifying intervals of diagenetic alteration of κ driven by organoclastic sulfate reduction (OSR) and the anaerobic oxidation of methane (AOM) in methane-bearing marine sediments offshore India. Using the correlation of a heavy mineral proxy from X-ray fluorescence data (Zr/Rb) and κ in unaltered sediments, we predict the primary detrital κ signal and identify intervals of decreased κ, which correspond to increased total sulfur content. Our approach is a rapid, high-resolution method that can identify overprinted κ resulting from pyritization of titanomagnetite due to H2S production in marine sediments. In addition, total organic carbon, total sulfur, and authigenic carbonate δ13C measurements indicate that both OSR and AOM can drive the observed κ loss, but AOM drives the greatest decreases in κ. Overall, our approach can enhance paleoenvironmental reconstructions and provide insight into paleo-positions of the sulfate-methane transition zone, past enhancements of OSR or paleo-methane seepage, and the role of detrital iron oxide minerals on the marine sediment sulfur sink, with consequences influencing the development of chemosynthetic biological communities at methane seeps.

Description: This research was supported by the American Chemical Society-Petroleum Research Fund Award #53006-ND8 and U.S. Department of Energy Grant #DE-FE0010120.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Goldstein, E. B., Buscombe, D., Lazarus, E. D., Mohanty, S. D., Rafique, S. N., Anarde, K. A., Ashton, A. D., Beuzen, T., Castagno, K. A., Cohn, N., Conlin, M. P., Ellenson, A., Gillen, M., Hovenga, P. A., Over, J.-S. R., Palermo, R., Ratliff, K. M., Reeves, I. R. B., Sanborn, L. H., Straub, J. A., Taylor, L. A., Wallace E. J., Warrick, J., Wernette, P., Williams, H. E. Labeling poststorm coastal imagery for machine learning: measurement of interrater agreement. Earth and Space Science, 8(9), (2021): e2021EA001896, https://doi.org/10.1029/2021EA001896.

Description: Classifying images using supervised machine learning (ML) relies on labeled training data—classes or text descriptions, for example, associated with each image. Data-driven models are only as good as the data used for training, and this points to the importance of high-quality labeled data for developing a ML model that has predictive skill. Labeling data is typically a time-consuming, manual process. Here, we investigate the process of labeling data, with a specific focus on coastal aerial imagery captured in the wake of hurricanes that affected the Atlantic and Gulf Coasts of the United States. The imagery data set is a rich observational record of storm impacts and coastal change, but the imagery requires labeling to render that information accessible. We created an online interface that served labelers a stream of images and a fixed set of questions. A total of 1,600 images were labeled by at least two or as many as seven coastal scientists. We used the resulting data set to investigate interrater agreement: the extent to which labelers labeled each image similarly. Interrater agreement scores, assessed with percent agreement and Krippendorff's alpha, are higher when the questions posed to labelers are relatively simple, when the labelers are provided with a user manual, and when images are smaller. Experiments in interrater agreement point toward the benefit of multiple labelers for understanding the uncertainty in labeling data for machine learning research.

Description: The authors gratefully acknowledge support from the U.S. Geological Survey (G20AC00403 to EBG and SDM), NSF (1953412 to EBG and SDM; 1939954 to EBG), Microsoft AI for Earth (to EBG and SDM), The Leverhulme Trust (RPG-2018-282 to EDL and EBG), and an Early Career Research Fellowship from the Gulf Research Program of the National Academies of Sciences, Engineering, and Medicine (to EBG). U.S. Geological Survey researchers (DB, J-SRO, JW, and PW) were supported by the U.S. Geological Survey Coastal and Marine Hazards and Resources Program as part of the response and recovery efforts under congressional appropriations through the Additional Supplemental Appropriations for Disaster Relief Act, 2019 (Public Law 116-20; 133 Stat. 871).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(8), (2021): e2021GL093127, https://doi.org/10.1029/2021GL093127.

Description: The seismometer deployed by the InSight lander measured the seismic velocity of the Martian crust. We use a rock physics model to interpret those velocities and constrain hydrogeological properties. The seismic velocity of the upper ∼10 km is too low to be ice-saturated. Hence there is no cryosphere that confines deeper aquifers and possibly no aquifers locally. An increase in seismic velocity at depths of ∼10 km could be explained by a few volume percent of mineral cement (1%–5%) in pore space and may document the past depth of aquifers.

Description: M. Manga was supported by NASA grant 80NSSC19K0545.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(18), (2021): e2021GL092621, https://doi.org/10.1029/2021GL092621.

Description: In the subsurface ocean, O2 depleted because of organic matter remineralization is generally estimated based on apparent oxygen utilization (AOU). However, AOU is an imperfect measure of oxygen utilization because of O2 air-sea disequilibrium at the site of deepwater formation. Recent methodological and instrumental advances have paved the way to further deconvolve the processes driving the O2 signature. Using numerical model simulations of the global ocean, we show that the measurements of the dissolved O2/Ar ratio, which so far have been confined to the ocean surface, can provide improved estimates of oxygen utilization, especially in regions where the disequilibrium at the site of deepwater formation is associated with physical processes. We discuss applications of this new approach and implications for the current tracers relying on O2 such as remineralization ratios, respiratory quotients, and preformed nutrients. Finally, we propose a new composite geochemical tracer, [O2]bio combining dissolved O2/Ar and phosphate concentration. Being insensitive to photosynthesis and respiration, the change in this new tracer reflects gas exchange at the air-sea interface at the sites of deepwater formation.

Description: Nicolas Cassar was supported by the “Laboratoire d'Excellence” LabexMER (ANR-10-LABX-19) and cofunded by a grant from the French government under the program “Investissements d'Avenir.” Samar Khatiwala was supported by UK NERC grant NE/T009357/1. Ellen Cliff acknowledges support from the Rhodes Trust.

Description: 2022-03-13

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bretschneider, L., Hathorne, E. C., Bolton, C. T., Gebregiorgis, D., Giosan, L., Gray, E., Huang, H., Holbourn, A., Kuhnt, W., & Frank, M. Enhanced late miocene chemical weathering and altered precipitation patterns in the watersheds of the Bay of Bengal recorded by detrital clay radiogenic isotopes. Paleoceanography and Paleoclimatology, 36(9), (2021): e2021PA004252, https://doi.org/10.1029/2021PA004252.

Description: The late Miocene was a period of declining CO2 levels and extensive environmental changes, which likely had a large impact on monsoon strength as well as on the weathering and erosion intensity in the South Asian Monsoon domain. To improve our understanding of these feedback systems, detrital clays from the southern Bay of Bengal (International Ocean Discovery Program Site U1443) were analyzed for the radiogenic isotope compositions of Sr, Nd, and Pb to reconstruct changes in sediment provenance and weathering regime related to South Asian Monsoon rainfall from 9 to 5 Ma. The 100 kyr resolution late Miocene to earliest Pliocene record suggests overall low variability in the provenance of clays deposited on the Ninetyeast Ridge. However, at 7.3 Ma, Nd and Pb isotope compositions indicate a switch to an increased relative contribution from the Irrawaddy River (by ∼10%). This shift occurred during the global benthic δ13C decline, and we suggest that global cooling and increasing aridity resulted in an eastward shift of precipitation patterns leading to a more focused erosion of the Indo-Burman Ranges. Sr isotope compositions were decoupled from Nd and Pb isotope signatures and became more radiogenic between 6 and 5 Ma. Grassland expansion generating thick, easily weatherable soils may have led to an environment supporting intense chemical weathering, which is likely responsible for the elevated detrital clay 87Sr/86Sr ratios during this time. This change in Sr isotope signatures may also have contributed to the late Miocene increase of the global seawater Sr isotope composition.

Description: This research used samples and data provided by the International Ocean Discovery Program and was funded by the German Research Foundation (DFG) (grants HA 5751/6-1 & -2). C. T. Bolton acknowledges funding from the French ANR project iMonsoon (ANR-16-CE01-0004-01) and IODP France. W. Kuhnt acknowledges funding from the DFG (grant Ku649/36-1).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Biogeosciences 126(7), (2021): e2020JG005977, https://doi.org/10.1029/2020JG005977.

Description: Increasing Arctic temperatures are thawing permafrost soils and liberating ancient organic matter, but the fate of this material remains unclear. Thawing of permafrost releases dissolved organic matter (DOM) into fluvial networks. Unfortunately, tracking this material in Arctic rivers such as the Kolyma River in Siberia has proven challenging due to its high biodegradability. Here, we evaluate late summer abruptly thawed yedoma permafrost dissolved organic carbon (DOC) inputs from Duvannyi Yar. We implemented ultrahigh-resolution mass spectrometry alongside ramped pyrolysis oxidation (RPO) and isotopic analyses. These approaches offer insight into DOM chemical composition and DOC radiocarbon values of thermochemical components for a permafrost thaw stream, the Kolyma River, and their biodegraded counterparts (n = 4). The highly aliphatic molecular formula found in undegraded permafrost DOM contrasted with the comparatively aliphatic-poor formula of Kolyma River DOM, represented by an 8.9% and 2.6% relative abundance, respectively, suggesting minimal inputs of undegraded permafrost DOM in the river. RPO radiocarbon fractions of Kolyma River DOC exhibited no “hidden” aged component indicative of permafrost influence. Thermostability analyses suggested that there was limited biodegraded permafrost DOC in the Kolyma River, in part determined by the formation of high-activation energy (thermally stable) biodegradation components in permafrost DOM that were lacking in the Kolyma River. A mixing model based on thermostability and radiocarbon allowed us to estimate a maximum input of between 0.8% and 7.7% of this Pleistocene-aged permafrost to the Kolyma River DOC. Ultimately, our findings highlight that export of modern terrestrial DOC currently overwhelms any permafrost DOC inputs in the Kolyma River.

Description: This work was funded by NSF grants ANT-1203885 and PLR-1500169 to R.G.M.S. The work was also supported by the National Science Foundation Division of Chemistry through DMR-1644779 and the State of Florida.

Description: 2022-01-09

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 22(4), (2021): e2020GC009481, https://doi.org/10.1029/2020GC009481.

Description: The impact of submarine hydrothermal systems on organic carbon in the ocean—one of the largest fixed carbon reservoirs on Earth—could be profound. Yet, different vent sites show diverse fluid chemical compositions and the subsequent biological responses. Observations from various vent sites are to evaluate hydrothermal systems' impact on the ocean carbon cycle. A response cruise in May 2009 to an on-going submarine eruption at West Mata Volcano, northeast Lau Basin, provided an opportunity to quantify the organic matter production in a back-arc spreading hydrothermal system. Hydrothermal vent fluids contained elevated dissolved organic carbon, particulate organic carbon (POC), and particulate nitrogen (PN) relative to background seawater. The δ13C-POC values for suspended particles in the diffuse vent fluids (−15.5‰ and −12.3‰) are distinct from those in background seawater (−23 ± 1‰), indicative of unique carbon synthesis pathways of the vent microbes from the seawater counterparts. The first dissolved organic nitrogen concentrations reported for diffuse vents were similar to or higher than those for background seawater. Enhanced nitrogen fixation and denitrification removed 37%–89% of the total dissolved nitrogen in the recharging background seawater in the hydrothermal vent flow paths. The hydrothermal plume samples were enriched in POC and PN, indicating enhanced biological production. The total “dark” organic carbon production within the plume matches the thermodynamic prediction based on available reducing chemical substances supplied to the plume. This research combines the measured organic carbon contents with thermodynamic modeled results and demonstrates the importance of hydrothermal activities on the water column carbon production in the deep ocean.

Description: This project was supported by N.S.F. (OCE0929881, J. P. Cowen and K. H. Rubin), the NOAA PMEL VENTS (now Earth-Ocean Interactions) Program and the Joint Institute for the Study of the Atmosphere and Ocean (JISAO) under NOAA Cooperative Agreement No. NA10OAR4320148, and the UH NASA Astrobiology Institute. The Ministry of Science and Technology of Taiwan award (MOST 107-2611-M-002-002, and MOST 108-2611-M-002-006 to H.-T. Lin). Ministry of Education (M.O.E.) Republic of China (Taiwan) 109L892601 to H.-T. Lin. SOEST contributions no. 11285, C-DEBI contribution no. 563. PMEL contribution no. 3996, JISAO contribution 2183.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Water Resources Research 57(4), (2021): e2020WR028430, https://doi.org/10.1029/2020WR028430.

Description: We use yearlong vertical temperature profile time-series (seven thermistors at evenly spaced depth intervals from 10 to 70 cm) from five sites in and around the Deep Hole thermal area, southeast of Stevenson Island, Yellowstone Lake, to investigate heat and mass fluxes across the lake floor. The records demonstrate that thermal gradients in surficial sediments are modulated by a rich spectrum of bottom water temperature variations generated by hydrodynamic processes, and that sites inside the thermal area also respond to hydrothermal variations. We develop and implement a new method for estimating the sediment effective thermal diffusivity and pore fluid vertical flow rate that exploits the full spectrum of observed temperature variations to generate the parameter estimates, uncertainties, and metrics to assess statistical significance. Sediments at sites outside thermal areas have gradients of ∼7.5°C/m, in situ thermal diffusivities of ∼1.6 × 10−7 m2/s consistent with highly porous (80–90%) siliceous sediments, and experience hypolentic flow in the upper ∼20 cm. Sites inside the Deep Hole thermal area exhibit considerable spatial and temporal variability, with gradients of 1–32°C/m, and higher thermal diffusivities of ∼2–12 × 10−7 m2/s, consistent with hydrothermal alteration of biogenic silica to clays, quartz, and pyrite. Upward pore fluid flow at these sites is observed across multiple depth intervals, with maximum values of ∼3 cm/day. The observed spatial and temporal variability within the thermal area is consistent with upward finger flow combined with short wavelength convection within the porous sediments above a steam reservoir.

Description: This research was supported by the National Science Foundation Grants EAR-1516361 to Robert A. Sohn and EAR-1515283 to Robert N. Harris, and by the Independent Research and Development Program at the Woods Hole Oceanographic Institution (Robert A. Sohn). All work in Yellowstone National Park was completed under an authorized Yellowstone research permit (YELL-2018-SCI-7018).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(5), (2021): e2020JC017091, https://doi.org/10.1029/2020JC017091.

Description: A region of exceptionally high macrofaunal benthic biomass exists in Barrow Canyon, implying a carbon export process that is locally concentrated. Here we offer an explanation for this benthic “hotspot” using shipboard data together with a set of dynamical equations. Repeat occupations of the Distributed Biological Observatory transect in Barrow Canyon reveal that when the northward flow is strong and the density front in the canyon is sharp, plumes of fluorescence and oxygen extend from the pycnocline to the seafloor in the vicinity of the hotspot. By solving the quasi-geostrophic omega equation with an analytical flow field fashioned after the observations, we diagnose the vertical velocity in the canyon. This reveals that, as the along stream flow converges into the canyon, it drives a secondary circulation cell with strong downwelling on the cyclonic side of the northward flow. The downwelling quickly advects material from the pycnocline to the seafloor in a vertical plume analogous to those seen in the observations. The plume occurs only when the phytoplankton reside in the pycnocline, since the near-surface vertical velocity is weak, also consistent with the observations. Using a wind-based proxy to represent the strength of the northward flow and hence the pumping, in conjunction with a satellite-derived phytoplankton source function, we construct a time series of carbon supply to the bottom of Barrow Canyon.

Description: This work was funded by National Science Foundation grants PLR-1504333 and OPP-1733564 (Robert S. Pickart, Frank Bahr), OPP-1822334 (Michael A. Spall), PLR-1304563 (Kevin R. Arrigo), OPP-1204082 and OPP-1702456 (Jacqueline M. Grebmeier); National Oceanic and Atmospheric Administration grants NA14OAR4320158 and NA19OAR4320074 (Robert S. Pickart, Peigen Lin, Leah T. McRaven), CINAR-22309.02 (Jacqueline M. Grebmeier).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(5), (2021): e2020JC016922, https://doi.org/10.1029/2020JC016922.

Description: Mesoscale eddies redistribute heat, salt, and nutrients in oceans. The South Atlantic Ocean (SA) is a basin that has active mesoscale eddies for which characteristics of the three-dimensional structure and its leading mechanism are complex but have yet been studied sufficiently. Here based on ocean reanalysis datasets we use a composite analysis approach to analyze the mixed layer anomalous heat budget and find distinct two types of spatial patterns: dipole and monopole – mainly present in the northern and southern regions of the SA, respectively. The dipole can be attributed to ocean horizontal advection, especially to the combined effect of eddy anomalous meridional current and meridional gradient of mean temperature. The monopole, on the other hand, is associated with complex contributions, for which zonal and meridional advections play opposite roles as cooling or heating around the eddies. At the eddy center, the vertical advection is non-negligible, especially the mean upwelling and vertical temperature gradient playing a vital role in the formation of a monopole. The analysis of eddy meridional heat transport shows that the stirring component is dominant, and poleward in most areas, especially at high latitudes. Such analysis on the leading mechanism of eddy-induced temperature anomaly could help improve our understanding on meso- and small-scale air-sea interactions and eddy-induced heat transport in the SA.

Description: This work is supported by the National Key R&D Program of China (2017YFC1404100 and 2017YFC1404104) and the National Natural Science Foundation of China (Grant No. 41775100, 41830964) as well as Shandong Province’s “Taishan” Scientist Program and Qingdao “Creative and Initiative” frontier Scientist Program. This research is also supported by the Center for High Performance Computing and System Simulation, Pilot National Laboratory for Marine Science and Technology (Qingdao).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(15), (2021): e2021GL093675, https://doi.org/10.1029/2021GL093675.

Description: Tide gauges provide a rich, long-term, record of the amplitude and spatiotemporal structure of interannual to multidecadal coastal sea-level variability, including that related to North American east coast sea level “hotspots.” Here, using wavelet analyses, we find evidence for multidecadal epochs of enhanced decadal (10–15 year period) sea-level variability at almost all long ( 70 years) east coast tide gauge records. Within this frequency band, large-scale spatial covariance is time-dependent; notably, coastal sectors north and south of Cape Hatteras exhibit multidecadal epochs of coherence ( 1960–1990) and incoherence ( 1990-present). Results suggest that previous interpretations of along coast covariance, and its underlying physical drivers, are clouded by time-dependence and frequency-dependence. Although further work is required to clarify the mechanisms driving sea-level variability in this frequency band, we highlight potential associations with the North Atlantic sea surface temperature tripole and Atlantic Multidecadal Variability.

Description: Christopher M. Little acknowledges funding support from NSF Grant OCE-1805029. CGP and RMP were funded through NASA Sea Level Change Team (CGP: Grant 80NSSC20K1241).

Description: 2022-01-15

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(4), (2021): e2020JB019395, https://doi.org/10.1029/2020JB019395.

Description: Improved understanding of the impact of crystal mush rheology on the response of magma chambers to magmatic events is critical for better understanding crustal igneous systems with abundant crystals. In this study, we extend an earlier model by Liao et al. (2018); https://doi.org/10.1029/2018jb015985 which considers the mechanical response of a magma chamber with poroelastic crystal mush, by including poroviscoelastic rheology of crystal mush. We find that the coexistence of the two mechanisms of poroelastic diffusion and viscoelastic relaxation causes the magma chamber to react to a magma injection event with more complex time-dependent behaviors. Specifically, we find that the system’s short-term evolution is dominated by the poroelastic diffusion process, while its long-term evolution is dominated by the viscoelastic relaxation process. We identify two post-injection timescales that represent these two stages and examine their relation to the material properties of the system. We find that better constraints on the poroelastic diffusion time are more important for the potential interpretation of surface deformation using the model.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Warner, J. C., Geyer, W. R., Ralston, D. K., & Kalra, T. Using tracer variance decay to quantify variability of salinity mixing in the Hudson River Estuary. Journal of Geophysical Research: Oceans, 125(12), (2020): e2020JC016096, https://doi.org/10.1029/2020JC016096.

Description: The salinity structure in an estuary is controlled by time‐dependent mixing processes. However, the locations and temporal variability of where significant mixing occurs is not well‐understood. Here we utilize a tracer variance approach to demonstrate the spatial and temporal structure of salinity mixing in the Hudson River Estuary. We run a 4‐month hydrodynamic simulation of the tides, currents, and salinity that captures the spring‐neap tidal variability as well as wind‐driven and freshwater flow events. On a spring‐neap time scale, salinity variance dissipation (mixing) occurs predominantly during the transition from neap to spring tides. On a tidal time scale, 60% of the salinity variance dissipation occurs during ebb tides and 40% during flood tides. Spatially, mixing during ebbs occurs primarily where lateral bottom salinity fronts intersect the bed at the transition from the main channel to adjacent shoals. During ebbs, these lateral fronts form seaward of constrictions located at multiple locations along the estuary. During floods, mixing is generated by a shear layer elevated in the water column at the top of the mixed bottom boundary layer, where variations in the along channel density gradients locally enhance the baroclinic pressure gradient leading to stronger vertical shear and more mixing. For both ebb and flood, the mixing occurs at the location of overlap of strong vertical stratification and eddy diffusivity, not at the maximum of either of those quantities. This understanding lends a new insight to the spatial and time dependence of the estuarine salinity structure.

Description: This study was funded through the Coastal Model Applications and Field Measurements Project and the Cross‐shore and Inlets Project, US Geological Survey Coastal Marine Hazards and Resources Program. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(4), (2021): e2020JC016757, https://doi.org/10.1029/2020JC016757.

Description: The along-shelf circulation in the Northwest Atlantic (NWA) Ocean is characterized by an equatorward flow from Greenland's south coast to Cape Hatters. The mean flow is considered to be primarily forced by freshwater discharges from rivers and glaciers while its variability is driven by both freshwater fluxes and wind stress. In this study, we hypothesize and test that the wind stress is important for the mean along-shelf flow. A two-layer model with realistic topography when forced by wind stress alone simulates a circulation system on the NWA shelves that is broadly consistent with that derived from observations, including an equatorward flow from Greenland coast to the Mid-Atlantic Bight (MAB). The along-shelf sea-level gradient is close to a previous estimate based on observations. The along-shelf flows exhibit strong seasonal variations with along-shelf transports being strong in fall/winter and weak in spring/summer, consistent with available observations. It is found that the NWA shelf circulation is affected by both wind-driven gyres through their western boundary currents and wind-stress forcing on the shelf especially along the coasts of Newfoundland and Labrador. The local wind stress forcing has more direct impacts on flows in shallower waters along the coast while the open-ocean gyres tend to affect the circulations along the outer shelf. Our conclusion is that wind stress is an important forcing of the main along-shelf flows in the NWA. One objective of this study is to motivate further examination of whether wind stress is as important as freshwater forcing for the mean flow.

Description: Both Yang and Chen are also supported by NOAA Climate Program Office's Climate Variability and Prediction Program under grant NA20OAR4310398. JY is supported by Woods Hole Oceanographic Institution (WHOI) W. V. A. Clark Chair for Excellence in Oceanography and NSF Ocean Science Division under grant OCE1634886. Chen is supported by WHOI Independent Research and Development award.

Description: 2021-09-30

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Donatelli, C., Kalra, T. S., Fagherazzi, S., Zhang, X., & Leonardi, N. Dynamics of marsh-derived sediments in lagoon-type estuaries. Journal of Geophysical Research: Earth Surface, 125(12), (2020): e2020JF005751, doi:10.1029/2020JF005751.

Description: Salt marshes are valuable ecosystems that must trap sediments and accrete in order to counteract the deleterious effect of sea level rise. Previous studies have shown that the capacity of marshes to build up vertically depends on both autogenous and exogenous processes including ecogeomorphic feedbacks and sediment supply from in‐land and coastal ocean. There have been numerous efforts to quantify the role played by the sediments coming from marsh edge erosion on the resistance of salt marshes to sea level rise. However, the majority of existing studies investigating the interplay between lateral and vertical dynamics use simplified modeling approaches, and they do not consider that marsh retreat can affect the regional‐scale hydrodynamics and sediment retention in back‐barrier basins. In this study, we evaluated the fate of the sediments originating from marsh lateral loss by using high‐resolution numerical model simulations of Jamaica Bay, a small lagoonal estuary located in New York City. Our findings show that up to 42% of the sediment released during marsh edge erosion deposits on the shallow areas of the basin and over the vegetated marsh platforms, contributing positively to the sediment budget of the remaining salt marshes. Furthermore, we demonstrate that with the present‐day sediment supply from the ocean, the system cannot keep pace with sea level rise even accounting for the sediment liberated in the bay through marsh degradation. Our study highlights the relevance of multiple sediment sources for the maintenance of the marsh complex.

Description: This study was supported by the Department of the Interior Hurricane Sandy Recovery program (ID G16AC00455, subaward to University of Liverpool). S. F. was partly supported by NSF awards 1637630 (PIE LTER) and 1832221 (VCR LTER). We thank Robert Chant from Rutgers University for sharing the hydrodynamic measurements in Jamaica Bay.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Biogeosciences 126(5), (2021): e2020JG006217, https://doi.org/10.1029/2020JG006217.

Description: It is assumed that to treat excess NO3− high soil organic matter content (%OM) is required to maintain high denitrification rates in natural or restored wetlands. However, this excess also represents a risk by increasing soil decomposition rates triggering peat collapse and wetland fragmentation. Here, we evaluated the role of %OM and temperature interactions controlling denitrification rates in eroding (Barataria Bay-BLC) and emerging (Wax Lake Delta-WLD) deltaic regions in coastal Louisiana using the isotope pairing (IPT) and N2:Ar techniques. We also assessed differences between total (direct denitrification + coupled nitrification-denitrification) and net (total denitrification minus nitrogen fixation) denitrification rates in benthic and wetland habitats with contrasting %OM and bulk density (BD). Sediment (benthic) and soil (wetland) cores were collected during summer, spring, and winter (2015–2016) and incubated at close to in-situ temperatures (30°C, 20°C, and 10°C, respectively). Denitrification rates were linearly correlated with temperature; maximum mean rates ranged from 40.1–124.1 μmol m−2 h−1 in the summer with lower rates (〈26.2 ± 5.3 μmol m−2 h−1) in the winter seasons. Direct denitrification was higher than coupled denitrification in all seasons. Denitrification rates were higher in WLD despite lower %OM, lower total N concentration, and higher BD in wetland soils. Therefore, in environments with low carbon availability, high denitrification rates can be sustained as long as NO3− concentrations are high (〉30 μM) and water temperature is 〉10°C. In coastal Louisiana, substrates under these regimes are represented by emergent supra-tidal flats or land created by sediment diversions under oligohaline conditions (〈1 ppt).

Description: This study was supported by the NOAA-Sea Grant Program-Louisiana (Grant 2013R/E-24) to Victor H. Rivera-Monroy and Kanchan Maiti. Victor H. Rivera-Monroy was also supported by the Department of the Interior South-Central Climate Adaptation Science Center (Cooperative Agreement #G12AC00002).

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(12),(2020): e2020JC016271, https://doi.org/10.1029/2020JC016271.

Description: Asian summer monsoon has a planetary‐scale, westward propagating “quasi‐biweekly” mode of variability with a 10–25 day period. Six years of moored observations at 18°N, 89.5°E in the north Bay of Bengal (BoB) reveal distinct quasi‐biweekly variability in sea surface salinity (SSS) during summer and autumn, with peak‐to‐peak amplitude of 3–8 psu. This large‐amplitude SSS variability is not due to variations of surface freshwater flux or river runoff. We show from the moored data, satellite SSS, and reanalyses that surface winds associated with the quasi‐biweekly monsoon mode and embedded weather‐scale systems, drive SSS and coastal sea level variability in 2015 summer monsoon. When winds are calm, geostrophic currents associated with mesoscale ocean eddies transport Ganga‐Brahmaputra‐Meghna river water southward to the mooring, salinity falls, and the ocean mixed layer shallows to 1–10 m. During active (cloudy, windy) spells of quasi‐biweekly monsoon mode, directly wind‐forced surface currents carry river water away to the east and north, leading to increased salinity at the moorings, and rise of sea level by 0.1–0.5 m along the eastern and northern boundary of the bay. During July–August 2015, a shallow pool of low‐salinity river water lies in the northeastern bay. The amplitude of a 20‐day oscillation of sea surface temperature (SST) is two times larger within the fresh pool than in the saltier ocean to the west, although surface heat flux is nearly identical in the two regions. This is direct evidence that spatial‐temporal variations of BoB salinity influences sub‐seasonal SST variations, and possibly SST‐mediated monsoon air‐sea interaction.

Description: The authors thank the Ministry of Earth Sciences (MoES) institutes NIOT and INCOIS, and the Upper Ocean Processes (UOP) group at WHOI for design, integration, and deployment of moorings in the BoB. The WHOI mooring was deployed from the ORV Sagar Nidhi and recovered from the ORV Sagar Kanya—we thank the officers, crew and science teams on the cruises for their support. Sengupta, Ravichandran and Sukhatme acknowledge MoES and the National Monsoon Mission, Indian Institute of Tropical Meteorology (IITM), Pune, for support; Lucas and Farrar acknowledge the US Office of Naval Research for support of ASIRI through grants N00014‐13‐1‐0489, N0001413‐100453, N0001417‐12880. We thank S. Shivaprasad, Dipanjan Chaudhuri and Jared Buckley for discussion on ocean currents and Ekman flow, and Fabien Durand for discussion on sea level. JSL would like to thank the Divecha Center for Climate Change, IISc., for support. DS acknowledges support from the Department of Science and Technology (DST), New Delhi, under the Indo‐Spanish Programme.

Description: 2021-05-16

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Biogeosciences 126(1), (2021): e2019JG005621, https://doi.org/10.1029/2019JG005621.

Description: Ongoing ocean warming can release methane (CH4) currently stored in ocean sediments as free gas and gas hydrates. Once dissolved in ocean waters, this CH4 can be oxidized to carbon dioxide (CO2). While it has been hypothesized that the CO2 produced from aerobic CH4 oxidation could enhance ocean acidification, a previous study conducted in Hudson Canyon shows that CH4 oxidation has a small short‐term influence on ocean pH and dissolved inorganic radiocarbon. Here we expand upon that investigation to assess the impact of widespread CH4 seepage on CO2 chemistry and possible accumulation of this carbon injection along 234 km of the U.S. Mid‐Atlantic Bight. Consistent with the estimates from Hudson Canyon, we demonstrate that a small fraction of ancient CH4‐derived carbon is being assimilated into the dissolved inorganic radiocarbon (mean fraction of 0.5 ± 0.4%). The areas with the highest fractions of ancient carbon coincide with elevated CH4 concentration and active gas seepage. This suggests that aerobic CH4 oxidation has a greater influence on the dissolved inorganic pool in areas where CH4 concentrations are locally elevated, instead of displaying a cumulative effect downcurrent from widespread groupings of CH4 seeps. A first‐order approximation of the input rate of ancient‐derived dissolved inorganic carbon (DIC) into the waters overlying the northern U.S. Mid‐Atlantic Bight further suggests that oxidation of ancient CH4‐derived carbon is not negligible on the global scale and could contribute to deepwater acidification over longer time scales.

Description: This study was sponsored by U.S. Department of Energy (DE‐FE0028980, awarded to J. D. K; DE‐FE0026195 interagency agreement with C. D. R.). We thank the crew of the R/V Hugh R. Sharp for their support, G. Hatcher, J. Borden, and M. Martini of the USGS for assistance with the LADCP, and Zach Bunnell, Lillian Henderson, and Allison Laubach for additional support at sea.

Description: 2021-06-23

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Fu, X., Waite, W. F., & Ruppel, C. D. Hydrate formation on marine seep bubbles and the implications for water column methane dissolution. Journal of Geophysical Research: Oceans, 126(9), (2021): e2021JC017363, https://doi.org/10.1029/2021JC017363.

Description: Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate-coated bubbles to deliver methane to shallower depths in the water column than hydrate-free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble-release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate-free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas-water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water-column observations that most bubbles become hydrate-coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

Description: X. Fu acknowledges support from the Miller Fellowship during her time at U.C. Berkeley. W. Waite and C. Ruppel are supported by the United States Geological Survey (USGS) Coastal/Marine Hazards and Resources Program and the Energy Resources Program, with research conducted under USGS-Department of Energy interagency agreements DE-FE0023495 and 89243320SFE000013.

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 47(22), (2020): e2020GL090431, doi:10.1029/2020GL090431.

Description: Vast quantities of solid CO2 reside in topographic basins of the south polar layered deposits (SPLD) on Mars and exhibit morphological features indicative of glacial flow. Previous experimental studies showed that CO2 ice is 1–2 orders of magnitude weaker than water ice under Martian polar conditions. Here we present data from deformation experiments on pure, fine‐grained CO2 ice, over a broader range of temperatures than previously explored (158–213 K). The experiments confirm previous observations of highly nonlinear power law creep at larger stresses, but also show a transition to a previously unseen linear‐viscous creep regime at lower stresses. We examine the viscosity of CO2 within the SPLD and predict that the CO2‐rich deposits are modestly stronger than previously thought. Nevertheless, CO2 ice flows much more readily than H2O ice, particularly on the steep flanks of SPLD topographic basins, allowing the CO2 to pond as observed.

Description: This work was funded by NASA grant NNH16ZDA001N‐SSW awarded to Smith and Goldsby. Additional salary support for Cross was provided by the WHOI Investment in Science Fund.

Description: 2021-04-29

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 35(4), (2021): e2020GB006887, https://doi.org/10.1029/2020GB006887.

Description: In this study we report full-depth water column profiles for nitrogen and oxygen isotopic composition (δ15N and δ18O) of nitrate (NO3−) during the GEOTRACES GA01 cruise (2014). This transect intersects the double gyre system of the subtropical and subpolar regions of the North Atlantic separated by a strong transition zone, the North Atlantic Current. The distribution of NO3− δ15N and δ18O shows that assimilation by phytoplankton is the main process controlling the NO3− isotopic composition in the upper 150 m, with values increasing in a NO3− δ18O versus δ15N space along a line with a slope of one toward the surface. In the subpolar gyre, a single relationship between the degree of NO3− consumption and residual NO3− δ15N supports the view that NO3− is supplied via Ekman upwelling and deep winter convection, and progressively consumed during the Ekman transport of surface water southward. The co-occurrence of partial NO3− assimilation and nitrification in the deep mixed layer of the subpolar gyre elevates subsurface NO3− δ18O in comparison to deep oceanic values. This signal propagates through isopycnal exchanges to greater depths at lower latitudes. With recirculation in the subtropical gyre, cycles of quantitative consumption-nitrification progressively decrease subsurface NO3− δ18O toward the δ18O of regenerated NO3−. The low NO3− δ15N observed south of the Subarctic Front is mostly explained by N2 fixation, although a contribution from the Mediterranean outflow is required to explain the lower NO3− δ15N signal observed between 600 and 1500 m depth close to the Iberian margin.

Description: The GEOVIDE project was co-funded by the French national program LEFE/INSU (GEOVIDE), ANR Blanc (GEOVIDE) and RPDOC, LabEX MER and IFREMER. F. Deman was supported by the Belgian Federal Science Policy Office (Belspo contract BL/12/C63) while writing the manuscript. This work was financed by Flanders Research Foundation (FWO contract G0715.12N) and Vrije Universiteit Brussel, R&D, Strategic Research Plan “Tracers of Past & Present Global Changes”. During the preparation of the manuscript, Debany Fonseca-Batista was supported by funding from the Canada First Research Excellence Fund, through an International Postdoctoral Fellowship of the Ocean Frontier Institute (OFI) at Dalhousie University.

Description: 2021-10-02

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Tectonics 39(11), (2020): e2020TC006409, doi:10.1029/2020TC006409.

Description: The dynamics of continental breakup at convergent margins has been described as the results of backarc opening caused by slab rollback or drag force induced by subduction direction reversal. Although the rollback hypothesis has been intensively studied, our understanding of the consequence of subduction direction reversal remains limited. Using thermo‐mechanical modeling based on constraints from the South China Sea (SCS) region, we investigate how subduction direction reversal controls the breakup of convergent margins. The numerical results show that two distinct breakup modes, namely, continental interior and edge breakup (“edge” refers to continent above the plate boundary interface), may develop depending on the “maturity” of the convergent margin and the age of the oceanic lithosphere. For a slab age of ~15 to ~45 Ma, increasing the duration of subduction promotes the continental interior breakup mode, where a large block of the continental material is separated from the overriding plate. In contrast, the continental edge breakup mode develops when the subduction is a short‐duration event, and in this mode, a wide zone of less continuous continental fragments and tearing of the subducted slab occur. These two modes are consistent with the interior (relic late Mesozoic arc) and edge (relic forearc) rifting characteristics in the western and eastern SCS margin, suggesting that variation in the northwest‐directed subduction duration of the Proto‐SCS might be a reason for the differential breakup locus along the strike of the SCS margin. Besides, a two‐segment trench associated with the northwest‐directed subduction is implied in the present‐day SCS region.

Description: This research was supported by the Guangdong NSF research team project (2017A030312002), the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0205), the K. C. Wong Education Foundation (GJTD‐2018‐13), the Strategic Priority Research Program of the Chinese Academy of Science (XDA13010303), the Chinese Academy of Sciences (Y4SL021001, QYZDY‐SSWDQC005, 133244KYSB20180029, and ISEE2019ZR01), the NSFC project (41606073, 41890813, and 41576070), the IODP‐China Foundation, the OMG Visiting Fellowship (OMG18‐15), and the Hong Kong Research Grant Council Grants (Nos. 14313816 and 14304820).

Description: 2021-04-06

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 125(12),(2020): e2020JB020040, https://doi.org/10.1029/2020JB020040.

Description: We model the magnetic signature of rift‐related volcanism to understand the distribution and volume of magmatic activity that occurred during the breakup of Pangaea and early Atlantic opening at the Eastern North American Margin (ENAM). Along‐strike variations in the amplitude and character of the prominent East Coast Magnetic Anomaly (ECMA) suggest that the emplacement of the volcanic layers producing this anomaly similarly varied along the margin. We use three‐dimensional magnetic forward modeling constrained by seismic interpretations to identify along‐margin variations in volcanic thickness and width that can explain the observed amplitude and character of the ECMA. Our model results suggest that the ECMA is produced by a combination of both first‐order (~600–1,000 km) and second‐order (~50–100 km) magmatic segmentation. The first‐order magmatic segmentation could have resulted from preexisting variations in crustal thickness and rheology developed during the tectonic amalgamation of Pangaea. The second‐order magmatic segmentation developed during continental breakup and likely influenced the segmentation and transform fault spacing of the initial, and modern, Mid‐Atlantic Ridge. These variations in magmatism show how extension and thermal weakening was distributed at the ENAM during continental breakup and how this breakup magmatism was related to both previous and subsequent Wilson cycle stages.

Description: Thanks to Anne Bécel, Dan Lizarralde, Collin Brandl, Brandon Shuck, and Mark Everett for beneficial discussion and assistance in compiling the archived data used in this study. We thank Debbie Hutchinson (USGS Woods Hole Coastal and Marine Science Center) for passing along her vast breadth of knowledge on the ENAM through numerous constructive suggestions to greatly strengthen our manuscript. We greatly appreciate the insightful comments from two reviewers, the Associate Editor, and the Editor that significantly improved the manuscript. Thanks to Maurice Tivey for providing codes that aided our magnetic modeling efforts. Project completed as part of J.A.G.'s Ph.D. dissertation at Texas A&M University.

Description: 2021-05-16

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(12), (2020): e2020JC016543, https://doi.org/10.1029/2020JC016543.

Description: On coral reefs, flow determines residence time of water influencing physical and chemical environments and creating observable microclimates within the reef structure. Understanding the physical mechanisms driving environmental variability on shallow reefs, which distinguishes them from the open ocean, is important for understanding what contributes to thermal resilience of coral communities and predicting their response to future anomalies. In June 2014, a field experiment conducted at Dongsha Atoll in the northern South China Sea investigated the physical forces that drive flow over a broad shallow reef flat. Instrumentation included current and pressure sensors and a distributed temperature sensing system, which resolved spatially and temporally continuous temperature measurements over a 3‐km cross‐reef section from the lagoon to reef crest. Spectral analysis shows that while diurnal variability was significant across the reef flat—a result expected from daily solar heating—temperature also varied at higher frequencies near the reef crest. These spatially variable temperature regimes, or thermal microclimates, are influenced by circulation on the wide reef flat, with spatially and temporally variable contributions from tides, wind, and waves. Through particle tracking simulations, we find the residence time of water is shorter near the reef crest (3.6 h) than near the lagoon (8.6 h). Tidal variability in flow direction on the reef flat leads to patterns in residence time that are different than what would be predicted from unidirectional flow. Circulation on the reef also determines the source (originating from offshore vs. the lagoon) of the water present on the reef flat.

Description: We thank S. Tyler, and J. Selker from the Center for Transformative Environmental Monitoring Programs (CTEMPs), funded by the National Science Foundation (EAR awards 1440596 and 1440506), for timely and effective provision of experimental design support, logistical support and equipment for the project. Support for S. Lentz is from NSF Grant No. OCE‐1558343. Support for A. Cohen from NSF Grant No. 1220529, by the Academia Sinica (Taiwan) through a thematic project grant to G. Wong and A. Cohen. Support for E. Reid and K. Davis is from National Science Foundation (NSF) Grant No. OCE‐1753317, and support to E. Reid from the Environmental Engineering Henry Samueli Endowed Fellowship and the UCI Oceans Graduate Fellowship.

Description: 2021-05-23

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Paleoceanography and Paleoclimatology 35(12), (2020): e2020PA003853, https://doi.org/10.1029/2020PA003853.

Description: During the middle Miocene, Earth's climate changed from a global warm period (Miocene Climatic Optimum) into a colder one with the expansion of the Antarctic ice sheet. This prominent climate transition was also a period of drastic changes in global atmospheric circulation. The development of the South Asian monsoon is not well understood and mainly derived from proxy records of wind strength. Data for middle Miocene changes in rainfall are virtually non‐existent for India and the Arabian Sea prior to 11 Ma. This study presents planktic foraminiferal trace element (Mg/Ca and Ba/Ca) and stable oxygen isotope records from NGHP‐01 Site 01A off the coast of West India in the Eastern Arabian Sea (EAS) to reconstruct the regional surface hydrography and hydroclimate in the South Asian monsoon (SAM) region during the middle Miocene. The Ba/Ca and local seawater δ18O (δ18Osw) changes reveal a notable gradual increase in SAM rainfall intensity during the middle Miocene. Additionally to this long‐term increase in precipitation, the seawater δ18O is punctuated by a prominent decrease, i.e. freshening, at ~14 Ma contemporary with Antarctic glaciation. This suggests that Southern Ocean Intermediate Waters (SOIW) transmitted Antarctic salinity changes into the Arabian Sea via an “oceanic tunnel” mechanism. The middle Miocene increase in SAM rainfall is consistent with climate model simulations of an overall strengthening Asian monsoon from the Eocene to the middle/late Miocene with a further acceleration after the middle Miocene climate transition.

Description: This study has been funded by the National Natural Science Foundation of China through a grant to S. Steinke (NSFC grant No. 41776055) and Z. Jian and S. Steinke (NSFC grant No. 919582080). We express our gratitude to H. Kuhnert (MARUM, University of Bremen) and his team for stable isotope analyses. We thank P. Qiao (Tongji University Shanghai) for technical and analytical support with the ICP‐MS analyses, A. Dolman (Alfred‐Wegener‐Institute, Helmholtz Centre for Polar and Marine Research, Potsdam, Germany) for statistical analyses, and B. Wang (State Key Laboratory of Marine Environmental Science, Xiamen University) and his team for the SEM‐EDAX Energy Dispersive X‐ray Spectroscopy (EDS) analyses. L. Giosan acknowledges funding from USSP and WHOI and thanks colleagues and crew from the NGHP‐01 expedition for intellectual interactions leading to long‐standing interests in the fluvial‐continental margin systems of Peninsular India. J. Groeneveld thanks the State Key Laboratory of Marine Environmental Science (Xiamen University) for a MEL Senior Visiting Fellowship (Project No. MELRS1915).

Description: 2021-05-27

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 21(12), (2020): e2020GC008914, https://doi.org/10.1029/2020GC008914.

Description: Rarely have small seamounts on the flanks of hotspot derived ocean‐island volcanoes been the targets of sampling, due to sparse high‐resolution mapping near ocean islands. In the Galápagos Archipelago, for instance, sampling has primarily targeted the subaerial volcanic edifices, with only a few studies focusing on large‐volume submarine features. Sampling restricted to these large volcanic features may present a selection bias, potentially resulting in a skewed view of magmatic and source processes because mature magmatic systems support mixing and volcanic accretion that overprints early magmatic stages. We demonstrate how finer‐scale sampling of satellite seamounts surrounding the volcanic islands in the Galápagos can be used to lessen this bias and thus, better constrain the evolution of these volcanoes. Seamounts were targeted in the vicinity of Floreana and Fernandina Islands, and between Santiago and Santa Cruz. In all regions, individual seamounts are typically monogenetic, but each seamount field requires multigenerational magmatic episodes to account for their geochemical variability. This study demonstrates that in the southern and eastern regions the seamounts are characterized by greater geochemical variability than the islands they surround but all three regions have (Sr‐Nd‐He) isotopic signatures that resemble neighboring islands. Variations in seamount chemistry from alkalic to tholeiitic near Fernandina support the concept that islands along the center of the hotspot track undergo greater mean depths of melting, as predicted by plume theory. Patterns of geochemical and isotopic enrichment of seamounts within each region support fine‐scale mantle heterogeneities in the mantle plume sourcing the Galápagos hotspot.

Description: This work was carried out with funding from National Science Foundation Division of Ocean Sciences (OCE‐1634952 to V. D. Wanless, OCE‐1634685 to S. A. Soule). The authors have no competing interests to declare. We thank Sally Gibson and three anonymous reviewers for providing detailed and critical feedback on this manuscript.

Description: 2021-05-06

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(1),(2021): e2020JB020982, https://doi.org/10.1029/2020JB020982.

Description: Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths 〉1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to 〈100‐m in width at depths 〉1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to 〈0.9 wt% H20. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.

Description: National Science Foundation. Grant Numbers: OCE-0625178, OCE-0841063

Description: 2021-06-15

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(10),(2020): e2020JC016507, https://doi.org/10.1029/2020JC016507.

Description: Survival of Gulf Stream (GS) warm core rings (WCRs) was investigated using a census consisting of a total of 961 rings formed during the period 1980–2017. Kaplan‐Meier survival probability and Cox hazard proportional models were used for the analysis. The survival analysis was performed for rings formed in four 5° zones between 75° W and 55° W. The radius, latitude, and distance from the shelf‐break of a WCR at formation all had a significant effect on the survival of WCRs. A pattern of higher survival was observed in WCRs formed in Zone 2 (70°–65° W) or Zone 3 (65°–60° W) and then demised in Zone 1 (75°–70° W). Survival probability of the WCRs increased to more than 70% for those formed within a latitude band from 39.5° to 41.5° N. Survival probability is reduced when the WCRs are formed near the New England Seamounts.

Description: We are grateful for financial supports from NOAA (NA11NOS0120038), NSF (OCE‐1851242), SMAST, and UMass Dartmouth. G. G. was supported by NSF under grant OCE‐1851261.

Description: 2021-04-14

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Anderson, M., Wanless, V. D., Perfit, M., Conrad, E., Gregg, P., Fornari, D., & Ridley, W. I. Extreme heterogeneity in mid-ocean ridge mantle revealed in lavas from the 8 degrees 20 ' N near-axis seamount chain. Geochemistry Geophysics Geosystems, 22(1), (2021): e2020GC009322, https://doi.org/10.1029/2020GC009322.

Description: Lavas that have erupted at near‐axis seamounts provide windows into mid‐ocean ridge mantle heterogeneity and melting systematics which are not easily observed on‐axis at fast‐spreading centers. Beneath ridges, most heterogeneity is obscured as magmas aggregate toward the ridge, where they efficiently mix and homogenize during transit and within shallow magma chambers prior to eruption. To understand the deeper magmatic processes contributing to oceanic crustal formation, we examine the compositions of lavas erupted along a chain of near‐axis seamounts and volcanic ridges perpendicular to the East Pacific Rise. We assess the chemistry of near‐ridge mantle using a ∼200 km‐long chain at ∼8°20′N. High‐resolution bathymetric maps are used with geochemical analyses of ∼300 basalts to evaluate the petrogenesis of lavas and the heterogeneity of mantle feeding these near‐axis eruptions. Major and trace element concentrations and radiogenic isotope ratios are highly variable on 〈1 km scales, and reveal a continuum of depleted, normal, and enriched basalts spanning the full range of ridge and seamount compositions in the northeast Pacific. There is no systematic compositional variability along the chain. Modeling suggests that depleted mid‐ocean ridge basalt (DMORB) lavas are produced by ∼5%–15% melting of a depleted mid‐ocean ridge (MOR) mantle. Normal mid‐ocean ridge basalts (NMORB) form from 5% to 15% melting of a slightly enriched MOR mantle. Enriched mid‐ocean ridge basalts (EMORB) range from 〈1% melting of 10% enriched mantle to 〉15% melting of 100% enriched mantle. The presence of all three lava types along the seamount chain, and on a single seamount closest to the ridge axis, confirms that the sub‐ridge mantle is much more heterogeneous than is commonly observed on‐axis and heterogeneity exists over small spatial scales.

Description: This work was supported by NSF OCE‐MGG 1356610 (Romano and Gregg), NSF OCE‐MGG 1356822 (Fornari), NSF OCE‐MGG 1357150 (Perfit), NSF OCE‐MGG 2001314 (Perfit and Wanless), the Burnham Research Grant at Boise State University, and the Graduate School Funding Fellowship at University of Florida.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Sutherland, K. M., Grabb, K. C., Karolewski, J. S., Plummer, S., Farfan, G. A., Wankel, S. D., Diaz, J. M., Lamborg, C. H., & Hansel, C. M. Spatial heterogeneity in particle-associated, light-independent superoxide production within productive coastal waters. Journal of Geophysical Research: Oceans, 125(10), (2020): e2020JC016747, https://doi.org/10.1029/2020JC016747.

Description: In the marine environment, the reactive oxygen species (ROS) superoxide is produced through a diverse array of light‐dependent and light‐independent reactions, the latter of which is thought to be primarily controlled by microorganisms. Marine superoxide production influences organic matter remineralization, metal redox cycling, and dissolved oxygen concentrations, yet the relative contributions of different sources to total superoxide production remain poorly constrained. Here we investigate the production, steady‐state concentration, and particle‐associated nature of light‐independent superoxide in productive waters off the northeast coast of North America. We find exceptionally high levels of light‐independent superoxide in the marine water column, with concentrations ranging from 10 pM to in excess of 2,000 pM. The highest superoxide concentrations were particle associated in surface seawater and in aphotic seawater collected meters off the seafloor. Filtration of seawater overlying the continental shelf lowered the light‐independent, steady‐state superoxide concentration by an average of 84%. We identify eukaryotic phytoplankton as the dominant particle‐associated source of superoxide to these coastal waters. We contrast these measurements with those collected at an off‐shelf station, where superoxide concentrations did not exceed 100 pM, and particles account for an average of 40% of the steady‐state superoxide concentration. This study demonstrates the primary role of particles in the production of superoxide in seawater overlying the continental shelf and highlights the importance of light‐independent, dissolved‐phase reactions in marine ROS production.

Description: This work was funded by grants from the Chemical Oceanography program of the National Science Foundation (OCE‐1355720 to C. M. H. and C. H. L.), NASA Earth and Space Science Fellowship (Grant NNX15AR62H to K. M. S.), Agouron Institute Postdoctoral Fellowship (K. M. S.), NSF GRFPs (2016230268 to K. C. G. and 2017250547 to S. P.), and a Sloan Research Fellowship (J. M. D.). The Guava flow cytometer was purchased through an NSF equipment improvement grant (1624593).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Paleoceanography and Paleoclimatology 36(7), (2021): e2020PA004088, https://doi.org/10.1029/2020PA004088.

Description: We reconstruct deep water-mass salinities and spatial distributions in the western North Atlantic during the Last Glacial Maximum (LGM, 19–26 ka), a period when atmospheric CO2 was significantly lower than it is today. A reversal in the LGM Atlantic meridional bottom water salinity gradient has been hypothesized for several LGM water-mass reconstructions. Such a reversal has the potential to influence climate, ocean circulation, and atmospheric CO2 by increasing the thermal energy and carbon storage capacity of the deep ocean. To test this hypothesis, we reconstructed LGM bottom water salinity based on sedimentary porewater chloride profiles in a north-south transect of piston cores collected from the deep western North Atlantic. LGM bottom water salinity in the deep western North Atlantic determined by the density-based method is 3.41–3.99 ± 0.15% higher than modern values at these sites. This increase is consistent with: (a) the 3.6% global average salinity change expected from eustatic sea level rise, (b) a northward expansion of southern sourced deep water, (c) shoaling of northern sourced deep water, and (d) a reversal of the Atlantic's north-south deep water salinity gradient during the LGM.

Description: This work was supported by the US National Science Foundation (grant numbers 1433150 and 1537485).

Description: 2021-10-24

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Paleoceanography and Paleoclimatology 36(9), (2021): e2021PA004226, https://doi.org/10.1029/2021PA004226.

Description: The deep ocean has long been recognized as the reservoir that stores the carbon dioxide (CO2) removed from the atmosphere during Pleistocene glacial periods. The removal of glacial atmospheric CO2 into the ocean is likely modulated by an increase in the degree of utilization of macronutrients at the sea surface and enhanced storage of respired CO2 in the deep ocean, known as enhanced efficiency of the biological pump. Enhanced biological pump efficiency during glacial periods is most easily documented in the deep ocean using proxies for oxygen concentrations, which are directly linked to respiratory CO2 levels. We document the enhanced storage of respired CO2 during the Last Glacial Maximum (LGM) in the Pacific Southern Ocean and deepest Equatorial Pacific using records of deglacial authigenic manganese, which form as relict peaks during increases in bottom water oxygen (BWO) concentration. These peaks are found at depths and regions where other oxygenation histories have been ambiguous, due to diagenetic alteration of authigenic uranium, another proxy for BWO. Our results require that the entirety of the abyssal Pacific below approximately 1,000 m was enriched in respired CO2 and depleted in oxygen during the LGM. The presence of authigenic Mn enrichment in the deep Equatorial Pacific for each of the last five deglaciations suggests that the storage of respired CO2 in the deep ocean is a ubiquitous feature of late-Pleistocene ice ages.

Description: This work was performed with support from the National Science Foundation (NSF) over about 30 years. The TT013 and NBP9802 cores were collected during the U.S. JGOFS program. Their collection and analyses were supported by NSF OCE-9022301 and OPP-95303398 to R. F. Anderson, and NSF OCE 9301097 to R. W. Murray. Coring and radiocarbon analyses on NBP1702 were funded by NSF OPP-1542962. XRF analysis on NBP9802 and NBP1702 cores, as well as additional radiocarbon measurements, was funded by an LDEO Climate Center Grant to F. J. Pavia.

Description: 2022-02-17

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(2), (2021): e2020JC016773, https://doi.org/10.1029/2020JC016773.

Description: A new modeling methodology for ripple dynamics driven by oscillatory flows using a Eulerian two‐phase flow approach is presented in order to bridge the research gap between near‐bed sediment transport via ripple migration and suspended load transport dictated by ripple induced vortices. Reynolds‐averaged Eulerian two‐phase equations for fluid phase and sediment phase are solved in a two‐dimensional vertical domain with a k‐ε closure for flow turbulence and particle stresses closures for short‐lived collision and enduring contact. The model can resolve full profiles of sediment transport without making conventional near‐bed load and suspended load assumptions. The model is validated with an oscillating tunnel experiment of orbital ripple driven by a Stokes second‐order (onshore velocity skewed) oscillatory flow with a good agreement in the flow velocity and sediment concentration. Although the suspended sediment concentration far from the ripple in the dilute region was underpredicted by the present model, the model predicts an onshore ripple migration rate that is in very good agreement with the measured value. Another orbital ripple case driven by symmetric sinusoidal oscillatory flow is also conducted to contrast the effect of velocity skewness. The model is able to capture a net offshore‐directed suspended load transport flux due to the asymmetric primary vortex consistent with laboratory observation. More importantly, the model can resolve the asymmetry of onshore‐directed near‐bed sediment flux associated with more intense boundary layer flow speed‐up during onshore flow cycle and sediment avalanching near the lee ripple flank which force the onshore ripple migration.

Description: This study is supported by National Science Foundation (Grant no. OCE‐1635151) and Strategic Environmental Research and Development Program (Grant no. MR20‐1478).

Description: 2021-06-29

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(2), (2021): e2020JC016856, https://doi.org/10.1029/2020JC016856.

Description: The genus Phaeocystis is distributed globally and has considerable ecological, biogeochemical, and societal impacts. Understanding its distribution, growth and ecological impacts has been limited by lack of extensive observations on appropriate scales. In 2018, we investigated the biological dynamics of the New England continental shelf and encountered a substantial bloom of Phaeocystis pouchetii. Based on satellite imagery during January through April, the bloom extended over broad expanses of the shelf; furthermore, our observations demonstrated that it reached high biomass levels, with maximum chlorophyll concentrations exceeding 16 µg L−1 and particulate organic carbon levels 〉 95 µmol L−1. Initially, the bloom was largely confined to waters with temperatures 〈6°C, which in turn were mostly restricted to shallow areas near the coast. As the bloom progressed, it appeared to sink into the bottom boundary layer; however, enough light and nutrients were available for growth. The bloom was highly productive (net community production integrated through the mixed layer from stations within the bloom averaged 1.16 g C m−2 d−1) and reduced nutrient concentrations considerably. Long‐term coastal observations suggest that Phaeocystis blooms occur sporadically in spring on Nantucket Shoals and presumably expand onto the continental shelf. Based on the distribution of Phaeocystis during our study, we suggest that it can have a significant impact on the overall productivity and ecology of the New England shelf during the winter/spring transition.

Description: This project was supported by the US National Science Foundation (Grants 1657855, 1657803, and 1657489). NES‐LTER contributions were supported by grants to HMS from NSF (Grant 1655686) and the Simons Foundation (Grant 561126). VPR operations were supported by the Dalio Explore Fund.

Description: 2021-07-15

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Water Resources Research 57(7), (2021): e2020WR028727, https://doi.org/10.1029/2020WR028727.

Description: Numerous wetlands in the prairies of Canada provide important ecosystem services, yet are threatened by climate and land-use changes. Understanding the impacts of climate change on prairie wetlands is critical to effective conservation planning. In this study, we construct a wetland model with surface water balance and ecoregions to project future distribution of wetlands. The climatic conditions downscaled from the Weather Research and Forecasting model were used to drive the Noah-MP land surface model to obtain surface water balance. The climate change perturbation is derived from an ensemble of general circulation models using the pseudo global warming method, under the RCP8.5 emission scenario by the end of 21st century. The results show that climate change impacts on wetland extent are spatiotemporally heterogenous. Future wetter climate in the western Prairies will favor increased wetland abundance in both spring and summer. In the eastern Prairies, particularly in the mixed grassland and mid-boreal upland, wetland areas will increase in spring but experience enhanced declines in summer due to strong evapotranspiration. When these effects of climate change are considered in light of historical drainage, they suggest a need for diverse conservation and restoration strategies. For the mixed grassland in the western Canadian Prairies, wetland restoration will be favorable, while the highly drained eastern Prairies will be challenged by the intensified hydrological cycle. The outcomes of this study will be useful to conservation agencies to ensure that current investments will continue to provide good conservation returns in the future.

Description: Z. Zhang was funded by a Mitacs Accelerate Fellowship funded by Ducks Unlimited Canada's Institute for Wetland and Waterfowl Research. Z. Zhang, Z. Li, and Y. Li acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant, and Global Water Futures Program, Canada First Research Excellence Fund. This project was supported by grants from Wildlife Habitat Canada, Bass Pro Shops Cabela’s Outdoor Fund, and the Alberta NAWMP Partnership.

Description: 2021-12-21

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(11), (2021): e2021GL093178, https://doi.org/10.1029/2021GL093178.

Description: The effects of heterogeneous reactions between river-borne particles and the carbonate system were studied in the plumes of the Mississippi and Brazos rivers. Measurements within these plumes revealed significant removal of dissolved inorganic carbon (DIC) and total alkalinity (TA). After accounting for all known DIC and TA sinks and sources, heterogeneous reactions (i.e., heterogeneous CaCO3 precipitation and cation exchange between adsorbed and dissolved ions) were found to be responsible for a significant fraction of DIC and TA removal, exceeding 10% and 90%, respectively, in the Mississippi and Brazos plume waters. This finding was corroborated by laboratory experiments, in which the seeding of seawater with the riverine particles induced the removal of the DIC and TA. The combined results demonstrate that heterogeneous reactions may represent an important controlling mechanism of the seawater carbonate system in particle-rich coastal areas and may significantly impact the coastal carbon cycle.

Description: This research was funded by the National Science Foundation (NSF) and the Bi-National Science Foundation U.S-Israel award number OCE-BSF 1635388.

Description: 2021-11-20

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(19), (2021): e2021GL094364, https://doi.org/10.1029/2021GL094364.

Description: The warm Gulf Stream sea surface temperatures strongly impact the evolution of winter clouds behind atmospheric cold fronts. Such cloud evolution remains challenging to model. The Gulf Stream is too wide within the ERA5 and MERRA2 reanalyses, affecting the turbulent surface fluxes. Known problems within the ERA5 boundary layer (too-dry and too-cool with too strong westerlies), ascertained primarily from ACTIVATE 2020 campaign aircraft dropsondes and secondarily from older buoy measurements, reinforce surface flux biases. In contrast, MERRA2 winter surface winds and air-sea temperature/humidity differences are slightly too weak, producing surface fluxes that are too low. Reanalyses boundary layer heights in the strongly forced winter cold-air-outbreak regime are realistic, whereas late-summer quiescent stable boundary layers are too shallow. Nevertheless, the reanalysis biases are small, and reanalyses adequately support their use for initializing higher-resolution cloud process modeling studies of cold-air outbreaks.

Description: This work was supported by NASA grant 80NSSC19K0390 to ACTIVATE, a NASA Earth Venture Suborbital-3 (EVS-3) investigation funded by NASA's Earth Science Division and managed through the Earth System Science Pathfinder Program Office. The Pacific Northwest National Laboratory (PNNL) is operated for the US Department of Energy (DOE) by Battelle Memorial Institute under Contract DE-AC06-76RLO 1830.

Description: 2022-03-08

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Loranger, S., & Weber, T. C. . Shipboard acoustic observations of flow rate from a seafloor-sourced oil spill. Journal of Geophysical Research: Oceans, 125(10), (2020): e2020JC016274, https://doi.org/10.1029/2020JC016274.

Description: In 2004 a debris flow generated by Hurricane Ivan toppled an oil production platform in Mississippi Canyon lease block 20 (MC20). Between 2004 and the installation of a containment system in 2019 MC20 became an in situ laboratory for a wide range of hydrocarbon in the sea‐related research, including different methods of assessing the volumetric flow rate of hydrocarbons spanning different temporal scales. In 2017 a shipboard acoustic Doppler current profiler (ADCP) and high‐frequency (90 to 154 kHz) broadband echosounder were deployed to assess the flow rate of liquid and gas phase hydrocarbons. Measurements of horizontal currents were combined with acoustic mapping to determine the rise velocity of the seep as it moved downstream. Models of the rise velocity for fluid particles were used to predict the size of oil droplets and gas bubbles in the seep. The amplitude and shape of the broadband acoustic backscatter were then used to differentiate between, and determine the flow rate of, hydrocarbons. Oil flow rate in the seep was estimated to be 56 to 86 barrels/day (mean urn:x-wiley:jgrc:media:jgrc24228:jgrc24228-math-0001 barrels/day) while the flow rate of gaseous hydrocarbons was estimated to be 98 to 359 m3/day (mean urn:x-wiley:jgrc:media:jgrc24228:jgrc24228-math-0002 m3/day).

Description: The work was supported by the National Oceanic and Atmospheric Administration (Grant NA15NOS4000200).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(1), (2021): e2020JC016771, https://doi.org/10.1029/2020JC016771.

Description: Estuaries worldwide have experienced modifications including channel deepening and intertidal reclamation over several centuries, resulting in altered fine sediment routing. Estuaries respond differently based on preexisting geometries, freshwater and sediment supplies, and extents and types of modification. The Coos Bay Estuary in Oregon is a relatively small estuary with complex geometry that has been extensively modified since 1865. A sediment transport model calibrated to modern conditions is used to assess the corresponding changes in sediment dynamics. Over ∼150 years, channel deepening (from ∼6.7 to 11 m), a 12% increase in area, and a 21% increase in volume have led to greater tidal amplitudes, salinity intrusion, and estuarine exchange flow. These changes have reduced current magnitudes, reduced bed stresses, and increased stratification, especially during rainy periods. Historically, fluvially derived sediment was dispersed across broad, deltaic‐style flats and through small tidal channels. Now, river water and sediments are diverted into a dredged navigation channel where an estuarine turbidity maximum (ETM) forms, with modeled concentrations 〉50 mg/L and measured concentrations 〉100 mg/L during discharge events. This “new” ETM supplies sediment to proximal embayments in the middle estuary and the shallow flats. Overall, sediment trapping during winter (and high river discharges) has increased more than two‐fold, owing to increased accommodation space, altered pathways of supply, and altered bed stresses and tidal asymmetries. In contrast to funnel‐shaped estuaries with simpler geometries and river‐channel transitions, these results highlight the importance of channel routing together with dredging in enhancing sediment retention and shifting pathways of sediment delivery.

Description: The Science Collaborative is funded by the National Oceanic and Atmospheric Administration and managed by the University of Michigan Water Center (NAI4NOS4190145).

Description: 2021-06-14

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Geyer, W. R., Ralston, D. K., & Chen, J. Mechanisms of exchange flow in an estuary with a narrow, deep channel and wide, shallow shoals. Journal of Geophysical Research: Oceans, 125(12), (2020): e2020JC016092, https://doi.org/10.1029/2020JC016092.

Description: Delaware Bay is a large estuary with a deep, relatively narrow channel and wide, shallow banks, providing a clear example of a “channel‐shoal” estuary. This numerical modeling study addresses the exchange flow in this channel‐shoal estuary, specifically to examine how the lateral geometry affects the strength and mechanisms of exchange flow. We find that the exchange flow is exclusively confined to the channel region during spring tides, when stratification is weak, and it broadens laterally over the shoals during the more stratified neap tides but still occupies a small fraction of the total width of the estuary. Exchange flow is relatively weak during spring tides, resulting from oscillatory shear dispersion in the channel augmented by weak Eulerian exchange flow. During neap tides, stratification and shear increase markedly, resulting in a strong Eulerian residual shear flow driven mainly by the along‐estuary density gradient, with a net exchange flow roughly 5 times that of the spring tide. During both spring and neap tides, lateral salinity gradients generated by differential advection at the edge of the channel drive a tidally oscillating cross‐channel flow, which strongly influences the stratification, along‐estuary salt balance, and momentum balance. The lateral flow also causes the phase variation in salinity that results in oscillatory shear dispersion and is an advective momentum source contributing to the residual circulation. Whereas the shoals make a negligible direct contribution to the exchange flow, they have an indirect influence due to the salinity gradients between the channel and the shoal.

Description: The ideas in this paper were influenced by discussions with Robert Chant. Funding was provided by National Science Foundation grants OCE‐1325136, OCE‐1634490, and OCE‐1736539.

Description: 2021-04-29

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(3), (2021): e2020GL091200, doi:10.1029/2020GL091200.

Description: Tabular calving events occur from Antarctica's large ice shelves only every few decades, and are preceded by rift propagation. We used high-resolution imagery and ICESat-2 data to determine the propagation rates for the three active rifts on Amery Ice Shelf (AIS; T1, T2, and E3) and observe the calving of D-28 on September 25, 2019 along T1. AIS front advance accelerated downstream of T1 in the years before calving, possibly increasing stress at the rift tip. T1 experienced significant acceleration for 12 days before calving, coinciding with a jump in propagation of E3. ICESat-2's high resolution and repeat acquisitions every 91 days allowed for analysis of the ice front before and after calving, and rift detection where it was not visible in imagery as a ∼1 m surface depression, suggesting that it propagates as a basal fracture. Our results show that ICESat-2 can provide process-scale information about iceberg calving.

Description: We received funding from the following sources: NASA NNX15AC80G and NSF grant 1443677 (Fricker and Becker) and NASA 80NSSC20K0960 (Walker). We are grateful to Mike Cloutier, Polar Geospatial Center for assistance with WorldView imagery. Geospatial support was provided by Polar Geospatial Center under NSF-OPP awards 1043681 and 1559691.

Description: 2021-07-12

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(7), (2021): e2020JC016899, https://doi.org/10.1029/2020JC016899.

Description: Circulation in the nearshore region, which is critical for material transport along the coast and between the surf zone and the inner shelf, includes strong vortical motions. The horizontal length scales and vertical structure associated with vortical motions are not well documented on alongshore-variable beaches. Here, a three-dimensional phase-resolving numerical model, Simulating WAves till SHore (SWASH), is compared with surfzone waves and flows on a barred beach, and is used to investigate surfzone eddies. Model simulations with measured bathymetry reproduce trends in the mean surfzone circulation patterns, including alongshore currents and rip current circulation cells observed for offshore wave heights from 0.5 to 2.0 m and incident wave directions from 0 to 15° relative to shore normal. The length scales of simulated eddies, quantified using the alongshore wavenumber spectra of vertical vorticity, suggest that increasing wave directional spread intensifies small-scale eddies ( (10) m). Simulations with bathymetric variability ranging from alongshore uniform to highly alongshore variable indicate that large-scale eddies ( (100) m) may be enhanced by surfzone bathymetric variability, whereas small-scale eddies ( (10) m) are less dependent on bathymetric variability. The simulated vertical dependence of the magnitude and mean length scale (centroid) of the alongshore wavenumber spectra of vertical vorticity and very low-frequency (f ≈ 0.005 Hz) currents is weak in the outer surf zone, and decreases toward the shoreline. The vertical dependence in the simulations may be affected by the vertical structure of turbulence, mean shear, and bottom boundary layer dynamics.

Description: Support was provided by the University of Washington Royalty Research Fund, the National Science Foundation, the Office of Naval Research, a National Defense Science and Engineering Graduate Fellowship, a Vannevar Bush Faculty Fellowship, the United States Army Corps of Engineers, the United States Coastal Research Program, Sea Grant, and the WHOI Investment in Science Program.

Description: 2021-12-26

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in O’Dea, A., Brodie, K., & Elgar, S. Field observations of the evolution of plunging-wave shapes. Geophysical Research Letters, 48(16), (2021): e2021GL093664, https://doi.org/10.1029/2021GL093664.

Description: There are few high-resolution field observations of the water surface during breaking owing to the difficulty of collecting spatially dense measurements in the surf zone, and thus the factors influencing breaking-wave shape in field conditions remain poorly understood. Here, the shape and evolution of plunging breakers is analyzed quantitatively using three-dimensional scans of the water surface collected at high spatial and temporal resolution with a multi-beam terrestrial lidar scanner. The observed internal void shapes in plunging breakers agree well with previously developed theoretical shapes at the onset of breaking, and become more elongated and less steep as breaking progresses. The normalized void area increases as the local bottom slope steepens and as the breaking depth decreases. The void shape becomes more circular as the local bottom slope and the ratio of breaking water depth to wavelength increase, as well as in conditions with opposing winds.

Description: Funding was provided by the U.S. Department of Defense (DoD) Laboratory University Collaboration Initiative program, the U.S. Army ERDC Military Engineering Basic Research Program from the Assistant Secretary of the Army for Acquisition, Logistics, and Technology, the Vannevar Bush Faculty Fellowship program, the National Science Foundation, and the U.S. Coastal Research Program. This project was supported in part by an appointment to the Research Participation Program at the DoD, administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and the DoD.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Earth Surface 126(7), (2021): e2021JF006132, https://doi.org/10.1029/2021JF006132.

Description: The evolution of ripple geometries and their equilibrium states due to different wave forcing parameters are investigated by a Reynolds-averaged two-phase model, SedFoam, in a two-dimensional domain. Modeled ripple geometries, for a given uniform grain diameter, show a good agreement with ripple predictors that include the wave period effect explicitly, in addition to the wave orbital excursion length (or wave orbital velocity amplitude). Furthermore, using a series of numerical experiments, the ripple's response to a step-change in the wave forcing is studied. The model is capable of simulating “splitting,” “sliding,” “merging,” and “protruding” as the ripples evolve to a new equilibrium state. The model can also simulate the transition to sheet flow in energetic wave conditions and ripple reformation from a nearly flat bed condition. Simulation results reveal that the equilibrium state is such that the “primary” vortices reach half of the ripple length. Furthermore, an analysis of the suspended load and near-bed load ratio in the equilibrium state indicates that in the orbital ripple regime, the near-bed load is dominant while the suspended load is conducive to the ripple decaying regime (suborbital ripples) and sheet flow condition.

Description: This study is supported by National Science Foundation (OCE-1635151 and OCE-1924532) and Strategic Environmental Research and Development Program (MR20-1478).

Description: 2021-12-24

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(1), (2021): e2020JC016802, https://doi.org/10.1029/2020JC016802.

Description: The neodymium isotopic composition of the detrital (lithogenic) fraction (εNd‐detrital) of surface sediments and sinking particles was examined to constrain transport trajectories associated with hemipelagic sedimentation on the northwest Atlantic margin. The provenance of resuspended sediments and modes of lateral transport in the water column were of particular interest given the energetic hydrodynamic regime that sustains bottom and intermediate nepheloid layers over the margin. A large across‐margin gradient of ∼5 εNd units was observed for surface sediments, implying strong contrasts in sediment provenance, with εNd‐detrital values on the lower slope similar to those of “upstream regions” (Scotian margin) under the influence of the Deep Western Boundary Current (DWBC). Sinking particles collected at three depths at a site (total water depth, ∼3,000 m) on the New England margin within the core of the DWBC exhibited a similarly large range in εNd‐detrital values. The εNd‐detrital values of particles intercepted at intermediate water depths (1,000 and 2,000 m) were similar to each other but significantly higher than those at 3,000 m (∼50 m above the seafloor). These observations suggest that lithogenic material accumulating in the upper two traps was primarily advected in intermediate nepheloid layers emanating from the adjacent shelf, while that at 3,000 m is strongly influenced by sediment resuspension and along‐margin, southward lateral transport within the bottom nepheloid layer via entrainment in the DWBC. Our results highlight the importance of both along‐ and across‐margin sediment transport as vectors for lithogenic material and associated organic carbon transport.

Description: This research was funded by the NSF Ocean Sciences Chemical Oceanography program (OCE‐0425677; OCE‐0851350). JH was partly supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (2020R1A2C1008378).

Description: 2021-06-04

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(2), (2021): e2020GL090287, https://doi.org/10.1029/2020GL090287.

Description: Salt marsh survival with sea‐level rise (SLR) increasingly relies on soil organic carbon (SOC) accumulation and preservation. Using a novel combination of geochemical approaches, we characterized fine SOC (≤1 mm) supporting marsh elevation maintenance. Overlaying thermal reactivity, source (δ13C), and age (F14C) information demonstrates several processes contributing to soil development: marsh grass production, redeposition of eroded material, and microbial reworking. Redeposition of old carbon, likely from creekbanks, represented ∼9%–17% of shallow SOC (≤26 cm). Soils stored marsh grass‐derived compounds with a range of reactivities that were reworked over centuries‐to‐millennia. Decomposition decreases SOC thermal reactivity throughout the soil column while the decades‐long disturbance of ponding accelerated this shift in surface horizons. Empirically derived estimates of SOC turnover based on geochemical composition spanned a wide range (640–9,951 years) and have the potential to inform predictions of marsh ecosystem evolution.

Description: This work was supported by NSF (OCE1233678) and NOAA (NA14OAR4170104 and NA14NOS4190145) grants to ACS, USGS Coastal & Marine Geology Program, and PIE‐LTER (NSF OCE1238212 and OCE1637630).

Description: 2021-06-11

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Schloesser, F., Thompson, P. R., & Piecuch, C. G. Meridional asymmetry in recent decadal sea-level trends in the subtropical Pacific Ocean. Geophysical Research Letters, 48(6), (2021): e2020GL091959, https://doi.org/10.1029/2020GL091959.

Description: Recent sea surface height (SSH) trends in the South Pacific are substantially greater than trends in the North Pacific. Here, we use the Estimating the Climate and Circulation of the Ocean Version 4 Release 4 ocean state estimate and the Ocean Reanalysis System 5 to identify the forcing and mechanisms underlying that meridional asymmetry during 2005–2015. Thermosteric contributions dominate the spatial structure in Pacific SSH trends, but contributions from local surface heat fluxes are small. Wind stress trends drive a spin-up of the South Pacific subtropical gyre and a northward shift of the North Pacific subtropical gyre. A reduced gravity model forced with reanalysis winds qualitatively reproduces the meridional seesaw in sea level, suggesting that asymmetric trends in subtropical wind stress drive a cross-equatorial heat transport. A reversal in forcing associated with this process could impact near-term rates of coastal sea-level change, particularly in Pacific Island communities.

Description: F. Schloesser and P. R. Thompson were supported by NASA grant 80NSSC17K0564 and NSF grant 1558980. C. G. Piecuch was supported by the NASA Sea Level Change Team (grant 80NSSC20K1241).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(2), (2021): e2020JB020237, https://doi.org/10.1029/2020JB020237.

Description: Alteration of mantle peridotite in the Samail ophiolite forms secondary minerals, mainly serpentine and Mg‐rich carbonates. Magnesium accounts for ∼25 – 30% of peridotite mass and its mobility can be used to trace this alteration. We report the first set of Mg isotope measurements from peridotites and their alteration products in Oman. Partially serpentinized peridotites have Mg isotope ratios that are indistinguishable from estimates for the average mantle and bulk silicate earth (δ26Mg = −0.25 ± 0.04‰). However, more extensively altered peridotite samples show large shifts in Mg isotopic composition. The range of δ26Mg values for our suite of alteration products from the mantle section is ∼4.5‰ (from −3.39‰ to 1.19‰), or 〉60% of the total range of terrestrial variability in δ26Mg values. Serpentine veins are typically enriched in 26Mg (max δ26Mg value = 0.96‰) whereas Mg‐carbonate veins are associated with low 26Mg/24Mg ratios (magnesite δ26Mg = −3.3‰, dolomite δ26Mg = −1.91‰). Our preferred explanation for the range in δ26Mg values involves coprecipitation of serpentine and carbonates at water‐to‐rock ratios 〉103. The coincidence of alteration products characterized by δ26Mg values that are both lower and higher than bulk silicate Earth and the finite 14C ages of the carbonates suggest that both serpentinization and carbonation are ongoing in Oman. Rates of calcite precipitation in travertines inferred from Δ26Mgcal‐fl suggest that travertine formation in Oman sequesters a total of 106–107 kg CO2/yr, consistent with previous estimates.

Description: This work was supported through the Sloan Foundation–Deep Carbon Observatory (Grant 2014‐3‐01, Kelemen PI), the U.S.‐National Science Foundation (NSF‐EAR‐1516300, Kelemen lead PI).

Description: 2021-06-04

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Brenner, S., Rainville, L., Thomson, J., Cole, S., & Lee, C. Comparing observations and parameterizations of ice-ocean drag through an annual cycle across the Beaufort Sea. Journal of Geophysical Research: Oceans, 126(4), (202): 1e2020JC016977, https://doi.org/10.1029/2020JC016977.

Description: Understanding and predicting sea ice dynamics and ice-ocean feedback processes requires accurate descriptions of momentum fluxes across the ice-ocean interface. In this study, we present observations from an array of moorings in the Beaufort Sea. Using a force-balance approach, we determine ice-ocean drag coefficient values over an annual cycle and a range of ice conditions. Statistics from high resolution ice draft measurements are used to calculate expected drag coefficient values from morphology-based parameterization schemes. With both approaches, drag coefficient values ranged from ∼1 to 10 × 10−3, with a minimum in fall and a maximum at the end of spring, consistent with previous observations. The parameterizations do a reasonable job of predicting the observed drag values if the under ice geometry is known, and reveal that keel drag is the primary contributor to the total ice-ocean drag coefficient. When translations of bulk model outputs to ice geometry are included in the parameterizations, they overpredict drag on floe edges, leading to the inverted seasonal cycle seen in prior models. Using these results to investigate the efficiency of total momentum flux across the atmosphere-ice-ocean interface suggests an inter-annual trend of increasing coupling between the atmosphere and the ocean.

Description: This work was supported by the Office of Naval Research as part of the Stratified Ocean Dynamics of the Arctic (SODA) research project. Funding was through grant numbers N00014-16-1-2349, N00014-14-1-2377, N00014-18-1-2687, and N00014-16-1-2381.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 35(4), (2021): e2020GB006895, https://doi.org/10.1029/2020GB006895.

Description: The Amazon River drains a diverse tropical landscape greater than 6 million km2, culminating in the world's largest export of freshwater and dissolved constituents to the ocean. Here, we present dissolved organic carbon (DOC), organic and inorganic nitrogen (DON, DIN), orthophosphate (PO43−), and major and trace ion concentrations and fluxes from the Amazon River using 26 samples collected over three annual hydrographs. Concentrations and fluxes were predominantly controlled by the annual wet season flood pulse. Average DOC, DON, DIN, and PO43− fluxes (±1 s.d.) were 25.5 (±1.0), 1.14 (±0.05), 0.82 (±0.03), and 0.063 (±0.003) Tg yr−1, respectively. Chromophoric dissolved organic matter absorption (at 350 nm) was strongly correlated with DOC concentrations, resulting in a flux of 74.8 × 106 m−2 yr−1. DOC and DON concentrations positively correlated with discharge while nitrate + nitrite concentrations negatively correlated, suggesting mobilization and dilution responses, respectively. Ammonium, PO43−, and silica concentrations displayed chemostatic responses to discharge. Major and trace ion concentrations displayed clockwise hysteresis (except for chloride, sodium, and rubidium) and exhibited either dilution or chemostatic responses. The sources of weathered cations also displayed seasonality, with the highest proportion of carbonate- and silicate-derived cations occurring during peak and baseflow, respectively. Finally, our seasonally resolved weathering model resulted in an average CO2 consumption yield of (3.55 ± 0.11) × 105 mol CO2 km−2 yr−1. These results represent an updated and temporally refined quantification of dissolved fluxes that highlight the strong seasonality of export from the world's largest river and set a robust baseline against which to gauge future change.

Description: This work was supported by a grant from the Harbourton Foundation to R. G. M. Spencer and R. M. Holmes. T. W. Drake was supported by ETH Zurich core funding to J. Six. R. G. M. Spencer was additionally supported by NSF OCE-1333157.

Description: 2021-09-15

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Paleoceanography and Paleoclimatology 36(5), (2021): e2021PA004219, https://doi.org/10.1029/2021PA004219.

Description: The Yucatán Peninsula (YP) has a complex hydroclimate with many proposed drivers of interannual and longer-term variability, ranging from coupled ocean–atmosphere processes to frequency of tropical cyclones. The mid-Holocene, a time of higher Northern Hemisphere summer insolation, provides an opportunity to test the relationship between YP precipitation and ocean temperature. Here, we present a new, ∼annually resolved speleothem record of stable isotope (δ18O and δ13C) and trace element (Mg/Ca and Sr/Ca) ratios for a section of the mid-Holocene (5.2–5.7 kyr BP), before extensive agriculture began in the region. A meter-long stalagmite from Río Secreto, a cave system in Playa del Carmen, Mexico, was dated using U–Th geochronology and layer counting, yielding multidecadal age uncertainty (median 2SD of ±70 years). New proxy data were compared to an existing late Holocene stalagmite record from the same cave system, allowing us to examine changes in hydrology over time and to paleoclimate records from the southern YP. The δ18O, δ13C, and Mg/Ca data consistently indicate higher mean precipitation and lower precipitation variability during the mid-Holocene compared to the late Holocene. Despite this reduced variability, multidecadal precipitation variations were persistent in regional hydroclimate during the mid-Holocene. We therefore conclude that higher summer insolation led to increased mean precipitation and decreased precipitation variability in the northern YP but that the region is susceptible to dry periods across climate mean states. Given projected decreases in wet season precipitation in the YP’s near future, we suggest that climate mitigation strategies emphasize drought preparation.

Description: This work was funded by US National Science Foundation grants AGS-1702848 (M. Medina-Elizalde) and AGS-1502877 (S. Burns). This material is based on work supported by the National Science Foundation Graduate Research Fellowship under grant 1122374 (G. Serrato Marks). Additional support was provided by the MIT EAPS Student Research Fund and the WHOI Ocean Ventures Fund.

Description: 2021-11-06

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Hayes, C. T., Costa, K. M., Anderson, R. F., Calvo, E., Chase, Z., Demina, L. L., Dutay, J., German, C. R., Heimburger-Boavida, L., Jaccard, S. L., Jacobel, A., Kohfeld, K. E., Kravchishina, M. D., Lippold, J., Mekik, F., Missiaen, L., Pavia, F. J., Paytan, A., Pedrosa-Pamies, R., Petrova, M., V., Rahman, S., Robinson, L. F., Roy-Barman, M., Sanchez-Vidal, A., Shiller, A., Tagliabue, A., Tessin, A. C., van Hulten, M., & Zhang, J. Global ocean sediment composition and burial flux in the deep sea. Global Biogeochemical Cycles, 35(4), (2021): e2020GB006769, https://doi.org/10.1029/2020GB006769.

Description: Quantitative knowledge about the burial of sedimentary components at the seafloor has wide-ranging implications in ocean science, from global climate to continental weathering. The use of 230Th-normalized fluxes reduces uncertainties that many prior studies faced by accounting for the effects of sediment redistribution by bottom currents and minimizing the impact of age model uncertainty. Here we employ a recently compiled global data set of 230Th-normalized fluxes with an updated database of seafloor surface sediment composition to derive atlases of the deep-sea burial flux of calcium carbonate, biogenic opal, total organic carbon (TOC), nonbiogenic material, iron, mercury, and excess barium (Baxs). The spatial patterns of major component burial are mainly consistent with prior work, but the new quantitative estimates allow evaluations of deep-sea budgets. Our integrated deep-sea burial fluxes are 136 Tg C/yr CaCO3, 153 Tg Si/yr opal, 20Tg C/yr TOC, 220 Mg Hg/yr, and 2.6 Tg Baxs/yr. This opal flux is roughly a factor of 2 increase over previous estimates, with important implications for the global Si cycle. Sedimentary Fe fluxes reflect a mixture of sources including lithogenic material, hydrothermal inputs and authigenic phases. The fluxes of some commonly used paleo-productivity proxies (TOC, biogenic opal, and Baxs) are not well-correlated geographically with satellite-based productivity estimates. Our new compilation of sedimentary fluxes provides detailed regional and global information, which will help refine the understanding of sediment preservation.

Description: This study was supported by the Past Global Changes (PAGES) project, which in turn received support from the Swiss Academy of Sciences and the US-NSF. The work grew out of a 2018 workshop in Aix-Marseille, France, funded by PAGES, GEOTRACES, SCOR, US-NSF, Aix Marseille Université, and John Cantle Scientific, and the authors would like to acknowledge all attendees of this meeting. The authors acknowledge the participants of the 68th cruise of RV Akademik Mstislav Keldysh for helping acquire samples. Christopher T. Hayes acknowledges support from US-NSF awards 1658445 and 1737023. Some data compilation on Arctic shelf seas was supported by the Russian Science Foundation, grant number 20-17-00157. This work was also supported through project CRESCENDO (grant no. 641816, European Commission). Zanna Chase acknowledges support from the Australian Research Council’s Discovery Projects funding scheme (project DP180102357). Christopher R. German acknowledges US-NSF awards 1235248 and 1234827. Some colorbars used in the figures were designed by Kristen Thyng et al. (2016) and Patrick Rafter.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(10), (2021): e2021GL092904, https://doi.org/10.1029/2021GL092904.

Description: We report marine dissolved organic carbon (DOC) concentrations, and DOC Δ14C and δ13C values in seawater collected from the Southern Ocean and eastern Pacific GOSHIP cruise P18 in 2016/2017. The aging of 14C in DOC in circumpolar deep water northward from 69°S to 20°N was similar to that measured in dissolved inorganic carbon in the same samples, indicating that the transport of deep waters northward is the primary control of 14C in DIC and DOC. Low DOC ∆14C and δ13C measurements between 1,200 and 3,400 m depth may be evidence of a source of DOC produced in nearby hydrothermal ridge systems (East Pacific Rise).

Description: This work was supported by NSF (OCE-1458941 and OCE-1951073 to Ellen R. M. Druffel), Fred Kavli Foundation, Keck Carbon Cycle AMS Laboratory, NSF/NOAA funded GO-SHIP Program, Canada Research Chairs program (to Brett D. Walker) and American Chemical Society Petroleum Research Fund New Directions (55,430-ND2 to Ellen R. M. Druffel and Brett D. Walker).

Description: 2021-11-24

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Chen, Z., Kwon, Y.-O., Chen, K., Fratantoni, P., Gawarkiewicz, G., Joyce, T. M., Miller, T. J., Nye, J. A., Saba, V. S., & Stock, B. C. Seasonal prediction of bottom temperature on the Northeast U.S. Continental Shelf. Journal of Geophysical Research: Oceans, 126(5), (2021): e2021JC017187, https://doi.org/10.1029/2021JC017187.

Description: The Northeast U.S. shelf (NES) is an oceanographically dynamic marine ecosystem and supports some of the most valuable demersal fisheries in the world. A reliable prediction of NES environmental variables, particularly ocean bottom temperature, could lead to a significant improvement in demersal fisheries management. However, the current generation of climate model-based seasonal-to-interannual predictions exhibits limited prediction skill in this continental shelf environment. Here, we have developed a hierarchy of statistical seasonal predictions for NES bottom temperatures using an eddy-resolving ocean reanalysis data set. A simple, damped local persistence prediction model produces significant skill for lead times up to ∼5 months in the Mid-Atlantic Bight and up to ∼10 months in the Gulf of Maine, although the prediction skill varies notably by season. Considering temperature from a nearby or upstream (i.e., more poleward) region as an additional predictor generally improves prediction skill, presumably as a result of advective processes. Large-scale atmospheric and oceanic indices, such as Gulf Stream path indices (GSIs) and the North Atlantic Oscillation Index, are also tested as predictors for NES bottom temperatures. Only the GSI constructed from temperature observed at 200 m depth significantly improves the prediction skill relative to local persistence. However, the prediction skill from this GSI is not larger than that gained using models incorporating nearby or upstream shelf/slope temperatures. Based on these results, a simplified statistical model has been developed, which can be tailored to fisheries management for the NES.

Description: This work was supported by NOAA's Climate Program Office's Modeling, Analysis, Predictions, and Projections (MAPP) Program (NA17OAR4310111, NA19OAR4320074), and Climate Program Office's Climate Variability and Predictability (CVP) Program (NA20OAR4310482). We acknowledge our participation in MAPP's Marine Prediction Task Force.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(5),(2021): e2020JB021098, https://doi.org/10.1029/2020JB021098.

Description: We report results from 149 heat flux measurements made over an ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 103 m2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m−2, with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m−2. We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s−1, a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.

Description: This work was funded by U.S. National Science Foundation (NSF) grants EAR-1515283 to R. N. Harris and J. E. Favorito, EAR-1516361 to R. A. Sohn, and EAR-1514865 to K. M. Luttrell All work in Yellowstone National Park was completed under research permit (YELL-2018-SCI-7018) and the authors thank Annie Carlson from the Yellowstone Center for Resources for logistical help.

Description: 2021-11-14

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Global Biogeochemical Cycles 35(5), (2021): e2020GB006706, https://doi.org/10.1029/2020GB006706.

Description: The Southern Ocean plays a critical role in regulating global uptake of atmospheric CO2. Trace elements like iron (Fe), cobalt (Co), and manganese (Mn) have been shown to modulate this primary productivity. Despite limited data, the vertical profiles for Mn, Fe, and Co in the Ross Sea show no evidence of scavenging, as typically observed in oceanic sites. This was previously attributed to low-particle abundance and/or by mixing rates exceeding scavenging rates. Scavenging of some trace metals such as cobalt (Co) is thought to be largely governed by Mn (oxyhydr)oxides, assumed to be the main component of particulate Mn (pMn). However, our data show that pMn has an average oxidation state below 3 and with nondetectable Mn oxides. In addition, soluble Co profiles show no evidence of scavenging and Co uptake measurements show little Co uptake in the euphotic zone and low/no scavenging at depth. Instead, high concentrations of dissolved Mn (dMn, up to 90 nM), which is primarily complexed as Mn(III)-L (up to 100%), are observed. Average dMn concentrations (10 ± 14 nM) are highest in bottom and surface waters. Manganese sources may include sediments and sea-ice melt, as elevated dMn was measured in sea ice (12 nM) compared to its surrounding waters (3 nM), and sea ice dMn was 97% Mn(III)-L. We contend that the lack of Co scavenging in the Ross Sea is due to a unique Mn redox cycle that favors the stabilization of Mn(III)-complexes at the expense of Mn oxide particle formation.

Description: The authors acknowledge support from the NSF 1643684 (MS), NSF 1644073 (GRD), NSF OCE-1355720 (CMH), and the Woods Hole Oceanographic Institution Post-Doctoral Scholarship (VEO). The Stanford Synchrotron Radiation Lightsource was utilized in this study. Use of the Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515.

Description: 2021-10-30

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Whitehead, J. A. Fluid flow with three upstream configurations in freezing tubes. Journal of Geophysical Research: Earth Surface, 126(6), (2021): e2020JF005969, https://doi.org/10.1029/2020JF005969.

Description: The accumulation of frozen liquid around a central passageway of melt as it flows through a freezing region can make calculations very challenging. To both illustrate and to quantify some of these challenges from freezing, a model equation is developed. It simplifies the solution of Holmes (2007, https://gfd.whoi.edu/wp-content/uploads/sites/18/2018/03/MHolmesGFDReport_30151.pdf) for low Reynolds number single component liquid flow through a long tube that has a wall kept at subfreezing temperature. This model equation is used in conjunction with three different upstream configurations, each with parameters expressing their behavior. Analytical and numerical results give the parameters that have criteria for: the freezing of a compressible upstream reservoir that includes oscillatory behavior; the freezing of flow fed through a constriction with a large upstream pressure, just like a dripping water faucet during winter; the evolution of flow in multiple tubes connected by an upstream manifold, where some tubes end up with full flow and others freeze shut. Numerical runs with 1,000 tubes give a formula for the spacing between actively flowing (non-frozen) tubes over wide ranges of the two upstream parameters (flow rate and manifold resistance). Results have implications in various areas in earth science. Some are: oscillatory and freezing shut criteria for flow of magma from a compressible region, a criterion for wintertime ice accumulation at natural springs, and the spacing between volcanos.

Description: Emeritus funds are provided by Woods Hole Oceanographic Institution.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Janout, M. A., Hellmer, H. H., Hattermann, T., Huhn, O., Sueltenfuss, J., Osterhus, S., Stulic, L., Ryan, S., Schroeder, M., & Kanzow, T. FRIS revisited in 2018: on the circulation and water masses at the Filchner and Ronne Ice Shelves in the Southern Weddell Sea. Journal of Geophysical Research: Oceans, 126(6), (2021): e2021JC017269, https://doi.org/10.1029/2021JC017269.

Description: The Filchner-Ronne Ice Shelf (FRIS) is characterized by moderate basal melt rates due to the near-freezing waters that dominate the wide southern Weddell Sea continental shelf. We revisited the region in austral summer 2018 with detailed hydrographic and noble gas surveys along FRIS. The FRIS front was characterized by High Salinity Shelf Water (HSSW) in Ronne Depression, Ice Shelf Water (ISW) on its eastern flank, and an inflow of modified Warm Deep Water (mWDW) entering through Central Trough. Filchner Trough was dominated by Ronne HSSW-sourced ISW, likely forced by a recently intensified circulation beneath FRIS due to enhanced sea ice production in the Ronne polynya since 2015. Glacial meltwater fractions and tracer-based water mass dating indicate two separate ISW outflow cores, one hugging the Berkner slope after a two-year travel time, and the other located in the central Filchner Trough following a ∼six year-long transit through the FRIS cavity. Historical measurements indicate the presence of two distinct modes, in which water masses in Filchner Trough were dominated by either Ronne HSSW-derived ISW (Ronne-mode) or more locally derived Berkner-HSSW (Berkner-mode). While the dominance of these modes has alternated on interannual time scales, ocean densities in Filchner Trough have remained remarkably stable since the first surveys in 1980. Indeed, geostrophic velocities indicated outflowing ISW-cores along the trough's western flank and onto Berkner Bank, which suggests that Ronne-ISW preconditions Berkner-HSSW production. The negligible density difference between Berkner- and Ronne-mode waters indicates that each contributes cold dense shelf waters to protect FRIS against inflowing mWDW.

Description: This study used samples and data provided by the Alfred Wegener Institute Helmholtz-Center for Polar- and Marine Research in Bremerhaven (Grant No. AWI-PS111_01). The authors thank Captain Schwarze and the crew of RV Polarstern for a very successful expedition. We acknowledge support from the EU Horizon 2020 grants 820575 (HHH, SØ) and 821001 (TK, SØ).

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(4), (2021): e2020JB021434, https://doi.org/10.1029/2020JB021434.

Description: Recent multi-channel seismic studies of fast spreading and hot-spot influenced mid-ocean ridges reveal magma bodies located beneath the mid-crustal Axial Magma Lens (AML), embedded within the underlying crustal mush zone. We here present new seismic images from the Juan de Fuca Ridge that show reflections interpreted to be from vertically stacked magma lenses in a number of locations beneath this intermediate-spreading ridge. The brightest reflections are beneath Northern Symmetric segment, from ∼46°42′-52′N and Split Seamount, where a small magma body at local Moho depths is also detected, inferred to be a source reservoir for the stacked magma lenses in the crust above. The imaged magma bodies are sub-horizontal, extend continuously for along-axis lengths of ∼1–8 km, with the shallowest located at depths of ∼100–1,200 m below the AML, and are similar to sub-AML bodies found at the East Pacific Rise. At both ridges, stacked sill-like lenses are detected beneath only a small fraction of the ridge length examined and are inferred to mark local sites of higher melt flux and active replenishment from depth. The imaged magma lenses are focused in the upper part of the lower crust, which coincides with the most melt rich part of the crystal mush zone detected in other geophysical studies and where sub-vertical fabrics are observed in geologic exposures of oceanic crust. We infer that the multi-level magma accumulations are ephemeral and may result from porous flow and mush compaction, and that they can be tapped and drained during dike intrusion and eruption events.

Description: This research was supported by NSF OCE 0002488 and 0648303 (LDEO), 0002551 (WHOI), 1658199 and 1357076 (UTIG). S. M. Carbotte was partially supported by Columbia University and J. P. Canales by the Independent Research & Development Program at WHOI.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(8), (2021): e2021JC017524, https://doi.org/10.1029/2021JC017524.

Description: Since 2014, an array of current meters deployed in the Iceland Basin as part of the Overturning in the Subpolar North Atlantic Program has provided new measurements of the southward flow of Iceland-Scotland Overflow Water (ISOW) along the eastern flank of the Reykjanes Ridge. The location of the array, near 58–59°N, captures the ISOW plume at the farthest downstream location in the Iceland Basin before significant amounts of ISOW can flow into the Irminger Basin through deep fractures in the Reykjanes Ridge. The net transport of the ISOW plume at this location—approximately 5.3 Sv based on the first 4 years of observations—is significantly larger than previous values obtained farther north in the Iceland Basin, suggesting that either previous measurements did not fully capture the plume transport or that additional entrainment into the ISOW plume occurs as it approaches the southern tip of the Reykjanes Ridge. A detailed water mass analysis of the plume from continuous temperature/salinity observations shows that about 50% of the plume transport (2.6 Sv) is derived from dense waters flowing over the Nordic Sea sills into the Iceland Basin, while the remainder is made up of nearly equal parts of entrained Atlantic thermocline water and modified Labrador Sea Water. The overall results from this study suggest that the ISOW plume approximately doubles its transport through entrainment, similar to that of the Denmark Strait overflow plume in the Irminger Sea that forms the other major overflow source of North Atlantic Deep Water.

Description: Financial support for this research was provided by the U.S. National Science Foundation under grants OCE-1259398 and OCE-1756231. S. Zou acknowledges support by the U.S. National Science Foundation under grant OCE-1756361.

Description: 2022-02-13

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 126(7), (2021): e2021JC017407, https://doi.org/10.1029/2021JC017407.

Description: The Arctic Ocean receives a large supply of dissolved organic matter (DOM) from its catchment and shelf sediments, which can be traced across much of the basin's upper waters. This signature can potentially be used as a tracer. On the shelf, the combination of river discharge and sea-ice formation, modifies water densities and mixing considerably. These waters are a source of the halocline layer that covers much of the Arctic Ocean, but also contain elevated levels of DOM. Here we demonstrate how this can be used as a supplementary tracer and contribute to evaluating ocean circulation in the Arctic. A fraction of the organic compounds that DOM consists of fluoresce and can be measured using in-situ fluorometers. When deployed on autonomous platforms these provide high temporal and spatial resolution measurements over long periods. The results of an analysis of data derived from several Ice Tethered Profilers (ITPs) offer a unique spatial coverage of the distribution of DOM in the surface 800 m below Arctic sea-ice. Water mass analysis using temperature, salinity and DOM fluorescence, can clearly distinguish between the contribution of Siberian terrestrial DOM and marine DOM from the Chukchi shelf to the waters of the halocline. The findings offer a new approach to trace the distribution of Pacific waters and its export from the Arctic Ocean. Our results indicate the potential to extend the approach to separate freshwater contributions from, sea-ice melt, riverine discharge and the Pacific Ocean.

Description: Danish Strategic Research Council for the NAACOS project (grant no. 10-093903), the Danish Center for Marine Research (grant no. 2012-01). C. A. S. has received funding from the Independent Research Fund Denmark Grant No. 9040-00266B. Funding for R.M.W.A. came from the US NSF, Arctic Natural Science program grant 1504469. RG-A has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Grant Agreement No. 839311. ITP93 and part of the work by MH and BR were a contribution to the Helmholtz society strategic investment Frontiers in Arctic Marine monitoring (FRAM). The work of BR is a contribution to the cooperative projects Regional Atlantic Circulation and global Change (RACE) grant #03F0824E funded by the German Ministry of Science and Education (BBMF) and Advective Pathways of nutrients and key Ecological substances in the Arctic (APEAR) grants NE/R012865/1, NE/R012865/2 and #03V01461, part of the Changing Arctic Ocean program, jointly funded by the UKRI Natural Environment Research Council (NERC) and the BMBF. Support for Krishfield was made possible by grants from the NSF Arctic Observing Network program (PLR-1303644 and OPP-1756100).

Description: 2021-12-27

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Dragone, N. B., Diaz, M. A., Hogg, I., Lyons, W. B., Jackson, W. A., Wall, D. H., Adams, B. J., & Fierer, N. Exploring the boundaries of microbial habitability in soil. Journal of Geophysical Research: Biogeosciences, 126(6), (2021): e2020JG006052, https://doi.org/10.1029/2020JG006052.

Description: Microbes are widely assumed to be capable of colonizing even the most challenging terrestrial surface environments on Earth given enough time. We would not expect to find surface soils uninhabited by microbes as soils typically harbor diverse microbial communities and viable microbes have been detected in soils exposed to even the most inhospitable conditions. However, if uninhabited soils do exist, we might expect to find them in Antarctica. We analyzed 204 ice-free soils collected from across a remote valley in the Transantarctic Mountains (84–85°S, 174–177°W) and were able to identify a potential limit of microbial habitability. While most of the soils we tested contained diverse microbial communities, with fungi being particularly ubiquitous, microbes could not be detected in many of the driest, higher elevation soils—results that were confirmed using cultivation-dependent, cultivation-independent, and metabolic assays. While we cannot confirm that this subset of soils is completely sterile and devoid of microbial life, our results suggest that microbial life is severely restricted in the coldest, driest, and saltiest Antarctic soils. Constant exposure to these conditions for thousands of years has limited microbial communities so that their presence and activity is below detectable limits using a variety of standard methods. Such soils are unlikely to be unique to the studied region with this work supporting previous hypotheses that microbial habitability is constrained by near-continuous exposure to cold, dry, and salty conditions, establishing the environmental conditions that limit microbial life in terrestrial surface soils.

Description: This work was supported by grants from the U.S. National Science Foundation (ANT 1341629 to B. J. Adams, N. Fierer, W. Berry Lyons, and D. H. Wall and OPP 1637708 to B. J. Adams) with additional support provided to N. B. Dragone from University Colorado Department of Ecology and Evolutionary Biology.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Andres, M. Spatil and temporal variability of the Gulf Stream near Cape Hatteras. Journal of Geophysical Research: Oceans, 126(9), (2021): e2021JC017579, https://doi.org/10.1029/2021JC017579.

Description: In situ observations from a 19-month deployment of current- and pressure-sensor equipped inverted echo sounders (CPIESs) along and across the Gulf Stream near Cape Hatteras capture spatial and temporal variability where this western boundary current separates from the continental margin. Regional hydrographic casts and two temperature cross-sections spanning the Gulf Stream southeast of Cape Hatteras are used with the CPIESs' records of acoustic travel time to infer changes in thermocline depth DT and Gulf Stream position. Wave-like Gulf Stream meanders are observed where the Stream approaches the separation location with periods less than 15 days, wavelengths less than 500-km, and phase speeds between 40 and 70 km d−1. Though meander amplitudes typically decrease by ∼30% on the final approach to Cape Hatteras, some signals are still coherent across the Gulf Stream separation location. Temporal variability in meander intensity may be related to the Loop Current ∼1,400 km upstream. Mesoscale variability is strongest downstream of the separation location where Gulf Stream position is no longer constrained by the steep continental slope. Low frequency transport changes in the Florida Straits are correlated with sea-surface height gradients along the entire South Atlantic Bight (SAB) and with DT inferred at the CPIES sites. The correlations with DT are likely due to coherent transport anomalies in the Gulf Stream approaching the separation location, which then drive Gulf Stream position changes downstream of the separation location. The patterns of coherent transport anomalies may reflect large-scale atmospheric forcing patterns or rapid equatorward propagation of barotropic signals along the SAB.

Description: This research was supported by the National Science Foundation through grant OCE-1558521.

Description: Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 48(8), (2021): e2020GL089471, https://doi.org/10.1029/2020GL089471.

Description: Major gaps exist in our understanding of the pathways between internal wave generation and breaking in the Southern Ocean, with important implications for the distribution of internal wave-driven mixing, the sensitivity of ocean mixing rates and patterns to changes in the ocean environment, and the necessary ingredients of mixing parameterizations. Here we assess the dominant processes in internal wave evolution by characterizing wave and mesoscale flow scales based on full-depth in situ measurements in a Southern Ocean mixing hot spot and a ray tracing calculation. The exercise highlights the importance of Antarctic Circumpolar Current jets as a dominant influence on internal wave life cycles through advection, the modification of wave characteristics via wave-mean flow interactions, and the set-up of critical layers for both upward- and downward-propagating waves. Our findings suggest that it is important to represent mesoscale flow impacts in parameterizations of internal wave-driven mixing in the Southern Ocean.

Description: The SOFine project was funded by the UK Natural Environmental Research Council (NERC) (grant NE/G001510/1). S. Waterman is currently supported by the National Science and Engineering Research Council of Canada (NSERC) Discovery Grant Program (NSERC-2020-05799). A. Meyer acknowledges current support from the ARC Centre of Excellence for Climate Extremes (CE170100023) and previous support from the joint CSIRO-University of Tasmania Quantitative Marine Science (QMS) program. A. N. Garabato acknowledges the support of the Royal Society and the Wolfson Foundation.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Kuo, Y.-N., Lo, M.-H., Liang, Y.-C., Tseng, Y.-H., & Hsu, C.-W. Terrestrial water storage anomalies emphasize interannual variations in global mean sea level during 1997-1998 and 2015-2016 El Nino Events. Geophysical Research Letters, 48(18), (2021): e2021GL094104, https://doi.org/10.1029/2021GL094104.

Description: Interannual variations in global mean sea level (GMSL) closely correlate with the evolution of El Niño-Southern Oscillation. However, GMSL differences occur in extreme El Niños; for example, in the 2015–2016 and 1997–1998 El Niños, the peak GMSL during the mature stage of the former (9.00 mm) is almost 2.5 times higher than the latter (3.72 mm). Analyses from satellite and reanalysis data sets show that the disparity in GMSL is primarily due to barystatic (ocean mass) changes. We find that the 2015–2016 event developed not purely as an Eastern Pacific El Niño event but with Central Pacific (CP) El Niño forcing. CP El Niños contribute to a stronger negative anomaly of global terrestrial water storage and subsequent higher barystatic heights. Our results suggest that the mechanism of hydrology-related interannual variations of GMSL should be further emphasized, as more CP El Niño events are projected to occur.

Description: This study was supported by a grant of MOST 106-2111-M-002-010-MY4 to National Taiwan University.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Casini, L., Maino, M., Sanfilippo, A., Ildefonse, B., & Dick, H. J. B. High-temperature strain localization and the nucleation of oceanic core complexes (16.5 degrees N, Mid-Atlantic Ridge). Journal of Geophysical Research: Solid Earth, 126(9), (2021): e2021JB022215, https://doi.org/10.1029/2021JB022215.

Description: Extension at slow to ultraslow midoceanic ridges is mostly accommodated by large detachment faults that expose mantle peridotite and/or lower-crustal rocks forming Oceanic Core Complexes (OCC). It is commonly accepted that OCC at slow spreading ridges form during the early stage of crystallization of the magmatic crust, when rocks are still close to their solidus temperature. This observation poses significant problems, as nucleation of detachment faults requires significant weakening, which instead is more easily obtained at low temperature. The RV Knorr cruise 210 Leg 5 on the 16.5°N OCC of the Mid-Atlantic Ridge recovered a narrow shear zone from the plutonic footwall of a mature detachment fault. Troctolites preserve a continuous transition from proto-mylonite to mylonite and ultra-mylonite equilibrated at temperature between 1100° and 900°C. EBSD analysis highlights increased phase mixing and weaker crystallographic fabrics in the ultra-mylonite with respect the mylonitic domains. While host troctolites were completely solidified at the deformation incoming, high-strain zones preserve evidences of syn-kinematic melt-related textures. Fabric patterns combined with plagioclase and olivine grain size piezometry and 1D rheological modeling indicate that the development of ultra-mylonite requires a switch from dislocation creep to melt-enhanced grain-boundary sliding. Activation of this mechanism was promoted by the occurrence of hydrous melt possibly produced by selective re-melting of plagioclase + Ti-pargasite microdomains in response to strain localization at subseismic strain rates. This study highlights the importance of hydrated magmatic phases to promote the onset of detachment faulting in OCC.

Description: L. Casini thanks to Regione Autonoma della Sardegna for partly supporting this research (RASSR14473), and Università di Sassari (FAR2019). Funding for H. Dick was provided by US National Science Foundation grant No. 1935837. Open access funding enabled and organized by Projekt DEAL.