Results for “Atmospheric correction” — Page 18

DOE OSTI2025

A Cloud-Tracking Data Set for the CSAPR2 Adaptive Scanning during TRACER

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility (Mather and Voyles 2013) deployed the first ARM Mobile Facility (AMF1; Miller et al. 2016) near LaPorte, Texas to support the Tracking Aerosol Convection Interactions Experiment (TRACER) (Jensen et al. 2025) near Houston, Texas. From October 2021 to September 2022, AMF1 was deployed to 29.67° N, 95.06° W near LaPorte, Texas and the 2nd Generation C-band Scanning ARM Precipitation Radar (CSAPR2) was deployed to a supplementary site at 29.53° N, 95.28° W (Figure 1). During an intensive operational period (IOP) from 1 June to 30 September 2022, the CSAPR2 sampled precipitation echoes in an adaptive scanning mode following the Multisensor Agile Adaptive Scanning (MAAS) framework (Kollias et al. 2020). MAAS helped optimize the CSAPR2 scan strategy to perform frequent plan position indicator (PPI) and range height indicator (RHI) scans (Lamer et al. 2023). Details of the CSAPR2 scanning, data processing, and calibration procedures used by the principal investigator (PI), and the PI data files are described by Oue et al. (2023). Details of the CSAPR2 operational performance, ARM data processing and correction procedures, and data quality masks are described by Feng et al. (2024a).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Convection-Permitting Ensembles of an Isolated Mountain Thunderstorm during RELAMPAGO/CACTI

Abstract The north–south-oriented Sierras de Córdoba (SDC) ridge in central Argentina is noted for initiating thunderstorms that may grow into intense mesoscale convective systems (MCSs). It also initiates more isolated, shorter-lived cells under weaker synoptic forcing. These cells are less impactful than MCSs but may be difficult to predict in convective-scale numerical weather prediction (NWP) due to their strong sensitivities to subgrid and partially resolved processes. To study the mechanisms and predictability of such cells, convection-permitting ensemble simulations were conducted of an isolated, diurnally forced SDC thunderstorm during Cloud, Aerosol, and Complex Terrain Interactions (CACTI)/Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO). The rich observational data facilitated detailed ensemble verification, where dry biases in the surface energy balance and soil moisture were identified. These biases promoted rapid removal of convective inhibition and an early onset of precipitating cells over the SDC that were shallower and weaker than the observed cell. Correction, and then overcorrection, of the soil moisture bias in two successive ensembles was required to rectify the surface energy balance and improve the representation of the SDC cell. Nevertheless, substantial ensemble variability in convective precipitation was found, with some members producing more widespread convection than observed and others producing no deep convection at all. This variability was largely explained by a combination of thermodynamic and dynamic mechanisms, dominated by a positive sensitivity of convective precipitation to preconvective moist instability over the ridge. Secondary sensitivities were found to low-level upward mass flux and midlevel cross-barrier winds, the latter of which caused gravity waves with elevated downdrafts that tended to suppress incipient clouds.

