Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data bias”

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 469 records · Page 26

Representing Global Soil Erosion and Sediment Flux in Earth System Models

Abstract Soil erosion produces enormous amounts of sediment, carbon, and nutrient fluxes from land to rivers, thus playing crucial roles in global biogeochemical cycles and food security. To predict soil erosion in the context of climate and land use changes, we explicitly parameterize cropland management actions (i.e., irrigation, conserved agriculture, and crop residue management) and geological factors (i.e., lithology and glacier) in the Energy Exascale Earth System Model (E3SM) soil erosion module. The erosion model is calibrated using a global‐scale regionalized parameter calibration method. The spatial variabilities of the modeled and the Revised Universal Soil Loss Equation (RUSLE) based benchmark soil erosion are consistent across vegetation, climate and soil properties. Compared with independent data, our model shows a bias reduction in 59% of the observations relative to the RUSLE‐based soil erosion, with 53% of the bias reduction exceeding 50%. This improvement is mainly due to a better representation of the topographic effect on soil erosion. Our results indicate that conserved agriculture practices have effectively reduced soil erosion in cropland by over 25% in the United States and Argentina. In contrast, irrigation has increased soil erosion in many Asian countries. For upland sediment flux, our model is consistent with the WBMsed benchmark data in inter‐basin variability but could be more skillful in simulating intra‐basin variability because it couples soil erosion and sediment flux explicitly. The developed model provides useful skills for more realistic predictions of soil erosion and river sediment dynamics under environmental changes.

58 GEOSCIENCES↗

How deep is your soil? Quantifying and spatially analyzing understudied deep soil in the United States

Deep soil is largely understudied and important in understanding biogeochemical processes in soil. Here, understudied soil is defined as the difference between soil studied to a known depth and the estimated bedrock depth. To understand more about deep soil, the understudied soil in the US was quantified and spatially analyzed using soil survey data and model estimates of bedrock depth. An equation was derived to find understudied soil using the dataset parameters “max lower depth studied”, “depth to bedrock”, and “likelihood of bedrock in the top 200 cm”. The survey data and bedrock model revealed that soil has been studied to an average depth of 1-2 meters, and the average depth to bedrock is 20 meters. Soil data density in the soil surveys was greatest in the West Coast, Midwest, and areas historically managed for agricultural, while the non-contiguous US and interior West were underrepresented. The soil had been studied deeper than the estimated soil depth in 455 out of 56,889 observation points concentrated in Alaska, California, Texas, Florida, Puerto Rico, and the US Virgin Islands. To understand the diversity and any taxonomic bias of the global soil data available, soil order was compared to US-based National Resource Conservation Service percentages and it was found that Oxisols, Alfisols, Ultisols, Andisols, and Histosols were overrepresented while Gelisols, Aridisols, Vertisols, Entisols, and Spodosols are underrepresented. Soil depth is important in exploring the complexity of biogeochemical processes that take place in soil.

Bedrock↗

New Methods for Retrieval of Chlorophyll Red Fluorescence from Hyperspectral Satellite Instruments: Simulations and Application to GOME-2 and SCIAMACHY

