Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Cloud Correspondance”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

Trifluoroacetic Acid from Degradation of HCFCs and HFCs: A Three-dimensional Modeling Study

Trifluoroacetic acid (TFA; CF3 COOH) is produced by the degradation of the halocarbon replacements HFC-134a, HCFC-124, and HCFC-123. The formation of TFA occurs by HFC/HCFC reacting with OH to yield CF3COX (X = F or CI), followed by in-cloud hydrolysis of CF3COX to form TFA. The TFA formed in the clouds may be reevaporated but is finally deposited onto the surface by washout or dry deposition. Concern has been expressed about the possible long-term accumulation of TFA in certain aquatic environments, pointing to the need to obtain information on the concentrations of TFA in rainwater over scales ranging from local to continental. Based on projected concentrations for HFC-134a, HCFC-124, and HCFC-123 of 80, 10, and 1 pptv in the year 2010, mass conservation arguments imply an annually averaged global concentration of 0.16 microg/L if washout were the only removal mechanism for TFA. We present 3-D simulations of the HFC/HCFC precursors of TFA that include the rates of formation and deposition of TFA based on assumed future emissions. An established (GISS/Harvard/ UCI) but coarse-resolution (8 deg latitude by 10 deg longitude) chemical transport model was used. The anually averaged rainwater concentration of 0.12 gg/L (global) was calculated for the year 2010, when both washout and dry deposition are included as the loss mechanism for TFA from the atmosphere. For some large regions in midnorthern latitudes, values are larger, 0.15-0.20 microg/L. The highest monthly averaged rainwater concentrations of TFA for northern midlatitudes were calculated for the month of July, corresponding to 0.3-0.45 microg/L in parts of North America and Europe. Recent laboratory experiments have suggested that a substantial amount of vibrationally excited CF3CHFO is produced in the degradation of HFC-134a, decreasing the yield of TFA from this compound by 60%. This decrease would reduce the calculated amounts of TFA in rainwater in the year 2010 by 26%, for the same projected concentrations of precursors.

Kotamarthi, V. R.↗

Trifluoroacetic Acid from Degradation of HCFCs and HFCs: A Three-Dimensional Modeling Study

Trifluoroacetic acid (TFA; CF3COOH) is produced by the degradation of the halocarbon replacements HFC-134a, HCFC-124, and HCFC-123. The formation of TFA occurs by HFC/HCFC reacting with OH to yield CF3COX (X = F or CI), followed by in-cloud hydrolysis of CF3COX to form TFA. The TFA formed in the clouds may be reevaporated but is finally deposited onto the surface by washout or dry deposition. Concern has been expressed about the possible long-term accumulation of TFA in certain aquatic environments, pointing to the need to obtain information on the concentrations of TFA in rainwater over scales ranging from local to continental. Based on projected concentrations for HFC-134a, HCFC-124, and HCFC-123 of 80, 10, and 1 pptv in the year 2010, mass conservation arguments imply an annually averaged global concentration of 0.16 micro g/L if washout were the only removal mechanism for TFA. We present 3-D simulations of the HFC/HCFC precursors of TFA that include the rates of formation and deposition of TFA based on assumed future emissions. An established (GISS/Harvard/ UCI) but coarse-resolution (8 deg latitude by 10 deg longitude) chemical transport model was used. The annually averaged rainwater concentration of 0.12 micro g/L (global) was calculated for the year 2010, when both washout and dry deposition are included as the loss mechanism for TFA from the atmosphere. For some large regions in midnorthern latitudes, values are larger, 0.15-0.20 micro g/L. The highest monthly averaged rainwater concentrations of TFA for northern midlatitudes were calculated for the month of July, corresponding to 0.3-0.45 micro g/L in parts of North America and Europe. Recent laboratory experiments have suggested that a substantial amount of vibrationally excited CF3CHFO is produced in the degradation of HFC-134a, decreasing the yield of TFA from this compound by 60%. This decrease would reduce the calculated amounts of TFA in rainwater in the year 2010 by 26%, for the same projected concentrations of precursors.

Kotamarthi, V. R.↗

Trifluoroacetic Acid from Degradation of HCFCs and HFCs: A Three-Dimensional Modeling Study

