For nearly a century, the functional response curves, which describe how predation rates vary with prey density, have been a mainstay of ecological modelling. While originally derived to describe terrestrial interactions, they have been adopted to characterize aquatic systems in marine biogeochemical, size-spectrum, and population models. However, marine ecological modellers disagree over the qualitative shape of the curve (e.g. Type II vs. III), whether its parameters should be mechanistically or empirically defined (e.g. disk vs. Michaelis-Menten scheme), and the most representative value of those parameters. As a case study, we focus on marine biogeochemical models, providing a comprehensive theoretical, empirical, and numerical road-map for interpreting, formulating, and parameterizing the functional response when used to prescribe zooplankton specific grazing rates on a single prey source. After providing a detailed derivation of each of the canonical functional response types explicitly for aquatic systems, we review the literature describing their parameterization. Empirical estimates of each parameter vary by over three orders of magnitude across 10 orders of magnitude in zooplankton size. However, the strength and direction of the allometric relationship between each parameter and size differs depending on the range of sizes being considered. In models, which must represent the mean state of different functional groups, size spectra or in many cases the entire ocean’s zooplankton population, the range of parameter values is smaller, but still varies by two to three orders of magnitude. Next, we conduct a suite of 0-D NPZ simulations to isolate the sensitivity of phytoplankton population size and stability to the grazing formulation. We find that the disk parameterizations scheme is much less sensitive to it parameterization than the Michaelis-Menten scheme, and quantify the range of parameters over which the Type II response, long known to have destabilizing properties, introduces dynamic instabilities. Finally, we use a simple theoretical model to show how the mean apparent functional response, averaged across sufficient sub-grid scale heterogeneity diverges from the local response. Collectively, we recommend using a type II disk response for models with smaller scales and finer resolutions but suggest that a type III Michaelis-Menten response may do a better job of capturing the complexity of all processes being averaged across in larger scale and coarser resolution modal, not just local consumption and capture rates. While we focus specifically on the grazing formulation in marine biogeochemical models, we believe these recommendations are robust across a much broader range of ecosystem models.

High-resolution underwater glider data collected in the Gulf of Oman (2015-16), combined with reanalysis datasets, describe the spatial and temporal variability of the mixed layer during winter and spring. We assess the effect of surface forcing and submesoscale processes on upper ocean buoyancy and their effects on mixed layer stratification. Episodic strong and dry wind events from the northwest (Shamals) drive rapid latent heat loss events which lead to intraseasonal deepening of the mixed layer. Comparatively, the prevailing southeasterly winds in the region are more humid, and do not lead to significant heat loss, thereby reducing intraseasonal upper ocean variability in stratification. We use this unique dataset to investigate the presence and strength of submesoscale flows, particularly in winter, during deep mixed layers. These submesoscale instabilities act mainly to restratify the upper ocean during winter through mixed layer eddies. The timing of the spring restratification differs by three weeks between 2015 and 2016 and matches the sign change of the net heat flux entering the ocean and the presence of restratifying submesoscale fluxes. These findings describe key high temporal and spatial resolution drivers of upper ocean variability, with downstream effects on phytoplankton bloom dynamics and ventilation of the oxygen minimum zone.

Idealized models can reveal insights into Earth’s climate system by reducing its complexities. However, their potential is undermined by the scarcity of fully coupled idealized models with components comparable to contemporary, comprehensive Earth System Models. To fill this gap, we compare and contrast the climates of two idealized planets which build on the Simpler Models initiative of the Community Earth System Model (CESM). Using the fully coupled CESM, the Aqua configuration is ocean-covered except for two polar land caps, and the Ridge configuration has an additional pole-to-pole grid-cell-wide continent. Contrary to most sea surface temperature profiles assumed for atmosphere-only aquaplanet experiments with the thermal maximum on the equator, the coupled Aqua configuration is characterized by a global cold belt of wind-driven equatorial upwelling, analogous to the eastern Pacific cold tongue. The presence of the meridional boundary on Ridge introduces zonal asymmetry in thermal and circulation features, similar to the contrast between western and eastern Pacific. This zonal asymmetry leads to a distinct climate state from Aqua, cooled by ~2{degree sign}C via the radiative feedback of clouds and water vapor. The meridional boundary of Ridge is also crucial for producing a more Earth-like climate state compared to Aqua, including features of atmospheric and ocean circulation, the seasonal cycle of the Intertropical Convergence Zone, and the meridional heat transport. The mean climates of these two basic configurations provide a baseline for exploring other idealized ocean geometries, and their application for investigating various features and scale interactions in the coupled climate system.

