Explore the vast range of titles available.

Although the surface eddy kinetic energy (EKE) has been well studied using satellite altimeter and surface drifter observations, our knowledge of EKE in the ocean interior is much more limited due to the sparsity of subsurface current measurements. Here we develop a new approach for estimating EKE over the full depth of the global ocean by combining 20 years of satellite altimeter and Argo float data to infer the vertical profile of eddies. The inferred eddy profiles are surface‐intensified at low latitudes and deep‐reaching at mid‐ and high latitudes. They compare favorably to the first empirical orthogonal function obtained from current meter velocities. The global‐integrated EKE estimated from the inferred profiles is about 3.1 × 1018 J, which is close to that estimated from the surface mode (3.0 × 1018 J) but about 30% smaller than that estimated from the traditional flat bottom modes (4.6 × 1018 J). Plain Language Summary The ocean is full of mesoscale eddies, analogous to weather systems in the atmosphere. Eddy kinetic energy in the surface ocean is generally well studied thanks to the availability of satellite and drifter data. The subsurface eddy energy, on the other hand, is not well known due to the relative lack of subsurface current observations. Using vertical eddy structures inferred from satellite altimeter and Argo float data, we provide the first observational estimate of eddy kinetic energy over the full depth of the global ocean. Our results have important implications for understanding the ocean energy budget and for representing the effects of mesoscale eddies in ocean and climate models. Key Points A new method is developed for estimating full‐depth eddy kinetic energy (EKE) from satellite altimeter and Argo float data Mesoscale eddy structures are surface‐intensified at low latitudes and deep‐reaching at high latitudes The total EKE in the global ocean is estimated to be about 3.1 × 1018 J

Recently introduced in oceanography to interpret the near-surface circulation, transition path theory (TPT) is a methodology that rigorously characterizes ensembles of trajectory pieces flowing out from a source last and into a target next, i.e., those that most productively contribute to transport. Here we use TPT to frame, in a statistically more robust fashion than earlier analysis, equatorward routes of North Atlantic Deep Water (NADW) in the subpolar North Atlantic. TPT is applied on all available RAFOS and Argo floats in the area by means of a discretization of the Lagrangian dynamics described by their trajectories. By considering floats at different depths, we investigate transition paths of NADW in its upper (UNADW) and lower (LNADW) layers. We find that the majority of UNADW transition paths sourced in the Labrador and southwestern Irminger Seas reach the western side of a target arranged zonally along the southern edge of the subpolar North Atlantic domain visited by the floats. This is accomplished in the form of a well-organized deep boundary current (DBC). LNADW transition paths sourced west of the Reykjanes Ridge reveal a similar pattern, while those sourced east of the ridge are found to hit the western side of the target via a DBC and also several other places along it in a less organized fashion, indicating southward flow along the eastern and western flanks of the Mid-Atlantic Ridge. Naked-eye inspection of trajectories suggest generally more diffused equatorward NADW routes. A source-independent dynamical decomposition of the flow domain into analogous backward-time basins of attraction, beyond the reach of direct inspection of trajectories, reveals a much wider influence of the western side of the target for UNADW than for LNADW. For UNADW, the average expected duration of the pathways from the Labrador and Irminger Seas was found to be of 2–3 years. For LNADW, the duration was found to be influenced by the Reykjanes Ridge, being as long as 8 years from the western side of the ridge and of about 3 years on average from its eastern side.

Four years (September 2012 to August 2016) of simultaneous current observations across the Yucatan Channel (~21.5°N) and the Straits of Florida (~81°W) have permitted us to investigate the characteristics of the flow through the Gulf of Mexico. The average transport in both channels is 27.6 Sv (1 Sv = 10 6 m 3 s −1 ), in accordance with previous estimates. At the Straits of Florida section, the transport related to the astronomical tide explains 55% of the observed variance with a mixed semidiurnal/diurnal character, while in the Yucatan Channel tides contribute 82% of the total variance and present a dominant diurnal character. At periods longer than a week the transports in the Yucatan and Florida sections have a correlation of 0.83 without any appreciable lag. The yearly running means of the transport time series in both channels are well correlated (0.98) and present a 3-Sv range variation in the 4 years analyzed. This long-term variability is well related to the convergence of the Sverdrup transport in the North Atlantic between 14.25° and 18.75°N. Using 2 years (July 2014–July 2016) of simultaneous currents observations in the Florida section, the Florida Cable section (~26.7°N), and a section across the Old Bahama Channel (~78.4°W), a mean northward transport of 28.4, 31.1, and 1.6 Sv, respectively, is obtained, implying that only 1.1 Sv is contributed by the Northwest Providence Channel to the mean transport observed at the Cable section during this 2-yr period.