Trifluoroacetic acid (TFA; CF3 COOH) is produced by the degradation of the halocarbon replacements HFC-134a, HCFC-124, and HCFC-123. The formation of TFA occurs by HFC/HCFC reacting with OH to yield CF3COX (X = F or CI), followed by in-cloud hydrolysis of CF to form TFA. The TFA formed in the clouds may be reevaporated but is finally deposited onto the surface by washout or dry deposition. Concern has been expressed about the possible long-term accumulation of TFA in certain aquatic environments, pointing to the need to obtain information on the concentrations of TFA in rainwater over scales ranging from local to continental. Based on projected concentrations for HFC-134a, HCFC-124, and HCFC-123 of 80, 10, and 1 pptv in the year 2010, mass conservation arguments imply an annually averaged global concentration of 0.16 micro g/L if washout were the only removal mechanism for TFA. We present 3-D simulations of the HFC/HCFC precursors of TFA that include the rates of formation and deposition of TFA based on assumed future emissions. An established (GISS[Harvard/ UCI) but coarse-resolution (8 deg latitude by 10 deg longitude) chemical transport model was used. The annually averaged rainwater concentration of 0.12 micro g/L (global) was calculated for the year 2010, when both washout and dry deposition are included as the loss mechanism for TFA from the atmosphere. For some large regions in midnorthern latitudes, values are larger. 0.15-0.20 micro g/L. The highest monthly averaged rainwater concentrations of TFA for northern midlatitudes were calculated for the month of July, corresponding to 0.3 - 0.45 micro g/L in parts of North America and Europe. Recent laboratory experiments have suggested that a substantial amount of vibrationally excited CF3CHFO is produced in the degradation of HFC-134a, decreasing the yield of TFA from this compound by 60%. This decrease would reduce the calculated amounts of TFA in rainwater in the year 2010 by 26%, for the same projected concentrations of precursors.

Kotamarthi, V. R.↗

Raindrop Size Spectrum in Deep Convective Regions of the Americas

This study compared drop size distribution (DSD) measurements on the surfaces, the corresponding properties, and the precipitation modes among three deep convective regions within the Americas. The measurement compilation corresponded to two sites in the midlatitudes: the U.S. Southern Great Plains and Córdoba Province in subtropical South America, as well as to one site in the tropics: Manacapuru in central Amazonia; these are all areas where intense rain-producing systems contribute to the majority of rainfall in the Americas’ largest river basins. This compilation included two types of disdrometers (Parsivel and 2D-Video Disdrometer) that were used at the midlatitude sites and one type of disdrometer (Parsivel) that was deployed at the tropical site. The distributions of physical parameters (such as rain rate R, mass-weighted mean diameter D m , and normalized droplet concentration N w ) for the raindrop spectra without rainfall mode classification seemed similar, except for the much broader N w distributions in Córdoba. The raindrop spectra were then classified into a light precipitation mode and a precipitation mode by using a cutoff at 0.5 mm h -1 based on previous studies that characterized the full drop size spectra. These segregated rain modes are potentially unique relative to previously studied terrain-influenced sites. In the light precipitation and precipitation modes, the dominant higher frequency observed in a broad distribution of N w in both types of disdrometers and the identification of shallow light precipitation in vertically pointing cloud radar data represent unique characteristics of the Córdoba site relative to the others. As a result, the co-variability between the physical parameters of the DSD indicates that the precipitation observed in Córdoba may confound existing methods of determining the rain type by using the drop size distribution.

54 ENVIRONMENTAL SCIENCES↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Variability of Water Vapor in the Tropical Middle Atmosphere Observed From Satellites and Interpreted Using SD-WACCM Simulations

Water vapor in the middle atmosphere plays an essential role in global warming, ozone depletion, and the formation of polar stratospheric and mesospheric clouds. We show that tropical middle atmospheric water vapor simulated with the specified-dynamics version of the Whole Atmosphere Community Climate Model (SD-WACCM) is consistent with changes observed in a merged satellite data set, which encompasses the period 1993–2020. Consistent with previous work, we find no significant trend in the stratosphere in either the observations or the simulation; in the mesosphere, we find a long-term trend of 0.1 ppmv per decade, but only in the observations. We also analyze an SD-WACCM simulation for the longer period 1980–2019 to quantify the contribution of various factors to the decadal variation of middle atmospheric water vapor. Over 1980–1995, the simulated water vapor in the upper stratosphere and mesosphere, averaged zonally and over ±30° latitude, increases by 0.30 ppmv per decade due to increasing methane emissions. After 1995, a significant abrupt decrease of water vapor of 0.37 ppmv per decade and then a gradual increase of 0.33 ppmv per decade result from changes in stratospheric cold point temperature. The cold-point temperature is strongly influenced by the strength of the Brewer-Dobson circulation. The acceleration of the Brewer-Dobson circulation before about 2003 leads to a cooler tropical tropopause and a decrease of water vapor, and the deceleration thereafter leads to corresponding warming of the tropopause and an increase in water vapor.

54 ENVIRONMENTAL SCIENCES↗

Variability of Water Vapor in the Tropical Middle Atmosphere Observed From Satellites and Interpreted Using SD-WACCM Simulations

Water vapor in the middle atmosphere plays an essential role in global warming, ozone depletion, and the formation of polar stratospheric and mesospheric clouds. We show that tropical middle atmospheric water vapor simulated with the specified-dynamics version of the Whole Atmosphere Community Climate Model (SD-WACCM) is consistent with changes observed in a merged satellite data set, which encompasses the period 1993–2020. Consistent with previous work, we find no significant trend in the stratosphere in either the observations or the simulation; in the mesosphere, we find a long-term trend of 0.1 ppmv per decade, but only in the observations. We also analyze an SD-WACCM simulation for the longer period 1980–2019 to quantify the contribution of various factors to the decadal variation of middle atmospheric water vapor. Over 1980–1995, the simulated water vapor in the upper stratosphere and mesosphere, averaged zonally and over ±30° latitude, increases by 0.30 ppmv per decade due to increasing methane emissions. After 1995, a significant abrupt decrease of water vapor of 0.37 ppmv per decade and then a gradual increase of 0.33 ppmv per decade result from changes in stratospheric cold point temperature. The cold-point temperature is strongly influenced by the strength of the Brewer-Dobson circulation. The acceleration of the Brewer-Dobson circulation before about 2003 leads to a cooler tropical tropopause and a decrease of water vapor, and the deceleration thereafter leads to corresponding warming of the tropopause and an increase in water vapor.

Wandi Yu↗

Cloud, dust, and ozone vertical profiles from solar occultation measurements: Implications for dynamics

An instrument was designed for solar occultation measurements of the martian atmosphere from the Phobos spacecraft. It was composed of three different dispersive systems working in the ultraviolet (UV: 0.22-0.32 microns) for the measurement of O3 and aerosols, in the near infrared (NIR: 0.76 microns, 0.94 microns) for the detection of O2 and H2O, and in the infrared (IR: 1.9 microns, 3.7 microns) where CO2 and H2O were measured. A detailed description of the instrument may be found in the special issue of Nature. Its principle objective is to measure from the Phobos orbit the spectrum of the Sun, modified by atmospheric extinction, during sunset. The UV-NIR spot has an angular diameter of 1 arcmin, or approximately 3 km vertical resolution, and is located near the center of the solar disk. The IR field is about twice as large and its line of sight is shifted by 8.5 arcmin, or approximately 20 km, relative to the previous one. It is therefore located near the edge of the solar disk. Sampling times are generally 0.5, 1, and 2 s for IR, UV, and NIR channels respectively, corresponding to vertical excursions of the line of sight of 1, 2, and 4 km respectively under nominal conditions. The instrument operated from February 8 to March 26 (the martian equinox occurred on February 17). The latitude of the intersection of the Sun-spacecraft axis with the surface of Mars varied from -11 to +20, the seasonal date L(sub s) being in the range 0 - 20. All measurements were therefore made near northern spring equinox in equatorial regions. Due to an error in the pointing system, only partial results were obtained, the region below approximately equals 30 km altitude being never sounded by the UV-NIR spectrometer. On the contrary, nine complete occultations were obtained in the IR channels, whose line of sight was fortunately approximately equals 20 km below the UV-NIR axis.

Chassefiere, E.↗

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗

Rocket Lab Venus - Enabling Low-Cost Interplanetary Missions

Rocket Lab’s mission to Venus [1], launching in January 2025, aims to demonstrate that small launch vehicles, such as Electron, and high-energy small spacecraft, such as Photon, can enable a new paradigm of regular, low-cost interplanetary missions. The Electron launch vehicle will deliver the Photon spacecraft carrying the entry probe into low earth orbit, at which point Photon will take over and put both itself and the attached probe on an Earth-Venus transfer trajectory, with assistance from a Lunar flyby. The primary science objective of the mission is to sample the Venusian cloud layers between atitudes of 45 and 60 km above the surface with an autofluorescing nephelometer. The ~1kg scientific instrument will be carried by a small ~17 kg direct entry probe which will separate from the Photon spacecraft approximately 30 minutes prior to atmospheric entry. The probe outer mould line is a scaled down version of the Deep Space 2 probes [2]. NASA’s newly developed 3D Woven Carbon Phenolic, a derivative of the Heatshield for Extreme Entry Environment Technology (HEEET) material [3], is an enabling technology for this mission, and it will comprise the probe forebody heat shield.The probe will communicate directly back to the Earth during its descent via a small on-board radio frequency (RF) antenna. The hyperbolic Venus encounter will see the probe enter the Venusia atmosphere at approximately 11 km/s with an entry flight path angle (EFPA) of about 10°. This relatively low entry angle was selected to increase the duration of the descent phase, reduce the peak deceleration, and reduce the sensitivity of the nominal entry environments to navigational errors. For reference, the Pioneer Venus probes entered at EFPAs varying from about 25° to 70°, resulting in peak decelerations between around 220 and 450 g’s and peak stagnation point heat fluxes ranging from about 40 to 75 MW/m2 [4]. In comparison, the Rocket Lab Venus probe will experience a lower preak deceleration of around 60 g’s and a lower peak stagnation point heat flux of about 15 MW/m2, at the cost of a greater total heat load and correspondingly higher thermal protection system (TPS) mass fraction. This paper recounts the development history of the TPS for the Rocket Lab Venus probe, with a focus on the unique requirements of the mission. In particular, the effect on the TPS design of the mission’s tight budget, short development timeline, slim mass and RF link margins, and mild but long duration entry heating. Rocket Lab’s initial investigations into a ‘hybrid’ TPS [5] based on a scaled version of the Pioneer Venus full-density carbon phenolic forebody shield, and the abandonment of this approach in favour of NASA’s 3D Woven Carbon Phenolic is presented. Additionally, the leveraging of previous testing under the HEEET project at higher conditions and the associated TPS margining strategy that was adopted to avoid the need for extensive mission-specific arc jet testing is discussed. Finally, this paper concludes with a brief consideration of the feasibility of a future low-cost aerocapture technology demonstrator mission based on the 2025 Rocket Lab Venus mission and entry probe.

Lyle Campbell↗

Improving our Understanding of Atlantic Tropical Cyclones through Knowledge of the Saharan Air Layer: Hope or Hype?

The existence of the Saharan air layer (SAL), a layer of warm, dry, dusty air that frequently moves westward off of the Saharan desert of Africa and over the tropical Atlantic Ocean, has long been appreciated. As air moves over the desert, it is strongly heated from below, producing a very hot air mass at low levels. Because there is no moisture source over the Sahara, the rise in temperature causes a sharp drop in relative humidity, thus drying the air. In addition, the warm air produces a very strong jet of easterly flow in the middle troposphere called the African easterly jet that is thought to play a critical role in hurricane formation. In recent years, there has been an increased focus on the impact that the SAL has on the formation and evolution of hurricanes in the Atlantic. However, the nature of its impact remains unclear, with some researchers arguing that the SAL amplifies hurricane development and with others arguing that it inhibits it. The argument for positively influencing hurricane development is based upon the fact that the African easterly jet produces the waves that eventually form hurricanes and that it leads to rising motion south of the jet that favors the development of deep thunderstorm clouds. The potential negative impacts of the SAL include 1) low-level vertical wind shear associated with the African easterly jet; 2) warm SAL air aloft, which increases thermodynamic stability and suppresses cloud development; and 3) dry air, which produces cold downdrafts in precipitating regions, thereby removing energy needed for storm development. As part of this recent focus on the SAL and hurricanes (which motivated a 2006 NASA field experiment), there has been little emphasis on the SAL s potential positive influences and almost complete emphasis on its possible negative influences, almost to the point of claims that the SAL is the major suppressing influence on hurricanes in the Atlantic. Multiple NASA satellite data sets (TRMM, MODIS, and AIRS/AMSU) and National Centers for Environmental Prediction global analyses are used to characterize the SAL s properties and evolution in relation to developing hurricanes. The results show that storms generally form on the southern side of the jet, where favorable background rotation is high. The jet often helps to form the northern side of the storms and rarely moves over their inner cores, so jet-induced vertical wind shear does not appear to be a negative influence on developing storms. Warm SAL air is confined to regions north of the jet and generally does not impact the tropical cyclone precipitation south of the jet. Of the three proposed negative influences, dry air appears to be the key influence; however, the presence of dry SAL air is not a good indicator of whether a storm will weaken since many examples of intensifying storms surrounded by such dry air can be found. In addition, a global view of relative humidity shows moisture distributions in other ocean basins that are almost identical to the Atlantic. The dry zones correspond to regions of descending air on the eastern and equatorward sides of semi-permanent oceanic high pressure systems. Thus, the dry air over the Atlantic appears to be primarily a product of the large-scale flow, but with enhanced drying at low levels associated with the Sahara. As a result, we conclude that the SAL is not a major negative influence on hurricanes. It is just one of many possible influences and can be both positive and negative.

