Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 163 records · Page 9

Using Host Galaxy Photometric Redshifts to Improve Cosmological Constraints with Type Ia Supernovae in the LSST Era

We perform a rigorous cosmology analysis on simulated Type Ia supernovae (SNe Ia) and evaluate the improvement from including photometric host galaxy redshifts compared to using only the "z spec " subset with spectroscopic redshifts from the host or SN. We use the Deep Drilling Fields (~50 deg 2 ) from the Photometric LSST Astronomical Time-Series Classification Challenge (PLAsTiCC) in combination with a low-z sample based on Data Challenge2. The analysis includes light-curve fitting to standardize the SN brightness, a high-statistics simulation to obtain a bias-corrected Hubble diagram, a statistical+systematics covariance matrix including calibration and photo-z uncertainties, and cosmology fitting with a prior from the cosmic microwave background. Compared to using the z spec subset, including events with SN+host photo-z results in (i) more precise distances for z > 0.5, (ii) a Hubble diagram that extends 0.3 further in redshift, and (iii) a 50% increase in the Dark Energy Task Force figure of merit (FoM) based on the w 0 w a CDM model. Analyzing 25 simulated data samples, the average bias on w 0 and w a is consistent with zero. The host photo-z systematic of 0.01 reduces FoM by only 2% because (i) most z < 0.5 events are in the z spec subset, (ii) the combined SN+host photo-z has ×2 smaller bias, and (iii) the anticorrelation between fitted redshift and color self-corrects distance errors. To prepare for analyzing real data, the next SN Ia cosmology analysis with photo-zs should include non–SN Ia contamination and host galaxy misassociations.

79 ASTRONOMY AND ASTROPHYSICS↗

Creating A Consistent Historical NASA POWER Solar Radiation Dataset to Support Renewable Energy, Building Energy Efficiency and Agro-Climatology Decisions

Prediction of Worldwide Energy Resources (POWER) project provides irradiance dataset to support renewable energy, building energy efficiency and agricultural needs. These datasets are derived from Global Energy and Water Cycle Experiment Surface Radiation Budget (GEWEX SRB) and Clouds and the Earth’s Radiant Energy System (CERES SYN1Deg). A systematic bias has been reported between these two datasets for the years with overlapping observations. For obtaining a consistent climate data record spanning the entire time record of observations, it is crucial to understand and remove the bias in the irradiance dataset. Inconsistency in solar radiation data can lead to inaccurate conclusions about solar energy potential and obscure real trends in solar radiation patterns that would impact energy availability assessments. In this study, we adapt quantile mapping approach to remove the systematic bias and to improve reliability of shortwave and longwave irradiance data. We present a validation of the bias corrected data against ground truth. For each 1° latitude and 1° longitude grid box across the globe, we match the CDFs of the reference dataset (CERES SYN1Deg) to that of the SRB dataset, thereby, adjusting the irradiance values to match the empirical distribution of two different measurements. The performance of quantile mapping is evaluated by using the metrics such as Mean Absolute Deviation (MAD). The results indicate that the quantile mapping significantly improves the accuracy and reliability of solar irradiance dataset especially for the weather conditions associated with high cloud cover and extreme irradiance values. The initial range of MAD for the studied sites for daily data was 4 to 19 Wm-2. After correction these reduced to 3 to 7 Wm-2. The findings from this study have important implications for solar energy system design, agricultural planning, and climate modeling community. Reducing the inconsistency and biases in solar irradiance dataset can enable better planning and operation of solar energy systems, leading to increased efficiency and cost-effectiveness. Additionally, this work also contributes to the statistical post-processing techniques in the renewable energy domain and highlights the potential of historical and near-real-time NASA POWER dataset as a valuable resource for solar energy research and applications.

POWER↗

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗

Developing A Continuous Ozone Record Through the SAGE and Aura Missions With NASA Reanalysis Products