Barotropic flow along depth contours is found in accordance with standard geostrophic theory. A numerical model is developed that studies the deviation from such a flow. The model gives a good approximation of the dynamical processes on the West Spitsbergen Shelf (WSS) and shows that the West Spitsbergen Current (WSC), the main gateway of Atlantic water (AW) toward the Arctic, connects more easily to the Isfjorden Trough than anywhere else along the shelf. The circulation of AW in the troughs along the WSS is here named the Spitsbergen Trough Current (STC). From hydrographical and ocean current observations it is evident that the STC is primarily barotropic and driven by the sea surface height. A connection between the along-coast wind stress and the STC is established, and it is demonstrated how the increased occurrence of winter cyclones in Fram Strait during January–February accelerates and widens the WSC. Ultimately, this results in a strengthened STC and dominance of AW on the WSS. The STC represents a slower route of AW toward the Arctic Ocean and a large heat transport toward the West Spitsbergen fjords during winter (0.2–0.4 TW toward Isfjorden). Heat flux estimates show that half of the AW heat loss in the Isfjorden Trough is due to heat loss to the surrounding water masses, while the rest is lost to the atmosphere. Sea ice production along West Spitsbergen has been reduced, or even nonexistent, in some fjords since 2006. Here, the authors argue that this is a consequence of the strong southerly wind periods along the WSS during winter.

Understanding carbon sources and sinks across the Earth's surface is fundamental in climate science and policy; thus, these topics have been extensively studied but have yet to be fully resolved and are associated with massive debate regarding the sign and magnitude of the carbon budget from global to regional scales. Developing new models and estimates based on state-of-the-art algorithms and data constraints can provide valuable knowledge and contribute to a final ensemble model in which various optimal carbon budget estimates are integrated, such as the annual global carbon budget paper. Here, we develop a new atmospheric inversion system based on the 4D local ensemble transform Kalman filter (4D-LETKF) coupled with the GEOS-Chem global transport model to infer surface-to-atmosphere net carbon fluxes from Orbiting Carbon Observatory-2 (OCO-2) V10r XCO2 retrievals. The 4D-LETKF algorithm is adapted to an OCO-2-based global carbon inversion system for the first time in this work. On average, the mean annual terrestrial and oceanic fluxes between 2015 and 2020 are estimated as − 2.02 and − 2.34 GtC yr−1, respectively, compensating for 21 % and 24 %, respectively, of global fossil carbon dioxide (CO2) emissions (9.80 GtC yr−1). Our inversion results agree with the CO2 atmospheric growth rates reported by the National Oceanic and Atmospheric Administration (NOAA) and reduce the modeled CO2 concentration biases relative to the prior fluxes against surface and aircraft measurements. Our inversion-based carbon fluxes are broadly consistent with those provided by other global atmospheric inversion models, although discrepancies still occur in the land–ocean flux partitioning schemes and seasonal flux amplitudes over boreal and tropical regions, possibly due to the sparse observational constraints of the OCO-2 satellite and the divergent prior fluxes used in different inversion models. Four sensitivity experiments are performed herein to vary the prior fluxes and uncertainties in our inversion system, suggesting that regions that lack OCO-2 coverage are sensitive to the priors, especially over the tropics and high latitudes. In the further development of our inversion system, we will optimize the data-assimilation configuration to fully utilize current observations and increase the spatial and seasonal representativeness of the prior fluxes over regions that lack observations.

A high-resolution reanalysis of the circulation in the Kattegat and Skagerrak is used to investigate the mechanisms that control the variability in the onset of the Norwegian Coastal Current. In the reanalysis, the authors have used all available in situ and remote sensing observations of salinity and temperature and use surface current observations from two coastal high-frequency radars that were ideally placed to monitor the exchange between the two basins. This study finds a strong correlation between the variability in the wind forcing in the Skagerrak and the transport in the Norwegian Coastal Current through the Torungen–Hirtshals section. Two cases with winds into and out of the Skagerrak are studied in more detail, and the results suggest asymmetries in the forcing mechanisms. For winds out of the Skagerrak, strong outflows of Baltic Sea Water associated with a deflection of the Kattegat–Skagerrak Front may disrupt local processes in the Skagerrak, which is not accounted for in previously published conceptual models for the variability of the coastal currents in this region.

Oceanic lithosphere, which forms two‐thirds of Earth's surface, is generated at mid‐ocean ridge spreading centers. Yet the internal structure of the lithosphere is not well characterized and often considered to be homogeneous relative to the structure of continental lithosphere. While geophysical observations clearly delineate the crust‐mantle boundary and the lithosphere‐asthenopshere boundary, other seismic anomalies known as mid‐lithosphere discontinuities (MLDs) have been challenging to detect and poorly constrained. Here we present melt transport models applied to the mid‐ocean ridge system that indicate MLDs are a widespread fundamental feature of oceanic lithosphere. In our models, some melt generated from decompression melting is frozen back into the lithosphere, forming a layered refertilization pattern. These refertilized layers are due to the stacked horizontal layering pattern of melt pooling beneath the freezing front. If the recrystallized melt is incorporated into the lithosphere as mafic lenses, the predicted seismic velocity is compatible with observations. Plain Language Summary Using state of the art numerical models of magma transport within Earth's interior at mid‐ocean ridges, we found that the layers observed in the oceanic lithosphere are due to the process of melt freezing close to the ridge axis and are likely ubiquitous beneath the ocean floor. We find that the expected seismic signature would match current observations. This supports the idea of a heterogeneous oceanic lithosphere that will eventually be subducted at a convergence margin. Considering both seismic and electromagnetic observations, the models suggest that the mantle is less permeable and therefore likely stronger than previously thought. Key Points Melt transport at mid‐ocean ridge models reproduce observed frozen melt layers in the oceanic lithosphere Frozen melt layers are present regardless of model parameters Predicted seismic signatures require further mechanisms to be compatible with observations