Braun, Scott A.↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Unique Observational Constraints on the Seasonal and Longitudinal Variability of the Earth’s Planetary Albedo and Cloud Distribution Inferred From EPIC Measurements

Thorough comparison to observations is key to developing a credible climate model forecasting capability. Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation provide a unique observational perspective that permits a more reliable model/data comparison than is possible with the otherwise available satellite data. The uniqueness is in the DSCOVR satellite’s viewing geometry, which enables continuous viewing of the Earth’s sunlit hemisphere from its Lissajous orbit around the Lagrangian L1 point. The key instrument is the Earth Polychromatic Imaging Camera (EPIC), which views the Earth’s sunlit hemisphere with 1024-by-1024-pixel imagery in 10 narrow spectral bands from 317 to 780 nm, acquiring up to 22 high spatial resolution images per day. The additional feature is that the frequency of EPIC image acquisition is nearly identical to that of the climate GCM data generation scheme where climate data for the entire globe are ‘instantaneously’ calculated at 1-h radiation time-step intervals. Implementation of the SHS (Sunlit Hemisphere Sampling) EPIC-view geometry for the in-line GCM output data sampling establishes a precise self-consistency in the space-time data sampling between EPIC observational and GCM output data generation and sampling. The remaining problem is that the GCM generated data are radiative fluxes, while the EPIC measurements are backscatter-dependent radiances. Radiance to flux conversion is a complex problem with no simple way to convert GCM radiative fluxes into spectral radiances. The more expedient approach is to convert the EPIC spectral radiances into broadband radiances by MODIS/CERES-based regression relationships and then into solar radiative fluxes using the CERES angular distribution models. Averaging over the sunlit hemisphere suppresses the meteorological weather noise, but preserves the intra-seasonal larger scale variability. Longitudinal slicing by the Earth’s rotation permits a self-consistent model/data comparison of the longitudinal model/data differences in the variability of the reflected solar radiation. Ancillary EPIC Composite data provide additional cloud property information for climate model diagnostics. Comparison of EPIC-derived seasonal and longitudinal variability of the Earth’s planetary albedo with the GISS ModelE2 results shows systematic overestimate of cloud reflectivity over the Pacific Ocean with corresponding underestimates over continental land areas.

