Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate model ensemble”

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 217 records · Page 12

Mesoscale Convective Systems During SCSMEX: Simulations with a Regional Climate Model and a Cloud-Resolving Model

The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convection and circulation changes associated with monsoons over the South China Sea region. SCSMEX also provided precipitation derived from atmospheric budgets and comparison to those obtained from the Tropical Rainfall Measuring Mission (TRMM). In this paper, a regional scale model (with grid size of 20 km) and Goddard Cumulus Ensemble (GCE) model (with 1 km grid size) are used to perform multi-day integration to understand the precipitation processes associated with the summer monsoon over Southeast Asia and southern China. The regional climate model is used to understand the soil-precipitation interaction and feedback associated with a flood event that occurred in and around China's Yantz River during SCSMEX Sensitivity tests on various land surface models, sea surface temperature (SST) variations, and cloud processes are performed to understand the precipitation processes associated with the onset of the monsoon over the S. China Sea during SCSMEX. These tests have indicated that the land surface model has a major impact on the circulation over the S. China Sea. Cloud processes can effect the precipitation pattern while SST variation can effect the precipitation amounts over both land and ocean. The exact location (region) of the flooding can be effected by the soil-rainfall feedback. The GCE-model results captured many observed precipitation characteristics because it used a fine grid size. For example, the model simulated rainfall temporal variation compared quite well to the sounding-estimated rainfall. The results show there are more latent heat fluxes prior to the onset of the monsoon. However, more rainfall was simulated after the onset of the monsoon. This modeling study indicates the latent heat fluxes (or evaporation) have more of an impact on precipitation processes and rainfall in the regional climate model simulations than in the cloud-resolving model simulations. Research is underway to determine if the difference in the grid sizes or the moist processes used in these two models is responsible for the differing influence of surface fluxes an precipitation processes.

Tao, W.-K.↗

Tropical Cyclones, Hurricanes, and Climate: NASA's Global Cloud-Scale Simulations and New Observations that Characterize the Lifecycle of Hurricanes

One of the primary interests of Global Change research is the impact of climate changes and climate variability on extreme weather events, such as intense tropical storms and hurricanes. Atmospheric climate models run at resolutions of global weather models have been used to study the impact of climate variability, as seen in sea surface temperatures, on the frequency and intensity of tropical cyclones. NASA's Goddard Earth Observing System Model, version 5 (GEOS-5) in ensembles run at 50 km resolution has been able to reproduce the interannual variations of tropical cyclone frequency seen in nature. This, and other global models, have found it much more difficult to reproduce the interannual changes in intensity, a result that reflects the inability of the models to simulate the intensities of the most extreme storms. Better representation of the structures of cyclones requires much higher resolution models. Such improved representation is also fundamental to making best use of satellite observations. In collaboration with NOAA's Geophysical Fluid Dynamics Laboratory, GEOS-5 now has the capability of running at much higher resolution to better represent cloud-scale resolutions. Global simulations at cloud-permitting resolutions (10- to 3.5-km) allows for the development of realistic tropical cyclones from tropical storm 119 km/hr winds) to category 5 (>249km1hr winds) intensities. GEOS-5 has produced realistic rain-band and eye-wall structures in tropical cyclones that can be directly analyzed against satellite observations. For the first time a global climate model is capable of representing realistic intensity and track variability on a seasonal scale across basins. GEOS-5 is also used in assimilation mode to test the impact of NASA's observations on tropical cyclone forecasts. One such test, for tropical cyclone Nargis in the Indian Ocean in May 2008, showed that observations from Atmospheric Infrared Sounder (AIRS) and the Advanced Microwave Sounding Unit (AMSU-A) on Aqua substantially reduced forecast track errors. Tropical cyclones in the northern Indian Ocean pose serious challenges to operational weather forecasting systems, partly due to their shorter lifespan and more erratic track, compared to those in the Atlantic and the Pacific. SA is also bringing several state of the art instruments in recent field campaigns to peer under the clouds and study the inner workings of the tropical storms. With the Genesis and Rapid Intensification Processes (GRIP) experiment, a NASA Earth science field experiment in 2010 that includes the Global Hawk Unmanned Airborne System (UAS) configured with a suite of in situ and remote sensing instruments that are observing and characterizing the lifecycle of hurricanes, we expect significant improvement in our understanding of the track and intensification processes with the assimilation of the satellite and field campaign observations of meteorological parameters in the numerical prediction models.