During the last quarter of the 20th century, the Stratospheric Aerosol and Gas Experiment (SAGE) missions were crucial in monitoring the loss and subsequent recovery of the stratospheric ozone layer. Due to the employed solar occultation and self-calibration method, the SAGE monitors have produced stable data throughout the lifetime of each instrument. However, over ten years passed between the end of the SAGE II and SAGE III/M3M missions in 2005 and the launch of SAGE III/ISS instrument in 2017, leaving a gap in the data that must be bridged in order to assess trends in the ozone record. The Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) reanalysis product, with output available starting in 1980, is an attractive candidate for trend analysis due to the statistically optimized combination of multiple observing systems and the regular temporal and spatial coverage. However, changes in the assimilated observation systems can introduce discontinuities within the MERRA-2 ozone record, such as in 2004 when the MERRA-2 system shifted from assimilating ozone retrievals collected by SBUV instruments to those collected by instruments onboard the Aura satellite. In this study, we explore using the SAGE II record as a transfer function to develop a stable reanalysis data product, suitable for trend analysis, from the start of the SAGE II record in 1984 through the present. We follow the procedure outlined by Wargan et al. (2018) to address discontinuities in the MERRA-2 ozone dataset at the 2004 transition and during the Aura record. SAGE II ozone profiles are used to correct discontinuities in upper stratospheric ozone associated with changes in the MERRA-2 meteorological observing system in 1998 and 1995. We will then assess the relative performance of the data from different SAGE sensors using the resulting bias-corrected MERRA-2 ozone fields.

Pamela Wales↗

Soil Moisture Data Assimilation to Estimate Irrigation Water Use

Knowledge of irrigation is essential to support food security, manage depleting water resources, and comprehensively understand the global water and energy cycles. Despite the importance of understanding irrigation, little consistent information exists on the amount of water that is applied for irrigation. In this study, we develop and evaluate a new method to predict daily to seasonal irrigation magnitude using a particle batch smoother data assimilation approach, where land surface model soil moisture is applied in different configurations to understand how characteristics of remotely sensed soil moisture may impact the performance of the method. The study employs a suite of synthetic data assimilation experiments, allowing for systematic diagnosis of known error sources. Assimilation of daily synthetic soil moisture observations with zero noise produces irrigation estimates with a seasonal bias of 0.66% and a correlation of 0.95 relative to a known truth irrigation. When synthetic observations were subjected to an irregular overpass interval and random noise similar to the Soil Moisture Active Passive satellite (0.04 cm(exp 3) cm(exp -3)), irrigation estimates produced a median seasonal bias of <1% and a correlation of 0.69. When systematic biases commensurate with those between NLDAS‐2 land surface models and Soil Moisture Active Passive are imposed, irrigation estimates show larger biases. In this application, the particle batch smoother outperformed the particle filter. The presented framework has the potential to provide new information into irrigation magnitude over spatially continuous domains, yet its broad applicability is contingent upon identifying new method(s) of determining irrigation schedule and correcting biases between observed and simulated soil moisture, as these errors markedly degraded performance.

R Abolafia-Rosenzweig↗

The dark energy survey supernova program: investigating beyond-ΛCDM

We report constraints on a variety of non-standard cosmological models using the full 5-yr photometrically classified type Ia supernova sample from the Dark Energy Survey (DES-SN5YR). Both Akaike Information Criterion (AIC) and Suspiciousness calculations find no strong evidence for or against any of the non-standard models we explore. When combined with external probes, the AIC and Suspiciousness agree that 11 of the 15 models are moderately preferred over Flat-|$\Lambda$|CDM suggesting additional flexibility in our cosmological models may be required beyond the cosmological constant. We also provide a detailed discussion of all cosmological assumptions that appear in the DES supernova cosmology analyses, evaluate their impact, and provide guidance on using the DES Hubble diagram to test non-standard models. An approximate cosmological model, used to perform bias corrections to the data holds the biggest potential for harbouring cosmological assumptions. We show that even if the approximate cosmological model is constructed with a matter density shifted by |$\Delta \Omega _{\rm m}\sim 0.2$| from the true matter density of a simulated data set the bias that arises is subdominant to statistical uncertainties. Nevertheless, we present and validate a methodology to reduce this bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of the Jason-2 Extension to the TOPEX/Poseidon, Jason-l Sea-Surface Height Time Series for Global Mean Sea Level Monitoring

The Jason-2 (OSTM) follow-on mission to Jason-I provides for the continuation of global and regional mean sea level estimates along the ground-track of the initial phase of the TOPEX/Poseidon mission. During the first several months, Jason-I and Jason-2 flew in formation separated by only 55 seconds, enabling the isolation of intermission instrument biases through direct collinear differencing of near simultaneous observations. The Jason-2 Ku-band range bias with respect to Jason-I is estimated to be -84 +/- 9 mm, based on the orbit altitudes provided on the Geophysical Data Records. Modest improved agreement is achieved with the GSFC replacement orbits, which further enables the isolation of subtle 1 cm) instrument-dependent range correction biases. Inter-mission bias estimates are confirmed with an independent assessment from comparisons to a 64-station tide-gauge network, also providing an estimate of the stability of the 17-year time series to be less than 0.1 mm/yr +/- 0.4 mm/yr. The global mean sea level derived from the multi-mission altimeter sea-surface height record from January 1993 through September 2009 is 3.3 +/- 0.4 mm/yr. Recent trends over the period from 2004 through 2008 are smaller and estimated to be 2.0 +/- 0.4 mm/yr.