Melt ponds have a strong impact on the Arctic surface energy balance and the ice-associated ecosystem because they transmit more solar radiation compared to bare ice. In the existing literature, melt ponds are considered as bright windows to the ocean, even during freeze-up in autumn. In the central Arctic during the summer-autumn transition in 2018, we encountered a situation where more snow accumulated on refrozen melt ponds compared to the adjacent bare ice, leading to a reduction in light transmittance of the ponds even below that of bare ice. Supporting results from a radiative transfer model suggest that melt ponds with a snow cover >0.04 m lead to lower light transmittance than adjacent bare ice. This scenario has not been described in the literature before, but has potentially strong implications for example on autumn ecosystem activity, oceanic heat budget and thermodynamic ice growth.

Alkalinity, the excess of proton acceptors over donors, plays a major role in ocean chemistry, in buffering and in calcium carbonate precipitation and dissolution. Understanding alkalinity dynamics is pivotal to quantify ocean carbon dioxide uptake during times of global change. Here we review ocean alkalinity and its role in ocean buffering as well as the biogeochemical processes governing alkalinity and pH in the ocean. We show that it is important to distinguish between measurable titration alkalinity and charge-balance alkalinity that is used to quantify calcification and carbonate dissolution and needed to understand the impact of biogeochemical processes on components of the carbon dioxide system. A general treatment of ocean buffering and quantification via sensitivity factors is presented and used to link existing buffer and sensitivity factors. The impact of individual biogeochemical processes on ocean alkalinity and pH is discussed and quantified using these sensitivity factors. Processes governing ocean alkalinity on longer time scales such as carbonate compensation, (reversed) silicate weathering and anaerobic mineralization are discussed and used to derive a close-to-balance ocean alkalinity budget for the modern ocean.

Observations of the spatio-temporal structure of turbulent mixing in a large, energetic strait were used to examine the interactions between wind- and tidally-forced boundary layers in a coastal environment. Te Moana-o-Raukawa (Cook Strait) of Aotearoa (New Zealand) is a relatively wide, energetic strait, known to experience substantial tidal currents and wind stress. A turbulence-enabled ocean glider mission measured O(40,000) turbulence samples that passed QAQC including the use of a vehicle-mounted speed through water sensor. The observations were compared to one-dimensional models of turbulence to understand the mechanisms that regulates the vertical structure of mixing. Tidal flows of O(1 m/s) and wind speeds of O(10 m/s) enhance dissipation to ε=O(10^{-5} W/kg) through boundary drag, shear-driven production of turbulent kinetic energy (P) and to a minor extent buoyancy flux (G). The benthic and wind-driven boundary layers behaved reasonably predictably when considering a 1D perspective. The interaction between the two boundary layers depended on mid-water column stratification which is to a large degree an externally-prescribed condition. Transient stratification can stabilize the mean flow (median Ri_g=0.6(>1/4)) and reduce both turbulence intensity (Re_b) and diapycnal diffusivity (K_z) by up to two orders of magnitude in the middle of the water column, insulating bottom and surface mixing-layers. Mid-water dissipation rate levels tend to be associated with marginal dynamical stability (median Ri_g=0.22(~1/4)) and canonical mixing efficiency (median R_f=0.17), while elevated levels are connected to unstable mean flow conditions (median Ri_g=0.14(<1/4)) and reduced mixing efficiency (median R_f=0.1(<0.17)) that promotes turbulence growth.

