Search NASASearch

Engineering topics

Benjamin T Johnson

Publications and source records attributed to Benjamin T Johnson.

The Community Radiative Transfer Model (CRTM): Community-Focused Collaborative Model Development Accelerating Research to Operations

The Joint Center for Satellite Data Assimilation (JCSDA) Community Radiative Transfer Model (CRTM) is a fast, 1-D radiative transfer model used in numerical weather prediction, calibration/validation, etc. across multiple federal agencies and universities. The key benefit of the CRTM is that it is a satellite simulator. It provides a highly accurate representation of satellite radiances by using the specific sensor response functions convolved with a line-by-line radiative transfer model (LBLRTM). CRTM covers the spectral ranges consistent with all present operational and most research satellites, from visible to microwave. The capability to simulate ultraviolet radiances and support space-based radar sensors is being added over the next two years in CRTM Version 3.0. In addition to simulated radiances, the CRTM also provides Jacobian outputs needed to interpret satellite observations for numerical weather prediction. The Jacobian estimates how changes in geophysical parameters affect simulated measurements from satellite sensors. Using the Jacobian in modeling and weather prediction improves the accuracy and efficiency of data analysis, leading to better weather predictions. The CRTM model's success and growth depend on community contributions and evaluation. To facilitate this, we have made the CRTM highly accessible through modular programming, clear documentation and tutorials, public domain licensing, unfettered public access via Github, and a clear path to operational implementation for innovative research. We encourage and welcome contributions from the community to help us continue to improve the CRTM.

Benjamin T Johnson

pyCRTM: A Python Interface for the Community Radiative Transfer Model

The Community Radiative Transfer Model (CRTM) is a powerful and versatile scalar radiative transfer model for satellite data assimilation and remote sensing applications. It is implemented as an object-oriented Fortran library, enabling flexible code development and optimal runtime performance on clusters. The downsides of the Fortran interface are a steep learning curve for students and the reduced productivity of users that is typical for static compiled languages, in contrast to dynamic interpreted languages like Python. pyCRTM is a new software framework that directly interfaces the CRTM Fortran data structures and procedures in Python, leveraging both the simplicity and ease of use of Python syntax as well as the flexibility arising from the vast contemporary Python ecosystem. The goal of pyCRTM is to lower the barrier of entry for university students to learn and use the CRTM and to boost the productivity of researchers seeking to create new methods in radiative transfer and data assimilation, or seeking to apply the CRTM to study atmospheric phenomena without having to go through the pre-existing complexity of the CRTM Fortran interface.

Python

Towards a Mass-Consistent Methodology for Realistic Melting Hydrometeor Retrieval

To address the acute challenge posed by the melting layer to accurate surface precipitation retrievals from space, we en-sure the compositional consistency in ice, liquid, and total masses of synthetic melting hydrometeors with a method of stochastic compensation. The method is applied to simulated melting hydrometeors prior to calculating their scattering properties using the discrete dipole approximation (DDA). We investigate the impact of this stochastic compensation to calculated scattering properties by contrasting it with a naïve approach and report our findings.

precipitation remote sensing

Microwave Properties of Ice-Phase Hydrometeors for Radar and Radiometers: Sensitivity to Model Assumptions

A simplied framework is presented for assessing the qualitative sensitivities of computed microwave properties, satellite brightness temperatures, and radar reflectivities to assumptions concerning the physical properties of ice-phase hydrometeors. Properties considered included the shape parameter of a gamma size distribution andthe melted-equivalent mass median diameter D0, the particle density, dielectric mixing formula, and the choice of complex index of refraction for ice. We examine these properties at selected radiometer frequencies of 18.7, 36.5, 89.0, and 150.0 GHz; and radar frequencies at 2.8, 13.4, 35.6, and 94.0 GHz consistent with existing and planned remote sensing instruments. Passive and active microwave observables of ice particles arefound to be extremely sensitive to the melted-equivalent mass median diameter D0 ofthe size distribution. Similar large sensitivities are found for variations in the ice vol-ume fraction whenever the geometric mass median diameter exceeds approximately 1/8th of the wavelength. At 94 GHz the two-way path integrated attenuation is potentially large for dense compact particles. The distribution parameter mu has a relatively weak effect on any observable: less than 1-2 K in brightness temperature and up to 2.7 dB difference in the effective radar reflectivity. Reversal of the roles of ice and air in the MaxwellGarnett dielectric mixing formula leads to a signicant change in both microwave brightness temperature (10 K) and radar reflectivity (2 dB). The choice of Warren (1984) or Warren and Brandt (2008) for the complex index of refraction of ice can produce a 3%-4% change in the brightness temperature depression.

Microwave properties

Detection Thresholds of Falling Snow from Satellite-Borne Active and Passive Sensors

