Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ensemble Impacts”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Numerical Study of Tropical Sea-Air Interactions Using a Cloud Resolving Model Coupled with an Ocean Mixed-Layer Model

Coupling a cloud resolving model (CRM) with an ocean mixed layer (OML) model can provide a powerful tool for better understanding impacts of atmospheric precipitation on sea surface temperature (SST) and salinity. The objective of this study is twofold. First, by using the three dimensional (3-D) CRM-simulated (the Goddard Cumulus Ensemble model, GCE) diabatic source terms, radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the OML model, the respective impact of individual component on upper ocean heat and salt budgets are investigated. Secondly, a two-way air-sea interaction between tropical atmospheric climates (involving atmospheric radiative-convective processes) and upper ocean boundary layer is also examined using a coupled two dimensional (2-D) GCE and OML model. Results presented here, however, only involve the first aspect. Complete results will be presented at the conference.

Shie, Chung-Lin↗
NASA NTRSDate not supplied

Coupled Aerosol-Chemistry Simulations of the January 2022 Eruptions of Hunga Tonga-Hunga Ha’apai in the NASA GEOS Earth System Model

The January 2022 eruptions of the underwater Hunga Tonga-Hunga Ha’apai volcano injected more than 100 Tg of water vapor and approximately 0.5 Tg of sulfur dioxide into the stratosphere. Injected materials reached as high as ~55 km in altitude, while the main plume from the eruption travelled west over Australia and the Indian Ocean between about 20 – 30 km altitude. We investigate the transport and the impact of the erupted materials in coupled aerosol-chemistry simulations performed with the NASA Goddard Earth Observing System (GEOS) Earth system model. A sulfur mechanism introduced into the Global Modeling Initiative (GMI) stratospheric-tropospheric chemistry package allows for an interactive simulation of the water vapor-chemistry impacts, and the large water vapor perturbation results in rapid conversion of sulfur dioxide to sulfate aerosol in our aerosol mechanism. We report on the results of a multi-year ensemble of simulations performed with this system that include a control ensemble (no eruption), a water vapor-only injection ensemble, and a water vapor and sulfur dioxide injection ensemble. Simulations of the near-field aerosol and chemistry transport are compared to available measurements from OMPS-LP, OMPS-NM, CALIOP, MLS, and other sensors. The extended impact of the volcanic materials on the stratospheric composition and chemistry over the next several years is further investigated in forecast simulations.

Peter Colarco↗
NASA NTRS2012

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low clean concentration and a high dirty concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated.

cloud-resolving model↗
NASA NTRS2010

Impact of Damping Uncertainty on SEA Model Response Variance

Statistical Energy Analysis (SEA) is commonly used to predict high-frequency vibroacoustic levels. This statistical approach provides the mean response over an ensemble of random subsystems that share the same gross system properties such as density, size, and damping. Recently, techniques have been developed to predict the ensemble variance as well as the mean response. However these techniques do not account for uncertainties in the system properties. In the present paper uncertainty in the damping loss factor is propagated through SEA to obtain more realistic prediction bounds that account for both ensemble and damping variance. The analysis is performed on a floor-equipped cylindrical test article that resembles an aircraft fuselage. Realistic bounds on the damping loss factor are determined from measurements acquired on the sidewall of the test article. The analysis demonstrates that uncertainties in damping have the potential to significantly impact the mean and variance of the predicted response.

Schiller, Noah↗
NASA NTRS2016

Similar Estimates of Temperature Impacts on Global Wheat Yield by Three Independent Methods

The potential impact of global temperature change on global crop yield has recently been assessed with different methods. Here we show that grid-based and point-based simulations and statistical regressions (from historic records), without deliberate adaptation or CO2 fertilization effects, produce similar estimates of temperature impact on wheat yields at global and national scales. With a 1 C global temperature increase, global wheat yield is projected to decline between 4.1% and 6.4%. Projected relative temperature impacts from different methods were similar for major wheat-producing countries China, India, USA and France, but less so for Russia. Point-based and grid-based simulations, and to some extent the statistical regressions, were consistent in projecting that warmer regions are likely to suffer more yield loss with increasing temperature than cooler regions. By forming a multi-method ensemble, it was possible to quantify 'method uncertainty' in addition to model uncertainty. This significantly improves confidence in estimates of climate impacts on global food security.