We investigate if the commonly neglected riverine detrital carbonate fluxes might balance several chemical mass balances of the global ocean. Particulate inorganic carbon (PIC) concentrations in riverine suspended sediments, i.e., carbon contained by these detrital carbonate minerals, was quantified at the basin and global scale. Our approach is based on globally representative datasets of riverine suspended sediment composition, catchment properties and a two-step regression procedure. The present day global riverine PIC flux is estimated at 3.1 ± 0.3 Tmol C/y (13% of total inorganic carbon export and 4 % of total carbon export), with a flux-weighted mean concentration of 0.26 ± 0.03 wt%. The flux prior to damming was 4.1 ± 0.5 Tmol C/y. PIC fluxes are concentrated in limestone-rich, rather dry and mountainous catchments of large rivers in Arabia, South East Asia and Europe with 2.2 Tmol C/y (67.6 %) discharged between 15 °N and 45 °N. Greenlandic and Antarctic meltwater discharge and ice-rafting additionally contribute 0.8 ± 0.3 Tmol C/y. This amount of detrital carbonate minerals annually discharged into the ocean implies a significant contribution of calcium (~ 4.75 Tmol Ca/y) and alkalinity fluxes (~ 10 Tmol(eq)/y) to marine mass balances and moderate inputs of strontium (~ 5 Gmol Sr/y), based on undisturbed riverine and cryospheric inputs and a dolomite/calcite ratio of 0.1. Magnesium fluxes (~ 0.25 Tmol Mg/y), mostly hosted by less-soluble dolomite, are rather negligible. These unaccounted fluxes help elucidating respective marine mass balances and potentially alter conclusions based on these budgets.

Studies on the molecular mechanisms of microbial adaptation in chaotropic and low water activity (aw) environments are poorly understood. Chaotropic environments are characterized as salt rich, MgCl2 and CaCl2, which lowers the availability of water for biological processes. PATRIC, an integrated genomic browsing tool containing vast libraries of sequenced genomes, can help us identify unique genetic markers in chaophilic and xerophilic microbes. Halophilic microbes are characterized as obligate hypersaline with the ability to tolerate exposure to chaotropic agents. Microbes with the greatest tolerance in these extreme environments must have advanced adaptive methods. Halobacterium salinarum and Haloquadratum walsbyi are chaotolerant and well adapted to low water activity. Haloquadratum walsbyi is unique among the halophilics as having the highest tolerance for chaotropes and its square shape. Performing comparative genomics using fully sequenced halophilic archaea such as Halobacterium salinarum NRC-1, a model halophile, and Haloquadratum walsbyi C23, we were able to identify genes that confer adaptation to chaotropic and low aw environments, as well as individual adaptations that may be responsible for the varying levels of tolerance in chaotropic environments . Characterizing genes associated with chaotolerance and low aw adaptations can help elucidate the cellular functions that make these microbes unique. Chaotropic brines may be used as analogs to study the origin of life and the possibility of suitable environments hosting extremophilic microbes on other planets like the Martian brines and the icy moons of Europa; therefore, studying the microbiome of chaotropic environments are essential in the field of astrobiology.

Poleward ocean heat transport is a key process in the earth system. We detail and review the northward Atlantic Water (AW) flow, Arctic Ocean heat transport, and heat loss to the atmosphere since 1900 in relation to sea ice cover. Our synthesis is largely based on a sea ice-ocean model forced by a reanalysis atmosphere (1900-2018) corroborated by a comprehensive hydrographic database (1950-), AW inflow observations (1996-), and other long-term time series of sea ice extent (1900-), glacier retreat (1984-) and Barents Sea hydrography (1900-). The Arctic Ocean, including the Nordic and Barents Seas, has warmed since the 1970s. This warming is congruent with increased ocean heat transport and sea ice loss and has contributed to the retreat of marine-terminating glaciers on Greenland. Heat loss to the atmosphere is largest in the Nordic Seas (60% of total) with large variability linked to the frequency of Cold Air Outbreaks and cyclones in the region, but there is no long-term statistically significant trend. Heat loss from the Barents Sea (~30%) and Arctic seas farther north (~10%) is overall smaller, but exhibit large positive trends. The AW inflow, total heat loss to the atmosphere, and dense outflow have all increased since 1900. These are consistently related through theoretical scaling, but the AW inflow increase is also wind-driven. The Arctic Ocean CO2 uptake has increased by ~30% over the last century - consistent with Arctic sea ice loss allowing stronger air-sea interaction and is ~8% of the global uptake.

