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: Reduction of Bias from Parameter Variance in Geophysical Data Estimation: Method and Application to Ice Water Content and Sedimentation Flux Estimated from Lidar
This paper addresses issues of statistical misrepresentation of the a priori parameters (henceforth called ancillary parameters) used in geophysical data estimation. Parameterizations using ancillary data are frequently needed to derive geophysical data of interest from remote measurements. Empirical fits to the ancillary data that do not preserve the distribution of such data may induce substantial bias. A semianalytical averaging approach based on Taylor expansion is presented to improve estimated cirrus ice water content and sedimentation flux for a range of volume extinction coefficients retrieved from spaceborne lidar observations by CALIOP combined with the estimated distribution of ancillary data from in situ aircraft measurements of ice particle microphysical parameters and temperature. It is shown that, given an idealized distribution of input parameters, the approach performs well against Monte Carlo benchmark predictions. Using examples with idealized distributions at the mean temperature for the tropics at 15 km, it is estimated that the commonly neglected variance observed in in situ measurements of effective diameters may produce a worst-case estimation bias spanning up to a factor of 2. For ice sedimentation flux, a similar variance in particle size distributions and extinctions produces a worst-case estimation bias of a factor of 9. The value of the bias is found to be mostly set by the correlation coefficient between extinction and ice effective diameter, which in this test ranged between all possible values. Systematic reporting of variances and covariances in the ancillary data and between data and observed quantities would allow for more accurate observational estimates.  more » « less
Award ID(s):
1660538 1743753
PAR ID:
10134944
Author(s) / Creator(s):
 ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Volume:
77
Issue:
3
ISSN:
0022-4928
Page Range / eLocation ID:
p. 835-857
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract As part of the analysis following the Seeded and Natural Orographic Wintertime Storms (SNOWIE) project, the ice water content (IWC) in ice and mixed-phase clouds is retrieved from airborne Wyoming Cloud Radar (WCR) measurements aboard the University of Wyoming King Air (UWKA), which has a suite of integrated in situ IWC, optical array probes, and remote sensing measurements, and it provides a unique dataset for this algorithm development and evaluation. A sensitivity study with different idealized ice particle habits shows that the retrieved IWC with aggregate ice particle habit agrees the best with the in situ measurement, especially in ice or ice-dominated mixed-phase clouds with a correlation coefficient (rr) of 0.91 and a bias of close to 0. For mixed-phase clouds with ice fraction ratio less than 0.8, the variances of IWC estimates increase (rr = 0.76) and the retrieved mean IWC is larger than in situ IWC by a factor of 2. This is found to be related to the uncertainty of in situ measurements, the large cloud inhomogeneity, and the retrieval assumption uncertainty. The simulated reflectivity Ze and IWC relationships assuming three idealized ice particle habits and measured particle size distributions show that hexagonal columns with the same Ze have a lower IWC than aggregates, whose Ze–IWC relation is more consistent with the observed WCR Ze and in situ IWC relation in those clouds. The 2D stereo probe (2DS) images also indicate that ice particle habit transition occurs in orographic mixed-phase clouds; hence, the retrieved IWC assuming modified gamma particle size distribution (PSD) of aggregate particles tends to have a greater bias in this kind of clouds. 
    more » « less
  2. Abstract The intrinsic uncertainty of radar-based retrievals in snow originates from a large diversity of snow growth habits, densities, and particle size distributions, all of which can make interpreting radar measurements of snow very challenging. The application of polarimetric radar for snow measurements can mitigate some of these issues. In this study, a novel polarimetric method for quantification of the extinction coefficient and visibility in snow, based on the joint use of radar reflectivity at horizontal polarization Z and specific differential phase K DP , is introduced. A large 2D-video-disdrometer snow dataset from central Oklahoma is used to derive a polarimetric bivariate power-law relation for the extinction coefficient, . The relation is derived for particle aspect ratios ranging from 0.5 to 0.8 and the width of the canting angle distribution ranging from 0° to 40°, values typical of aggregated snow, and validated via theoretical and analytical derivations/simulations. The multiplier of the relation is sensitive to variations in particles’ densities, the width of the canting angle distribution, and particles’ aspect ratios, whereas the relation’s exponents are practically invariant to changes in the latter two parameters. This novel approach is applied to polarimetric S-band WSR-88D data and verified against previous studies and in situ measurements of the extinction coefficient for four snow events in the eastern United States. The polarimetric radar estimates of the extinction coefficient exhibit smaller biases in comparison to previous studies concerning the ground measurements. The results indicate that there is good potential for reliable radar estimates of visibility from polarimetric weather radars, a parameter inversely proportional to the extinction coefficient. 
    more » « less
  3. Abstract Properties of frozen hydrometeors in clouds remain difficult to sense remotely. Estimates of number concentration, distribution shape, ice particle density, and ice water content are essential for connecting cloud processes to surface precipitation. Progress has been made with dual-frequency radars, but validation has been difficult because of lack of particle imaging and sizing observations collocated with the radar measurements. Here, data are used from two airborne profiling (up and down) radars, the W-band Wyoming Cloud Radar and the Ka-band Profiling Radar, allowing for Ka–W-band dual-wavelength ratio (DWR) profiles. The aircraft (the University of Wyoming King Air) also carried a suite of in situ cloud and precipitation probes. This arrangement is optimal for relating the “flight-level” DWR (an average from radar gates below and above flight level) to ice particle size distributions measured by in situ optical array probes, as well as bulk properties such as minimum snow particle density and ice water content. This comparison reveals a strong relationship between DWR and the ice particle median-volume diameter. An optimal range of DWR values ensures the highest retrieval confidence, bounded by the radars’ relative calibration and DWR saturation, found here to be about 2.5–7.5 dB. The DWR-defined size distribution shape is used with a Mie scattering model and an experimental mass–diameter relationship to test retrievals of ice particle concentration and ice water content. Comparison with flight-level cloud-probe data indicate good performance, allowing microphysical interpretations for the rest of the vertical radar transects. 
    more » « less
  4. Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting. 
    more » « less
  5. This is the first paper to approach the problem of bias in the output of a stochastic simulation due to us- ing input distributions whose parameters were estimated from real-world data. We consider, in particular, the bias in simulation-based estimators of the expected value (long-run average) of the real-world system performance; this bias will be present even if one employs unbiased estimators of the input distribution parameters due to the (typically) nonlinear relationship between these parameters and the output response. To date this bias has been assumed to be negligible because it decreases rapidly as the quantity of real-world input data increases. While true asymptotically, this property does not imply that the bias is actually small when, as is always the case, data are finite. We present a delta-method approach to bias estimation that evaluates the nonlinearity of the expected-value performance surface as a function of the input-model parameters. Since this response surface is unknown, we propose an innovative experimental design to fit a response-surface model that facilitates a test for detecting a bias of a relevant size with specified power. We evaluate the method using controlled experiments, and demonstrate it through a realistic case study concerning a healthcare call centre. 
    more » « less