DSCOVR↗

Frozen Hydrometeor Terminal Fall Velocity Dependence on Particle Habit and Riming as Observed by Vertically Pointing Radars

Vertically pointing Ka-band radar measurements are used to derive fall velocity–reflectivity factor ($V$ t = $aZ$$^{b}_{e}$) relations for frozen hydrometeor populations of different habits during snowfall events observed at Oliktok Point, Alaska, and at the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC). Case study events range from snowfall with highly rimed particles observed during periods with large amounts of supercooled liquid water path (LWP > 320 g m –2 ) to unrimed snowflakes including instances when pristine planar crystals were the dominant frozen hydrometeor habit. The prefactor a and the exponent b in the observed $V$ t –$Z$ e relations scaled to the sea level vary in the approximate ranges 0.5–1.4 and 0.03–0.13, respectively (reflectivities are in mm 6 m –3 and velocities are in m s –1 ). The coefficient a values are the smallest for planar crystals (a ~ 0.5) and the largest (a > 1.2) for particles under severe riming conditions with high LWP. There is no clear distinction between b values for high and low LWP conditions. The range of the observed $V$ t –$Z$ e relation coefficients is in general agreement with results of modeling using fall velocity–size (υ t = αD β ) relations for individual particles found in literature for hydrometeors of different habits, though there is significant variability in α and β coefficients from different studies even for a same particle habit. Correspondences among coefficients in the $V$ t –$Z$ e relations for particle populations and in the individual particle υ t –$D$ relations are analyzed. Furthermore, these correspondences and the observed $V$ t –$Z$ e relations can be used for evaluating different frozen hydrometeor fall velocity parameterizations in models.