Large Antarctic ice volume changes characterized the middle to Late Oligocene and the first million years of climate evolution during the Miocene. However, the sea surface temperature (SST) evolution over this period remains poorly constrained, as only a few records from contrasting proxies are available. In this study, we present a long-term alkenone-derived SST record from sediments drilled by the Ocean Drilling Program (ODP) at Site 1168 in the west Tasmanian Sea spanning 29.8 Ma to 16.7 Ma. The SST record highlight that the long-term warming in the Late Oligocene linked to the end of the Middle Oligocene Glacial Interval can be recognized also at mid-to-high latitudes of the Southern Hemisphere. Warmer average temperatures (25.5°C) characterize the period from 24.6 to 22 Ma; average temperatures then decrease by 1 to 2°C into the Miocene and stabilize by 20.1 Ma. The reconstructed temperatures are highly variable in the warm Late Oligocene waters, and more stable and slightly colder in the Early to Middle Miocene. We confirm that this temperature trend is not an artefact of the latitudinal drift of the site, as the temperature anomaly relative to the modern water temperature at the paleolocation confirms the SST trends of the Oligocene. This is the first alkenone-derived record to reproduce the long-term Oligocene climate trend previously interpreted from the benthic δ18O, which recorded a warming and/or reduction in ice volume from the Middle Oligocene Glacial Interval through the latest Oligocene.

Marine sedimentary records are a key archive when reconstructing past climate; however, mixing at the seabed (bioturbation) can strongly influence climate records, especially when sedimentation rates are low. By commingling the climate signal from different time periods, bioturbation both smooths climate records, by damping fast climate variations, and creates noise when measurements are made on samples containing small numbers of individual proxy carriers, such as foraminifera. Bioturbation also influences radiocarbon-based age-depth models, as sample ages may not represent the true ages of the sediment layers from which they were picked. While these effects were first described several decades ago, the advent of ultra-small-sample 14C dating now allows samples containing very small numbers of foraminifera to be measured, thus enabling us to directly measure the age-heterogeneity of sediment for the first time. Here, we use radiocarbon dates measured on replicated samples of 3-30 foraminifera to estimate age-heterogeneity for five marine sediment cores with sedimentation rates ranging from 2-30 cm / kyr. From their age-heterogeneities and sedimentation rates we infer mixing depths of 10-20 cm for our core sites. Our results show that when accounting for age-heterogeneity, the true error of radiocarbon dating can be several times larger than the reported measurement. We present estimates of this uncertainty as a function of sedimentation rate and the number of individuals per radiocarbon date. A better understanding of this uncertainty will help us to optimise radiocarbon measurements, construct age models with appropriate uncertainties and better interpret marine paleo records.

Ocean bathymetry exerts a strong control on ice sheet-ocean interactions within Antarctic ice-shelf cavities, where it can limit the access of warm, dense water at depth to the underside of floating ice shelves. However, ocean bathymetry is challenging to measure within or close to ice-shelf cavities. It remains unclear how uncertainty in existing bathymetry datasets affect simulated sub-ice shelf melt rates. Here we infer linear sensitivities of ice shelf melt rates to bathymetric shape with grid-scale detail by means of the adjoint of an ocean general circulation model. Both idealised and realistic-geometry experiments of sub-ice shelf cavities in West Antarctica reveal that bathymetry has a strong impact on melt in localised regions such as topographic obstacles to flow. Moreover, response of melt to bathymetric perturbation is found to be non-monotonic, with deepening leading to either increased or decreased melt depending on location. Our computational approach provides a comprehensive way of identifying regions where refined knowledge of bathymetry is most impactful, and also where bathymetric errors have relatively little effect on modelled ice sheet-ocean interactions.

Diurnal tidal currents are the dominate contributors to diapycnal mixing in many regions along the pathways for warm Atlantic Water (AW) circulating within the Arctic Ocean along the continental slope. This mixing diffuses AW heat and salt into the cooler and fresher surroundings, including the upper ocean where ocean heat fluxes play a role in the stability of the ice pack. The strongest diurnal currents are associated with topographically-trapped vorticity waves, which are sensitive to stratification and mean flow. In models, these waves are also sensitive to choices for forcing and geometry. Sensitivity to background conditions implies that tidal currents and mixing will change as the Arctic evolves towards a new climate state. Here, as a first step towards understanding how diurnal tidal currents might change in a future Arctic Ocean, we describe results from a suite of high-resolution (dx=2 km) 2-D and 3-D models for Arctic diurnal tides, focusing on their currents at locations along the AW pathways. We first demonstrate that accurate representation of barotropic diurnal tides requires forcing with both open boundary conditions and the direct potential tide. Next, we use 3-D models with realistic, ocean background stratification and mean flow to describe the annual cycle of depth-averaged diurnal tidal currents. Finally, we investigate the baroclinic structure of diurnally forced waves including the generation of harmonics (semidiurnal and higher) that can contribute to mixing within the water column. Our results show that tides should be explicitly included in ocean and coupled predictive models for the Arctic to represent the feedbacks between tidal energetics and ocean mean state via mixing.

