Search Results
You are looking at 1 - 10 of 11 items for
- Author or Editor: Maximilian Maahn x
- Refine by Access: All Content x
Abstract
Retrievals of ice-cloud properties from cloud-radar observations are challenging because the retrieval methods are typically underdetermined. Here, the authors investigate whether additional information can be obtained from higher-order moments and the slopes of the radar Doppler spectrum such as skewness and kurtosis as well as the slopes of the Doppler peak. To estimate quantitatively the additional information content, a generalized Bayesian retrieval framework that is based on optimal estimation is developed. Real and synthetic cloud-radar observations of the Indirect and Semi-Direct Aerosol Campaign (ISDAC) dataset obtained around Barrow, Alaska, are used in this study. The state vector consists of the microphysical (particle-size distribution, mass–size relation, and cross section–area relation) and kinematic (vertical wind and turbulence) quantities required to forward model the moments and slopes of the radar Doppler spectrum. It is found that, for a single radar frequency, more information can be retrieved when including higher-order moments and slopes than when using only reflectivity and mean Doppler velocity but two radar frequencies. When using all moments and slopes with two or even three frequencies, the uncertainties of all state variables, including the mass–size relation, can be considerably reduced with respect to the prior knowledge.
Abstract
Retrievals of ice-cloud properties from cloud-radar observations are challenging because the retrieval methods are typically underdetermined. Here, the authors investigate whether additional information can be obtained from higher-order moments and the slopes of the radar Doppler spectrum such as skewness and kurtosis as well as the slopes of the Doppler peak. To estimate quantitatively the additional information content, a generalized Bayesian retrieval framework that is based on optimal estimation is developed. Real and synthetic cloud-radar observations of the Indirect and Semi-Direct Aerosol Campaign (ISDAC) dataset obtained around Barrow, Alaska, are used in this study. The state vector consists of the microphysical (particle-size distribution, mass–size relation, and cross section–area relation) and kinematic (vertical wind and turbulence) quantities required to forward model the moments and slopes of the radar Doppler spectrum. It is found that, for a single radar frequency, more information can be retrieved when including higher-order moments and slopes than when using only reflectivity and mean Doppler velocity but two radar frequencies. When using all moments and slopes with two or even three frequencies, the uncertainties of all state variables, including the mass–size relation, can be considerably reduced with respect to the prior knowledge.
A New Criterion to Improve Operational Drizzle Detection with Ground-Based Remote SensingAbstract
Light shallow precipitation in the form of drizzle is one of the mechanisms for liquid water removal, affecting cloud lifetime and boundary layer dynamics and thermodynamics. The early formation of drizzle drops is of particular interest for quantifying aerosol–cloud–precipitation interactions. In models, drizzle initiation is represented by the autoconversion, that is, the conversion of liquid water from a cloud liquid water category (where particle sedimentation is ignored) into a precipitating liquid water category. Various autoconversion parameterizations have been proposed in recent years, but their evaluation is challenging due to the lack of proper observations of drizzle development in the cloud. This work presents a new algorithm for Classification of Drizzle Stages (CLADS). CLADS is based on the skewness of the Ka-band radar Doppler spectrum. Skewness is sensitive to the drizzle growth in the cloud: the observed Gaussian Doppler spectrum has skewness zero when only cloud droplets are present without any significant fall velocity. Defining downward velocities positive, skewness turns positive when embryonic drizzle forms and becomes negative when drizzle starts to dominate the spectrum. CLADS identifies spatially coherent structures of positive, zero, and negative skewness in space and time corresponding to drizzle seeding, drizzle growth/nondrizzle, and drizzle mature, respectively. We test CLADS on case studies from the Jülich Observatory for Cloud Evolution Core Facility (JOYCE-CF) and the Barbados Cloud Observatory (BCO) to quantitatively estimate the benefits of CLADS compared to the standard Cloudnet target categorization algorithm. We suggest that CLADS can provide additional observational constraints for understanding the processes related to drizzle formation better.
Abstract
Light shallow precipitation in the form of drizzle is one of the mechanisms for liquid water removal, affecting cloud lifetime and boundary layer dynamics and thermodynamics. The early formation of drizzle drops is of particular interest for quantifying aerosol–cloud–precipitation interactions. In models, drizzle initiation is represented by the autoconversion, that is, the conversion of liquid water from a cloud liquid water category (where particle sedimentation is ignored) into a precipitating liquid water category. Various autoconversion parameterizations have been proposed in recent years, but their evaluation is challenging due to the lack of proper observations of drizzle development in the cloud. This work presents a new algorithm for Classification of Drizzle Stages (CLADS). CLADS is based on the skewness of the Ka-band radar Doppler spectrum. Skewness is sensitive to the drizzle growth in the cloud: the observed Gaussian Doppler spectrum has skewness zero when only cloud droplets are present without any significant fall velocity. Defining downward velocities positive, skewness turns positive when embryonic drizzle forms and becomes negative when drizzle starts to dominate the spectrum. CLADS identifies spatially coherent structures of positive, zero, and negative skewness in space and time corresponding to drizzle seeding, drizzle growth/nondrizzle, and drizzle mature, respectively. We test CLADS on case studies from the Jülich Observatory for Cloud Evolution Core Facility (JOYCE-CF) and the Barbados Cloud Observatory (BCO) to quantitatively estimate the benefits of CLADS compared to the standard Cloudnet target categorization algorithm. We suggest that CLADS can provide additional observational constraints for understanding the processes related to drizzle formation better.
Abstract
The influence of ice hydrometeor shape on the dual-wavelength ratio (DWR) of radar reflectivities at millimeter-wavelength frequencies is studied theoretically and on the basis of observations. Data from dual-frequency (Ka–W bands) radar show that, for vertically pointing measurements, DWR increasing trends with reflectivity Z e are very pronounced when Ka-band Z e is greater than about 0 dBZ and that DWR and Z e values are usually well correlated. This correlation is explained by strong relations between hydrometeor characteristic size and both of these radar variables. The observed DWR variability for a given level of reflectivity is as large as 8 dB, which is in part due to changes in mean hydrometeor shape as expressed in terms of the particle aspect ratio. Hydrometeors with a higher degree of nonsphericity exhibit lower DWR values when compared with quasi-spherical particles because of near-zenith reflectivity enhancements for particles outside the Rayleigh-scattering regime. When particle mass–size relations do not change significantly (e.g., for low-rime conditions), DWR can be used to differentiate between quasi-spherical and highly nonspherical hydrometeors because (for a given reflectivity value) DWR tends to increase as particles become more spherical. Another approach for differentiating among different degrees of nonsphericity for larger scatterers is based on analyzing DWR changes as a function of radar elevation angle. These changes are more pronounced for highly nonspherical particles and can exceed 10 dB. Measurements of snowfall spatiotemporally collocated with spaceborne CloudSat W-band radar and ground-based S-band operational weather radars also indicate that DWR values are generally smaller for ice hydrometeors with higher degrees of nonsphericity, which, for the same level of S-band reflectivity, exhibit greater differential reflectivity values.
Abstract
The influence of ice hydrometeor shape on the dual-wavelength ratio (DWR) of radar reflectivities at millimeter-wavelength frequencies is studied theoretically and on the basis of observations. Data from dual-frequency (Ka–W bands) radar show that, for vertically pointing measurements, DWR increasing trends with reflectivity Z e are very pronounced when Ka-band Z e is greater than about 0 dBZ and that DWR and Z e values are usually well correlated. This correlation is explained by strong relations between hydrometeor characteristic size and both of these radar variables. The observed DWR variability for a given level of reflectivity is as large as 8 dB, which is in part due to changes in mean hydrometeor shape as expressed in terms of the particle aspect ratio. Hydrometeors with a higher degree of nonsphericity exhibit lower DWR values when compared with quasi-spherical particles because of near-zenith reflectivity enhancements for particles outside the Rayleigh-scattering regime. When particle mass–size relations do not change significantly (e.g., for low-rime conditions), DWR can be used to differentiate between quasi-spherical and highly nonspherical hydrometeors because (for a given reflectivity value) DWR tends to increase as particles become more spherical. Another approach for differentiating among different degrees of nonsphericity for larger scatterers is based on analyzing DWR changes as a function of radar elevation angle. These changes are more pronounced for highly nonspherical particles and can exceed 10 dB. Measurements of snowfall spatiotemporally collocated with spaceborne CloudSat W-band radar and ground-based S-band operational weather radars also indicate that DWR values are generally smaller for ice hydrometeors with higher degrees of nonsphericity, which, for the same level of S-band reflectivity, exhibit greater differential reflectivity values.
Abstract
A polarimetric radar–based method for retrieving atmospheric ice particle shapes is applied to snowfall measurements by a scanning Ka-band radar deployed at Oliktok Point, Alaska (70.495°N, 149.883°W). The mean aspect ratio, which is defined by the hydrometeor minor-to-major dimension ratio for a spheroidal particle model, is retrieved as a particle shape parameter. The radar variables used for aspect ratio profile retrievals include reflectivity, differential reflectivity, and the copolar correlation coefficient. The retrievals indicate that hydrometeors with mean aspect ratios below 0.2–0.3 are usually present in regions with air temperatures warmer than approximately from −17° to −15°C, corresponding to a regime that has been shown to be favorable for growth of pristine ice crystals of planar habits. Radar reflectivities corresponding to the lowest mean aspect ratios are generally between −10 and 10 dBZ. For colder temperatures, mean aspect ratios are typically in a range between 0.3 and 0.8. There is a tendency for hydrometeor aspect ratios to increase as particles transition from altitudes in the temperature range from −17° to −15°C toward the ground. This increase is believed to result from aggregation and riming processes that cause particles to become more spherical and is associated with areas demonstrating differential reflectivity decreases with increasing reflectivity. Aspect ratio retrievals at the lowest altitudes are consistent with in situ measurements obtained using a surface-based multiangle snowflake camera. Pronounced gradients in particle aspect ratio profiles are observed at altitudes at which there is a change in the dominant hydrometeor species, as inferred by spectral measurements from a vertically pointing Doppler radar.
Abstract
A polarimetric radar–based method for retrieving atmospheric ice particle shapes is applied to snowfall measurements by a scanning Ka-band radar deployed at Oliktok Point, Alaska (70.495°N, 149.883°W). The mean aspect ratio, which is defined by the hydrometeor minor-to-major dimension ratio for a spheroidal particle model, is retrieved as a particle shape parameter. The radar variables used for aspect ratio profile retrievals include reflectivity, differential reflectivity, and the copolar correlation coefficient. The retrievals indicate that hydrometeors with mean aspect ratios below 0.2–0.3 are usually present in regions with air temperatures warmer than approximately from −17° to −15°C, corresponding to a regime that has been shown to be favorable for growth of pristine ice crystals of planar habits. Radar reflectivities corresponding to the lowest mean aspect ratios are generally between −10 and 10 dBZ. For colder temperatures, mean aspect ratios are typically in a range between 0.3 and 0.8. There is a tendency for hydrometeor aspect ratios to increase as particles transition from altitudes in the temperature range from −17° to −15°C toward the ground. This increase is believed to result from aggregation and riming processes that cause particles to become more spherical and is associated with areas demonstrating differential reflectivity decreases with increasing reflectivity. Aspect ratio retrievals at the lowest altitudes are consistent with in situ measurements obtained using a surface-based multiangle snowflake camera. Pronounced gradients in particle aspect ratio profiles are observed at altitudes at which there is a change in the dominant hydrometeor species, as inferred by spectral measurements from a vertically pointing Doppler radar.
Abstract
Observing ice clouds using zenith pointing millimeter cloud radars is challenging because the transfer functions relating the observables to meteorological quantities are not uniquely defined. Here, the authors use a spectral radar simulator to develop a consistent dataset containing particle mass, area, and size distribution as functions of size. This is an essential prerequisite for radar sensitivity studies and retrieval development. The data are obtained from aircraft in situ and ground-based radar observations during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) campaign in Alaska. The two main results of this study are as follows: 1) An improved method to estimate the particle mass–size relation as a function of temperature is developed and successfully evaluated by combining aircraft in situ and radar observations. The method relies on a functional relation between reflectivity and Doppler velocity. 2) The impact on the Doppler spectrum by replacing measurements of particle area and size distribution by recent analytical expressions is investigated. For this, higher-order moments such as skewness and kurtosis as well as the slopes of the Doppler spectrum are also used as a proxy for the Doppler spectrum. For the area–size relation, it is found that a power law is not sufficient to describe particle area and small deviations from a power law are essential for obtaining consistent higher moments. For particle size distributions, the normalization approach for the gamma distribution of Testud et al., adapted to maximum diameter as size descriptor, is preferred.
Abstract
Observing ice clouds using zenith pointing millimeter cloud radars is challenging because the transfer functions relating the observables to meteorological quantities are not uniquely defined. Here, the authors use a spectral radar simulator to develop a consistent dataset containing particle mass, area, and size distribution as functions of size. This is an essential prerequisite for radar sensitivity studies and retrieval development. The data are obtained from aircraft in situ and ground-based radar observations during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) campaign in Alaska. The two main results of this study are as follows: 1) An improved method to estimate the particle mass–size relation as a function of temperature is developed and successfully evaluated by combining aircraft in situ and radar observations. The method relies on a functional relation between reflectivity and Doppler velocity. 2) The impact on the Doppler spectrum by replacing measurements of particle area and size distribution by recent analytical expressions is investigated. For this, higher-order moments such as skewness and kurtosis as well as the slopes of the Doppler spectrum are also used as a proxy for the Doppler spectrum. For the area–size relation, it is found that a power law is not sufficient to describe particle area and small deviations from a power law are essential for obtaining consistent higher moments. For particle size distributions, the normalization approach for the gamma distribution of Testud et al., adapted to maximum diameter as size descriptor, is preferred.
Atmospheric Ice Particle Shape Estimates from Polarimetric Radar Measurements and In Situ ObservationsAbstract
A remote sensing approach to retrieve the degree of nonsphericity of ice hydrometeors using scanning polarimetric Ka-band radar measurements from a U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Program cloud radar operated in an alternate transmission–simultaneous reception mode is introduced. Nonsphericity is characterized by aspect ratios representing the ratios of particle minor-to-major dimensions. The approach is based on the use of a circular depolarization ratio (CDR) proxy reconstructed from differential reflectivity Z DR and copolar correlation coefficient ρ hυ linear polarization measurements. Essentially combining information contained in Z DR and ρ hυ , CDR-based retrievals of aspect ratios are fairly insensitive to hydrometeor orientation if measurements are performed at elevation angles of around 40°–50°. The suggested approach is applied to data collected using the third ARM Mobile Facility (AMF3), deployed to Oliktok Point, Alaska. Aspect ratio retrievals were also performed using Z DR measurements that are more strongly (compared to CDR) influenced by hydrometeor orientation. The results of radar-based retrievals are compared with in situ measurements from the tethered balloon system (TBS)-based video ice particle sampler and the ground-based multiangle snowflake camera. The observed ice hydrometeors were predominantly irregular-shaped ice crystals and aggregates, with aspect ratios varying between approximately 0.3 and 0.8. The retrievals assume that particle bulk density influencing (besides the particle shape) observed polarimetric variables can be deduced from the estimates of particle characteristic size. Uncertainties of CDR-based aspect ratio retrievals are estimated at about 0.1–0.15. Given these uncertainties, radar-based retrievals generally agreed with in situ measurements. The advantages of using the CDR proxy compared to the linear depolarization ratio are discussed.
Abstract
A remote sensing approach to retrieve the degree of nonsphericity of ice hydrometeors using scanning polarimetric Ka-band radar measurements from a U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Program cloud radar operated in an alternate transmission–simultaneous reception mode is introduced. Nonsphericity is characterized by aspect ratios representing the ratios of particle minor-to-major dimensions. The approach is based on the use of a circular depolarization ratio (CDR) proxy reconstructed from differential reflectivity Z DR and copolar correlation coefficient ρ hυ linear polarization measurements. Essentially combining information contained in Z DR and ρ hυ , CDR-based retrievals of aspect ratios are fairly insensitive to hydrometeor orientation if measurements are performed at elevation angles of around 40°–50°. The suggested approach is applied to data collected using the third ARM Mobile Facility (AMF3), deployed to Oliktok Point, Alaska. Aspect ratio retrievals were also performed using Z DR measurements that are more strongly (compared to CDR) influenced by hydrometeor orientation. The results of radar-based retrievals are compared with in situ measurements from the tethered balloon system (TBS)-based video ice particle sampler and the ground-based multiangle snowflake camera. The observed ice hydrometeors were predominantly irregular-shaped ice crystals and aggregates, with aspect ratios varying between approximately 0.3 and 0.8. The retrievals assume that particle bulk density influencing (besides the particle shape) observed polarimetric variables can be deduced from the estimates of particle characteristic size. Uncertainties of CDR-based aspect ratio retrievals are estimated at about 0.1–0.15. Given these uncertainties, radar-based retrievals generally agreed with in situ measurements. The advantages of using the CDR proxy compared to the linear depolarization ratio are discussed.
Abstract
Remote sensing instruments are heavily used to provide observations for both the operational and research communities. These sensors do not provide direct observations of the desired atmospheric variables, but instead, retrieval algorithms are necessary to convert the indirect observations into the variable of interest. It is critical to be aware of the underlying assumptions made by many retrieval algorithms, including that the retrieval problem is often ill posed and that there are various sources of uncertainty that need to be treated properly. In short, the retrieval challenge is to invert a set of noisy observations to obtain estimates of atmospheric quantities. The problem is often complicated by imperfect forward models, by imperfect prior knowledge, and by the existence of nonunique solutions. Optimal estimation (OE) is a widely used physical retrieval method that combines measurements, prior information, and the corresponding uncertainties based on Bayes’s theorem to find an optimal solution for the atmospheric state. Furthermore, OE also allows the relative contributions of the different sources of error to the uncertainty in the final retrieved atmospheric state to be understood. Here, we provide a novel Python library to illustrate the use of OE for inverse problems in the atmospheric sciences. We introduce two example problems: how to retrieve drop size distribution parameters from radar observations and how to retrieve the temperature profile from ground-based microwave sensors. Using these examples, we discuss common pitfalls, how the various error sources impact the retrieval, and how the quality of the retrieval results can be quantified.
Abstract
Remote sensing instruments are heavily used to provide observations for both the operational and research communities. These sensors do not provide direct observations of the desired atmospheric variables, but instead, retrieval algorithms are necessary to convert the indirect observations into the variable of interest. It is critical to be aware of the underlying assumptions made by many retrieval algorithms, including that the retrieval problem is often ill posed and that there are various sources of uncertainty that need to be treated properly. In short, the retrieval challenge is to invert a set of noisy observations to obtain estimates of atmospheric quantities. The problem is often complicated by imperfect forward models, by imperfect prior knowledge, and by the existence of nonunique solutions. Optimal estimation (OE) is a widely used physical retrieval method that combines measurements, prior information, and the corresponding uncertainties based on Bayes’s theorem to find an optimal solution for the atmospheric state. Furthermore, OE also allows the relative contributions of the different sources of error to the uncertainty in the final retrieved atmospheric state to be understood. Here, we provide a novel Python library to illustrate the use of OE for inverse problems in the atmospheric sciences. We introduce two example problems: how to retrieve drop size distribution parameters from radar observations and how to retrieve the temperature profile from ground-based microwave sensors. Using these examples, we discuss common pitfalls, how the various error sources impact the retrieval, and how the quality of the retrieval results can be quantified.
Abstract
Thorough understanding of aerosols, clouds, boundary layer structure, and radiation is required to improve the representation of the Arctic atmosphere in weather forecasting and climate models. To develop such understanding, new perspectives are needed to provide details on the vertical structure and spatial variability of key atmospheric properties, along with information over difficult-to-reach surfaces such as newly forming sea ice. Over the last three years, the U.S. Department of Energy (DOE) has supported various flight campaigns using unmanned aircraft systems [UASs, also known as unmanned aerial vehicles (UAVs) and drones] and tethered balloon systems (TBSs) at Oliktok Point, Alaska. These activities have featured in situ measurements of the thermodynamic state, turbulence, radiation, aerosol properties, cloud microphysics, and turbulent fluxes to provide a detailed characterization of the lower atmosphere. Alongside a suite of active and passive ground-based sensors and radiosondes deployed by the DOE Atmospheric Radiation Measurement (ARM) program through the third ARM Mobile Facility (AMF-3), these flight activities demonstrate the ability of such platforms to provide critically needed information. In addition to providing new and unique datasets, lessons learned during initial campaigns have assisted in the development of an exciting new community resource.
Abstract
Thorough understanding of aerosols, clouds, boundary layer structure, and radiation is required to improve the representation of the Arctic atmosphere in weather forecasting and climate models. To develop such understanding, new perspectives are needed to provide details on the vertical structure and spatial variability of key atmospheric properties, along with information over difficult-to-reach surfaces such as newly forming sea ice. Over the last three years, the U.S. Department of Energy (DOE) has supported various flight campaigns using unmanned aircraft systems [UASs, also known as unmanned aerial vehicles (UAVs) and drones] and tethered balloon systems (TBSs) at Oliktok Point, Alaska. These activities have featured in situ measurements of the thermodynamic state, turbulence, radiation, aerosol properties, cloud microphysics, and turbulent fluxes to provide a detailed characterization of the lower atmosphere. Alongside a suite of active and passive ground-based sensors and radiosondes deployed by the DOE Atmospheric Radiation Measurement (ARM) program through the third ARM Mobile Facility (AMF-3), these flight activities demonstrate the ability of such platforms to provide critically needed information. In addition to providing new and unique datasets, lessons learned during initial campaigns have assisted in the development of an exciting new community resource.
