Atlantic meridional overturning circulation increases flood risk along the United States southeast coast.

The system of oceanic flows constituting the Atlantic Meridional Overturning Circulation (AMOC) moves heat and other properties to the subpolar North Atlantic, controlling regional climate, weather, sea levels, and ecosystems. Climate models suggest a potential AMOC slowdown towards the end of this century due to anthropogenic forcing, accelerating coastal sea level rise along the western boundary and dramatically increasing flood risk. While direct observations of the AMOC are still too short to infer long-term trends, we show here that the AMOC-induced changes in gyre-scale heat content, superimposed on the global mean sea level rise, are already influencing the frequency of floods along the United States southeastern seaboard. We find that ocean heat convergence, being the primary driver for interannual sea level changes in the subtropical North Atlantic, has led to an exceptional gyre-scale warming and associated dynamic sea level rise since 2010, accounting for 30-50% of flood days in 2015-2020.

Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean-atmosphere models.

High-frequency precipitation variance is calculated in 12 different free-running (non-data-assimilative) coupled high resolution atmosphere-ocean model simulations, an assimilative coupled atmosphere-ocean weather forecast model, and an assimilative reanalysis. The results are compared with results from satellite estimates of precipitation and rain gauge observations. An analysis of irregular sub-daily fluctuations, which was applied by Covey et al. (Geophys Res Lett 45:12514-12522, 2018. 10.1029/2018GL078926) to satellite products and low-resolution climate models, is applied here to rain gauges and higher-resolution models. In contrast to lower-resolution climate simulations, which Covey et al. (2018) found to be lacking with respect to variance in irregular sub-daily fluctuations, the highest-resolution simulations examined here display an irregular sub-daily fluctuation variance that lies closer to that found in satellite products. Most of the simulations used here cannot be analyzed via the Covey et al. (2018) technique, because they do not output precipitation at sub-daily intervals. Thus the remainder of the paper focuses on frequency power spectral density of precipitation and on cumulative distribution functions over time scales (2-100 days) that are still relatively "high-frequency" in the context of climate modeling. Refined atmospheric or oceanic model grid spacing is generally found to increase high-frequency precipitation variance in simulations, approaching the values derived from observations. Mesoscale-eddy-rich ocean simulations significantly increase precipitation variance only when the atmosphere grid spacing is sufficiently fine (< 0.5°). Despite the improvements noted above, all of the simulations examined here suffer from the "drizzle effect", in which precipitation is not temporally intermittent to the extent found in observations.

Attribution of Space-Time Variability in Global-Ocean Dissolved Inorganic Carbon.

The inventory and variability of oceanic dissolved inorganic carbon (DIC) is driven by the interplay of physical, chemical, and biological processes. Quantifying the spatiotemporal variability of these drivers is crucial for a mechanistic understanding of the ocean carbon sink and its future trajectory. Here, we use the Estimating the Circulation and Climate of the Ocean-Darwin ocean biogeochemistry state estimate to generate a global-ocean, data-constrained DIC budget and investigate how spatial and seasonal-to-interannual variability in three-dimensional circulation, air-sea CO2 flux, and biological processes have modulated the ocean sink for 1995-2018. Our results demonstrate substantial compensation between budget terms, resulting in distinct upper-ocean carbon regimes. For example, boundary current regions have strong contributions from vertical diffusion while equatorial regions exhibit compensation between upwelling and biological processes. When integrated across the full ocean depth, the 24-year DIC mass increase of 64 Pg C (2.7 Pg C year-1) primarily tracks the anthropogenic CO2 growth rate, with biological processes providing a small contribution of 2% (1.4 Pg C). In the upper 100 m, which stores roughly 13% (8.1 Pg C) of the global increase, we find that circulation provides the largest DIC gain (6.3 Pg C year-1) and biological processes are the largest loss (8.6 Pg C year-1). Interannual variability is dominated by vertical advection in equatorial regions, with the 1997-1998 El Niño-Southern Oscillation causing the largest year-to-year change in upper-ocean DIC (2.1 Pg C). Our results provide a novel, data-constrained framework for an improved mechanistic understanding of natural and anthropogenic perturbations to the ocean sink.

Local and Remote Forcing of Interannual Sea-Level Variability at Nantucket Island.