Harmful algal blooms (HABs) caused by the dinoflagellate Karenia brevis on the West Florida Shelf have become a nearly annual occurrence causing widespread ecological and economic harm. Effects range from minor respiratory irritation and localized fish kills to large-scale and long-term events causing massive mortalities to marine organisms. Reports of hypoxia on the shelf have been infrequent; however, there have been some indications that some HABs have been associated with localized hypoxia. We examined oceanographic data from 2004 to 2019 across the West Florida Shelf to determine the frequency of hypoxia and to assess its association with known HABs. Hypoxia was present in 5 of the 16 years examined and was always found shoreward of the 50-meter bathymetry line. There were 2 clusters of recurrent hypoxia: midshelf off the Big Bend coast and near the southwest Florida coast. We identified 3 hypoxic events that were characterized by multiple conductivity, temperature, and depth (CTD) casts and occurred concurrently with extreme HABs in 2005, 2014, and 2018. These HAB-hypoxia events occurred when K. brevis blooms initiated in early summer months and persisted into the fall likely driven by increased biological oxygen demand from decaying algal biomass and reduced water column ventilation due to stratification. There were also four years, 2011, 2013, 2015, and 2017, with low dissolved oxygen located near the shelf break that were likely associated with upwelling of deeper Gulf of Mexico water onto the shelf. We had difficulty in assessing the spatiotemporal extent of these events due to limited data availability and potentially unobserved hypoxia due to the inconsistent difference between the bottom of the CTD cast and the seafloor. While we cannot unequivocally explain the association between extreme HABs and hypoxia on the West Florida Shelf, there is sufficient evidence to suggest a causal linkage between them.

Barrier islands are especially vulnerable to hurricanes and other large storms, owing to their mobile composition, low elevations, and detachment from the mainland. Conceptual models of barrier-island evolution emphasize ocean-side processes that drive landward migration through overwash, inlet migration, and aeolian transport. In contrast, we found that the impact of Hurricane Dorian (2019) on North Core Banks, a 36-km barrier island on the Outer Banks of North Carolina, was primarily driven by inundation of the island from Pamlico Sound, as evidenced by storm-surge model results and observations of high-water marks and wrack lines. Analysis of photogrammetry products from aerial imagery collected before and after the storm indicate the loss of about 18% of the subaerial volume of the island through the formation of over 80 erosional washout channels extending from the marsh and washover platform, through gaps in the foredunes, to the shoreline. The washout channels were largely co-located with washover fans deposited by earlier events. Net seaward export of sediment resulted in the formation of deltaic bars offshore of the channels, which became part of the post-storm berm recovery by onshore bar migration and partial filling of the washouts with washover deposits within two months. The partially filled features have created new ponds and lowland habitats that will likely persist for years. We conclude that this event represents a setback in the overwash/rollover behavior required for barrier transgression.

Understanding the droplet cloud and spray dynamics is important on the study of the ocean surface and marine boundary layer.The role that the wave energy and the type of wave breaking plays in the resulting distribution and dynamics of droplets is yet to be understood. The aim of this work was to generate violent plunging breakers in the laboratory and analyze the spray production by the crest of the wave when it impacts in the free surface. The droplet sizes and their dynamics were measured and the effect of different wind speeds on the droplet production was also considered. It was found that the mean radius increases with the wave energy and the presence of larger droplets (radius > 2 mm) in the vertical direction increases with the presence of wind. Furthermore, the normalized distribution of droplet sizes is consistent with the distribution of ligament-mediated spray formation. Also, indications of turbulence affecting the droplet dynamics at wind speeds of 5 m/s were found. The amount of large droplets (radius >1 mm) found in this work was larger than reported in the literature. An improved estimation of the initial distribution of large droplets can largely affect the evolution of the Sea Spray Generation Function, and therefore the estimation of energy and mass transport in the marine boundary layer.

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 set the scales of the subducted water masses in the pycnocline. As a result, the tracer-derived vertical transport rate spectrum is flatter than the vertical velocity spectrum. 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 pathways that are unraveled in this study along with the Lagrangian evolution of properties on water parcels, emphasize the role of submesoscale dynamics coupled with mesoscale frontogenesis.