Putman, William M.↗

Quantifying the role of ocean coupling in Arctic amplification and sea-ice loss over the 21st century

The enhanced warming of the Arctic, relative to other parts of the Earth, a phenomenon known as Arctic amplification, is one of the most striking features of climate change, and has important climatic impacts for the entire Northern Hemisphere. Several mechanisms are believed to be responsible for Arctic amplification; however, a quantitative understanding of their relative importance is still missing. Here, using ensembles of model integrations, we quantify the contribution of ocean coupling, both its thermodynamic and dynamic components, to Arctic amplification over the 20th and 21st centuries. We show that ocean coupling accounts for ~80% of the amplification by 2100. In particular, we show that thermodynamic coupling is responsible for future amplification and sea-ice loss as it overcomes the effect of dynamic coupling which reduces the amplification and sea-ice loss by ~35%. Our results demonstrate the utility of targeted numerical experiments to quantify the role of specific mechanisms in Arctic amplification, for better constraining climate projections.

Ocean heat transport↗

A Note on the Relationship between Temperature and Water Vapor in Quasi-Equilibrium and Climate States

An ideal and simple formulation is successfully derived that well represents a quasi-linear relationship found between the domain-averaged water vapor, q (mm), and temperature, T (K), fields obtained from a series of quasi-equilibrium (long-term) simulations for the Tropics using the two-dimensional Goddard Cumulus Ensemble (GCE) model. Earlier model work showed that the forced maintenance of two different wind profiles in the Tropics leads to two different equilibrium states. Investigating this finding required investigation of the slope of the moisture-temperature relations, which turns out to be linear in the Tropics. The extra-tropical climate equilibriums become more complex, but insight on modeling sensitivity can be obtained by linear stepwise regression of the integrated temperature and humidity. A globally curvilinear moisture-temperature distribution, similar to the famous Clausius-Clapeyron curve (i.e., saturated water vapor pressure versus temperature), is then found in this study. Such a genuine finding clarifies that the dynamics are crucial to the climate (shown in the earlier work) but the thermodynamics adjust. The range of validity of this result is further examined herein. The GCE-modeled tropical domain-averaged q and T fields form a linearly-regressed "q-T" slope that genuinely resides within an ideal range of slopes obtained from the aforementioned formulation. A quantity (denoted as dC2/dC1) representing the derivative between the static energy densities due to temperature (C2) and water vapor (C1) for various quasi-equilibrium states can also be obtained. A dC2/dC1 value near unity obtained for the GCE-modeled tropical simulations implies that the static energy densities due to moisture and temperature only differ by a pure constant for various equilibrium states. An overall q-T relation also including extra-tropical regions is, however, found to have a curvilinear relationship. Accordingly, warm/moist regions favor change in water vapor faster than temperature, while cold/dry regions favor an increase in temperature quicker than water vapor.

Shie, C.-L.↗

The Subseasonal Experiment (SubX): A Multi-Model Subseasonal Prediction Experiment

SubX is a multi-model subseasonal prediction experiment designed around operational requirements with the goal of improving subseasonal forecasts. Seven global models have produced seventeen years of retrospective (re-) forecasts and more than a year of weekly real-time forecasts. The re-forecasts and forecasts are archived at the Data Library of the International Research Institute for Climate and Society, Columbia University, providing a comprehensive database for research on subseasonal to seasonal predictability and predictions. The SubX models show skill for temperature and precipitation three weeks ahead of time in specific regions. The SubX multi-model ensemble mean is more skillful than any individual model overall. Skill in simulating the Madden-Julian Oscillation (MJO) and the North Atlantic Oscillation (NAO), two sources of subseasonal predictability, is also evaluated with skillful predictions of the MJO four weeks in advance and of the NAO 2 weeks in advance. SubX is also able to make useful contributions to operational forecast guidance at the Climate Prediction Center. Additionally, SubX provides information on the potential for extreme precipitation associated with tropical cyclones which can help emergency management and aid organizations to plan for disasters. (Capsule Summary) A research to operations project in service of developing better operational subseasonal forecasts.

