skip to main content


Title: Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean–atmosphere models
Abstract

High-frequency precipitation variance is calculated in 12 different free-running (non-data-assimilative) coupled high resolution atmosphere–ocean model simulations, an assimilative coupled atmosphere–ocean weather forecast model, and an assimilative reanalysis. The results are compared with results from satellite estimates of precipitation and rain gauge observations. An analysis of irregular sub-daily fluctuations, which was applied by Covey et al. (Geophys Res Lett 45:12514–12522, 2018.https://doi.org/10.1029/2018GL078926) to satellite products and low-resolution climate models, is applied here to rain gauges and higher-resolution models. In contrast to lower-resolution climate simulations, which Covey et al. (2018) found to be lacking with respect to variance in irregular sub-daily fluctuations, the highest-resolution simulations examined here display an irregular sub-daily fluctuation variance that lies closer to that found in satellite products. Most of the simulations used here cannot be analyzed via the Covey et al. (2018) technique, because they do not output precipitation at sub-daily intervals. Thus the remainder of the paper focuses on frequency power spectral density of precipitation and on cumulative distribution functions over time scales (2–100 days) that are still relatively “high-frequency” in the context of climate modeling. Refined atmospheric or oceanic model grid spacing is generally found to increase high-frequency precipitation variance in simulations, approaching the values derived from observations. Mesoscale-eddy-rich ocean simulations significantly increase precipitation variance only when the atmosphere grid spacing is sufficiently fine (< 0.5°). Despite the improvements noted above, all of the simulations examined here suffer from the “drizzle effect”, in which precipitation is not temporally intermittent to the extent found in observations.

 
more » « less
Award ID(s):
1851164
NSF-PAR ID:
10372600
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Climate Dynamics
Volume:
59
Issue:
9-10
ISSN:
0930-7575
Page Range / eLocation ID:
p. 2887-2913
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    New estimates ofpCO2from profiling floats deployed by the Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) project have demonstrated the importance of wintertime outgassing south of the Polar Front, challenging the accepted magnitude of Southern Ocean carbon uptake (Gray et al., 2018,https://doi:10.1029/2018GL078013). Here, we put 3.5 years of SOCCOM observations into broader context with the global surface carbon dioxide database (Surface Ocean CO2Atlas, SOCAT) by using the two interpolation methods currently used to assess the ocean models in the Global Carbon Budget (Le Quéré et al., 2018,https://doi:10.5194/essd‐10‐2141‐2018) to create a ship‐only, a float‐weighted, and a combined estimate of Southern Ocean carbon fluxes (<35°S). In our ship‐only estimate, we calculate a mean uptake of −1.14 ± 0.19 Pg C/yr for 2015–2017, consistent with prior studies. The float‐weighted estimate yields a significantly lower Southern Ocean uptake of −0.35 ± 0.19 Pg C/yr. Subsampling of high‐resolution ocean biogeochemical process models indicates that some of the differences between float and ship‐only estimates of the Southern Ocean carbon flux can be explained by spatial and temporal sampling differences. The combined ship and float estimate minimizes the root‐mean‐squarepCO2difference between the mapped product and both data sets, giving a new Southern Ocean uptake of −0.75 ± 0.22 Pg C/yr, though with uncertainties that overlap the ship‐only estimate. An atmospheric inversion reveals that a shift of this magnitude in the contemporary Southern Ocean carbon flux must be compensated for by ocean or land sinks within the Southern Hemisphere.

     
    more » « less
  2. Abstract

    A recent laboratory study suggests that water vapor displays structured absorption features over the 290–350 nm region, with maximum and minimum cross‐sections of 8.4 × 10−25and 1.4 × 10−25 cm2/molecule at room temperature (Pei et al. 2019,https://doi.org/10.1029/2019jd030724; Du et al., 2013,https://doi.org/10.1002/grl.50935). To observe water vapor absorption features in the ultraviolet (UV) region in the atmosphere, a United States Department of Agriculture reference spectroradiometer was upgraded with a new fore‐optical module, enabling it to measure direct solar beam and sky radiance at given azimuth and elevation angles. This double Czerny‐Turner spectroradiometer enables wavelength scanning from 290 to 410 nm, with a nominal bandwidth of 0.1 nm. It can operate with a step‐size of 0.0005 nm and a full width at half maximum of 0.1 nm. It has an out‐of‐band rejection ratio of approximately 10−10. This high resolution spectroradiometer can be used as a reference instrument for UV radiation measurements and for monitoring atmospheric gases such as O3, SO2, and NO2. A series of field observations were made using this spectroradiometer in the University at Albany campus. A residual analysis method is developed to analyze absorption by atmospheric components and to retrieve atmospheric optical depth. The residual optical depth was calculated by subtracting the optical depths of Rayleigh scattering, aerosol extinction, and absorption of typical atmospheric gases such as O3, SO2, and NO2from the retrieved total optical depth. Multiple case studies show that residual optical depth from the observed UV spectra is sensitive to the atmospheric water vapor amount. The greater the water vapor path, the larger the magnitude of residual optical depth. The ozone amount was inferred from the residual analysis; it is comparable to the satellite measurements. For example, in a case with water vapor path of 13 mm on October 24, 2019, the inferred ozone amount from residual analysis is 2.7% lower than retrievals from the Ozone Monitoring Instrument‐Total Ozone Mapping Spectrometer.

     
    more » « less
  3. 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
  4. 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
  5. Abstract

    On the dayside of August 25–26, 2018 (main phase, MP of the storm), we unveiled the storm time effects on the latitudinal distribution of ionospheric total electron content (TEC). We used 17 and 19 Global Positioning System receivers in American and Asian‐Australian sectors, respectively. Also, we employed a pair of magnetometers in each sector to unveil storm time effects on verticalE×Bupward directed inferred drift velocity in the F region ionosphere. Also used is NASA Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite airglow instrument to investigate storm time changes in neutral composition, O/N2ratio. In this investigation, we corrected the latitudinal offset found in the works of Younas et al. (2020,https://doi.org/10.1029/2020JA027981). Interestingly, we observed that a double‐humped increase (DHI) seen at a middle latitude station (MGUE, ∼22°S) after the MP on the dayside in American sector (Younas et al., 2020,https://doi.org/10.1029/2020JA027981) did straddle ∼23.58°N and ∼22°S. On August 25, 2018, storm commencement was evident in Sym‐H (∼−8 nT) around 18:00 UT. It later became intensified (∼−174 nT) on August 26 around 08:00 UT. During storm's MP (after the MP), fountain effect operation was significantly enhanced (inhibited) in Asian‐Australian (American) sector. Middle latitude TEC during MP got reduced in American sector (13:00 LT–15:40 LT) compared to those seen in Asian‐Australian sector (13:00 LT–15:40 LT). The northern equatorial peak (∼25 TECU) seen at IHYO (14:00 LT) after MP in the American sector is higher when compared with that (∼21 TECU) seen at PPPC (11:40 LT) during MP in Asian‐Australian sector.

     
    more » « less