Fundamental process understanding and description of heat, mass, and momentum exchanges across the land‐atmosphere interface in model boundary forcing parameterizations is critical to the simulation of near‐surface soil moisture dynamics (e.g., bare‐soil evaporation). This study explores the sensitivity of a continuum‐scale porous media heat and mass transfer model to the spatial‐discretization length‐scales (i.e., spatial‐resolution) of near‐surface atmospheric data; the goal is to determine how much data are needed to force the model and adequately capture evaporative water losses and subsurface state variable distributions. The requisite atmospheric forcing data were taken from the high‐resolution, precision bare‐soil evaporation experiments of Trautz et al. (2018,
This content will become publicly available on August 16, 2024
Given the key role that atmospheric heat transport plays in Earth's climate system, efforts to document its changes over the satellite era are valuable. Clark et al. (2022,
- Award ID(s):
- 2019647
- NSF-PAR ID:
- 10469812
- Publisher / Repository:
- Geophysical Research Letters
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 50
- Issue:
- 15
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract https://doi.org/10.1029/2018WR023102 ). Simulation results demonstrated that shallow subsurface mass and heat transfer dynamics can be adequately captured with forcing data averaged over large length‐scales, or a minimal number of measurements, provided that soil conditions are properly described. The soil moisture spatial distributions were found to be insensitive to horizontal variations in the forcing data. The model failed to capture small‐scale trends observed experimentally; this did not impact the accuracy of total evaporative water loss estimates however. These results indicate that in future physical experimental efforts conducted at 1–10‐m length‐scales, there is no need to focus on the generation of high‐spatial resolution atmospheric measurements—time and effort would be better spent in characterizing soil conditions and properties. Even though a theoretical foundation was not provided to directly extrapolate this work to the field scale, these findings have practical value in designing field data collection strategies. -
more » « less
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
Abstract Following the reanalysis of individual experimental runs of some widely cited studies (Jain et al., 2018,
https://doi.org/10.1002/2017JB014847 ), we revisit the global data analysis of Korenaga and Karato (2008,https://doi.org/10.1029/2007JB005100 ) with a significantly improved version of their Markov chain Monte Carlo inversion. Their algorithm, previously corrected by Mullet et al. () to minimize potential parameter bias, is further modified here to estimate more efficiently interrun biases in global data sets. Using the refined Markov chain Monte Carlo inversion technique, we simultaneously analyze experimental data on the deformation of olivine aggregates compiled from different studies. Realistic composite rheological models, including both diffusion and dislocation creep, are adopted, and the role of dislocation‐accommodated grain boundary sliding is also investigated. Furthermore, the influence of interrun biases on inversion results is studied using experimental and synthetic data. Our analysis shows that existing data can tightly constrain the grain‐size exponent for diffusion creep at ∼2, which is different from the value commonly assumed ( p = 3). Different data sets and model assumptions, however, yield nonoverlapping estimates on other flow‐law parameters, and the flow‐law parameters for grain boundary sliding are poorly resolved in most cases. We thus provide a few plausible candidate flow‐law models for olivine rheology to facilitate future geodynamic modeling. The availability of more data that explore a wider range of experimental conditions, especially higher pressures, is essential to improve our understanding of upper mantle rheology. -
more » « less
Abstract An extraordinary strong wind‐driven upwelling event occurred in Lake Superior in summer of 2010 when the lake was strongly stratified. In this paper, a detailed three‐dimensional (3‐D) investigation of the current and thermal structures during the upwelling event was conducted using in situ observations, remote sensing products, and the results of a long‐term numerical simulation. A 3‐D finite‐volume coupled ice–ocean model tailored for the Laurentian Great Lakes was employed for this purpose. The model was validated with temperature observations at National Oceanic and Atmospheric Administration buoys and mooring data from 2010. The upwelling event observed in satellite imagery and at a mooring station was reproduced by the model, showing a cooling of as much as 10°C in August 2010 along the northwestern coast. The relationship between the alongshore wind and the offshore thermocline displacement (upwelling front) derived in theoretical work (Csanady, 1977,
https://doi.org/10.1029/JC082i003p00397 ) was used to calculate upwelling front movement offshore and found to be in in close agreement with model prediction. A significant correlation between alongshore wind stress and lake temperature change in the upwelling zone was found with a correlation coefficient of −0.87. A simple linear heat balance model explained most of variability in temperature. -
more » « less
Abstract We examine the behavior of natural basaltic and trachytic samples during paleointensity experiments on both the original and laboratory‐acquired thermal remanences and characterize the samples using proxies for domain state including curvature (
k ) and the bulk domain stability parameters of Paterson (2011,https://doi.org/10.1029/2011JB008369 ) and Paterson et al. (2017,https://doi.org/10.1073/pnas.1714047114 ), respectively. A curvature value of 0.164 (suggested by Paterson, 2011,https://doi.org/10.1029/2011JB008369 ) as a critical threshold that separates single‐domain‐like remanences from multidomain‐like remanances on the original paleointensity data was used to separate samples into “straight” (single‐domain‐like) and “curved” (multidomain‐like) groups. Specimens from the two sample sets were given a “fresh” thermal remanent magnetization in a 70 μT field and subjected to an infield‐zerofield, zerofield‐infield (IZZI)‐type (Yu et al., 2004,https://doi.org/10.1029/2003GC000630 ) paleointensity experiment. The straight sample set recovered the laboratory field with high precision while the curved set had much more scattered results (70.5 ± 1.5 and 71.9 ± 5.2 μT, respectively). The average intensity of both sets for straight and curved was quite close to the laboratory field of 70 μT, however, suggesting that if experiments contain a sufficient number of specimens, there does not seem to be a large bias in the field estimate. We found that the dependence of the laboratory thermal remanent magnetization on cooling rate was significant in most samples and did not depend on domain states inferred from proxies based on hysteresis measurements and should be estimated for all samples whose cooling rates differ from that used in the laboratory.