Global satellite measurements of solar-induced fluorescence (SIF) from chlorophyll over land and ocean have proven useful for a number of different applications related to physiology, phenology, and productivity of plants and phytoplankton. Terrestrial chlorophyll fluorescence is emitted throughout the red and far-red spectrum, producing two broad peaks near 683 and 736nm. From ocean surfaces, phytoplankton fluorescence emissions are entirely from the red region (683nm peak). Studies using satellite-derived SIF over land have focused almost exclusively on measurements in the far red (wavelengths greater than 712nm), since those are the most easily obtained with existing instrumentation. Here, we examine new ways to use existing hyperspectral satellite data sets to retrieve red SIF (wavelengths less than 712nm) over both land and ocean. Red SIF is thought to provide complementary information to that from the far red for terrestrial vegetation. The satellite instruments that we use were designed to make atmospheric trace-gas measurements and are therefore not optimal for observing SIF; they have coarse spatial resolution and only moderate spectral resolution (0.5nm). Nevertheless, these instruments, the Global Ozone Monitoring Instrument 2 (GOME-2) and the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY), offer a unique opportunity to compare red and far-red terrestrial SIF at regional spatial scales. Terrestrial SIF has been estimated with ground-, aircraft-, or satellite-based instruments by measuring the filling-in of atmospheric andor solar absorption spectral features by SIF. Our approach makes use of the oxygen (O2) gamma band that is not affected by SIF. The SIF-free O2 gamma band helps to estimate absorption within the spectrally variable O2 B band, which is filled in by red SIF. SIF also fills in the spectrally stable solar Fraunhofer lines (SFLs) at wavelengths both inside and just outside the O2 B band, which further helps to estimate red SIF emission. Our approach is then an extension of previous approaches applied to satellite data that utilized only the filling-in of SFLs by red SIF. We conducted retrievals of red SIF using an extensive database of simulated radiances covering a wide range of conditions. Our new algorithm produces good agreement between the simulated truth and retrievals and shows the potential of the O2 bands for noise reduction in red SIF retrievals as compared with approaches that rely solely on SFL filling. Biases seen with existing satellite data, most likely due to instrumental artifacts that vary in time, space, and with instrument, must be addressed in order to obtain reasonable results. Our 8-year record of red SIF observations over land with the GOME-2 allows for the first time reliable global mapping of monthly anomalies. These anomalies are shown to have similar spatiotemporal structure as those in the far red, particularly for drought-prone regions. There is a somewhat larger percentage response in the red as compared with the far red for these areas that are drought sensitive. We also demonstrate that good-quality ocean fluorescence line height retrievals can be achieved with GOME-2, SCIAMACHY, and similar instruments by utilizing the full complement of radiance measurements that span the red SIF emission feature.

spatiotemporal structure↗

Skill Assessment of a Spectral Ocean-Atmosphere Radiative Model