The Southern California Bight (SCB), an eastern boundary upwelling system, is impacted by global warming, acidification and deoxygetation, and receives anthropogenic nutrients from a coastal population of 20 million people. We describe the configuration, forcing, and validation of a realistic, submesoscale resolving ocean model as a tool to investigate coastal eutrophication. This modeling system represents an important achievement because it strikes a balance of capturing the forcing by U.S. Pacific Coast-wide phenomena, while representing the bathymetric features and submesoscale circulation that affect the vertical and horizontal transport of nutrients from natural and human sources. Moreover, the model allows to run simulations at timescales that approach the interannual frequencies of ocean variability, making the grand challenge of disentangling natural variability, climate change, and local anthropogenic forcing a tractable task in the near-term. The model simulation is evaluated against a broad suite of observational data throughout the SCB, showing realistic depiction of mean state and its variability with remote sensing and in situ physical-biogeochemical measurements of state variables and biogeochemical rates. The simulation reproduces the main structure of the seasonal upwelling front, the mean current patterns, the dispersion of plumes, as well as their seasonal variability. It reproduces the mean distributions of key biogeochemical and ecosystem properties. Biogeochemical rates reproduced by the model, such as primary productivity and nitrification, are also consistent with measured rates. Results of this validation exercise demonstrate the utility of fine-scale resolution modeling in support of management decisions on local anthropogenic nutrient discharges to coastal zones.

Organic matter (OM) sulfurization can enhance the preservation and sequestration of carbon in anoxic sediments, and it has been observed in sinking marine particles from marine O2-deficient zones. The magnitude of this effect on carbon burial remains unclear, however, because the transformations that occur when sinking particles encounter sulfidic conditions remain undescribed. Here, we briefly expose sinking marine particles from the eastern tropical North Pacific O2-deficient zone to environmentally relevant sulfidic conditions (20C, 0.5 mM [poly]sulfide, two days) and then characterize the resulting solid-phase organic and inorganic products in detail. During these experiments, the abundance of organic sulfur in both hydrolyzable and hydrolysis-resistant solids roughly triples, indicating extensive OM sulfurization. Lipids also sulfurize on this timescale, albeit less extensively. In all three pools, OM sulfurization produces organic monosulfides, thiols, and disulfides. Hydrolyzable sulfurization products appear within ≤ 200-m regions of relatively homogenous composition that are suggestive of sulfurized extracellular polymeric substances (EPS). Concurrently, reactions with particulate iron oxyhydroxides generate low and fairly uniform concentrations of iron sulfide (FeS) within these same EPS-like materials. Iron oxyhydroxides were not fully consumed during the experiment, which demonstrates that organic materials can be competitive with reactive iron for sulfide. These experiments support the hypothesis that sinking, OM- and EPS-rich particles in a sulfidic water mass can sulfurize within days, potentially contributing to enhanced sedimentary carbon sequestration. Additionally, sulfur-isotope and chemical records of organic S and iron sulfides in sediments have the potential to incorporate signals from water column processes.

Prediction of rapid intensification in tropical cyclones prior to landfall is a major societal issue. While air-sea interactions are clearly linked to storm intensity, the connections between the underlying thermal conditions over continental shelves and rapid intensification are limited. Here, an exceptional set of in-situ and satellite data are used to identify spatial heterogeneity in sea surface temperatures across the inner core of Hurricane Sally (2020), a storm that rapidly intensified over the shelf. A leftward shift in the region of maximum cooling was observed as the hurricane transited from the open gulf to the shelf. This shift was generated, in part, by the surface heat flux in conjunction the along and across-shelf transport of heat from storm-generated coastal circulation. The spatial differences in the sea surface temperatures were large enough to potentially influence rapid intensification processes suggesting that coastal thermal features need to be accounted for to improve storm forecasting as well as to better understand how climate change will modify interactions between tropical cyclones and the coastal ocean.