Beckley, B. D.↗

Evaluation of the Near-Surface Variables in the HRRR Weather Model Using Observations from the ARM SGP Site

Abstract The performance of version 4 of the NOAA High-Resolution Rapid Refresh (HRRR) numerical weather prediction model for near-surface variables, including wind, humidity, temperature, surface latent and sensible fluxes, and longwave and shortwave radiative fluxes, is examined over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) region. The study evaluated the model’s bias and bias-corrected mean absolute error relative to the observations on different time scales. Forecasts of near-surface geophysical variables at five SGP sites (HRRR at 3-km scale) were found to agree well with observations, but some consistent observation–forecast differences also occurred. Sensible and latent heat fluxes are the most challenging variables to be reproduced. The diurnal cycle is the main temporal scale affecting observation–forecast differences of the near-surface variables, and almost all of the variables showed different biases throughout the diurnal cycle. Results show that the overestimation of downward shortwave and the underestimation of downward longwave radiative flux are the two major biases found in this study. The timing and magnitude of downward longwave flux, wind speed, and sensible and latent heat fluxes are also different with contributions from model representations, data assimilation limitations, and differences in scales between HRRR and SGP sites. The positive bias in downward shortwave and negative bias in longwave radiation suggests that the model is underestimating cloud fraction in the study domain. The study concludes by showing a brief comparison with version 3 of the HRRR and shows that version 4 has better performance in almost all near-surface variables. Significance Statement A correct representation of the near-surface variables is important for numerical weather prediction models. This study investigates the capability of the latest NOAA High-Resolution Rapid Refresh (HRRRv4) model in simulating the near-surface variables by comparing against the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) in situ observations. Among others, we find that the surface heat fluxes, such as sensible and latent heat fluxes, are the most difficult variables to be reproduced. This study also shows that the diurnal cycle has the dominant impact on the model’s performance, which means the majority of the outputted near-surface variables have the strong diurnal cycle in their bias errors.

54 ENVIRONMENTAL SCIENCES↗

Toward a Framework for Systematic Error Modeling of NASA Spaceborne Radar with NOAA/NSSL Ground Radar-Based National Mosaic QPE

Characterization of the error associated to satellite rainfall estimates is a necessary component of deterministic and probabilistic frameworks involving space-born passive and active microwave measurement") for applications ranging from water budget studies to forecasting natural hazards related to extreme rainfall events. We focus here on the error structure of NASA's Tropical Rainfall Measurement Mission (TRMM) Precipitation Radar (PR) quantitative precipitation estimation (QPE) at ground. The problem is addressed by comparison of PR QPEs with reference values derived from ground-based measurements using NOAA/NSSL ground radar-based National Mosaic and QPE system (NMQ/Q2). A preliminary investigation of this subject has been carried out at the PR estimation scale (instantaneous and 5 km) using a three-month data sample in the southern part of US. The primary contribution of this study is the presentation of the detailed steps required to derive trustworthy reference rainfall dataset from Q2 at the PR pixel resolution. It relics on a bias correction and a radar quality index, both of which provide a basis to filter out the less trustworthy Q2 values. Several aspects of PR errors arc revealed and quantified including sensitivity to the processing steps with the reference rainfall, comparisons of rainfall detectability and rainfall rate distributions, spatial representativeness of error, and separation of systematic biases and random errors. The methodology and framework developed herein applies more generally to rainfall rate estimates from other sensors onboard low-earth orbiting satellites such as microwave imagers and dual-wavelength radars such as with the Global Precipitation Measurement (GPM) mission.

Kirstettier, Pierre-Emmanual↗

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗

GEWEX Water and Energy Budget Study