Ocean phytoplankton, detrital material, and water absorb and scatter light spectrally. The Ocean- Atmosphere Spectral Irradiance Model (OASIM) is intended to provide surface irradiance over the oceans with sufficient spectral resolution to support ocean ecology, biogeochemistry, and heat exchange investigations, and of sufficient duration to support inter-annual and decadal investigations. OASIM total surface irradiance (integrated 200 nm to 4 microns) was compared to in situ data and three publicly available global data products at monthly 1-degree resolution. OASIM spectrally-integrated surface irradiance had root mean square (RMS) difference= 20.1 W/sq m (about 11%), bias=1.6 W/sq m (about 0.8%), regression slope= 1.01 and correlation coefficient= 0.89, when compared to 2322 in situ observations. OASIM had the lowest bias of any of the global data products evaluated (ISCCP-FD, NCEP, and ISLSCP 11), and the best slope (nearest to unity). It had the second best RMS, and the third best correlation coefficient. OASIM total surface irradiance compared well with ISCCP-FD (RMS= 20.7 W/sq m; bias=-11.4 W/sq m, r=0.98) and ISLSCP II (RMS =25.2 W/sq m; bias= -13.8 W/sq m; r=0.97), but less well with NCEP (RMS =43.0 W/sq m ;bias=-22.6 W/sq m; x=0.91). Comparisons of OASIM photosynthetically available radiation (PAR) with PAR derived from SeaWiFS showed low bias (-1.8 mol photons /sq m/d, or about 5%), RMS (4.25 mol photons /sq m/d ' or about 12%), near unity slope (1.03) and high correlation coefficient (0.97). Coupled with previous estimates of clear sky spectral irradiance in OASIM (6.6% RMS at 1 nm resolution), these results suggest that OASIM provides reasonable estimates of surface broadband and spectral irradiance in the oceans, and can support studies on ocean ecosystems, carbon cycling, and heat exchange.

Gregg, Watson, W.↗

Sea Surface Salinity Distribution in the Southern Ocean as Observed from Space

Large-scale spatial and temporal variabilities of sea surface salinity (SSS) in the Southern Ocean from 2011 to 2017 were studied using products derived from microwave sensors on board Aquarius, Soil Moisture and Ocean Salinity (SMOS), and Soil Moisture Active and Passive (SMAP) satellites. Four products, three from Aquarius and one from SMOS, were evaluated and shown to be generally consistent within 0.3 to 0.6 psu and agree favorably with in situ measurements. However, although the Aquarius products show consistent seasonality of SSS with high values of 34.45 psu in October and low values of 33.40 psu in May, theSMOS and SMAP products lack such seasonal variations. This may be caused by larger uncertainties in the SMOS and SMAP data due in part to the lack of concurrent scatterometer measurements that is used to correct for roughness effects. The four products provide similar spatial distributions of SSS with root mean square difference from 0.25 to 0.58 psu. Differences among Aquarius products are mainly due to varying salinity retrieval algorithms, smoothing, and masking of sea ice, while the SMOS product showed the highest SSS deviation that is likely due to the bias adjustment done on the data set. Our analyses show that SSS in the Southern Ocean region has significant meridional variations with the lowest SSS near the ice edge and highest at lower latitudes. The SSS is also lowest in summer indicating the predominant influence of sea ice and glacial melt, but it stays low near ice edges even during the growth season. Salinity together with temperature is the key parameters that control the state and circulation of the World's oceans. The large‐scale distribution of the sea surface salinity of the Southern Ocean is quantified for the first time using Aquarius, Soil Moisture and Ocean Salinity, and Soil Moisture Active and Passive satellite microwave data. Seasonal and interannual changes in the distribution as affected by the sea ice cover, surface temperature, and precipitation are evaluated. Comparative studies were done, using four different products that are publicly available, in conjunction with in situ observations to gain insights into the true nature of the distribution and how consistently the sea surface salinity is depicted by the different products. There are general consistencies in the products, and discrepancies are attributed to different algorithms, smoothing techniques, and sea ice masking. Aquarius data are shown to have higher accuracy than Soil Moisture and Ocean Salinity and Soil Moisture Active and Passive products in part because of having concurrent scatterometer that provides accurate correction to roughness effects.

Garcia-Eidell, Cynthia↗

Comparisons of the v11.1 Orbiting Carbon Observatory‐2 (OCO‐2) X CO2 Measurements With GGG2020 TCCON

The Orbiting Carbon Observatory 2 (OCO-2) is NASA's first Earth observation satellite mission dedicated to studying the sources and sinks of carbon dioxide (CO 2 ) on a global scale. The observations of reflected sunlight are inverted in a retrieval algorithm to produce estimates of the dry air mole-fractions of CO 2 (X CO2 ). The OCO-2 Level 2 data release, version 11.1 (v11.1) retrievals from the Atmospheric Carbon Observations from Space (ACOS) algorithm, includes significant improvements in the X CO2 data product compared to older OCO-2 data versions. This work compares the v11.1 X CO2 from OCO-2 against X CO2 estimates collected from a global ground-based network known as the Total Carbon Column Observing Network (TCCON), OCO-2's primary validation source. The OCO-2 project provides a version of the Level 2 data product, called “lite” files that include calibrated and bias-corrected XCO2 values, accessible together with all OCO-2 data products through the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). This work shows that OCO-2 X CO2 observations made between September 2014 and December 2023, after quality filtering and the application of an averaging kernel correction, agree well with coincident TCCON data for all OCO-2 observational modes of land (nadir, glint, target) and ocean (glint). The aggregated, bias-corrected, and quality-filtered absolute average bias values are less than or equal to 0.20 parts per million (ppm) globally for all OCO-2 observation modes, where the biases do not indicate a statistically significant time dependence. The land nadir/glint mode has the lowest bias value of −0.03 ± 0.85 ppm.

54 ENVIRONMENTAL SCIENCES↗

Repository of HydroSMADE: Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites

This repository presents HydroSMADE—Hydropower Site-level Monthly Availability Data Ensemble, a new open dataset that provides monthly hydropower availability for 1,593 existing and 124,333 potential sites worldwide over the period 1950–2100. The dataset is generated by using a global hydrologic model (Xanthos) with explicit representation of hydropower operation. Specifically, HydroSMADE distinguishes between storage and diversion sites, applies optimized operating rules, and incorporates site-specific characteristics such as generation capacity, maximum turbine flow, and reservoir storage. Driven by bias-corrected meteorological inputs, the data is provided for 30 alternative future scenarios. The scenarios consist of the full factorial combination of three standard CMIP6 atmospheric forcing pathways (SSP1-2.6, SSP3-7.0, and SSP5-8.5) and ten CMIP6 General Circulation Models (GCMs): GFDL-ESM4, IPSL-CM6A-LR, MPI-ESM1-2-HR, MRI-ESM2-0, EC-Earth3, CanESM5, MIROC6, CNRM-ESM2-1, UKESM1-0-LL, and CNRM-CM6-1. The repository contains a total of 122 files: a text file (readme.txt) containing a brief description of the included data, a CSV file containing site attributes, and the remaining 120 files (in CSV) containing site-level monthly hydropower availability. Example Jupyter Notebooks to explore the HydroSMADE dataset are available on GitHub at https://github.com/kamal0013/HydroSMADE More details on the methods and technical validation of HydroSMADE are available in the following paper by the same authors: Chowdhury, A. K., Abeshu, G. W., Zhao, M., Wild, T. B., Hassan, N., Ying, Z., Kim, G. J., Matthew, B., Jonathan, L., & Li, H.-Y. (Submitted). Hydropower Site-level Monthly Availability Data Ensemble for 1950-2100 at Existing and Potential Global Sites.