54 ENVIRONMENTAL SCIENCES↗

The Transition of High-Resolution NASA MODIS Sea Surface Temperatures into the WRF Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composite at 2-km resolution that has been implemented in version 3 of the National Weather Service (NWS) Weather Research and Forecasting (WRF) Environmental Modeling System (EMS). The WRF EMS is a complete, full physics numerical weather prediction package that incorporates dynamical cores from both the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). The installation, configuration, and execution of either the ARW or NMM models is greatly simplified by the WRF EMS to encourage its use by NWS Weather Forecast Offices (WFOs) and the university community. The WRF EMS is easy to run on most Linux workstations and clusters without the need for compilers. Version 3 of the WRF EMS contains the most recent public release of the WRF-NMM and ARW modeling system (version 3 of the ARW is described in Skamarock et al. 2008), the WRF Pre-processing System (WPS) utilities, and the WRF Post-Processing program. The system is developed and maintained by the NWS National Science Operations Officer Science and Training Resource Coordinator. To initialize the WRF EMS with high-resolution MODIS SSTs, SPoRT developed the composite product consisting of MODIS SSTs over oceans and large lakes with the NCEP Real-Time Global (RTG) filling data over land points. Filling the land points is required due to minor inconsistencies between the WRF land-sea mask and that used to generate the MODIS SST composites. This methodology ensures a continuous field that adequately initializes all appropriate arrays in WRF. MODIS composites covering the Gulf of Mexico, western Atlantic Ocean and the Caribbean are generated daily at 0400, 0700, 1600, and 1900 UTC corresponding to overpass times of the NASA Aqua and Terra polar orbiting satellites. The MODIS SST product is output in gridded binary-1 (GRIB-1) data format for a seamless incorporation into WRF via the WPS utilities. The full-resolution, 1-km MODIS product is sub-sampled to 2-km grid spacing due to limitations in handling very large dimensions in the GRIB-1 data format. The GRIB-1 files are posted online at ftp://ftp.nsstc.org/sstcomp/WRF/, which is directly accessed by the WRF EMS scripts. The MODIS SST composites are also downloaded to the EMS data server, which is accessible by the WRF EMS users and NWS WFOs. The SPoRT MODIS SST composite provides the model with superior detail of the ocean gradients around Florida and surrounding waters, whereas the operational RTG SST typically depicts a relatively smooth field and is not able to capture sharp horizontal gradients in SST. Differences of 2-3 C are common over small horizontal distances, leading to enhanced SST gradients on either side of the Gulf Stream and along the edges of the cooler shelf waters. These sharper gradients can in turn produce atmospheric responses in simulated temperature and wind fields as depicted in LaCasse et al. Differences in atmospheric verification statistics over a several month study were generally small in the vicinity of south Florida; however, the validation of SSTs at specific buoy locations revealed important improvements in the biases and RMS errors, especially in the vicinity of the cooler shelf waters off the east-central Florida coast. A current weakness in the MODIS SST product is the occurrence of occasional discontinuities caused by high latency in SST coverage due to persistent cloud cover. An enhanced method developed by Jedlovec et al. (2009, GHRSST User Symposium) reduces the occurrence of these problems by adding Advanced Microwave Scanning Radiometer -- EOS (AMSR-E) SST data to the compositing process. Enhanced SST composites are produced over the ocean regions surrounding the Continental U.S. at four times each day corresponding to Terra and Aqua equator crossing times. For a given day and overpass time, both MODInd AMSR-E data from the previous seven days form a collection used in the compositing. At each MODIS pixel, cloud-free SST values from the collection are used to form a weighted average based on their latency (number of days from the current day). In this way, recent SST data are given more weight than older data. One of the primary issues involved in incorporating the AMSR-E microwave data in the composites is the tradeoff between the decreased spatial resolution of the AMSR-E data (25 km) and the increased coverage due to its near all-weather capability. Currently, the AMSR-E is given a weight of 20% compared to MODIS data, thereby preserving the spatial structure observed in the MODIS data. Day-time (night-time) AMSR-E SST data from Aqua are used with both Terra and Aqua MODIS day-time (night-time) SST data sets.