Despite its importance in controlling the abundance of methane (CH4) and a myriad of other tropospheric species, the hydroxyl radical (OH) is poorly constrained due to its large spatial heterogeneity and the inability to measure tropospheric OH with satellites. Here, we present a methodology to infer tropospheric column OH (TCOH) in the tropics over the open oceans using a combination of a machine learning model, output from a simulation of the GEOS model, and satellite observations. Our overall goals are to assess the feasibility of our methodology, to identify potential limitations, and to suggest areas of improvement in the current observational network. The methodology reproduces the variability of TCOH from independent 3D model output and of observations from the Atmospheric Tomography mission (ATom). While the methodology also reproduces the magnitude of the 3D model validation set, the accuracy of the magnitude when applied to observations is uncertain because current observations are insufficient to fully evaluate the machine learning model. Despite large uncertainties in some of the satellite retrievals necessary to infer OH, particularly for NO2 and formaldehyde (HCHO), current satellite observations are of sufficient quality to apply the machine learning methodology, resulting in an error comparable to that of in situ OH observations. Finally, the methodology is not limited to a specific suite of satellite retrievals. Comparison of TCOH determined from two sets of retrievals does show, however, that systematic biases in NO2, resulting both from retrieval algorithm and instrumental differences, lead to relative biases in the calculated TCOH. Further evaluation of NO2 retrievals in the remote atmosphere is needed to determine their accuracy. With slight modifications, a similar methodology could likely be expanded to the extratropics and over land, with the benefits of increasing our understanding of the atmospheric oxidation capacity and, for instance, informing understanding of recent CH4 trends.

Observations from Coastal Data Information Program (CDIP) moored buoys off the coast of Florida reveal tidally driven wave–current interactions that modify significant wave heights by up to 25% and shift peak periods by up to a second. A case study at Fernandina Beach, Florida, shows surface waves steepening on following tidal currents and becoming less steep on opposing tidal currents, with the largest modulations occurring in the long-period swell band. To better understand tidal modulations as a function of the phase of the tide, we use simplified analytical and numerical solutions to the equations of geometrical optics and conservation of wave action under the assumption of a one-dimensional tide acting as a progressive shallow-water wave. The theoretical frameworks allow us to identify parameters that characterize the magnitude of variation in surface waves due to tidally induced currents and changes in water depth. We compute modulations to the omnidirectional and directional wave spectrum (between 0.05 and 0.15 Hz), as well as characteristic bulk parameters such as significant wave height and peak period. The theory is corroborated using directional wave and surface current observations from the Fernandina Beach CDIP station (located in water of average depth of 16 m). We find that the numerical results reproduce the observed wave modulations due to tidal currents and changes in water depth. Specifically, surface waves traveling in the direction of the tide are strongly modulated, and the relative speeds between the tide and surface waves set the sign and magnitude of these modulations. Given knowledge of tidal currents, water-depth variations, and wave climatology, theoretical and numerical predictions may be used to provide both statistical and instantaneous estimates of wave-height variations due to tides. Because operational forecasts and nowcasts do not routinely include tides or currents, these findings can help to accurately represent nearshore surface wave variability.

The drag coefficient under tropical cyclones and its dependence on sea states are investigated by combining upper-ocean current observations [using electromagnetic autonomous profiling explorer (EM-APEX) floats deployed under five tropical cyclones] and a coupled ocean–wave (Modular Ocean Model 6–WAVEWATCH III) model. The estimated drag coefficient averaged over all storms is around 2–3 × 10 −3 for wind speeds of 25–55 m s −1 . While the drag coefficient weakly depends on wind speed in this wind speed range, it shows stronger dependence on sea states. In particular, it is significantly reduced when the misalignment angle between the dominant wave direction and the wind direction exceeds about 45°, a feature that is underestimated by current models of sea state–dependent drag coefficient. Since the misaligned swell is more common in the far front and in the left-front quadrant of the storm (in the Northern Hemisphere), the drag coefficient also tends to be lower in these areas and shows a distinct spatial distribution. Our results therefore support ongoing efforts to develop and implement sea state–dependent parameterizations of the drag coefficient in tropical cyclone conditions.