skip to main content


Title: Retrieval of Aerosol Size Distributions From In Situ Particle Counter Measurements: Instrument Counting Efficiency and Comparisons With Satellite Measurements
Abstract

The method to derive aerosol size distributions from in situ stratospheric measurements from the University of Wyoming is modified to include an explicit counting efficiency function (CEF) to describe the channel‐dependent instrument counting efficiency. This is motivated by Kovilakam and Deshler's (2015,https://doi.org/10.1002/2015JD023303) discovery of an error in the calibration method applied to the optical particle counter (OPC40) developed in the late 1980s and used from 1991 to 2012. The method can be applied to other optical aerosol instruments for which counting efficiencies have been measured. The CEF employed is the integral of the Gaussian distribution representing the instrument response at any one aerosol channel, the aerosol counting efficiency. Results using the CEF are compared to previous derivations of aerosol size distributions (Deshler et al., 2003,https://doi.org/10.1029/2002JD002514) applied to the measurements before and after Kovilakam and Deshler's correction of number concentration for the OPC40 calibration error. The CEF method is found, without any tuning parameter, to reproduce or improve upon the Kovilakam and Deshler's results, thus accounting for the calibration error without any external comparisons other than the laboratory determined counting efficiency at each aerosol channel. Moments of the new aerosol size distributions compare well with aerosol extinctions measured by Stratospheric Aerosol and Gas Experiment II and Halogen Occultation Experiment in the volcanic period 1991–1996, generally within ±40%, the precision of OPC40 moments, and in the nonvolcanic period after 1996, generally within ±20%. Stratospheric Aerosol and Gas Experiment II and Halogen Occultation Experiment estimates of aerosol surface area are generally in agreement with those derived using the new CEF method.

 
more » « less
NSF-PAR ID:
10460667
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Atmospheres
Volume:
124
Issue:
9
ISSN:
2169-897X
Page Range / eLocation ID:
p. 5058-5087
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    We present results and analysis of finite‐difference time‐domain (FDTD) simulations of electromagnetic waves scattering off meteor head plasma using an analytical model and a simulation‐derived model of the head plasma distribution. The analytical model was developed by (Dimant & Oppenheim, 2017b,https://doi.org/10.1002/2017JA023963) and the simulation‐derived model is based on particle‐in‐cell (PIC) simulations presented in (Sugar et al., 2019,https://doi.org/10.1029/2018JA026434). Both of these head plasma distribution models show the meteor head plasma is significantly different than the spherically symmetric distributions used in previous studies of meteor head plasma. We use the FDTD simulation results to fit a power law model that relates the meteoroid ablation rate to the head echo radar cross section (RCS), and show that the RCS of plasma distributions derived from the Dimant‐Oppenheim analytical model and the PIC simulations agree to within 4 dBsm. The power law model yields more accurate meteoroid mass estimates than previous methods based on spherically symmetric plasma distributions.

     
    more » « less
  4. Abstract

    Using Magnetospheric Multiscale (MMS) observations and combined MHD/test particle simulations, we further explore characteristic ion velocity distributions in the plasma sheet boundary layer. The observations are characterized by earthward beams, which at a slightly later time are accompanied by weaker but faster tailward beams. Two events are presented showing different histories. The first event happens at entry from the lobe into the plasma sheet. Energy‐time dispersion indicates a source region about 25 tailward of the satellite. The second event follows the passage of a dipolarization front closer to Earth. In contrast to earlier MHD simulations, but in better qualitative agreement with the first observation, reconnection in the present simulation was initiated near. Simulated distributions right at the boundary are characterized by a single crescent‐shaped earthward beam, as discussed earlier (Birn, Hesse, et al., 2015,https://doi.org/10.1002/2015JA021573). Farther inside, or at a later time, the distributions now also show a simple reflected beam, evolving toward a more ring‐like distribution. The simulations provide insight into the acceleration sites: The innermost edges of the direct and reflected beams consist of ions accelerated in the vicinity of the reconnection site. This supports the validity of estimating the acceleration location based on a time‐of‐flight analysis (after Onsager et al., 1990,https://doi.org/10.1029/GL017i011p01837). However, this assumption becomes invalid at later times when the acceleration becomes dominated by the earthward propagating dipolarization electric field, such that earthward and tailward reflected beams are no longer accelerated at the same location and the same time.

     
    more » « less
  5. Abstract

    We present modeling results of tube and knot (T&K) dynamics accompanying thermospheric Kelvin Helmholtz Instabilities (KHI) in an event captured by the 2018 Super Soaker campaign (R. L. Mesquita et al., 2020,https://doi.org/10.1029/2020JA027972). Chemical tracers released by a rocketsonde on 26 January 2018 showed coherent KHI in the lower thermosphere that rapidly deteriorated within 45–90 s. Using wind and temperature data from the event, we conducted high resolution direct numerical simulations (DNS) employing both wide and narrow spanwise domains to facilitate (wide domain case) and prohibit (narrow domain case) the axial deformation of KH billows that allows tubes and knots to form. KHI T&K dynamics are shown to produce accelerated instability evolution consistent with the observations, achieving peak dissipation rates nearly two times larger and 1.8 buoyancy periods faster than axially uniform KHI generated by the same initial conditions. Rapidly evolving twist waves are revealed to drive the transition to turbulence; their evolution precludes the formation of secondary convective instabilities and secondary KHI seen to dominate the turbulence evolution in artificially constrained laboratory and simulation environments. T&K dynamics extract more kinetic energy from the background environment and yield greater irreversible energy exchange and entropy production, yet they do so with weaker mixing efficiency due to greater energy dissipation. The results suggest that enhanced mixing from thermospheric KHI T&K events could account for the discrepancy between modeled and observed mixing in the lower thermosphere (Garcia et al., 2014,https://doi.org/10.1002/2013JD021208; Liu, 2021,https://doi.org/10.1029/2020GL091474) and merits further study.

     
    more » « less