Existing and Potential Sites↗

A noise source identification technique using an inverse Helmholtz integral equation method

A technique is developed which utilizes numerical models and field pressure information to characterize acoustic fields and identify acoustic sources. The numerical models are based on boundary element numerical procedures. Either pressure, velocity, or passive boundary conditions, in the form of impedance boundary conditions, may be imposed on the numerical model. Alternatively, if no boundary information is known, a boundary condition can be left unspecified. Field pressure data may be specified to overdetermine the numerical problem. The problem is solved numerically for the complete sound field from which the acoustic sources may be determined. The model can then be used to idenfify acoustic intensity paths in the field. The solution can be modified and the model used to evaluate design alternatives. In this investigation the method is tested analytically and verified. In addition, the sensitivity of the method to random and bias error in the input data is demonstrated.

Gardner, B. K.↗

Identification of Sea Surface Temperature and Sea Surface Salinity Fronts along the California Coast: Application Using Saildrone and Satellite Derived Products

Coastal upwelling regions are one of the most dynamic areas of the world’s oceans. The California and Baja California Coasts are impacted by both coastal upwelling and the California Current, leading to frontal activity that is captured by gradients in both Sea Surface Temperature (SST) and Sea Surface Salinity (SSS). Satellite data are a great source of spatial data to study fronts. However, biases near coastal areas and coarse resolutions can impair its usefulness in upwelling areas. In this work gradients in SST from NASA Multi-Scale Ultra-High Resolution (MUR) and in two SSS products derived from the Soil Moisture Active Passive (SMAP) NASA mission are compared directly with gradients derived from the Saildrone uncrewed vehicles to validate the gradients as well as to assess their ability to detect known frontal features. The three remotely sensed data sets (MURSST/JPL, SMAP SSS/RSS, SMAP SSS) were co-located with the Saildrone data prior to the calculation of the gradients. Wavelet analysis is used to determine how well the satellite derived SST and SSS products are reproducing the Saildrone derived gradients. Overall results indicate the remote sensing products are reproducing features of known areas of coastal upwelling. Differences between the SST and SSS gradients are mainly associated with the limitations of the microwave derived SSS coverage near land and its reduced spatial resolution. The results are promising for using remote sensing data sets to monitor frontal structure along the California Coast and the application to long term changes in coastal upwelling and dynamics.

Sea Surface Salinity↗

Bias Reduction in Short Records of Satellite Soil Moisture

Although surface soil moisture data from different sources (satellite retrievals, ground measurements, and land model integrations of observed meteorological forcing data) have been shown to contain consistent and useful information in their seasonal cycle and anomaly signals, they typically exhibit very different mean values and variability. These biases pose a severe obstacle to exploiting the useful information contained in satellite retrievals through data assimilation. A simple method of bias removal is to match the cumulative distribution functions (cdf) of the satellite and model data. However, accurate cdf estimation typically requires a long record of satellite data. We demonstrate here that by wing spatial sampling with a 2 degree moving window we can obtain local statistics based on a one-year satellite record that are a good approximation to those that would be derived from a much longer time series. This result should increase the usefulness of relatively short satellite data records.

Reichle, Rolf H.↗

Ionospheric Signatures in Radio Occultation Data

We can extend robustly the radio occultation data record by 6 years (+60%) by developing a singlefrequency processing method for GPS/MET data. We will produce a calibrated data set with profile-byprofile data characterization to determine robust upper bounds on ionospheric bias. Part of an effort to produce a calibrated RO data set addressing other key error sources such as upper boundary initialization. Planned: AIRS-GPS water vapor cross validation (water vapor climatology and trends).

