Engineering topics
Williams, Christopher
Publications and source records attributed to Williams, Christopher.
Ka-Band ARM Zenith RADAR (KAZR): precipitation mode, spectrally-decomposed
kazrspectrascopepr.c0
Ka-Band ARM Zenith RADAR (KAZR): moderate sensitivity mode, spectrally-decomposed
kazrspectrascopemd.c0
Advanced Precipitation and Boundary Layer Data Products Derived from ARM Radar Wind Profilers
This research project was successful in delivering on four main objectives. First, software was developed to accurately calculate 915-MHz radar wind profiler (RWP) spectrum moments from the recorded Doppler velocity power spectra. Second, software was developed to calibrate the RWP reflectivity factor using collocated surface disdrometer observations. Third, the Python processing code was documented and given to the ARM Infrastructure to produce ARM ‘b level’ calibrated RWP products. Fourth, calibrated RWP products were uploaded to the ARM Archive as PI Products for 10 years of SGP RWP observations and for GoAmazon and TRACER field campaign RWP observations. In addition to working with RWP observations, this research project also worked with KAZR observations to distinguish insects from boundary layer clouds to help improve the ARSCL cloud mask product. The PI worked with senior and early career ARM funded scientists at BNL exploring how to include calibrated RWP moments into future versions of the ARSCL product.
ARM Cloud and Precipitation Measurements and Science Group (CPMSG) 2024 Workshop Report
The mission of the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility is to improve the understanding and representation of cloud and aerosol processes and their interaction with the Earth's surface in Earth system models (ESMs) by providing comprehensive field observations and supporting advanced data analytics. The ARM Cloud and Precipitation Measurements and Science Group (CPMSG) was chartered in March 2019 to help improve the performance and scientific impact of ARM measurements of clouds and precipitation. The group aims to identify and address gaps in measurement capabilities, maximize the scientific impact of ARM data, and effectively serve the scientific community. To achieve these goals, the group includes experts in cloud and precipitation science, as well as representatives from ARM infrastructure, including instrument mentors, engineers, data quality officers, and data product translators. Prior to CPMSG, early discussions on cloud and precipitation measurements primarily focused on improving radar systems, but have since evolved to include a broader scope involving radiometers and other instruments. Since its formation, the CPMSG has gathered feedback using science traceability matrices. CPMSG aims to keep these as living documents to show the measurement needs, scientific drivers, roadblocks, maturity of measurements and retrievals, and pathways to model improvements. The group meets quarterly to discuss and prioritize measurement and operational improvements.
915rwpprecipcor
Radar Wind Profiler (RWP) has operational modes: short duration high resolution (referred to as 'high' or 'high mode') and long duration low resolution (referred to as 'low' or 'low mode'). The motivation for this is to create a VAP that quality controls (noise, sig echo filtering), merges, co-grids, and performs initial analyses on these RWP modes. This is done to streamline downstream VAP development on the RWP, as well as enable incorporation of echo properties into the ARSCL chain. The proposed VAP will use the RWPPRECIP b0/b1 level data stream as input. Additional consideration for the incorporation of other ARM data streams to improve this VAP was performed, but not included at this time.
Evidence and attribution of the enhanced land carbon sink
Climate change has been partially mitigated by an increasing net land carbon sink in the terrestrial biosphere; understanding the processes that drive this sink is thus essential for protecting, managing and projecting this important ecosystem service. In this Review, we examine evidence for an enhanced land carbon sink and attribute the observed response to drivers and processes. This sink has doubled from 1.2 ± 0.5 PgC yr -1 in the 1960s to 3.1 ± 0.6 PgC yr -1 in the 2010s. This trend results largely from carbon dioxide fertilization increasing photosynthesis (driving an increase in the annual land carbon sink of >2 PgC globally since 1900), mainly in tropical forest regions, and elevated temperatures reducing cold limitation, mainly at higher latitudes. Continued long-term land carbon sequestration is possible through the end of this century under multiple emissions scenarios, especially if nature-based climate solutions and appropriate ecosystem management are used. In conclusion, a new generation of globally distributed field experiments is needed to improve understanding of future carbon sink potential by measuring belowground carbon release, the response to carbon dioxide enrichment, and long-term shifts in carbon allocation and turnover.
Calibrated Radar Wind Profiler (RWP) Moments
The SGP Central Facility (C1) radar wind profiler (RWP) was calibrated using nearby surface disdrometer observations. Between 2011 and 2019, the SGP C1 RWP operated in two modes. The vertically pointing mode (named the precipitation mode) transmitted a short and long pulse length to have two different range resolutions and the beam-swinging mode (named the wind mode) transmitted one pulse length into three different beam directions. The precipitation-mode observations were available and calibrated from April 2011 through mid-August 2019. The wind-mode observations were available and calibrated between April 2014 and March 2019. The RWP spectra were processed to account for Nyquist velocity aliasing and coherent integration filtering effects before calculating the spectrum moments. During intense precipitation events, the calculated signal-to-noise ratio (SNR) is biased low due to signal power being distributed across the velocity spectrum such that some signal power is erroneously included in the noise level estimate, causing the noise level to be biased high. To correct for the low SNR bias, a new noise level is estimated using observations without precipitation and the SNR is increased accordingly. The adjusted SNR was converted to radar reflectivity factor and then calibrated against a nearby surface disdrometer. The calibration methodology is fully described in: Williams, CR, J Barrio, PE Johnston, P Muradyan, and S Giangrande. 2023. “Calibrating radar wind profiler reflectivity factor using surface disdrometer observations.” Atmospheric Measurement Techniques, https://egusphere.copernicus.org/preprints/2023/egusphere-2022-1405/
Across‐model spread and shrinking in predicting peatland carbon dynamics under global change
Abstract Large across‐model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across‐model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO 2 levels and five temperature levels. We developed a matrix‐based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter‐model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter‐model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO 2 (eCO 2 ) concentrations in the measurements, simulated NEP responded positively to eCO 2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across‐model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.
Multi-atom Pt and PtRu catalysts for high performance AEMFCs with ultra-low PGM content
To reduce the platinum group metal (PGM) loading in anion exchange membrane fuel cells (AEMFCs), it is important to transition to catalysts with very low PGM content, and eventually to catalysts that are completely PGM-free. Here, in this work, four supported low-PGM Pt and PtRu catalysts were prepared using a new, simple, scalable technique: Controlled Surface Tension (CST) method. CST allows for a high density of very small multi-atom clusters. Catalysts were physically characterized using a wide array of techniques and tested for their ORR and HOR activity both ex-situ and integrated into operating AEMFCs. The PGM loading was reduced by a factor of 14 while achieving comparable performance to commercial catalysts. AEMFCs were also assembled with ultralow PGM loading (0.05 mgPGM cm -2 ), where PtRu anodes were paired with Fe–N–C cathodes to achieve a specific power of 25 W/mg PGM (40 W/mg Pt ).
Evaluation of Cloud-Resolving and Limited Area Model Intercomparison Simulations Using TWP-ICE Observations: Precipitation Microphysics - Part 2
Ten 3-D cloud-resolving model (CRM) simulations and four 3-D limited area model (LAM) simulations of an intense mesoscale convective system observed on 23-24 January 2006 during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) are compared with each other and with observations and retrievals from a scanning polarimetric radar, colocated UHF and VHF vertical profilers, and a Joss-Waldvogel disdrometer in an attempt to explain a low bias in simulated stratiform rainfall. Despite different forcing methodologies, similar precipitation microphysics errors appear in CRMs and LAMs with differences that depend on the details of the bulk microphysics scheme used. One-moment schemes produce too many small raindrops, which biases Doppler velocities low, but produces rainwater contents (RWCs) that are similar to observed. Two-moment rain schemes with a gamma shape parameter (mu) of 0 produce excessive size sorting, which leads to larger Doppler velocities than those produced in one-moment schemes but lower RWCs. Two-moment schemes also produce a convective median volume diameter distribution that is too broad relative to observations and, thus, may have issues balancing raindrop formation, collision-coalescence, and raindrop breakup. Assuming a mu of 2.5 rather than 0 for the raindrop size distribution improves one-moment scheme biases, and allowing mu to have values greater than 0 may improve excessive size sorting in two-moment schemes. Underpredicted stratiform rain rates are associated with underpredicted ice water contents at the melting level rather than excessive rain evaporation, in turn likely associated with convective detrainment that is too high in the troposphere and mesoscale circulations that are too weak. A limited domain size also prevents a large, well-developed stratiform region like the one observed from developing in CRMs, although LAMs also fail to produce such a region.
Drop Size Distribution - Based Separation of Stratiform and Convective Rain
For applications in hydrology and meteorology, it is often desirable to separate regions of stratiform and convective rain from meteorological radar observations, both from ground-based polarimetric radars and from space-based dual frequency radars. In a previous study by Bringi et al. (2009), dual frequency profiler and dual polarization radar (C-POL) observations in Darwin, Australia, had shown that stratiform and convective rain could be separated in the log10(Nw) versus Do domain, where Do is the mean volume diameter and Nw is the scaling parameter which is proportional to the ratio of water content to the mass weighted mean diameter. Note, Nw and Do are two of the main drop size distribution (DSD) parameters. In a later study, Thurai et al (2010) confirmed that both the dual-frequency profiler based stratiform-convective rain separation and the C-POL radar based separation were consistent with each other. In this paper, we test this separation method using DSD measurements from a ground based 2D video disdrometer (2DVD), along with simultaneous observations from a collocated, vertically-pointing, X-band profiling radar (XPR). The measurements were made in Huntsville, Alabama. One-minute DSDs from 2DVD are used as input to an appropriate gamma fitting procedure to determine Nw and Do. The fitted parameters - after averaging over 3-minutes - are plotted against each other and compared with a predefined separation line. An index is used to determine how far the points lie from the separation line (as described in Thurai et al. 2010). Negative index values indicate stratiform rain and positive index indicate convective rain, and, moreover, points which lie somewhat close to the separation line are considered 'mixed' or 'transition' type precipitation. The XPR observations are used to evaluate/test the 2DVD data-based classification. A 'bright-band' detection algorithm was used to classify each vertical reflectivity profile as either stratiform or convective, depending on whether or not a clearly-defined melting layer is present at an expected height, and if present, maximum reflectivity within the melting layer as well as the corresponding height are determined. We will present results of quantitative comparisons between the XPR observations-based classifications and the simultaneous 2DVD data-based classifications. Time series comparisons will be presented for thirteen events in Huntsville.
A Framework for the Generation and Dissemination of Drop Size Distribution (DSD) Characteristics Using Multiple Platforms
Proper characterization of the precipitation drop size distribution (DSD) is integral to providing realistic and accurate space- and ground-based precipitation retrievals. Current technology allows for the development of DSD products from a variety of platforms, including disdrometers, vertical profilers and dual-polarization radars. Up to now, however, the dissemination or availability of such products has been limited to individual sites and/or field campaigns, in a variety of formats, often using inconsistent algorithms for computing the integral DSD parameters, such as the median- and mass-weighted drop diameter, total number concentration, liquid water content, rain rate, etc. We propose to develop a framework for the generation and dissemination of DSD characteristic products using a unified structure, capable of handling the myriad collection of disdrometers, profilers, and dual-polarization radar data currently available and to be collected during several upcoming GPM Ground Validation field campaigns. This DSD super-structure paradigm is an adaptation of the radar super-structure developed for NASA s Radar Software Library (RSL) and RSL_in_IDL. The goal is to provide the DSD products in a well-documented format, most likely NetCDF, along with tools to ingest and analyze the products. In so doing, we can develop a robust archive of DSD products from multiple sites and platforms, which should greatly benefit the development and validation of precipitation retrieval algorithms for GPM and other precipitation missions. An outline of this proposed framework will be provided as well as a discussion of the algorithms used to calculate the DSD parameters.
Intercomparison of Various Algorithms in Determining Convective and Stratiform Precipitation
The paper presents the intercomparison of various algorithms in determining the convective and stratiform precipitation in tropical convection. The latent heat released from precipitation is the driving mechanism for the general circulation of the atmosphere and exhibits different vertical profile in convective and stratiform regimes of the tropical convection. The primary goal of this paper is to evaluate the performance of the operational convective/stratiform algorithm that is employed by the NASA Tropical Rainfall Measuring Mission (TRMM) ground validation program. The operational algorithm is based on the texture of radar reflectivity field at constant altitude. The TRMM ground validation program continuously monitors the tropical convection through ground based radar and rain gauge network at selected sites around the tropics and generates the rainfall products including the convective/stratiform rainfall map that can readily be used by the TRMM satellite program. An evaluation and improvement of these rainfall products requires additional measurements that can be available through field campaigns. The data used in this study was collected during the second phase of the Texas Florida Underflights (TEFLUN-B) field campaign in East Central Florida, a primary site of the TRMM ground validation program. The TEFLUN-B field campaign was conducted in August and September 1998 in support of the NASA Tropical Rainfall Measuring Mission (TRMM) satellite observations. An objective of the TEFLUN-B field campaign was to obtain the ground based precipitation measurements of Florida convection at different temporal and spatial scales. The field campaign had various components of precipitation measurements including dual wavelength precipitation profiler, various types of disdrometers and rain gauges, collocated about 40 km southwest of Melbourne NEXRAD radar. In this study, we examine the profiler based reflectivity and Doppler velocity and the disdrometer based drop size distribution observations to evaluate the operational scanning radar (i.e. NEXRAD) based convective/stratiform algorithm.