Lopez, Andres [Department of Atmospheric and Ocean↗

DOE OSTI2020

Longitudinal coherence and short-term wind speed prediction based on a nacelle-mounted Doppler lidar

The spatial structure of turbulence in atmospheric boundary layer flows is highly relevant to wind energy. In particular, wind turbine control strategies based on inflow preview measurements require knowledge of the longitudinal evolution of turbulent flow as it approaches the rotor. These upstream measurements are usually obtained with nacelle-mounted wind lidars. In contrast to traditional in situ anemometry, lidars collect measurements within a probe volume which varies in size depending on the technology of the commercial system being used. Here, we address two issues related to the use of wind lidar to measure the incoming flow to a wind turbine: (i) whether existing longitudinal coherence models can be used to predict flow at the rotor, based on measurements performed at a distance away from the rotor; and (ii) what effect probe-volume averaging has on the inflow predictions. These two questions are critical to the design and implementation of robust wind turbine control strategies. To address these questions, we perform field measurements and large-eddy simulations to determine which incoming flow structures can be readily predicted with existing coherence models, and which require additional corrections to account for lidar volumetric averaging effects. Results reveal that the wind turbine induction zone has a negligible impact on the longitudinal coherence and first-order turbulence quantities, such as the standard deviation of velocity fluctuations. However, the phase of the signal, from which advection time periods of the turbulent structures are derived, is affected by the rotor blockage effect.

17 WIND ENERGY↗

DOE OSTI2023

Mini-Micropulse Lidar data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (mini-MPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. The mini-MPL deployed with ROAM-V is a 532 nm elastic and polarization lidar. The mini-MPL outputs normalized relative backscatter (NRB) derived from raw signal after after-pulse, overlap, and dead-time correction calibration. The depolarization ratio is calculated from the co-polarized and cross-polarized NRB (Flynna et al., 2007). The NRB and depolarization ratio data are resampled from the original data at 1-minute intervals. The vertical resolution of the mini-MPL data is 15 meters. The mini-MPL data can be used to determine the boundary layer, cloud top, and cloud bottom height and can be used to retrieve aerosol type and concentration profile. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu . Flynna, C. J., Mendozaa, A., Zhengb, Y., & Mathurb, S. (2007). Novel polarization-sensitive micropulse lidar measurement technique. Optics express, 15(6), 2785-2790.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Unveiling four decades of intensifying precipitation from tropical cyclones using satellite measurements

Increases in precipitation rates and volumes from tropical cyclones (TCs) caused by anthropogenic warming are predicted by climate modeling studies and have been identified in several high intensity storms occurring over the last half decade. However, it has been difficult to detect historical trends in TC precipitation at time scales long enough to overcome natural climate variability because of limitations in existing precipitation observations. We introduce an experimental global high-resolution climate data record of precipitation produced using infrared satellite imagery and corrected at the monthly scale by a gauge-derived product that shows generally good performance during two hurricane case studies but estimates higher mean precipitation rates in the tropics than the evaluation datasets. General increases in mean and extreme rainfall rates during the study period of 1980–2019 are identified, culminating in a 12–18%/40-year increase in global rainfall rates. Overall, all basins have experienced intensification in precipitation rates. Increases in rainfall rates have boosted the mean precipitation volume of global TCs by 7–15% over 40 years, with the starkest rises seen in the North Atlantic, South Indian, and South Pacific basins (maximum 59–64% over 40 years). In terms of inland rainfall totals, year-by-year trends are generally positive due to increasing TC frequency, slower decay over land, and more intense rainfall, with an alarming increase of 81–85% seen from the strongest global TCs. As the global trend in precipitation rates follows expectations from warming sea surface temperatures (11.1%/°C), we hypothesize that the observed trends could be a result of anthropogenic warming creating greater concentrations of water vapor in the atmosphere, though retrospective studies of TC dynamics over the period are needed to confirm.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

DOE OSTI2023

CPC data from TAMU TRACER campaign in the Houston, TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.3 LPM is directed to the CPC. To correct the data for particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate total particle loss through the inlet for every SMPS scan at a single deployment location. Second, the CPC data was corrected for particle losses using the average of the total losses per scan. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. An additional flatline correction of 14% was applied to all CPC data based on laboratory calibration prior to and after the campaign. In laboratory calibrations, we identified that for the same sample of air the field CPC (GRIMM Model 5.403 CPC) undercounted the total concentration of particles by 14% compared to the ground truth laboratory CPC (TSI Model 3750). This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

CACTI CSAPR2 Taranis Retrievals

Taranis is an end-to-end processing chain for radar data written in Python with C extension for computation performance. Features include: masking for quality control, specific differential phase (Kdp), attenuation correction for reflectivity factor (Z) and differential reflectivity (Zdr) in rain, and additional geophysical retrievals. Retrievals are mostly drawn from literature or open-source software when appropriate, and have been tested, tuned, and modified to work with one another cohesively rather than using isolated off-the-shelf algorithms. Incorporated algorithms include hydrometeor (echo) identification, rain water content, raindrop mass-weighted mean diameter (gamma size distribution assumption), and rainfall rate (QPE). Taranis data sets exist for CSAPR2 PPI, HSRHI, and sector RHI scans. Cartesian-gridded data sets were also produced as well as a near-surface rain rate retrieval. More details can be found in the README.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

When Will MISR Detect Rising High Clouds?

It is predicted by both theory and models that high-altitude clouds will occur higher in the atmosphere as a result of climate warming. This produces a positive longwave feedback and has a substantial impact on the Earth's response to warming. This effect is well established by theory, but is poorly constrained by observations, and there is large spread in the feedback strength between climate models. For this study, we use the NASA Multi-angle Imaging SpectroRadiometer (MISR) to examine changes in Cloud-Top-Height (CTH). MISR uses a stereo-imaging technique to determine CTH. This approach is geometric in nature and insensitive to instrument calibration and therefore is well suited for trend analysis and studies of variability on long time scales. In this article we show that the current MISR record does have an increase in CTH for high-altitude cloud over Southern Hemisphere (SH) oceans but not over Tropical or the Northern Hemisphere (NH) oceans. We use climate model simulations to estimate when MISR might be expected to detect trends in CTH, that include the NH. The analysis suggests that according to the models used in this study MISR should detect changes over the SH ocean earlier than the NH, and if the model predictions are correct should be capable of detecting a trend over the Tropics and NH very soon (3–10 years). This result highlights the potential value of a follow-on mission to MISR, which no longer maintains a fixed equator crossing time and is unlikely to be making observations for another 10 years.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Shockwave Arrival Times from Operation Redwing and Operation Upshot-Knothole

Back in the 1940s, 50s, and 60s, 210 nuclear tests were conducted in the atmosphere. As part of the testing program, EG&G performed shockwave arrival-time measurements for different shots in different operations using pressure transducers at different distances from ground zero. They then developed correlations unique to each set of data. They did not attempt to correlate all data onto one curve. As a fireball moves, material is picked up and pulled in, which is called entrainment. This extra material increases the mass of the fireball, and as the energy is spread out over a larger and larger mass, the specific internal energy of the fireball gases drops, which causes the temperature to drop as well. When the energy in the solid angle of a spherical shockwave hits the surface, it is reflected. Shortly thereafter, it catches up with the outward-moving shockwave and re-combines with it, causing it to expand at a slightly faster rate. As the shockwave grows, it approaches a hemispherical shape. The purpose of this project is to demonstrate that the shockwave arrival-time data can, in fact, be correlated onto a single function when the shockwave time-dependent solution is modified to include a geometric factor, θ, and is corrected for an entrainment coefficient.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2022

Understanding the Physics Representation of Deep Learning Models in Environmental Applications

Deep learning (DL) models have been popular in earth and environmental modeling and analysis, which exhibit huge potential in capturing and reconstructing the non-linearity of relevant environmental processes. They are extensively used as analytical tools or emulators for multiple domains (atmosphere, land surface, ocean, and biogeochemistry). Despite their success, their internal working mechanism remains largely unknown. Such a lack of knowledge hinders the identification of physically consistent models that are fully adaptive to non-stationary climate, as well as the development of physics-informed machine learning such as physics-informed neural network (PINN). To establish preliminary knowledge and framework of such physics representation evaluation, this project focuses on an improved understanding of DL models in the environmental applications. DL models are increasingly applied to environmental modeling and prediction. However, they have been evaluated mostly from a performance perspective, and there is a gap in understanding how they represent the known physics internally. Such knowledge is especially critical when applying DL models under climate change conditions, where new inputs are likely outside the ranges of the training datasets. In this project, we reveal how the known physical processes are represented within DL models from both statistical and mechanistic perspectives. Leveraging the traditional model evaluations that focus more on the accuracies of predictions, we establish a framework that examines both the accuracy and physics representation of DL models. This analysis framework can identify DL models that make the correct predictions based on correct physics, thus enhancing the existing explainable artificial intelligence (explainable-AI) portfolio. It lays a foundation for developing novel metrics to evaluate the emerging DL models in environmental applications. This knowledge also informs the development of physics-informed DL models by revealing the direct connections between the known physical processes and specific model components or structures.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

xsacrcfrspccopol (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

xsacrcfrspcxpol (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

xsacrcfrspccross (a1)

ARM's scanning cloud radars are fully coherent dual-frequency, dual-polarization Doppler radars mounted on a common scanning pedestal. Each pedestal includes a Ka-band radar (2kW peak power) and the deployment location determines whether the second radar is a W-band (WSACR; 1.7 kW peak power) or X-band (XSACR; 20 kW peak power). At ARM's tropical sites, X-band radars are paired with the Ka-band because they are better suited for the atmospheric attenuation in this region. Beamwidth for the X-band is approximately 1 degree, and the Ka-band beamwidth is roughly 0.3 degrees. Due to the narrow antenna beamwidth, ARM’s scanning cloud radars use scanning strategies unlike typical weather radars. Rather than focusing on plan position indicator, or PPI, scans, the XSACR uses range height indicator, or RHI, scans at numerous azimuths to obtain cloud volume data. For the second ARM Mobile Facility (AMF2), the XSACR can be dismounted from the pedestal it shares with the KASACR and mounted on a separate pedestal. This allows the XSACR to operate more like a weather (precipitation) radar in deployments where such local coverage is lacking. Measurements collected with the XSACR are copolar and cross-polar radar reflectivity, Doppler velocity, spectra width and spectra when not scanning, differential Reflectivity (Zdr), correlation coefficient (rho-hv), and specific differential phase (phi-dp). XSACR data from the 2018–2019 Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina are now available as b1-level products. Building on the original CACTI operational data, the b1-level products feature improved data quality resulting from extensive analyses and corrections. The data are cross-calibrated to a common point, datastreams are corrected for operational issues that occurred during the campaign, and several data quality masks and basic derived products are incorporated. For more information, read the CACTI radar b1-level processing report.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part II: Idealized Markov Chain Monte Carlo Experiments

Observationally informed development of a new framework for bulk rain microphysics, the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS; described in Part I of this study), is demonstrated. This scheme’s development is motivated by large uncertainties in cloud and weather simulations associated with approximations and assumptions in existing microphysics schemes. Here, a proof-of-concept study is introduced using a Markov chain Monte Carlo sampling algorithm with BOSS to probabilistically estimate microphysical process rates and parameters directly from a set of synthetically generated rain observations. The framework utilized is an idealized steady-state one-dimensional column rainshaft model with specified column-top rain properties and a fixed thermodynamical profile. Different configurations of BOSS—flexibility being a key feature of this approach—are constrained via synthetic observations generated from a traditional three-moment bulk microphysics scheme. The ability to retrieve correct parameter values when the true parameter values are known is illustrated. For cases when there is no set of true parameter values, the accuracy of configurations of BOSS that have different levels of complexity is compared. It is found that addition of the sixth moment as a prognostic variable improves prediction of the third moment (proportional to bulk rain mass) and rain rate. In contrast, increasing process rate formulation complexity by adding more power terms has little benefit—a result that is explained using further-idealized experiments. BOSS rainshaft simulations are shown to well estimate the true process rates from constraint by bulk rain observations, with the further benefit of rigorously quantified uncertainty of these estimates.

58 GEOSCIENCES↗

DOE OSTI2024

Estimating scalar turbulent fluxes with slow-response sensors in the stable atmospheric boundary layer

Conventional and recently developed approaches for estimating turbulent scalar fluxes under stable atmospheric conditions are evaluated, with a focus on gases for which fast sensors are not readily available. First, the relaxed eddy accumulation (REA) classical approach and a recently proposed mixing length parameterization, labeled A22, are tested against eddy-covariance computations. Using high-frequency measurements collected from two contrasting sites (the frozen tundra near Utqiaġvik, Alaska, and a sparsely vegetated grassland in Wendell, Idaho, during winter), it is shown that the REA and A22 models outperform the conventional Monin–Obukhov similarity theory (MOST) utilized widely to infer fluxes from mean gradients. Second, scenarios where slow trace gas sensors are the only viable option in field measurements are investigated using digital filtering applied to fast-response sensors to simulate their slow-response counterparts. With a filtered scalar signal, the observed filtered eddy-covariance fluxes are referred to here as large-eddy-covariance (LEC) fluxes. A virtual eddy accumulation (VEA) approach, akin to the REA model but not requiring a mechanical apparatus to separate the gas flows, is also formulated and tested. A22 outperforms VEA and LEC in predicting the observed unfiltered (total) eddy-covariance (EC) fluxes; however, VEA can still capture the LEC fluxes well. This finding motivates the introduction of a sensor response time correction into the VEA formulation to offset the effect of sensor filtering on the underestimated net averaged fluxes. The only needed parameter for this correction is the mean velocity at the instrument height, a surrogate of the advective timescale. The VEA approach is very suitable and simple to use with gas sensors of intermediate speed (∼ 0.5 to 1 Hz) and with conventional open- or closed-path setups.

58 GEOSCIENCES↗

DOE OSTI2024

Recalibration of missing low-frequency variability and trends in the North Atlantic Oscillation

Abstract Multi-decadal trends in the wintertime North Atlantic Oscillation (NAO) are under-represented by coupled general circulation models (CGCMs), consistent with a lack of autocorrelation in their NAO index series. This study proposes and tests two simple “reddening” approaches for correcting this problem in simulated indices based on simple one parameter short-term (AR; Auto-Regressive order 1) and long-term (FD; Fractional-Difference) time series filters. Using CGCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6), the FD filter successfully improves the autocorrelation structure of the NAO, and in turn the simulation of extreme trends, while the AR filter is less successful. The 1963–1993 NAO trend is the maximum 31-year trend in the historical period. Raw CGCMs underestimate the likelihood of this trend by a factor of ten but this discrepancy is corrected after reddening. CMIP6 future projections show that long-term (2024–2094) NAO ensemble mean trends systematically increase with the magnitude of radiative forcing: -2.4 to 3.5 hPa/century for low-to-high forcing after reddening (more than double the range using raw output). The related likelihood of future maximum 31year trends comparable to 1963–1993 ranges from 3 to 7% whereas none of these CMIP6 projections simulate this without reddening. Near-term projections of the next 31 years (2024–2054) are less sensitive than long term trends to the future scenario, showing weak-to-no forced trend. However, reddening increases the ensemble range by 74% (to +/-1 standard deviation/decade), which could increase/decrease regional climate change signals in the Northern Hemisphere by magnitudes that are underestimated when using raw CGCM output.

Meteorology & Atmospheric Sciences↗

DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