The relative contributions of local and remote wind stress and air-sea buoyancy forcing to sea-level variations along the East Coast of the United States are not well quantified, hindering the understanding of sea-level predictability there. Here, we use an adjoint sensitivity analysis together with an Estimating the Circulation and Climate of the Ocean (ECCO) ocean state estimate to establish the causality of interannual variations in Nantucket dynamic sea level. Wind forcing explains 67% of the Nantucket interannual sea-level variance, while wind and buoyancy forcing together explain 97% of the variance. Wind stress contribution is near-local, primarily from the New England shelf northeast of Nantucket. We disprove a previous hypothesis about Labrador Sea wind stress being an important driver of Nantucket sea-level variations. Buoyancy forcing, as important as wind stress in some years, includes local contributions as well as remote contributions from the subpolar North Atlantic that influence Nantucket sea level a few years later. Our rigorous adjoint-based analysis corroborates previous correlation-based studies indicating that sea-level variations in the subpolar gyre and along the United States northeast coast can both be influenced by subpolar buoyancy forcing. Forward perturbation experiments further indicate remote buoyancy forcing affects Nantucket sea level mostly through slow advective processes, although coastally trapped waves can cause rapid Nantucket sea level response within a few weeks.

Ocean forcing drives glacier retreat in Greenland.

The retreat and acceleration of Greenland glaciers since the mid-1990s have been attributed to the enhanced intrusion of warm Atlantic Waters (AW) into fjords, but this assertion has not been quantitatively tested on a Greenland-wide basis or included in models. Here, we investigate how AW influenced retreat at 226 marine-terminating glaciers using ocean modeling, remote sensing, and in situ observations. We identify 74 glaciers in deep fjords with AW controlling 49% of the mass loss that retreated when warming increased undercutting by 48%. Conversely, 27 glaciers calving on shallow ridges and 24 in cold, shallow waters retreated little, contributing 15% of the loss, while 10 glaciers retreated substantially following the collapse of several ice shelves. The retreat mechanisms remain undiagnosed at 87 glaciers without ocean and bathymetry data, which controlled 19% of the loss. Ice sheet projections that exclude ocean-induced undercutting may underestimate mass loss by at least a factor of 2.

Long-Term Earth-Moon Evolution With High-Level Orbit and Ocean Tide Models.

Tides and Earth-Moon system evolution are coupled over geological time. Tidal energy dissipation on Earth slows E a r t h ' s rotation rate, increases obliquity, lunar orbit semi-major axis and eccentricity, and decreases lunar inclination. Tidal and core-mantle boundary dissipation within the Moon decrease inclination, eccentricity and semi-major axis. Here we integrate the Earth-Moon system backwards for 4.5 Ga with orbital dynamics and explicit ocean tide models that are "high-level" (i.e., not idealized). To account for uncertain plate tectonic histories, we employ Monte Carlo simulations, with tidal energy dissipation rates (normalized relative to astronomical forcing parameters) randomly selected from ocean tide simulations with modern ocean basin geometry and with 55, 116, and 252 Ma reconstructed basin paleogeometries. The normalized dissipation rates depend upon basin geometry and E a r t h ' s rotation rate. Faster Earth rotation generally yields lower normalized dissipation rates. The Monte Carlo results provide a spread of possible early values for the Earth-Moon system parameters. Of consequence for ocean circulation and climate, absolute (un-normalized) ocean tidal energy dissipation rates on the early Earth may have exceeded t o d a y ' s rate due to a closer Moon. Prior to ∼ 3 Ga , evolution of inclination and eccentricity is dominated by tidal and core-mantle boundary dissipation within the Moon, which yield high lunar orbit inclinations in the early Earth-Moon system. A drawback for our results is that the semi-major axis does not collapse to near-zero values at 4.5 Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation.

Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future.

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change.

Global Estimates of the Energy Transfer From the Wind to the Ocean, With Emphasis on Near-Inertial Oscillations.