Case, Jonathan L.↗

Monitoring Mars with the Hubble Space Telescope: 1990-1991 observations

Images of Mars obtained using the Planetary Camera on the Hubble Space Telescope during 1990 and 1991 are described and analyzed. Multispectral images of Mars record the martian season between L(sub s) = 348 deg and 60 deg corresponding to late winter and spring in the northern hemisphere. The wavelengths of these observations varied between 230 nm in the ultraviolet and 890 nm in the near infrared. We use these images to discuss atmospheric and polar phenomena, to constrain the aerosols and ozone in the martian atmosphere, and to compare surface albedo features within this data set as well as with previous observations. Two methods are used to constrain the opacity of the martian atmosphere: comparisons of ultraviolet images with atmospheric scattering models and examination of contrast changes of surface features observed at visible wavelengths. Our observations are consistent with a Mars which was quite different from the planet when it was viewed by Viking; we estimate an upper limit of 0.1 for equatorial dust opacity near vernal equinox, compared to values of 0.4 to 0.5 measured in Chryse by Viking Lander 1. Later, at L(sub s) = 60 deg, and upper limit of 0.2 was found for dust in the southern hemisphere. Except for suggestions of local dust activity near the north polar cap in late spring, there was no visual evidence of dust storms on the planet. The Hartley band of ozone is within the F230W filter bandpass, and the 230W/336W ratio is therefore sensitive to ozone. High latitude ozone is observed at a concentration similar to that observed by Mariner 9; a seasonal dependence in ozone concentration was observed with concentration decreasing at northern arctic latitudes and increasing at southern midlatitudes as L(sub s) increases from 350 deg to 60 deg. Despite the difference in dust loading, the behavior of condensate clouds and the recession of the north polar cap are consistent with those observed during Viking years. Late winter images of the north polar region reveal an active north polar hood with a maximum (UV) opacity of about 0.25 near 50 deg latitude; the hood seems to have an annular form, with the central, polar regions having smaller opacity. Clouds observed in the northern hemisphere to the south of the hood boundary are suggestive of fronts similar to those observed during the Mariner 9 mission. A diurnal variation in the hood clouds which is (anti) correlated with the amount of water vapor in the atmosphere is noted. Images were deconvolved using the Richardson-Lucy iterative algorithm to remove the defocussing effects of the spherical aberration of the HST primary. Even when Mars subtended less than 5 arcsec, in May 1991, the images were scientifically useful and comparable to photographic images obtained at opposition. Although the data are subject to the solar pointing constraint for HST, the relatively slow movement of the allowed observation window, which is more than one half martian year, with respect to the martian seasons will permit us to address the question of major variations in the martian atmosphere.