Climate impacts↗
NASA NTRS2010

Shorebird Migration Patterns in Response to Climate Change: A Modeling Approach

The availability of satellite remote sensing observations at multiple spatial and temporal scales, coupled with advances in climate modeling and information technologies offer new opportunities for the application of mechanistic models to predict how continental scale bird migration patterns may change in response to environmental change. In earlier studies, we explored the phenotypic plasticity of a migratory population of Pectoral sandpipers by simulating the movement patterns of an ensemble of 10,000 individual birds in response to changes in stopover locations as an indicator of the impacts of wetland loss and inter-annual variability on the fitness of migratory shorebirds. We used an individual based, biophysical migration model, driven by remotely sensed land surface data, climate data, and biological field data. Mean stop-over durations and stop-over frequency with latitude predicted from our model for nominal cases were consistent with results reported in the literature and available field data. In this study, we take advantage of new computing capabilities enabled by recent GP-GPU computing paradigms and commodity hardware (general purchase computing on graphics processing units). Several aspects of our individual based (agent modeling) approach lend themselves well to GP-GPU computing. We have been able to allocate compute-intensive tasks to the graphics processing units, and now simulate ensembles of 400,000 birds at varying spatial resolutions along the central North American flyway. We are incorporating additional, species specific, mechanistic processes to better reflect the processes underlying bird phenotypic plasticity responses to different climate change scenarios in the central U.S.

Smith, James A.↗
NASA NTRS2016

Multi-Wheat-Model Ensemble Responses to Interannual Climate Variability

We compare 27 wheat models' yield responses to interannual climate variability, analyzed at locations in Argentina, Australia, India, and The Netherlands as part of the Agricultural Model Intercomparison and Improvement Project (AgMIP) Wheat Pilot. Each model simulated 1981e2010 grain yield, and we evaluate results against the interannual variability of growing season temperature, precipitation, and solar radiation. The amount of information used for calibration has only a minor effect on most models' climate response, and even small multi-model ensembles prove beneficial. Wheat model clusters reveal common characteristics of yield response to climate; however models rarely share the same cluster at all four sites indicating substantial independence. Only a weak relationship (R2 0.24) was found between the models' sensitivities to interannual temperature variability and their response to long-termwarming, suggesting that additional processes differentiate climate change impacts from observed climate variability analogs and motivating continuing analysis and model development efforts.

uncertainty↗
NASA NTRS2015

Combined Observational and Modeling Efforts to Better Understand Aerosol-Cloud-Precipitation Interactions Over Land: Preliminary Results from 7-SEAS/BASELInE 2013

This talk presents some of the detailed observations of low-level stratocumulus over northern Vietnam during 7-SEASBASELInE 2013 by SMARTLabs' ACHIEVE W-band cloud radar and other remote sensing instruments. These observations are the first of their kind for this region and will aid in ongoing studies of biomass-burning aerosol impacts on local and regional weather and climate. Preliminary results from simulations using the Goddard Cumulus Ensemble (GCE) with recently implemented triple-moment bulk microphysics to examine the sensitivity of low-level stratocumulus over land to aerosols are also presented. Recommendations for future observational activities in the 7-SEAS northern region in collaboration with international partners will also be discussed.

satellite remote sensing↗
NASA NTRS2023

The Impact of ENSO and NAO Initial Conditions and Anomalies on the Modeled Response to Pinatubo-Sized Volcanic Forcing

