skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Sea State Dependent Gas Transfer Velocity for CO 2 Unifying Theory, Model, and Field Data
Abstract Wave breaking induced bubbles contribute a significant part of air‐sea gas fluxes. Recent modeling of the sea state dependent CO2flux found that bubbles contribute up to ∼40% of the total CO2air‐sea fluxes (Reichl & Deike, 2020,https://doi.org/10.1029/2020gl087267). In this study, we implement the sea state dependent bubble gas transfer formulation of Deike and Melville (2018,https://doi.org/10.1029/2018gl078758) into a spectral wave model (WAVEWATCH III) incorporating the spectral modeling of the wave breaking distribution from Romero (2019,https://doi.org/10.1029/2019gl083408). We evaluate the accuracy of the sea state dependent gas transfer parameterization against available measurements of CO2gas transfer velocity from 9 data sets (11 research cruises, see Yang et al. (2022,https://doi.org/10.3389/fmars.2022.826421)). The sea state dependent parameterization for CO2gas transfer velocity is consistent with observations, while the traditional wind‐only parameterization used in most global models slightly underestimates the observations of gas transfer velocity. We produce a climatology of the sea state dependent gas transfer velocity using reanalysis wind and wave data spanning 1980–2017. The climatology shows that the enhanced gas transfer velocity occurs frequently in regions with developed sea states (with strong wave breaking and high significant wave height). The present study provides a general sea state dependent parameterization for gas transfer, which can be implemented in global coupled models.  more » « less
Award ID(s):
2121646 2122042
PAR ID:
10477517
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Earth and Space Science
Volume:
10
Issue:
11
ISSN:
2333-5084
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract This letter compares the predictions of two expressions proposed for the porosity evolution in the context of rate and state friction. One (Segall & Rice, 1995,https://doi.org/10.1029/95jb02403) depends only on the sliding velocity; the other (Sleep, 1995,https://doi.org/10.1029/94jb03340) depends only on the state variable. Simulations of both are similar for velocity stepping and slide‐hold‐slide experiments. They differ significantly for normal effective stress jumps at constant sliding velocity. Segall and Rice (1995,https://doi.org/10.1029/95jb02403) predicts no change in the porosity; Sleep (1995,https://doi.org/10.1029/94jb03340) does. Simulation with a spring‐block model indicates that the magnitude of rapid slip events is essentially the same for the two formulations. Variations of porosity and induced pore pressure near rapid slip events are similar and consistent with experimental observations. Predicted porosity variations during slow slip intervals and the time at which rapid slip events occur are significantly different. The simulation indicates that changes in friction stress due to pore pressure changes exceed those due to rate and state effects. 
    more » « less
  2. Abstract We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML‐AIM). ML‐AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020,https://doi.org/10.1029/2020JA027908), the FAC‐derived auroral conductance model of Robinson et al. (2020,https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993,https://doi.org/10.1029/92gl02109). The ML‐AIM inputs are 60‐min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML‐AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML‐AIM produces physically accurate ionospheric potential patterns such as the two‐cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML‐AIM, the Weimer (2005,https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data‐assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML‐AIM than others. ML‐AIM is unique and innovative because it predicts ionospheric responses to the time‐varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005,https://doi.org/10.1029/2004ja010884) designed to provide a quasi‐static ionospheric condition under quasi‐steady solar wind/IMF conditions. Plans are underway to improve ML‐AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times. 
    more » « less
  3. Abstract We analyze three substorms that occur on (1) 9 March 2008 05:14 UT, (2) 26 February 2008 04:05 UT, and (3) 26 February 2008 04:55 UT. Using ACE solar wind velocity and interplanetary magnetic fieldBzvalues, we calculate the rectified (southwardBz) solar wind voltage propagated to the magnetosphere. The solar wind conditions for the two events were vastly different, 300 kV for 9 March 2008 substorm, compared to 50 kV for 26 February 2008. The voltage is input to a nonlinear physics‐based model of the magnetosphere called WINDMI. The output is the westward auroral electrojet current which is proportional to the auroral electrojet (AL) index from World Data Center for Geomagnetism Kyoto and the SuperMAG auroral electrojet index (SML). Substorm onset times are obtained from the superMAG substorm database, Pu et al. (2010,https://doi.org/10.1029/2009JA014217), Lui (2011,https://doi.org/10.1029/2010JA016078) and synchronized to Time History of Events and Macroscale Interactions during Substorms satellite data. The timing of onset, model parameters, and intermediate state space variables are analyzed. The model onsets occurred about 5 to 10 min earlier than the reported onsets. Onsets occurred when the geotail current in the WINDMI model reached a critical threshold of 6.2 MA for the 9 March 2008 event, while, in contrast, a critical threshold of 2.1 MA was obtained for the two 26 February 2008 events. The model estimates 1.99 PJ of total energy transfer during the 9 March 2008 event, with 0.95 PJ deposited in the ionosphere. The smaller events on 26 February 2008 resulted in a total energy transfer of 0.37 PJ according to the model, with 0.095 PJ deposited in the ionosphere. 
    more » « less
  4. Abstract We compared the performance of DREAM3D simulations in reproducing the long‐term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low‐energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics insideL* = 6, with a better model performance at lowerμthan higherμ, whereμis the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data‐driven outer boundary condition is critical to the model performance, while adding in the data‐driven seed population doesn't further improve the performance. (b) The model shows comparable performance withDLLfrom Brautigam and Albert (2000,https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014,https://doi.org/10.1002/2013ja019204), and Liu et al. (2016,https://doi.org/10.1002/2015gl067398), while withDLLfrom Ali et al. (2016,https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data‐based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higherμusing the Wang et al. (2019,https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under‐reproduced in the current model. 
    more » « less
  5. Abstract Long‐term efforts have sought to extend global model resolution to smaller scales enabling more accurate descriptions of gravity wave (GW) sources and responses, given their major roles in coupling and variability throughout the atmosphere. Such studies reveal significant improvements accompanying increasing resolution, but no guidance on what is sufficient to approximate reality. We take the opposite approach, using a finite‐volume model solving the Navier‐Stokes equations exactly. The reference simulation addresses mountain wave (MW) generation and responses over the Southern Andes described using isotropic 500 m, central resolution by Fritts et al. (2021),https://doi.org/10.1175/JAS-D-20-0207.1and Lund et al. (2020),https://doi.org/10.1175/JAS-D-19-0356.1. Reductions of horizontal resolution to 1 and 2 km result in (a) systematic increases in initial MW breaking altitudes, (b) weaker, larger‐scale generation of secondary GWs and acoustic waves accompanying these dynamics, and (c) significantly weaker and less extended responses in the mesosphere in latitude and longitude. Horizontal resolution of 4 km largely suppresses instabilities, but allows weak, sustained mean‐flow interactions. Responses for 8 km resolution are very weak and fail to capture any aspects of the high‐resolution responses. The chosen mean winds allow efficient MW penetration into the mesosphere and lower thermosphere, hence only exhibit strong pseudo‐momentum deposition and mean wind decelerations at higher altitudes. A companion paper by Fritts et al. (2022),https://doi.org/10.1029/2021JD036035explores the impacts of decreasing resolution on responses in the thermosphere. 
    more » « less