Closing the global water and energy budgets has been an elusive Global Energy and Water-cycle Experiment (GEWEX) goal. It has been difficult to gather many of the needed global water and energy variables and processes, although, because of GEWEX, we now have globally gridded observational estimates for precipitation and radiation and many other relevant variables such as clouds and aerosols. Still, constrained models are required to fill in many of the process and variable gaps. At least there are now several atmospheric reanalyses ranging from the early National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and NCEP/Department of Energy (DOE) reanalyses to the more recent ERA40 and JRA-25 reanalyses. Atmospheric constraints include requirements that the models state variables remain close to in situ observations or observed satellite radiances. This is usually done by making short-term forecasts from an analyzed initial state; these short-term forecasts provide the next guess, which is corrected by comparison to available observations. While this analysis procedure is likely to result in useful global descriptions of atmospheric temperature, wind and humidity, there is no guarantee that relevant hydroclimate processes like precipitation, which we can observe and evaluate, and evaporation over land, which we cannot, have similar verisimilitude. Alternatively, the Global Land Data Assimilation System (GLDAS), drives uncoupled land surface models with precipitation, surface solar radiation, and surface meteorology (from bias-corrected reanalyses during the study period) to simulate terrestrial states and surface fluxes. Further constraints are made when a tuned water balance model is used to characterize the global runoff observational estimates. We use this disparate mix of observational estimates, reanalyses, GLDAS and calibrated water balance simulations to try to characterize and close global and terrestrial atmospheric and surface water and energy budgets to within 10-20% for long term (1986-1995), large-scale global to regional annual means.

Roads, J.↗

Pooling Data Improves Multimodel IDF Estimates over Median-Based IDF Estimates: Analysis over the Susquehanna and Florida

Traditional multimodel methods for estimating future changes in precipitation intensity, duration, and frequency (IDF) curves rely on mean or median of models’ IDF estimates. Such multimodel estimates are impaired by large estimation uncertainty, shadowing their efficacy in planning efforts. Here, assuming that each climate model is one representation of the underlying data generating process, i.e., the Earth system, we propose a novel extension of current methods through pooling model data: (i) evaluate performance of climate models in simulating the spatial and temporal variability of the observed annual maximum precipitation (AMP), (ii) bias-correct and pool historical and future AMP data of reasonably performing models, and (iii) compute IDF estimates in a nonstationary framework from pooled historical and future model data. Pooling enhances fitting of the extreme value distribution to the data and assumes that data from reasonably performing models represent samples from the “true” underlying data generating distribution. Through Monte Carlo simulations with synthetic data, we show that return periods derived from pooled data have smaller biases and lesser uncertainty than those derived from ensembles of individual model data. We apply this method to NA-CORDEX models to estimate changes in 24-h precipitation intensity–frequency (PIF) estimates over the Susquehanna watershed and Florida peninsula. Our approach identifies significant future changes at more stations compared to median-based PIF estimates. The analysis suggests that almost all stations over the Susquehanna and at least two-thirds of the stations over the Florida peninsula will observe significant increases in 24-h precipitation for 2–100-yr return periods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain Corrected Satellite Sea Surface Salinity Using the GMAO S2S Forecast System

During the past years, we have seen that the La Nina to El Nino transition has had a significant global socio-economic impact and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help to improve the initialization of the thermocline, while satellite Sea Surface Temperature (SST) aids in constraining surface heat-fluxes, leading to improved subseasonal to seasonal forecasts of the coupled system. However, few studies have focused on improving the fresh-water flux and near-surface density and mixing through assimilation of satellite sea surface salinity (SSS). For expediency, the few ocean models that do assimilate SSS bias-correct the satellite SSS data to normalize towards the near-surface Argo data. However, in rainy regions, where buoyant water sits as a fresh lens at the surface, this assumption is likely inadequate. In previous work, we have shown that adjusting SSS data to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 has improved the near-surface density and mixed layer depth, leading to deeper thermocline and improved the NINO3.4 SST forecasts. Now we utilize the Soil Moisture/Ocean Salinity, Rain Corrected (SMOS_RC) SSS product provided by the Centre Aval de Traitement des données SMOS (CATDS CPDC) to represent the Sb more accurately at first model layer (in our case 5 m). Rather than using a diffusivity model as with RIM, SMOS_RC relies on an observed relationship between the spatial heterogeneity of SMOS SSS and instantaneous rain rate (RR) (Supply et al., 2020). In order to test the impact of SMOS_RC versus SMOS, we compare two reanalyses over the period 2014 to 2021. For both reanalysis experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is (i.e., with the fresh bias), and a separate reanalysis is performed assimilating the SMOS_RC data. We assess the impact for near-surface and subsurface dynamics within ocean reanalyses by validating against observations and explore how SSS assimilation (SMOS versus SMOS_RC) impacts dynamical ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (GEOS S2S-3, Molod et al., 2020, Hackert et al., 2023). We will show that improved SSS estimates and near-surface density and mixing led to more accurate coupled air/sea interaction and better ENSO forecasts.