radio occultation↗

Estimating Subhourly Inverter Clipping Loss From Satellite-Derived Irradiance Data

Photovoltaic system production simulations are conventionally run using hourly weather datasets. Hourly simulations are sufficiently accurate to predict the majority of long-term system behavior but cannot resolve high-frequency effects like inverter clipping caused by short-duration irradiance variability. Direct modeling of this subhourly clipping error is only possible for the few locations with high-resolution irradiance datasets. This paper describes a method of predicting the magnitude of this error using a machine learning regressor ensemble model, comprised of a random forest and an XGBoost model, and 30-minute satellite irradiance data. The method predicts a correction for each 30-minute interval with the potential to roll up into 60-minute corrections to match an hourly energy model. The model is trained and validated at locations where the error can be directly simulated from 1-minute ground data. The validation shows low bias at most ground station locations. The model is also applied to gridded satellite irradiance to produce a heatmap of the estimated clipping error across the United States. Finally, the relative importance of each predictor satellite variable is retrieved from the model and discussed.

41 EE - Solar Energy Technologies Office (EE-4S)↗

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↗

Propagating Uncertainties in the SALT3 Model-training Process to Cosmological Constraints

Type Ia supernovae (SNe Ia) are standardizable candles that must be modeled empirically to yield cosmological constraints. To understand the robustness of this modeling to variations in the model-training procedure, we build an end-to-end pipeline to test the recently developed SALT3 model. We explore the consequences of removing pre-2000s low-z or poorly calibrated U-band data, adjusting the amount and fidelity of SN Ia spectra, and using a model-independent framework to simulate the training data. We find that the SALT3 model surfaces are improved by having additional spectra and U-band data, and can be shifted by ~5% if host-galaxy contamination is not sufficiently removed from SN spectra. We find that resulting measurements of w are consistent to within 2.5% for all of the training variants explored in this work, with the largest shifts coming from variants that add color-dependent calibration offsets or host-galaxy contamination to the training spectra and those that remove pre-2000s low-z data. These results demonstrate that the SALT3 model-training procedure is largely robust to reasonable variations in the training data, but that additional attention must be paid to the treatment of spectroscopic data in the training process. We also find that the training procedure is sensitive to the color distributions of the input data—the resulting w measurement can be biased by ~2% if the color distribution is not sufficiently wide. Future low-z data, particularly u-band observations and high signal-to-noise ratio SN Ia spectra, will help to significantly improve SN Ia modeling in the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗

Single-Receiver GPS Phase Bias Resolution

Existing software has been modified to yield the benefits of integer fixed double-differenced GPS-phased ambiguities when processing data from a single GPS receiver with no access to any other GPS receiver data. When the double-differenced combination of phase biases can be fixed reliably, a significant improvement in solution accuracy is obtained. This innovation uses a large global set of GPS receivers (40 to 80 receivers) to solve for the GPS satellite orbits and clocks (along with any other parameters). In this process, integer ambiguities are fixed and information on the ambiguity constraints is saved. For each GPS transmitter/receiver pair, the process saves the arc start and stop times, the wide-lane average value for the arc, the standard deviation of the wide lane, and the dual-frequency phase bias after bias fixing for the arc. The second step of the process uses the orbit and clock information, the bias information from the global solution, and only data from the single receiver to resolve double-differenced phase combinations. It is called "resolved" instead of "fixed" because constraints are introduced into the problem with a finite data weight to better account for possible errors. A receiver in orbit has much shorter continuous passes of data than a receiver fixed to the Earth. The method has parameters to account for this. In particular, differences in drifting wide-lane values must be handled differently. The first step of the process is automated, using two JPL software sets, Longarc and Gipsy-Oasis. The resulting orbit/clock and bias information files are posted on anonymous ftp for use by any licensed Gipsy-Oasis user. The second step is implemented in the Gipsy-Oasis executable, gd2p.pl, which automates the entire process, including fetching the information from anonymous ftp

Bertiger, William I.↗

Evaluation of Hanford 200 West Area Tank Farms (241-S/241-SX-/241-U tank farms) Physical Properties Data for Use in Development of West Area Tank Treatment (WATT) Processing