The ocean’s “biological pump” significantly modulates atmospheric carbon dioxide levels. However, the complexity and variability of processes involved introduces uncertainty in interpretation of transient observations and future climate projections. Much work has focused on “parametric uncertainty”, particularly determining the exponent(s) of a power-law relationship of sinking particle flux with depth. Varying this relationship’s functional form introduces additional “structural uncertainty”. We use an ocean biogeochemistry model substituting six alternative remineralization profiles fit to a reference power-law curve, to characterize structural uncertainty, which, in atmospheric pCO2 terms, is roughly 50% of the parametric uncertainty associated with varying the power-law exponent within its plausible global range, and similar to uncertainty associated with regional variation in power-law exponents. The substantial contribution of structural uncertainty to total uncertainty highlights the need to improve characterization of biological pump processes, and compare the performance of different profiles within Earth System Models to obtain better constrained climate projections.

Melt ponds play an important role in the seasonal evolution of Arctic sea ice. During the melt season, snow atop the sea ice begins to metamorphose and melt, forming ponds on the ice. These ponds reduce the albedo of the surface, allowing for increased solar energy absorption and thus further melting of snow and ice. Analyzing the spatial distribution and temporal evolution of melt ponds helps us understand the sea ice processes that occur during the summer melt season. It has been shown that the inclusion of melt pond parameters in sea ice models increases the skill of predicting the summer sea ice minimum extent. Previous studies have used remote sensing imagery to characterize surface features and calculate melt pond statistics. Here we use new observations of melt ponds obtained by the Digital Mapping System (DMS) flown onboard NASA Operation IceBridge (OIB) during two Arctic summer melt campaigns which surveyed thousands of kilometers of sea ice and resulted in more than 45,000 images. One campaign was conducted in the Beaufort Sea (July 2016), and one in the Lincoln Sea and the Arctic Ocean north of Greenland (July 2017). Using these data we expect to advance our understanding of the differences and similarities between melt pond features on young, thin sea ice seen in the Beaufort Sea versus those on multi-year ice. We have developed a pixel-based classification scheme by considering the different RGB spectral values associated with each surface type. We identify four sea ice surface types (level ice, rubbled ice, open water, and melt ponds). The classification scheme enables the calculation of parameters including melt pond fraction, ice concentration, melt pond area, and melt pond dimensions. We compare results with data from the Airborne Topographic Mapper (ATM), a laser altimeter also operated during these OIB missions. Given the extent over which the OIB data are available, regional information may be derived. Leveraging existing satellite data products, we examine whether the high-resolution airborne statistics are representative of the region and can be scaled up for comparison against satellite-derived parameters such as ice concentration and extent.

We decompose the monthly global Ocean Bottom Pressure (OBP) from GRACE(-FO) mass concentration solutions, with trends and seasonal harmonics removed from the signal, to extract 23 significant regional modes of variability. The 23 modes are analyzed and discussed considering Sea-Level Anomalies (SLA), Wind Stress Curl (WSC), and major climate indices. Two-thirds of the patterns correspond to extratropical regions and are substantially documented in other global or regional studies. Over the equatorial band, the identified modes are unprecedented, with an amplitude ranging between 0.5 and 1 centimeter. With smaller amplitude than extratropical patterns, they appear to be less correlated with the local SLA or WSC; yet, they present significantly coherent dynamics. The Pacific Ocean modes show significant correlations with the Pacific Decadal Oscillation (PDO) and El Niño Southern Oscillation (ENSO).

Plate motion is a remarkable Earth process and is widely ascribed to two primary driving forces: slab pull and ridge push. With the release of the first- and second-order stress fields since 1989, a few features of tectonic stresses provide strong constrain on these forces. The observed stresses are mainly distributed on the uppermost brittle part of the lithosphere. A modeling analysis, however, reveals that the stress produced by ridge push is dominantly distributed in the lower part of the lithosphere; Doglioni and Panza recently made an in-depth investigation on slab pull and found this force cannot be in accordance with observations. These findings of ridge push and slab pull suggest that there needs other force to be responsible for plate motion and tectonic stress. Here, we propose that the pressure of deep ocean water against the wall of continent yields enormous force (i.e., ocean-generated force) on the continent. The continent is fixed on the top of the lithosphere, this attachment allows ocean-generated force to be laterally transferred to the lithospheric plate. We show that this force may combine other forces to form force balances for the lithospheric plate, consequently, the African, Indian, South American, Australian, and Pacific plates obtain a movement of 4.52, 6.09, 2.11, 3.52, and 6.62 cm/yr, respectively. A torque balance modelling shows that the error between the movements calculated for 121 sample locations and the movements extracted from GSRM v.2.1 is less than 0.8 mm/yr in speed and 0.3o in azimuth.