Eric Hackert↗

Using Reanalysis in Crop Monitoring and Forecasting Systems

Weather observations are essential for crop monitoring and forecasting but they are not always available and in some cases they have limited spatial representativeness. Thus, reanalyses represent an alternative source of information to be explored. In this study, we assess the feasibility of reanalysis-based crop monitoring and forecasting by using the system developed and maintained by the European Commission- Joint Research Centre, its gridded daily meteorological observations, the biased-corrected reanalysis AgMERRA and the ERA-Interim reanalysis. We focus on Europe and on two crops, wheat and maize, in the period 1980-2010 under potential and water-imited conditions. In terms of inter-annual yield correlation at the country scale, the reanalysis-driven systems show a very good performance for both wheat and maize (with correlation values higher than 0.6 in almost all EU28 countries) when compared to the observations-driven system. However, significant yield biases affect both crops. All simulations show similar correlations with respect to the FAO reported yield time series. These findings support the integration of reanalyses in current crop monitoring and forecasting systems and point to the emerging opportunities linked to the coming availability of higher-resolution reanalysis updated at near real time.

Reanalysis↗

Reconstructing Hanford worker external doses from photons for epidemiology

The accurate reconstruction of external photon doses is essential for credible radiation epidemiology. This article presents the methodology used to derive dose estimates for 37 012 Hanford Site workers included in the Million Person Study. The approach employs historical dose records from the Hanford Radiation Exposure database and a previous epidemiology study. Bias correction factors specific to dosimeter type and period of use were applied and missing annual doses were estimated using a hierarchical nearby method to estimate deep dose equivalent for each worker. For early years with limited detection sensitivity, missed doses were quantified based on expected time-period-specific, low-dose statistical distributions. The revised dose estimates resulted in lower median and mean career doses than unadjusted data, while increasing the number of person-years with nonzero dose. Sensitivity analyses assessed the influence of bias in dosimetry measurements, missed doses and gap years on dose estimates. Differences in cumulative dose estimates between unadjusted and revised annual estimates are most prominent in the early operational years due to the highest bias during that time period.

dose reconstruction↗

Evaluation of the Surface PM2.5 in Version 1 of the NASA MERRA Aerosol Reanalysis over the United States

We use surface fine particulate matter (PM2.5) measurements collected by the United States Environmental Protection Agency (US EPA) and the Interagency Monitoring of Protected Visual Environments (IMPROVE) networks as independent validation for Version 1 of the Modern Era Retrospective analysis for Research and Applications Aerosol Reanalysis (MERRAero) developed by the Global Modeling Assimilation Office (GMAO). MERRAero is based on a version of the GEOS-5 model that is radiatively coupled to the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) aerosol module and includes assimilation of bias corrected Aerosol Optical Depth (AOD) from Moderate Resolution Imaging Spectroradiometer (MODIS) sensors on both Terra and Aqua satellites. By combining the spatial and temporal coverage of GEOS-5 with observational constraints on AOD, MERRAero has the potential to provide improved estimates of PM2.5 compared to the model alone and with greater coverage than available observations.Importantly, assimilation of AOD data constrains the total column aerosol mass in MERRAero subject to assumptions about optical properties for each of the species represented in GOGART. However, single visible wavelength AOD data does not contain sufficient information content to correct errors in either aerosol vertical placement or composition, critical elements for a proper characterization of surface PM2.5. Despite this, we find that the data-assimilation equipped version of GEOS-5 better represents observed PM2.5 between 2003 and 2012 compared to the same version of the model without AOD assimilation. Compared to measurements from the EPA-AQS network, MERRAero shows better PM2.5 agreement with the IMPROVE network measurements, which are composed essentially of rural stations. Regardless the data network, MERRAero PM2.5 are closer to observation values during the summer while larger discrepancies are observed during the winter. Comparing MERRAero to PM2.5 data collected by the Chemical Speciation Network (CSN) offers greater insight on the species MERRAero predicts well and those for which there are biases relative to the EPA observations. Analysis of this speciated data indicates that the lack of nitrate emissions in MERRAero and an underestimation of carbonaceous emissions in the Western US explains much of the reanalysis bias during the winter. To further understand discrepancies between the reanalysis and observations, we use complimentary data to assess two important aspects of MERRAero that are of relevance to the diagnosis of PM2.5, in particular AOD and vertical structure

Air Pollution↗