Precipitation↗

Assimilation of Terrestrial Water Storage from GRACE in a Snow-Dominated Basin

Terrestrial water storage (TWS) information derived from Gravity Recovery and Climate Experiment (GRACE) measurements is assimilated into a land surface model over the Mackenzie River basin located in northwest Canada. Assimilation is conducted using an ensemble Kalman smoother (EnKS). Model estimates with and without assimilation are compared against independent observational data sets of snow water equivalent (SWE) and runoff. For SWE, modest improvements in mean difference (MD) and root mean squared difference (RMSD) are achieved as a result of the assimilation. No significant differences in temporal correlations of SWE resulted. Runoff statistics of MD remain relatively unchanged while RMSD statistics, in general, are improved in most of the sub-basins. Temporal correlations are degraded within the most upstream sub-basin, but are, in general, improved at the downstream locations, which are more representative of an integrated basin response. GRACE assimilation using an EnKS offers improvements in hydrologic state/flux estimation, though comparisons with observed runoff would be enhanced by the use of river routing and lake storage routines within the prognostic land surface model. Further, GRACE hydrology products would benefit from the inclusion of better constrained models of post-glacial rebound, which significantly affects GRACE estimates of interannual hydrologic variability in the Mackenzie River basin.

Forman, Barton A.↗

ISMIP6 Antarctica: A Multi-Model Ensemble of the Antarctic Ice Sheet Evolution Over the 21st Century

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between -7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica ass change varies between -6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

Hélène Seroussi↗

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

Comparative Analysis of Upper Ocean Heat Content Variability from Ensemble Operational Ocean Analyses

Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyzed

Xue, Yan↗

Processing Tomato Production Is Expected to Decrease By 2050 Due to the Projected Increase in Temperature

The global processing tomato production is concentrated in a small number of regions where climate change will have a significant impact on the future supply. Process-based tomato models project that the production in the three main producing countries (the United States, Italy, and China, representing 65% of global production) will decrease 6% by 2050, compared to the baseline period of 1980-2009. The predicted reduction in processing tomato production is due to a simulated increase in air temperature. Under an ensemble of projected climate scenarios, California and Italy might not be able to sustain the current levels of processing tomato production due to water resources constraints. Cooler producing regions, such as China and the northern parts of California, stand to improve their competitive advantage. The projected environmental changes indicate that the main growing regions of processing tomatoes might change in the coming decades.

Tomato production↗

Cold Season Performance of the NU-WRF Regional Climate Model in the Great Lakes Region

As Earth’s largest collection of freshwater, the Laurentian Great Lakes have enormous ecological and socio-economic value. Their basin has become a regional hotspot of climatic and limnological change, potentially threatening its vital natural resources. Consequentially, there is a need to assess the current state of climate models regarding their performance across the Great Lakes region and develop the next generation of high-resolution regional climate models to address complex limnological processes and lake-atmosphere interactions. In response to this need, the current paper focuses on the generation and analysis of a 20-member ensemble of 3-km National Aeronautics and Space Administration (NASA)-Unified Weather Research and Forecasting (NU-WRF) simulations for the 2014-2015 cold season. The study aims to identify the model’s strengths and weaknesses; optimal configuration for the region; and the impacts of different physics parameterizations, coupling to a 1D lake model, time-variant lake-surface temperatures, and spectral nudging. Several key biases are identified in the cold-season simulations for the Great Lakes region, including an atmospheric cold bias that is amplified by coupling to a 1D lake model but diminished by applying the Community Atmosphere Model radiation scheme and Morrison microphysics scheme; an excess precipitation bias; anomalously early initiation of fall lake turnover and subsequent cold lake bias; excessive and overly persistent lake ice cover; and insufficient evaporation over Lakes Superior and Huron. The research team is currently addressing these key limitations by coupling NU-WRF to a 3D lake model in support of the next generation of regional climate models for the critical Great Lakes Basin.