Precipitation, including rain and snow, is a critical part of the Earth's energy and hydrology cycles. Precipitation impacts latent heating profiles locally while global circulation patterns distribute precipitation and energy from the equator to the poles. For the hydrological cycle, falling snow is a primary contributor in northern latitudes during the winter seasons. Falling snow is the source of snow pack accumulations that provide fresh water resources for many communities in the world. Furthermore, falling snow impacts society by causing transportation disruptions during severe snow events. In order to collect information on the complete global precipitation cycle, both liquid and frozen precipitation must be collected. The challenges of estimating falling snow from space still exist though progress is being made. These challenges include weak falling snow signatures with respect to background (surface, water vapor) signatures for passive sensors over land surfaces, unknowns about the spherical and non-spherical shapes of the snowflakes, their particle size distributions (PSDs) and how the assumptions about the unknowns impact observed brightness temperatures or radar reflectivities, differences in near surface snowfall and total column snow amounts, and limited ground truth to validate against. While these challenges remain, knowledge of their impact on expected retrieval results is an important key for understanding falling snow retrieval estimations. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types. For example, can a lake effect snow system with low (approx 2.5 km) cloud tops having an ice water content (IWC) at the surface of 0.25 g / cubic m and dendrite snowflakes be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. Here, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases since simulations provide all the information to determine the measurements from space and the ground truth. Sensitivity analyses were performed to better ascertain the relationships between multifrequency microwave and millimeter-wave sensor observations and the falling snow/underlying field of view. In addition, thresholds of detection for various sensor channel configurations, snow event system characteristics, snowflake particle assumptions, and surface types were studied. Results will be presented for active radar at Ku, Ka, and W-band and for passive radiometer channels from 10 to 183 GHz.

Gail Skofronick-Jackson

Surface and Atmospheric Contributions to Passive Microwave Brightness Temperatures

Physically-based passive microwave precipitation retrieval algorithms require a set of relationships between satellite observed brightness temperatures (TB) and the physical state of the underlying atmosphere and surface. These relationships are typically non-linear, such that inversions are ill-posed especially over variable land surfaces. In order to better understand these relationships, this work presents a theoretical analysis using brightness temperature weighting functions to quantify the percentage of the TB resulting from absorption/emission/reflection from the surface, absorption/emission/scattering by liquid and frozen hydrometeors in the cloud, the emission from atmospheric water vapor, and other contributors. The results are presented for frequencies from 10 to 874 GHz and for several individual precipitation profiles as well as for three cloud resolving model simulations of falling snow. As expected, low frequency channels (<89 GHz) respond to liquid hydrometeors and the surface, while the higher frequency channels become increasingly sensitive to ice hydrometeors and the water vapor sounding channels react to water vapor in the atmosphere. Low emissivity surfaces (water and snow-covered land) permit energy downwelling from clouds to be reflected at the surface thereby increasing the percentage of the TB resulting from the hydrometeors. The slant path at a 53deg viewing angle increases the hydrometeor contributions relative to nadir viewing channels and show sensitivity to surface polarization effects. The TB percentage information presented in this paper answers questions about the relative contributions to the brightness temperatures and provides a key piece of information required to develop and improve precipitation retrievals over land surfaces.

Gail Skofronick-Jackson

Geophysical Retrievals and Cloud Analyses from Merged Airborne Radiometer Datasets Covering 10–684 GHz

Airborne microwave radiometers provide insight about numerous aspects of Earth’s atmosphere and yield critical validation datasets for spaceborne radiometers. Three radiometers that are important to NASA’s airborne remote-sensing arsenal include: the Advanced Microwave Precipitation Radiometer (AMPR), covering 10–85 GHz; the Conical Scanning Millimeter-wave Imaging Radiometer (CoSMIR), covering 50–183 GHz; and the Compact Scanning Submillimeter-wave Imaging Radiometer (CoSSIR), covering 170–684 GHz. The NASA field campaigns of interest to this study include: the Integrated Precipitation and Hydrology Experiment (IPHEx) in 2014, the Olympic Mountains Experiment and Radar Definition Experiment (OLYMPEX/RADEX) in 2015–2016, the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) in 2020–2023, and the Airborne Lightning Observatory for FEGS and TGFs (ALOFT) in 2023. To provide a more comprehensive perspective on clouds and precipitation observed during these airborne field campaigns, AMPR data were merged spatiotemporally with CoSMIR data for IPHEx, OLYMPEX/RADEX, and IMPACTS (2020 and 2022), while considering differences in instrument characteristics and operations, providing brightness temperature (Tb) values from 10–183 GHz in a common background grid throughout each flight. AMPR and CoSSIR data were similarly merged for IMPACTS (2023) and ALOFT, providing a common background grid with Tb values covering 10–684 GHz throughout each flight. These merged Tb data were employed in geophysical retrievals using the Community Radiative Transfer Model (CRTM), an Eddington radiative transfer model, and a one-dimensional variational (1DVAR) inversion method. Retrievals of cloud liquid water path were of primary interest. This presentation will include an overview of the methods for the radiometer data mergers, the radiative transfer methods, the geophysical retrievals, and detailed results from examining trends in Tb and cloud liquid water path in clouds, precipitation, and cloud-to-precipitation transition zones.

Corey G Amiot