Estimates of the kinetic energy transfer from the wind to the ocean are often limited by the spatial and temporal resolution of surface currents and surface winds. Here we examine the wind work in a pair of global, very high-resolution (1/48° and 1/24°) MIT general circulation model simulations in Latitude-Longitude-polar Cap (LLC) configuration that provide hourly output at spatial resolutions of a few kilometers and include tidal forcing. A cospectrum analysis of wind stress and ocean surface currents shows positive contribution at large scales (>300 km) and near-inertial frequency and negative contribution from mesoscales, tidal frequencies, and internal gravity waves. Larger surface kinetic energy fluxes are in the Kuroshio in winter at large scales (40 mW/m2) and mesoscales (-30 mW/m2). The Kerguelen region is dominated by large scale (∼20 mW/m2), followed by inertial oscillations in summer (13 mW/m2) and mesoscale in winter (-12 mW/m2). Kinetic energy fluxes from internal gravity waves (-0.1 to -9.9 mW/m2) are generally stronger in summer. Surface kinetic energy fluxes in the LLC simulations are 4.71 TW, which is 25-85% higher than previous global estimates from coarser (1/6-1/10°) general ocean circulation models; this is likely due to improved representation of wind variability (6-hourly, 0.14°, operational European Center for Medium-Range Weather Forecasts). However, the low wind power input to the near-inertial frequency band obtained with LLC (0.16 TW) compared to global slab models suggests that wind variability on time scales less than 6 hr and spatial scales less than 15 km are critical to better representing the wind power input in ocean circulation models.

Pathways of ocean heat towards Pine Island and Thwaites grounding lines.

In the Amundsen Sea, modified Circumpolar Deep Water (mCDW) intrudes into ice shelf cavities, causing high ice shelf melting near the ice sheet grounding lines, accelerating ice flow, and controlling the pace of future Antarctic contributions to global sea level. The pathways of mCDW towards grounding lines are crucial as they directly control the heat reaching the ice. A realistic representation of mCDW circulation, however, remains challenging due to the sparsity of in-situ observations and the difficulty of ocean models to reproduce the available observations. In this study, we use an unprecedentedly high-resolution (200 m horizontal and 10 m vertical grid spacing) ocean model that resolves shelf-sea and sub-ice-shelf environments in qualitative agreement with existing observations during austral summer conditions. We demonstrate that the waters reaching the Pine Island and Thwaites grounding lines follow specific, topographically-constrained routes, all passing through a relatively small area located around 104°W and 74.3°S. The temporal and spatial variabilities of ice shelf melt rates are dominantly controlled by the sub-ice shelf ocean current. Our findings highlight the importance of accurate and high-resolution ocean bathymetry and subglacial topography for determining mCDW pathways and ice shelf melt rates.

Global-scale dispersal and connectivity in mangroves.

Dispersal provides a key mechanism for geographical range shifts in response to changing environmental conditions. For mangroves, which are highly susceptible to climate change, the spatial scale of dispersal remains largely unknown. Here we use a high-resolution, eddy- and tide-resolving numerical ocean model to simulate mangrove propagule dispersal across the global ocean and generate connectivity matrices between mangrove habitats using a range of floating periods. We find high rates of along-coast transport and transoceanic dispersal across the Atlantic, Pacific, and Indian Oceans. No connectivity is observed between populations on either side of the American and African continents. Archipelagos, such as the Galapagos and those found in Polynesia, Micronesia, and Melanesia, act as critical stepping-stones for dispersal across the Pacific Ocean. Direct and reciprocal dispersal routes across the Indian Ocean via the South Equatorial Current and seasonally reversing monsoon currents, respectively, allow connectivity between western Indian Ocean and Indo-West Pacific sites. We demonstrate the isolation of the Hawaii Islands and help explain the presence of mangroves on the latitudinal outlier Bermuda. Finally, we find that dispersal distance and connectivity are highly sensitive to the minimum and maximum floating periods. We anticipate that our findings will guide future research agendas to quantify biophysical factors that determine mangrove dispersal and connectivity, including the influence of ocean surface water properties on metabolic processes and buoyancy behavior, which may determine the potential of viably reaching a suitable habitat. Ultimately, this will lead to a better understanding of global mangrove species distributions and their response to changing climate conditions.

Response to Comment on "Contrasting carbon cycle responses of the tropical continents to the 2015-2016 El Niño".