Strong, strato-volcanic eruptions are a substantial, intermittent source of natural climate variability. Initial atmospheric and oceanic conditions, such as El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO), also naturally impact climate on interannual timescales. We examine how initial conditions of ENSO and NAO contribute to the evolution of climate in the period following a Pinatubo-type eruption using a large (81-member) ensemble of model simulations in GISS model E2.1-G. Simulations are initialized from sampled conditions of ENSO and NAO using the protocol of the coordinated CMIP6 Volcanic Model Intercomparison Project (VolMIP) – where aerosols are forced with respect to time, latitude, and height. We analyze paired anomalous variations (perturbed – control) to understand changes in global and regional climate responses under positive, negative, and neutral ENSO and NAO conditions. In particular, we find that for paired anomalies there is a high probability of strong (∼1.5 °C) warming of northern Eurasia surface air temperature in the first winter after the volcanic eruption for negative NAO ensembles coincident with decreased lower stratospheric temperature at the poles, decreased geopotential height, and strengthening of the stratospheric polar vortex. Climate anomalies (relative to average conditions across the control period), however, show no mean warming and suggest that the strength of this response is impacted by conditions present in the selected period of the control run. Again using paired anomalies, we also observe that under both +ENSO and −ENSO ensembles sea surface temperature decreases in the first post-eruptive boreal winter coinciding with surface cooling from volcanic aerosols. Neutral ENSO ensembles, on the other hand, show variability in their response with no clear trend in post-eruptive warming or cooling. In general, paired anomalies from unperturbed simulations give insight into the evolution of the climate response to volcanic forcing; however, when compared with anomalies from climatological conditions, it is clear that paired anomalies are significantly affected by sampled initial conditions occurring at the time of the volcanic eruption.

volcanic eruptions↗
NASA NTRS2020

An Examination of the Impact of Grid Spacing on WRF Simulations of Wintertime Precipitation in the Mid-Atlantic United States

A large set of deterministic and ensemble forecasts was produced to identify the optimal spacing for forecasting U.S. East Coast snowstorms. WRF forecasts were produced on cloud-allowing (~1-km grid spacing) and convection-allowing (3–4 km) grids, and compared against forecasts with parameterized convection (>~10 km). Performance diagrams were used to evaluate 19 deterministic forecasts from the winter of 2013–14. Ensemble forecasts of five disruptive snowstorms spanning the years 2015–18 were evaluated using various methods to evaluate probabilistic forecasts. While deterministic forecasts using cloud-allowing grids were not better than convection-allowing forecasts, both had lower bias and higher success ratios than forecasts with parameterized convection. All forecasts were underdispersive. Nevertheless, forecasts on the higher-resolution grids were more reliable than those with parameterized convection. Forecasts on the cloud-allowing grid were best able to discriminate areas that received heavy snow and those that did not, while the forecasts with parameterized convection were least able to do so. It is recommended to use convection-resolving and (if computationally possible) to use cloud-allowing forecast grids when predicting East Coast winter storms.

Snowstorms↗
NASA NTRS2023

Seasonal Soil Freeze/Thaw Variability Across North America via Ensemble Land Surface Modeling

Land surface modeling provides the opportunity to investigate winter soil processes such as soil freezing and thawing over large domains. However, the variability in simulated winter soil characteristics among land surface models and forcing datasets is not well understood. In this study, a nine-member ensemble was employed to characterize the spatial and inter-annual variability of three winter soil characteristics, annual number of frozen days, annual minimum temperature and annual number of freeze-thaw cycles over North America during the water years 2010–2016. The ensemble, including three land surface model (JULES, Noah2.7.1 and Noah-MP) and three forcing datasets (ECMWF, GDAS and MERRA2), was developed through the Snow Ensemble Uncertainty Project (SEUP). In many regions, there was remarkably good agreement across the ensemble for the winter soil temperatures. However, the variability among the ensemble's annual number of frozen days, as quantified with standard deviation, exceeded 150 days at the northern Pacific coastline. While the differences among the ensemble members' annual minimum temperature were typically <3 °C, the differences exceeded 6 °C north of 50°N. The ensemble members generally agreed within one to three freeze-thaw cycles in mid latitude regions, Alaska, and the south and west coasts of the USA. High variability among the ensemble, more than six freeze-thaw cycles, occurred in the Great Plains, northern Pacific coastline, and along the Appalachian Mountains. Differences in winter soil temperature characteristics were more apparent among the LSMs rather than the meteorological forcing datasets. Except for maritime regions, the Noah2.7.1 members had the fewest annual number of frozen days and the least number of freeze-thaw cycles. Noah2.7.1 also had the coldest annual minimum temperatures except for ephemeral regions. Comparisons between in-situ observations and the SEUP estimates of winter soil characteristics revealed that the modeled frozen period was much longer than observed, that the modeled annual minimum temperatures were much colder than observed, and the modeled freeze-thaw cycles occurred more frequently than observed. Excluding the high latitude sites, the observed frozen period is less than two months with the minimum temperature above − 5 °C at most of the studied sites, while the ensemble members simulated, on average, a four month frozen period with − 10 °C minimum temperature. Errors in capturing snow during the accumulation period appear to impact differences in modeled versus observed soil temperature throughout the entire winter.