With the recent acceptance of West Area Tank Treatment disposition alternative for 200 West Area tanks at the Hanford Site by the State of Washington and the U.S. Department of Energy, a review was initiated to identify the physical properties data available in the literature for the Hanford 241-S, 241-SX, and 241-U tank farms. The literature reviewed indicated that there is a relatively small set of useful data on physical properties of 200 West Area tanks, and the data that do exist are biased around a narrow range of tank samples. Much of the testing between the 1990s and mid-2010s was intended to support either enhanced sludge washing or feed delivery to the Pretreatment Facility at the Hanford Waste Treatment and Immobilization Plant. As such, some physical properties of 200 West Area samples were measured under conditions that are no longer relevant. Because of the distinctly different nature of many past processes at the 200 West Area compared to the 200 East Area, insight from waste testing in the 200 East Area waste should be used with caution, as there may be significantly different qualities in the physical properties data between these two areas (both in situ and as measured in laboratory analyses). Based on this assessment, there is a need to collect additional physical property data to support planning for 200 West Area retrievals.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Real-time Sub-cm Differential Orbit Determination of two Low-Earth Orbiters with GPS Bias Fixing

An effective technique for real-time differential orbit determination with GPS bias fixing is formulated. With this technique, only real-time GPS orbits and clocks are needed (available from the NASA Global Differential GPS System with 10-20 cm accuracy). The onboard, realtime orbital states of user satellites (few meters in accuracy) are used for orbit initialization and integration. An extended Kalman filter is constructed for the estimation of the differential orbit between the two satellites as well as a reference orbit, together with their associating dynamics parameters. Due to close proximity of the two satellites and of similar body shapes, the differential dynamics are highly common and can be tightly constrained which, in turn, strengthens the orbit estimation. Without explicit differencing of GPS data, double-differenced phase biases are formed by a transformation matrix. Integer-valued fixing of these biases are then performed which greatly strengthens the orbit estimation. A 9-day demonstration between GRACE orbits with baselines of approx.200 km indicates that approx.80% of the double-differenced phase biases can successfully be fixed and the differential orbit can be determined to approx.7 mm as compared to the results of onboard K-band ranging.

global positioning system (GPS)↗

Performance Before and After Irradiation of Pixelated 3D Silicon Sensors for the HL-LHC CMS Tracker

The Large Hadron Collider (LHC) particle accelerator at European Center for Nuclear Research (CERN) will be shut down starting in 2026 to achieve the High Luminosity Large Hadron Collider (HL-LHC) upgrade. The upgrade will allow for higher fluences, in order to increase the probability of detecting increasingly rare particles, and to obtain higher precision measurements of known particles. To accommodate the new accelerator conditions, many aspects of the Compact Muon Solenoid (CMS) detector will be upgraded; of particular interest for this thesis are the silicon pixel detectors located in the inner tracker. These will be replaced and upgraded to accommodate the higher fluences of the HL-LHC upgrade, as well as to replace existing sensors which have sustained radiation damage. In order for the new sensors to operate under high luminosity conditions, they must be increasingly radiation hard, and in order to detect rare particles, they must be increasingly more precise. The performance of one Centro Nacional de Microelectronica (CNM) 3D silicon sensor before and after undergoing irradiation at fluences similar to those which will be observed at the HL-LHC was investigated to determine radiation hardness and precision. Data was collected at Fermi National Laboratory (Fermilab), in the Fermi National Laboratory Test Beam Facility (FTBF) silicon tracker telescope, which can be used to determine the number of particles, and tracks made by high energy protons passing through. The sensor was also irradiated at Fermilab in the Irradiation Test Area (ITA). Prior to data collection a tuning procedure is carried out to determine ideal bias voltage operating conditions, mask noisy and dead pixels, adjust to the ideal threshold, and map sensor gain. Data is then collected at the FTBF, where the sensor is installed in the center of the FTBF silicon telescope. Variables, including angle and bias voltage, are varied throughout data collection. Data is then processed using an alignment software to determine the exact telescope geometry, along with the tracks which were observed passing through the sensor and telescope. Sensor performance was found to be comparable before and after irradiation, with irradiated results showing slightly lower efficiencies and cluster sizes. Position resolution is comparable both before and after irradiation, and similar distributions of cluster shape are observed. After irradiation, the sensor shows increasing collected charge with bias, an indication of increased width of the depletion region. Peak charge pre-irradiation is higher than post-irradiation peak charge, indicating the irradiated results are not taken under fully-depleted conditions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