Michael Notaro↗

Assimilation of Satellite-Derived Skin Temperature Observations into Land Surface Models

Land surface (or "skin") temperature (LST) lies at the heart of the surface energy balance and is a key variable in weather and climate models. Here we assimilate LST retrievals from the International Satellite Cloud Climatology Project (ISCCP) into the Noah and Catchment (CLSM) land surface models using an ensemble-based, off-line land data assimilation system. LST is described very differently in the two models. A priori scaling and dynamic bias estimation approaches are applied because satellite and model LST typically exhibit different mean values and variability. Performance is measured against 27 months of in situ measurements from the Coordinated Energy and Water Cycle Observations Project at 48 stations. LST estimates from Noah and CLSM without data assimilation ("open loop") are comparable to each other and superior to that of ISCCP retrievals. For LST, RMSE values are 4.9 K (CLSM), 5.6 K (Noah), and 7.6 K (ISCCP), and anomaly correlation coefficients (R) are 0.62 (CLSM), 0.61 (Noah), and 0.52 (ISCCP). Assimilation of ISCCP retrievals provides modest yet statistically significant improvements (over open loop) of up to 0.7 K in RMSE and 0.05 in anomaly R. The skill of surface turbulent flux estimates from the assimilation integrations is essentially identical to the corresponding open loop skill. Noah assimilation estimates of ground heat flux, however, can be significantly worse than open loop estimates. Provided the assimilation system is properly adapted to each land model, the benefits from the assimilation of LST retrievals are comparable for both models.

Reichle, Rolf H.↗

Solar Rotational Modulations of Spectral Irradiance and Correlations with the Variability of Total Solar Irradiance

Aims: We characterize the solar rotational modulations of spectral solar irradiance (SSI) and compare them with the corresponding changes of total solar irradiance (TSI). Solar rotational modulations of TSI and SSI at wavelengths between 120 and 1600 nm are identified over one hundred Carrington rotational cycles during 2003-2013. Methods: The SORCE (Solar Radiation and Climate Experiment) and TIMED (Thermosphere Ionosphere Mesosphere Energetics and Dynamics)/SEE (Solar EUV Experiment) measured and SATIRE-S modeled solar irradiances are analyzed using the EEMD (Ensemble Empirical Mode Decomposition) method to determine the phase and amplitude of 27-day solar rotational variation in TSI and SSI. Results: The mode decomposition clearly identifies 27-day solar rotational variations in SSI between 120 and 1600 nm, and there is a robust wavelength dependence in the phase of the rotational mode relative to that of TSI. The rotational modes of visible (VIS) and near infrared (NIR) are in phase with the mode of TSI, but the phase of the rotational mode of ultraviolet (UV) exhibits differences from that of TSI. While it is questionable that the VIS to NIR portion of the solar spectrum has yet been observed with sufficient accuracy and precision to determine the 11-year solar cycle variations, the temporal variations over one hundred cycles of 27-day solar rotation, independent of the two solar cycles in which they are embedded, show distinct solar rotational modulations at each wavelength.

Spectral Solar Irradiance↗

Scaling Observation Error for Optimal Assimilation of CCI SST Data into a Regional HYCOM EnOI System

South Africa currently possesses no operational ocean forecasting system for the purpose of predicting ocean state variables including temperature,salinity and velocity. Substantial initial efforts towards this goal have been made and resulted in a system using a regional Hybrid Coordinate Ocean Model (HYCOM) along with the Ensemble Optimal Interpolation (EnOI)assimilation scheme. Assimilating only sea surface temperature (SST) observations from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) product into the system resulted in a degraded forecast. Aiming to address this, Climate Change Initiative (CCI) SSTs are assimilated into the system in an effort to improve the forecast skill. Observation errors in the assimilated product are used in the EnOI to determine whether more confidence should be placed in the model or observations in producing the analysis, but overconfidence in observations can shock the model and result in failure. To tweak the impact of the assimilation, a scaling factor is applied in the assimilation code. A scaling factor of 25 was found to produce a favourable result with lowest mean root mean square error (RMSE;1.098°C) between the model and observations over time. Postulating the error to be overconfident, a floor value is introduced in order to set a minimum value for the observation error thereby reducing confidence in the observations. These experiments fared less favourably with a floor value of 0.5 and a scaling factor of 15 producing the best mean RMSE (1.118°C).