Chevallier showed a column CO2 ([Formula: see text]) anomaly of ±0.5 parts per million forced by a uniform net biosphere exchange (NBE) anomaly of 2.5 gigatonnes of carbon over the tropical continents within a year, so he claimed that the inferred NBE uncertainties should be larger than presented in Liu et al We show that a much concentrated NBE anomaly led to much larger [Formula: see text] perturbations.

Origin of Circumpolar Deep Water intruding onto the Amundsen and Bellingshausen Sea continental shelves.

Melting of West Antarctic ice shelves is enhanced by Circumpolar Deep Water (CDW) intruding onto the Amundsen and Bellingshausen Seas (ABS) continental shelves. Despite existing studies of cross-shelf and on-shelf CDW transports, CDW pathways onto the ABS originating from further offshore have never been investigated. Here, we investigate CDW pathways onto the ABS using a regional ocean model. Simulated CDW tracers from a zonal section across 67°S (S04P) circulate along the Antarctic Circumpolar Current (ACC) and Ross Gyre (RG) and travel into ABS continental shelf after 3-5 years, but source locations are shifted westward by ~900 km along S04P in 2001-2006 compared to 2009-2014. We find that simulated on- and off-shelf CDW is ~0.1-0.2 °C warmer in the 2009-2014 case than in the 2001-2006 case together with changes in simulated ocean circulation. These differences are primarily caused by lateral, rather than surface, boundary conditions, implying that large-scale atmospheric and ocean circulations are able to control CDW pathways and thus off- and on-shelf CDW properties.

Inferring Connectivity Range in Submerged Aquatic Populations (Ruppia L.) Along European Coastal Lagoons From Genetic Imprint and Simulated Dispersal Trajectories.

Coastal salt- and brackish water lagoons are unique shallow habitats characterized by beds of submerged seagrasses and salt-tolerant Ruppia species. Established long-term and large-scale patterns of connectivity in lagoon systems can be strongly determined by patterns of nearshore and coastal currents next to local bird-mediated seed dispersal. Despite the importance of dispersal in landscape ecology, characterizing patterns of connectivity remains challenging in aquatic systems. Here, we aimed at inferring connectivity distances of Ruppia cirrhosa along European coastal lagoons using a population genetic imprint and modeled dispersal trajectories using an eddy-resolving numerical ocean model that includes tidal forcing. We investigated 1,303 individuals of 46 populations alongside subbasins of the Mediterranean (Balearic, Tyrrhenian, Ionian) and the Atlantic to Baltic Sea coastline over maximum distances of 563-2,684 km. Ten microsatellite loci under an autotetraploid condition revealed a mixed sexual and vegetative reproduction mode. A pairwise FST permutation test of populations revealed high levels of historical connectivity only for distance classes up to 104-280 km. Since full range analysis was not fully explanatory, we assessed connectivity in more detail at coastline and subbasin level using four approaches. Firstly, a regression over restricted geographical distances (300 km) was done though remained comparable to full range analysis. Secondly, piecewise linear regression analyses yielded much better explained variance but the obtained breakpoints were shifted toward greater geographical distances due to a flat slope of regression lines that most likely reflect genetic drift. Thirdly, classification and regression tree analyses revealed threshold values of 47-179 km. Finally, simulated ocean surface dispersal trajectories for propagules with floating periods of 1-4 weeks, were congruent with inferred distances, a spatial Bayesian admixed gene pool clustering and a barrier detection method. A kinship based spatial autocorrelation showed a contemporary within-lagoon connectivity up to 20 km. Our findings indicate that strong differentiation or admixtures shaped historical connectivity and that a pre- and post LGM genetic imprint of R. cirrhosa along the European coasts was maintained from their occurrence in primary habitats. Additionally, this study demonstrates the importance of unraveling thresholds of genetic breaks in combination with ocean dispersal modeling to infer patterns of connectivity.

Scaling Properties of Arctic Sea Ice Deformation in a High-Resolution Viscous-Plastic Sea Ice Model and in Satellite Observations.