James, Philip B.↗

Combined Mesonet and Tracker

Title: Combined Mesonet and Trackers (UNL Mobile Mesonets) Authors University of Nebraska PI: Adam Houston, UNL Professor (ahouston2@unl.edu) Mailing Address: 126 Bessey Hall P.O. Box 880340 Lincoln, NE 68588-0340 CoMeT Overview The University of Nebraska-Lincoln operates three Combined Mesonet and Tracker (CoMeTs). CoMeTs are Ford Explorers (model years 2015, 2017, and 2019) with forward-mounted suites of meteorological sensors and dual moonroofs, combining the capability of a mobile mesonet to collect near-surface observations with the capability of an unmanned aircraft systems (UAS) tracker vehicle, which enables an observer in the second row of seats to see the aircraft and maintain compliance with Federal Aviation Administration policies on UAS operation. The CoMeTs collect observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155A, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, and vehicle heading using a KVH Industries C-100 fluxgate compass (Barbieri et al. 2019). The HMP155A and 109SS-L thermistor are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). This list of sensors is also included in the CoMeT data file metadata. Manufacturer specifications for these instruments are given in Table 1 of Hanft and Houston (2018). The reported measured quantities are summarized below.CoMeT-3 was funded through an equipment allocation included in the NSF TORUS award (AGS-1824649). Instrument Description The specific sensors included on each CoMeT are summarized in the table at the end of this section. In general each CoMeT collects observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210 barometer with a Gill pressure port, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, position using a Garmin 19x HVS receiver, and vehicle heading using a KVH Industries C-100 fluxgate compass. The HMP155 and 109SS are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). Fast temperature and corrected RH measurements (using sensors housed within the U-tube) have a time constant of 10-12 s based on data collected across a temperature and RH shock during the CLOUD-MAP 2017 calibration/validation tests on June 26, 2017. Vehicle speed was < 10 kts for this test. CoMeT-1 CoMeT-2 CoMeT-3 Slow Temperature Slow RH Vaisala HMP155A-L20-PT Part #: 22280-7 Vaisala HMP155E Part #: E1AA11A0B1A1A0A Vaisala HMP155E Part #: E1AA11A0B1A1A0A Fast temperature Campbell Scientific 109SS-L20-PT Part #: 21448-3 Campbell Scientific 109SS-L12-PW Part #: 21448-109 Campbell Scientific 109SS-L12-PW Part #: 21448-150 Pressure Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Wind RM Young 05103-L20-PT Part #: 18435-310 RM Young 05103-L20-PW Part #: 18435-244 RM Young 05103-L20-PW Part #: 18435-244 GPS Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Compass KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 Logger Campbell Scientific CR6-NA-XT-SW Part #: 28385-9 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Data Collection and Real-Time Processing The reported measured quantities are summarized in the table below. Quantity Units Source Epoch time Seconds GPS Latitude and longitude Degrees GPS Altitude m GPS Pressure hPa PTB210 Temperature (fast) deg C 109SS-L Temperature (slow) deg C HMP155 RH (slow) % HMP155 Vehicle speed m/s GPS Vehicle heading deg C-100 and GPS In addition to the measured variables, several derived variables are calculated. Corrected/fast relative humidity (%) Relative humidity is adjusted to the fast temperature following Richardson et al. (1998) and Houston et al. (2016). Water vapor mixing ratio (g/kg) Dew point temperature (&deg;C) Potential temperature (Kelvin) Virtual potential temperature (Kelvin) Equivalent potential temperature (Kelvin) Regular intercomparisons between all three CoMeTs were performed during TORUS 2019. Comparisons were also conducted between CoMeT-1 and CoMeT-2 during LAPSE-RATE (2018) on 14 July. In these intercomparisons, the vehicles were parked adjacent to each other aligned perpendicular to (and facing into) the wind. To minimize engine heating effects, intercomparisons were only conducted when the wind speed was >10 kts. Data Format Original data files for each deployment are saved as text files and then converted to NetCDF. NetCDF versions have units that are CF compliant and may not match the original units in the txt files. The naming convention for the NetCDF files is as follows: UNL.CoMeT3.{deployment date YYYYMMDD}.{start time of observation collection in UTC HHMM}.L2.{post-processing codes}.cdf example: UNL.CoMeT3.20190627.1931.L2.g1.f1.cdf Post-processing codes are included to track modifications to the raw data. These codes are closely connected to error flags associated with each record. Each letter corresponds to a particular instrument: g: GPS p: Barometer tf: Fast temperature ts: Slow temperature rh: Relative humidity f: Compass w: Wind monitor a: All instruments Each number corresponds to a particular post-processing action described more below. Measured and derived variables are included in the following table. Variable Heading Standard Name Units time Time seconds since 00:00:00, 01-01-1970 Alt Altitude meters lat Latitude degrees north lon Longitude degrees east fast_temp Air Temperature kelvin slow_temp Air Temperature kelvin pressure Air Pressure pascals logger_RH Relative Humidity percent calc_corr_RH Relative Humidity percent wind_speed Wind Speed meters per second wind_dir Wind From Direction degrees vehicle_dir Vehicle Direction degrees dewpoint Dew Point Temperature kelvin mixing_ratio Humidity Mixing Ratio g/g theta Air Potential Temperature kelvin theta_v Virtual Potential Temperature kelvin theta_e Equivalent Potential Temperature Kelvin error_flag The error_flag variable is a string that matches the post-processing codes listed above. All instruments will have an associated code, but will have a &ldquo;0&rdquo; if the datum is unchanged from the initial processed value. Error Codes The following table summarizes the error codes for data collected before 2020: Error Code Relevant CoMeT Description g1 1,2,3 Exact correction. GPS position and time reprocessed from raw data g2 1 As far as we can tell this is an exact correction to an error in the GPS time. During the correct time periods the time suddenly went backwards ~250s and stayed at this offset for 750s when it corrected itself. The offset was applied to the &ldquo;time warp&rdquo; period. p1 2 Approximate correction. Hole in the pressure tube connecting the pressure port to the barometer. Resulted in erroneously low air pressure measurements when the vehicle was in motion. Derived variables recalculated (dew point temperature [e depends on qv and p], water vapor mixing ratio, potential temperature, virtual potential temperature, equivalent potential temperature) a1 3 Exact correction. Missing data reprocessed from raw data a2 1 Bug fix to bias correction for ts1, ts2, and rh1: water vapor mixing ratio was off by a factor of 10 and virtual potential temperature was wrong because of this. f1 3 No correction, missing data. Fluxgate compass inoperable. Wind speed and direction calculated using GPS-derived vehicle heading instead. rh1 1 Approximate correction. Constant bias of +1.7% removed from relative humidity. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts1 1 Approximate correction. Constant bias of +0.6 K removed from slow temperature. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts2 1 Approximate correction. Constant bias of +1.0 K removed. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) References Bolton, D., 1980: The Computation of Equivalent Potential Temperature. Mon. Wea. Rev., 108, 1046&ndash;1053, https://doi.org/10.1175/1520-0493(1980)108<1046:TCOEPT>2.0.CO;2. Hanft, W., and A. L. Houston, 2018: An Observational and Modeling Study of Mesoscale Air Masses with High Theta-E. Mon. Wea. Rev., 146, 2503&ndash;2524, https://doi.org/10.1175/MWR-D-17-0389.1.Wexler Houston, A. L., R. J. Laurence III, T. W. Nichols, S. Waugh, B. Argrow, and C. L. Ziegler, 2016: Intercomparison of unmanned aircraft-borne and mobile mesonet atmospheric sensors. Journal of Atmospheric and Oceanic Technology. 33, 1569-1582, doi: 10.1175/JTECH-D-15-0178.1. Lowe, P. R., 1977: An Approximating Polynomial for the Computation of Saturation Vapor Pressure. J. Applied Meteorology, 16, 100&ndash;103. Richardson, S. J., S. E. Frederickson, F. V. Brock, and J. A. Brotzge, 1998: Combination temperature and relative humidity probes: Avoiding large air temperature errors and associated relative humidity errors. Preprints, 10th Symp. On Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 278&ndash;283. Waugh, S., and S. E. Frederickson, 2010: An improved aspirated temperature system for mobile meteorological observations, especially in severe weather. 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., P5.2. [Available online at https://ams.confex.com/ams/25SLS/techprogram/paper_176205.htm.]

54 ENVIRONMENTAL SCIENCES↗