Luyt, Hermann↗

Evaluating the Global Water and Energy Cycles in Reanalyses Using Ensemble Spread

Reanalyses provide continuous observation based global data over weather and climate time scales through observational assimilation. Ensemble forecasts are a component of modern data assimilation systems, providing information to the analysis process. Storing this spread over the duration of a climate reanalysis will provide insight to the variations of the model that underlies the reanalysis data. Global water and energy cycles (WECs) are a critical aspect of both weather and climate, yet there are some large gaps in the observation of key components as well as their representation in models and reanalyses. NASA’s Water and Energy cycle Studies (NEWS) program is developing a corrected and balanced budget data set that is closed at monthly and regional scales during the 21st century, so far. In this paper, we will evaluate reanalysis WEC components against the initial NEWS collection of observation, but integrate the analysis spread from the reanalysis data assimilation as a new perspective on the evaluation. The initial evaluation will cover a few more recent years using ERA5 10 member spread. As the GMAO’s production of the Modern-Era Retrospective analysis for Research and Applications 21st Century (MERRA-21C) progresses, the effort will include the MERRA-21C 32 member spread. Ultimately, we will explore the potential of the ensemble spread as a measure of uncertainty in the reanalysis.

Michael Bosilovich↗

Assimilation of GRACE Terrestrial Water Storage into a Land Surface Model: Evaluation 1 and Potential Value for Drought Monitoring in Western and Central Europe

A land surface model s ability to simulate states (e.g., soil moisture) and fluxes (e.g., runoff) is limited by uncertainties in meteorological forcing and parameter inputs as well as inadequacies in model physics. In this study, anomalies of terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission were assimilated into the NASA Catchment land surface model in western and central Europe for a 7-year period, using a previously developed ensemble Kalman smoother. GRACE data assimilation led to improved runoff correlations with gauge data in 17 out of 18 hydrological basins, even in basins smaller than the effective resolution of GRACE. Improvements in root zone soil moisture were less conclusive, partly due to the shortness of the in situ data record. In addition to improving temporal correlations, GRACE data assimilation also reduced increasing trends in simulated monthly TWS and runoff associated with increasing rates of precipitation. GRACE assimilated root zone soil moisture and TWS fields exhibited significant changes in their dryness rankings relative to those without data assimilation, suggesting that GRACE data assimilation could have a substantial impact on drought monitoring. Signals of drought in GRACE TWS correlated well with MODIS Normalized Difference Vegetation Index (NDVI) data in most areas. Although they detected the same droughts during warm seasons, drought signatures in GRACE derived TWS exhibited greater persistence than those in NDVI throughout all seasons, in part due to limitations associated with the seasonality of vegetation.

Li, Bailing↗

Uncertainties in Evapotranspiration Estimates over West Africa

An evapotranspiration (ET) ensemble composed of 36 land surface model (LSM) experiments and four diagnostic datasets (GLEAM, ALEXI, MOD16, and FLUXNET) is used to investigate uncertainties in ET estimate over five climate regions in West Africa. Diagnostic ET datasets show lower uncertainty estimates and smaller seasonal variations than the LSM-based ET values, particularly in the humid climate regions. Overall, the impact of the choice of LSMs and meteorological forcing datasets on the modeled ET rates increases from north to south. The LSM formulations and parameters have the largest impact on ET in humid regions, contributing to 90% of the ET uncertainty estimates. Precipitation contributes to the ET uncertainty primarily in arid regions. The LSM-based ET estimates are sensitive to the uncertainty of net radiation in arid region and precipitation in humid region. This study serves as support for better determining water availability for agriculture and livelihoods in Africa with earth observations and land surface models.

land surface model↗