Sea ice models with the traditional viscous-plastic (VP) rheology and very small horizontal grid spacing can resolve leads and deformation rates localized along Linear Kinematic Features (LKF). In a 1 km pan-Arctic sea ice-ocean simulation, the small-scale sea ice deformations are evaluated with a scaling analysis in relation to satellite observations of the Envisat Geophysical Processor System (EGPS) in the Central Arctic. A new coupled scaling analysis for data on Eulerian grids is used to determine the spatial and temporal scaling and the coupling between temporal and spatial scales. The spatial scaling of the modeled sea ice deformation implies multifractality. It is also coupled to temporal scales and varies realistically by region and season. The agreement of the spatial scaling with satellite observations challenges previous results with VP models at coarser resolution, which did not reproduce the observed scaling. The temporal scaling analysis shows that the VP model, as configured in this 1 km simulation, does not fully resolve the intermittency of sea ice deformation that is observed in satellite data.

Ocean submesoscales as a key component of the global heat budget.

Recent studies highlight that oceanic motions associated with horizontal scales smaller than 50 km, defined here as submesoscales, lead to anomalous vertical heat fluxes from colder to warmer waters. This unique transport property is not captured in climate models that have insufficient resolution to simulate these submesoscale dynamics. Here, we use an ocean model with an unprecedented resolution that, for the first time, globally resolves submesoscale heat transport. Upper-ocean submesoscale turbulence produces a systematically-upward heat transport that is five times larger than mesoscale heat transport, with winter-time averages up to 100 W/m2 for mid-latitudes. Compared to a lower-resolution model, submesoscale heat transport warms the sea surface up to 0.3 °C and produces an upward annual-mean air-sea heat flux anomaly of 4-10 W/m2 at mid-latitudes. These results indicate that submesoscale dynamics are critical to the transport of heat between the ocean interior and the atmosphere, and are thus a key component of the Earth's climate.

Contrasting carbon cycle responses of the tropical continents to the 2015-2016 El Niño.

The 2015-2016 El Niño led to historically high temperatures and low precipitation over the tropics, while the growth rate of atmospheric carbon dioxide (CO2) was the largest on record. Here we quantify the response of tropical net biosphere exchange, gross primary production, biomass burning, and respiration to these climate anomalies by assimilating column CO2, solar-induced chlorophyll fluorescence, and carbon monoxide observations from multiple satellites. Relative to the 2011 La Niña, the pantropical biosphere released 2.5 ± 0.34 gigatons more carbon into the atmosphere in 2015, consisting of approximately even contributions from three tropical continents but dominated by diverse carbon exchange processes. The heterogeneity of the carbon-exchange processes indicated here challenges previous studies that suggested that a single dominant process determines carbon cycle interannual variability.

Contrasting Effects of Historical Sea Level Rise and Contemporary Ocean Currents on Regional Gene Flow of Rhizophora racemosa in Eastern Atlantic Mangroves.

Mangroves are seafaring taxa through their hydrochorous propagules that have the potential to disperse over long distances. Therefore, investigating their patterns of gene flow provides insights on the processes involved in the spatial genetic structuring of populations. The coastline of Cameroon has a particular geomorphological history and coastal hydrology with complex contemporary patterns of ocean currents, which we hypothesize to have effects on the spatial configuration and composition of present-day mangroves within its spans. A total of 982 trees were sampled from 33 transects (11 sites) in 4 estuaries. Using 11 polymorphic SSR markers, we investigated genetic diversity and structure of Rhizophora racemosa, a widespread species in the region. Genetic diversity was low to moderate and genetic differentiation between nearly all population pairs was significant. Bayesian clustering analysis, PCoA, estimates of contemporary migration rates and identification of barriers to gene flow were used and complemented with estimated dispersal trajectories of hourly released virtual propagules, using high-resolution surface current from a mesoscale and tide-resolving ocean simulation. These indicate that the Cameroon Volcanic Line (CVL) is not a present-day barrier to gene flow. Rather, the Inter-Bioko-Cameroon (IBC) corridor, formed due to sea level rise, allows for connectivity between two mangrove areas that were isolated during glacial times by the CVL. Genetic data and numerical ocean simulations indicated that an oceanic convergence zone near the Cameroon Estuary complex (CEC) presents a strong barrier to gene flow, resulting in genetic discontinuities between the mangrove areas on either side. This convergence did not result in higher genetic diversity at the CEC as we had hypothesized. In conclusion, the genetic structure of Rhizophora racemosa is maintained by the contrasting effects of the contemporary oceanic convergence and historical climate change-induced sea level rise.