skip to main content


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
NSF-PAR ID:
10477516
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

    Canopy‐snow unloading is the complex physical process of snow unloading from the canopy through meltwater drip, sublimation to the atmosphere, or solid snow unloading to the snowpack below. This process is difficult to parameterize due to limited observations. Time‐lapse photographs of snow in the canopy were characterized by citizen scientists to create a data set of snow interception observations at multiple locations across the western United States. This novel interception data set was used to evaluate three snow unloading parameterizations in the Structure for Unifying Multiple Modeling Alternatives (SUMMA) modular hydrologic modeling framework. SUMMA was modified to include a third snow unloading parameterization, termed Wind‐Temperature (Roesch et al., 2001,https://doi.org/10.1007/s003820100153), which includes wind‐dependent and temperature‐dependent unloading functions. It was compared to a meltwater drip unloading parameterization, termed Melt (Andreadis et al., 2009,https://doi.org/10.1029/2008wr007042), and a time‐dependent unloading parameterization, termed Exponential‐Decay (Hedstrom & Pomeroy, 1998,https://doi.org/10.1002/(SICI)1099-1085(199808/09)12:10/11<1611::AID-HYP684>3.0.CO;2-4). Wind‐Temperature performed well without calibration across sites, specifically in cold climates, where wind dominates unloading and rime accretion is low. At rime prone sites, Wind‐Temperature should be calibrated to account for longer interception events with less sensitivity to wind, otherwise Melt can be used without calibration. The absence of model physics in Exponential‐Decay requires local calibration that can only be transferred to sites with similar unloading patterns. The choice of unloading parameterization can result in 20% difference in SWE on the ground below the canopy and 10% difference in estimated average winter canopy albedo. These novel observations shed light on processes that are often overlooked in hydrology.

     
    more » « less
  2. Abstract

    The Mackenzie Mountains (MMs) in the Yukon and Northwest Territories, Canada, are an enigmatic mountain range. They are currently uplifting (Leonard et al., 2008,https//doi.org/10.1029/2007JB005456), yet are about 700 km from the nearest plate boundary. Their arcuate shape is distinct and extends over 100 km eastward from the general trend of the Northern Canadian Cordillera. To better assess the cause and conditions of the current uplift, we processed ambient seismic noise data from a linear array of broadband seismographs crossing the mountains, along with other regional seismic stations, to estimate Rayleigh wave phase velocities between 6 and 40 s periods. From this, we estimated phase velocity dispersion and performed a tomographic inversion to estimateVS. Tomography reveals a low‐velocity structure that extends upward from the base of the ∼50–66 km thick lithosphere to the upper crust, and we hypothesize that inferred low density and low rigidity associated with theVSanomaly localizes the ongoing uplift and thrust‐dominated seismicity of the MMs. Additionally, we find relatively low crustal velocities that extend to the west of the MMs, suggesting that strain transfer from the Gulf of Alaska plate boundary plays a driving role as the crust translates to the northeast and buckles up against the craton consistent with the orogenic float hypothesis of Mazzotti and Hyndman (2002,https//doi.org/10.1130/0091-7613(2002)030〈0495:YCASTA〉2.0.CO;2). Finally, we observe lithospheric azimuthal anisotropy with an NW‐SE fast direction. This is nearly orthogonal to teleseismic shear wave splitting measurements in the central MMs, and suggests that asthenosphere flow and lithospheric strain are not aligned in this region.

     
    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

    Atomic oxygen (O) in the mesosphere and lower thermosphere (MLT) results from a balance between production via photo‐dissociation in the lower thermosphere and chemical loss by recombination in the upper mesosphere. The transport of O downward from the lower thermosphere into the mesosphere is preferentially driven by the eddy diffusion process that results from dissipating gravity waves and instabilities. The motivation here is to probe the intra‐annual variability of the eddy diffusion coefficient (kzz) and eddy velocity in the MLT based on the climatology of the region, initially accomplished by Garcia and Solomon (1985,https://doi.org/10.1029/JD090iD02p03850). In the current study, the intra‐annual cycle was divided into 26 two‐week periods for each of three zones: the northern hemisphere (NH), southern hemisphere (SH), and equatorial (EQ). Both 16 years of SABER (2002–2018) and 10 years of SCIAMACHY (2002–2012) O density measurements, along with NRLMSIS®2.0 were used for calculation of atomic oxygen eddy diffusion velocities and fluxes. Our prominent findings include a dominant annual oscillation below 87 km in the NH and SH zones, with a factor of 3–4 variation between winter and summer at 83 km, and a dominant semiannual oscillation at all altitudes in the EQ zone. The measured global average kzzat 96 km lacks the intra‐annual variability of upper atmosphere density data deduced by Qian et al. (2009,https://doi.org/10.1029/2008JA013643). The very large seasonal (and hemispherical) variations in kzzand O densities are important to separate and isolate in satellite analysis and to incorporate in MLT models.

     
    more » « less
  5. Abstract

    All else equal, if the ocean's “biological [carbon] pump” strengthens, the dissolved oxygen (O2) content of the ocean interior declines. Confidence is now high that the ocean interior as a whole contained less oxygen during the ice ages. This is strong evidence that the ocean's biological pump stored more carbon in the ocean interior during the ice ages, providing the core of an explanation for the lower atmospheric carbon dioxide (CO2) concentrations of the ice ages. Vollmer et al. (2022,https://doi.org/10.1029/2021PA004339) combine proxies for the oxygen and nutrient content of bottom waters to show that the ocean nutrient reservoir was more completely harnessed by the biological pump during the Last Glacial Maximum, with an increase in the proportion of dissolved nutrients in the ocean interior that were “regenerated” (transported as sinking organic matter from the ocean surface to the interior) rather than “preformed” (transported to the interior as dissolved nutrients by ocean circulation). This points to changes in the Southern Ocean, the dominant source of preformed nutrients in the modern ocean, with an apparent additional contribution from a decline in the preformed nutrient content of North Atlantic‐formed interior water. Vollmer et al. also find a lack of LGM‐to‐Holocene difference in the preformed13C/12C ratio of dissolved inorganic carbon. This finding may allow future studies to resolve which of the proposed Southern Ocean mechanisms was most responsible for enhanced ocean CO2storage during the ice ages: (a) coupled changes in ocean circulation and biological productivity, or (b) physical limitations on air‐sea gas exchange.

     
    more » « less