Soil freeze/thaw↗
NASA NTRS2015

The Impact of Assimilating Precipitation-affected Radiance on Cloud and Precipitation in Goddard WRF-EDAS Analyses

High-frequency TMI and AMSR-E radiances, which are sensitive to precipitation over land, are assimilated into the Goddard Weather Research and Forecasting Model- Ensemble Data Assimilation System (WRF-EDAS) for a few heavy rain events over the continental US. Independent observations from surface rainfall, satellite IR brightness temperatures, as well as ground-radar reflectivity profiles are used to evaluate the impact of assimilating rain-sensitive radiances on cloud and precipitation within WRF-EDAS. The evaluations go beyond comparisons of forecast skills and domain-mean statistics, and focus on studying the cloud and precipitation features in the jointed rainradiance and rain-cloud space, with particular attentions on vertical distributions of height-dependent cloud types and collective effect of cloud hydrometers. Such a methodology is very helpful to understand limitations and sources of errors in rainaffected radiance assimilations. It is found that the assimilation of rain-sensitive radiances can reduce the mismatch between model analyses and observations by reasonably enhancing/reducing convective intensity over areas where the observation indicates precipitation, and suppressing convection over areas where the model forecast indicates rain but the observation does not. It is also noted that instead of generating sufficient low-level warmrain clouds as in observations, the model analysis tends to produce many spurious upperlevel clouds containing small amount of ice water content. This discrepancy is associated with insufficient information in ice-water-sensitive radiances to address the vertical distribution of clouds with small amount of ice water content. Such a problem will likely be mitigated when multi-channel multi-frequency radiances/reflectivity are assimilated over land along with sufficiently accurate surface emissivity information to better constrain the vertical distribution of cloud hydrometers.

WRF-EDAS Analyses↗
NASA NTRSDate not supplied

Oxygen Deficiency in Spaceflight & its Impact on Plants’ Adaptive Changes

The goal of this study was to investigate the effects of hypoxic conditions in spaceflight. The distribution of genes involved with hypoxia in Arabidopsis thaliana and Brassica rapa were analyzed with the results from past spaceflight experiments to evaluate genes for future studies. Transcriptomes data of two different spaceflight studies of Arabidopsis thaliana from the NASA GeneLab database, GLDS-7 and GLDS-17, were compared. DNA microarrays were utilized for transcription profiling to conduct these studies. For GLDS-7, the response in spaceflight was studied with approaches that collected gene expression data. Leaves, hypocotyls, and root tissues were compared to the whole plant. For GLDS-17, seedlings and undifferentiated cultured cells were placed in the Biological Research in Canisters (BRIC), specifically BRIC-16. The genes related to hypoxia in Arabidopsis thaliana from these two studies were compared to genes in Brassica rapa with the TOAST database to evaluate similarities. When transcriptomes were analyzed for GLDS-7 and 17, genes that were considered significant had p-values ≤ 0.05 and log fold change values ≤ -1 or ≥1. Sixteen genes fulfilled the criteria. The genes related to hypoxia were alcohol dehydrogenase, elongation factor, ethylene-responsive factor, GUS, heat-shock proteins, NAP, RAP2.12, and RD20. The genes most impacted by spaceflight were heat-shock proteins. These genes were compared with Brassica rapa through Arabidopsis Ensemble Orthology from the TOAST Database. Similarities were seen in alcohol dehydrogenase, elongation factor, ethylene-responsive factor, heat-shock proteins, NAP, and RAP2.12. Overall, transcription profiling indicates that plants’ survival in spaceflight is dependent on adaptive changes with gene expression. This study also indicates that there are similarities in gene expression between Arabidopsis thaliana and Brassica rapa with comparable gene expression. Future studies could include analyzing additional species to understand which genes could be modified to ensure better yield of space crops amid hypoxic conditions.

hypoxia↗
NASA NTRS2021

Dynamic Ensemble Prediction of Cognitive Performance in Space

Astronauts are exposed to a unique set of stressors in spaceflight. Microgravity, isolation, confinement, and environmental and operational hazards: all of these can impact sleep, vigilant attention, and alertness, which are critical to mission success. In this paper, we seek to understand the most important predictors of alertness over the course of a space mission, using self-reported, cognitive, and environmental data collected from 24 astronauts on 6-month missions to the International Space Station (ISS). Alertness was repeatedly and objectively assessed on the ISS with a brief 3-minute Psychomotor Vigilance Test (PVT) that is highly sensitive to sleep deprivation. To relate PVT performance to time-varying and sparsely-measured environmental, operational, and psychological covariates, we propose a n ensemble prediction model comprising of linear mixed effects regression, random forest, and functional concurrent regression models. An extensive cross-validation procedure reveals that this ensemble outperforms any one of its components alone. We also discover that a participant’s past performance, reported fatigue and stress, and temperature and radiation exposure were among the most important variables associated with alertness. This method is broadly applicable to environmental studies where the main goal is accurate, individualized prediction involving a mixture of person-level traits and irregularly measured time series.

Danni Tu↗
NASA NTRS2025

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗
NASA NTRS2018

Toward Integrated Seasonal Predictions of Land and Ocean Carbon Flux: Lessons from the 2015-16 El Nino

Seasonal forecasts made by coupled atmosphere-ocean general circulation models (GCMs) are increasingly able to provide skillful forecasts of climate anomalies. At some centers, the capabilities of these models are being expanded to represent carbon-climate feedbacks including ocean biogeochemistry (OB), terrestrial biosphere (TB) interactions, and fires. These advances raise the question of whether such models can support skillful forecasts of carbon fluxes.Here, we examine whether land and ocean carbon flux anomalies associated with the 2015-16 El Nino could have been predicted months in advance. This El Nino was noteworthy for the magnitude of the ocean temperature perturbation, the skill with which this perturbation was predicted, and the extensive satellite observations that can be used to track its impact. We explore this topic using NASA's Goddard Earth Observing System (GEOS) model, which routinely produces an ensemble of seasonal climate forecasts, and a suite of offline dynamical and statistical models that estimate carbon flux processes. Using GEOS forecast fields from 2015-16 to force flux model hindcasts shows that these models are able to reproduce significant features observed by satellites. Specifically, OB hindcasts are able to predict anomalies in chlorophyll distributions with lead times of 3-4 months. The ability of TB hindcasts to reproduce NDVI anomalies is driven by the skill of the climate forecast, which is greatest at short lead times over tropical landmasses. Statistical fire forecasts driven by ocean climate indices are able to predict burned area in the tropics with lead times of 3-12 months. We also integrate the ocean and land hindcast fluxes into the GEOS GCM to examine the magnitude of the atmospheric carbon dioxide anomaly and compare with satellite and ground-based observations.While seasonal forecasting remains an active area of research, these results demonstrate that forecasts of carbon flux processes can support a variety of applications, potentially allowing scientists to understand carbon-climate feedbacks as they happen and to capitalize on more flexible satellite technologies that allow areas of interest to be targeted with lead times of weeks to months. We also provide a first glimpse at the spring 2019 carbon forecast using the GEOS-based forecasting system.

Ott, L.↗
NASA NTRS2022

Ultraviolet Spectropolarimetry Conservative and Nonconservative Mass Transfer in OB Interacting Binaries

The current consensus is that at least half of the OB stars are formed in binary or multiple star systems. The evolution of OB stars is greatly influenced by whether the stars begin as close binaries, and the evolution of the binary systems depend on whether the mass transfer is conservative or nonconservative. FUV/NUV spectropolarimetry is poised to answer the latter question. This paper discusses how the Polstar spectropolarimetry mission can characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase, and quantify its amount. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 10 −3 precision or better. The spectroscopy will identify absorption by mass streams and other plasmas seen in projection against the stellar disk as a function of orbital phase, as well as scattering from extended splash structures, including jets. The polarimetry tracks the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Special attention will be given to the epochs of enhanced systemic mass loss that have been identified from IUE observations (pre-mass reversal and tangential gas stream impact). We show how the history of systemic mass and angular momentum loss/gain episodes can be inferred via ensemble evolution through the r–q diagram. Combining the above elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass and ultimate supernova of a large number of massive stars found in binaries at close enough separation to undergo interaction.

Geraldine J. Peters↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