Engineering topics
Publications and source records attributed to Koster, Randal.
This is an overview of forecast skill related to land surface initialization in prediction systems. It is provided as "general interest" to the USCLIVAR PPAI (Predictability, Predictions, and Applications Interface) community.
Recently NASA's Global Modeling and Assimilation Office (GMAO) has developed a new Subseasonal to Seasonal Prediction system Version 3 (GEOS-S2S-3). This upgrade replaces the GEOS-S2S-2 which is NASA's current contribution to the North American Multi-Model Experiment seasonal prediction project (Kirtman et al., 2014). The main improvements for our S2S-3 system include 1) a higher resolution MOM5 (Griffies et al., 2005) ocean model (now 0.25o x 0.25o x 50 layers), 2) an improved atmospheric/ocean interface layer (Akella and Suarez, 2018), and 3) assimilation of a long-track satellite salinity into the ocean model (Hackert et al, 2019). Atmospheric forcing is provided by the NASA MERRA-2 reanalysis (Gelaro et al., 2017). Initialization for the ocean relies on the GMAO ocean reanalysis system which assimilates all available in situ temperature and salinity, satellite sea surface salinity, and sea level using the Local Ensemble Transform Kalman Filter (LETKF) implementation of (Penny et al., 2013) on a 5 day assimilation cycle with 20 fixed ensemble members.In this presentation, we will authenticate our new S2S-3 ocean reanalysis using standard GODAE validation metrics. For example, we will compare gridded fields of mean and standard deviation of the ocean reanalysis versus observed fields. We will show correlation/RMS of model versus observations and temperature and salinity mean profiles for the various basins and latitude bands. Basin-scale volume transports, such as the Atlantic Meridional Overturning Circulation and the Indonesian Throughflow will be validated. Equatorial ocean waves will be compared by decomposing sea level into Kelvin and Rossby components. For each of these metrics, we plan to validate the results and then compare our new S2S-3 against the current production version, S2S-2. Finally, we will compare 9-month seasonal forecasts initialized from these two systems for the tropical Pacific NINO3.4 region over the period 1981-present.
The shape of the nonlinear relationship between evapotranspiration and soil moisture (the "ET-W relationship") helps control the evolution of soil moisture with time. Together, the shape of the relationship and the magnitude of the soil moisture anomaly at the beginning of a subseasonal forecast help determine whether a given anomaly will still be present at subseasonal leads, allowing it to contribute to skill in subseasonal temperature and precipitation prediction at those leads. In this study we examine subseasonal prediction in the context of soil moisture initialization using a suite of forecasts performed with the NASA GEOS seasonal forecast system. Large soil moisture anomalies are in fact found to be harbingers of increased skill in the subseasonal forecasts. Furthermore, accounting explicitly for the nonlinear shape of the ET-W relationship improves our ability to quantity the increase in forecast reliability associated with soil moisture initialization.
The moisture held within the top meter or two of soil is a very tiny fraction (less than 0.01) of the Earthstotal water (Eagleson, 1970). Nevertheless, its presence at the interface of the land and atmosphere gives it inordinateimportance in the context of climate variability. Simply put, soil moisture variations can help determine meteorologicalvariations. Consider, for example, an anomalously high evapotranspiration rate induced by a high soil moisture content.The high evapotranspiration can lead to an anomalously cooled land surface and thus cooler air temperatures(Seneviratne et al. 2010), and it can also lead to modifications in the evolution of the boundary layer, with concomitantimpacts on the generation of convective rainfall (Betts et al., 1994).Given this potential for feedback on the atmosphere,soil moisture is particularly important in the context of prediction. Atmospheric physics and dynamics are fast andcontribute little to the lifetime of an atmospheric anomaly, as reflected in the well-known time scale of typical weatherforecaststhese rely heavily on atmospheric initialization and are valid for only a week or so. Soil moisture processes, incontrast, are relatively slow. A soil moisture anomaly may persist for a month or more and is thus predictable at such timescales. As a result, and because (as noted above) the atmosphere may respond in a predictable way to a given soilmoisture anomaly, aspects of the atmosphere may also be predictable at the monthly time scale (NRC, 2010).
Spatial correlation structures can describe the degree to which soil moisture at a specified location co-varies in time with that at other, remote locations. Using four years of warm season SMAP Level 2 near-surface soil moisture data, we compute these spatial correlation structures for points across North America. The character of these structures is seen to differ geographically; the structures found for the west-central US, for example, are significantly more spatially extensive. We then demonstrate how these structures can potentially be used to reconstruct soil moisture fields during the pre-SMAP era. In this exercise, we consider as "truth" the soil moistures produced in a long-term offline land surface model simulation (1980-2014) that utilizes precipitation forcing based on a high density of precipitation gauges. Then, for a given location within the continent, we construct an"estimated" soil moisture time series based solely on historical soil moisture information simulated at least 300 km distant from the location, using the SMAP-based spatial correlation structures to determine how to make best use of the remote information. The reconstructed soil moistures are found to have significant skill relative to the assumed truth, suggesting that the same approach, when applied in areas of low rain gauge density (i.e., in areas for which historically simulated soil moistures are necessarily inaccurate), could provide useful historical soil moisture estimates through the SMAP-guided extraction of relevant information from neighboring gauged regions.
Flash droughts uncharacteristically rapid dryings of the land system are naturally associated with extreme precipitation deficits. Such precipitation deficits, however, do not tell the whole story. Land surface drying can be exacerbated by anomalously high evapotranspiration (ET) rates driven, for example, by anomalously high temperatures (e.g., during heatwaves) and/or anomalously high incoming radiation (e.g., from reduced cloudiness). In this study the relative contributions of precipitation and ET anomalies to flash drought generation are quantified through the analysis of diagnostic fields contained within the MERRA-2 reanalysis product. Unique to the overall approach is the explicit treatment of soil moisture impacts on ET an ET anomaly that is negative relative to the local long-term climatological mean is still considered positive in terms of its contribution to a flash drought, if it is high for the concurrent value of soil moisture. Maps produced in the analysis show the fraction of flash drought production stemming specifically from ET anomalies and the degree to which these ET anomalies are related to temperature and radiation anomalies. ET is found to have a large impact, for example, on flash drought production in the east-central US and in parts of Russia known from past studies to be prone to heatwave-related drought.
A new version of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was recently released by the NASA/Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution than the previous, (approximately 1/2 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter, and now includes the assimilation of along-track sea surface height. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a series of retrospective seasonal forecasts. Results show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. Analysis of the ensemble spread shows improvements relative to the previous system, including generally better reliability. The precipitation forecast skill is comparable to previous S2S systems, and the only tradeoff is an increased "double ITCZ", which is expected as we go to higher atmospheric resolution.
In this talk, I presented my current work at GMAO about the carbon cycle research. This includes (1) summary of our recently published paper about the impact of atmospheric CO2 variability on the global land carbon fluxes, and (2) our ongoing AGCM study with fully coupled carbon-water-energy cycles between the land and the atmosphere.
Tropical-extratropical teleconnections are considered key to advancing subseasonal prediction. The Madden Julian oscillation (MJO), characterized by large scale convective envelopes propagating along the tropical Indo-Pacific sector, is known to modulate midlatitude circulation and associated weather patterns. Although there is a general consensus on the MJO's influence on the midlatitude circulation, which is thought to be due to modulations of the North Atlantic Oscillation (NAO) and the Pacific North America (PNA) pattern, relatively less is known about the predictability of these teleconnection signals in dynamical forecast models. The composite evolution of the midlatitude circulation anomalies and associated wave train structure as the delayed response to tropical heating are reported in many studies that have examined reanalyses and long model simulations. However, it is yet to be determined whether they lend any beneficial subseasonal forecast skill, especially to weekly mean surface temperature and precipitation over North America. Investigating useful predictable signals from the MJO teleconnections is also complicated by the fact that the MJO is a moving heat source with an approximate periodicity of 30-60 days, and that the structure and amplitude of the midlatitude response can be sensitive to the longitudinal positioning of the heating anomaly as well as the propagation speed of the MJO. The objective of this study is to investigate the impact of MJO teleconnections on forecast accuracy at 2-3 week lead over North America, with an emphasis on the above-mentioned lesser known aspects of these teleconnections. To this end, we utilize a suite of subseasonal reforecasts performed with the latest NASA GEOS-5 seasonal-to-subseasonal (S2S) system. These reforecasts were performed as part of the NOAA SubX project, wherein the NASA GEOS-5 atmosphere-ocean coupled model was run at degree horizontal resolution, initialized every 5 days for the period 1999-2016. The GEOS-5 model shows skillful predictions of the MJO, with the correlation coefficient based on the real-time multivariate MJO (RMM) index staying at or above 0.5 up to forecast lead 26-36 days. The system is thus a useful tool for investigating MJO teleconnection processes.
By providing global measurements of near-surface soil moisture (down to about 5 cm) with unprecedented accuracy, the Soil Moisture Active/Passive (SMAP) satellite mission has opened the door to new and (in my opinion) exciting hydrological science. In this seminar, I present the results of a recent series of analyses performed with SMAP soil moisture data, covering a wide range of topics: (a) the characterization of the dynamics of near-surface soil moisture, with implications for forecasting soil moisture days into the future; (b) the multi-faceted character of the SMAP data, in the sense that different, established analysis approaches can extract information from the data that is largely (and perhaps unexpectedly) complementary; and (c) the interpretation of the data in the context of large-scale water fluxes. This final analysis is particularly exciting to me because it shows that, once the relevant algorithms are calibrated, precipitation and streamflow rates in hydrological basins can be estimated from the SMAP data alone - a reflection of the fact that the near-surface soil is a critical gateway between the atmospheric and subsurface branches of the hydrological cycle.
The SMAP estimates of rainfall and streamflow are not perfect, but they do contain relevant information. At the very least, they should prove useful for constraining, or otherwise contributing to, rainfall and streamflow estimates obtained with more conventional approaches.
In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System, current users and products, and methods for validation and evaluation of the system. Methods for evaluation include baseline evaluations metrics, the ability to simulate key modes of variability, and evaluation of new development areas.
In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System with a focus on the computing time and resources and actual time it takes to complete a full set of hindcasts. The goal is to come up with some solutions to allow us to run more ensemble members for the next version of the system which will be higher resolution and take many more resources.
In this study, we examine the prediction skill and predictability of the Madden Julian Oscillation (MJO) in a recent version of the NASA GEOS-5 atmosphere-ocean coupled model run at at 1/2 degree horizontal resolution. The results are based on a suite of hindcasts produced as part of the NOAA SubX project, consisting of seven ensemble members initialized every 5 days for the period 1999-2015. The atmospheric initial conditions were taken from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and the ocean and the sea ice were taken from a GMAO ocean analysis. The land states were initialized from the MERRA-2 land output, which is based on observation-corrected precipitation fields. We investigated the MJO prediction skill in terms of the bivariate correlation coefficient for the real-time multivariate MJO (RMM) indices. The correlation coefficient stays at or above 0.5 out to forecast lead times of 26-36 days, with a pronounced increase in skill for forecasts initialized from phase 3, when the MJO convective anomaly is located in the central tropical Indian Ocean. A corresponding estimate of the upper limit of the predictability is calculated by considering a single ensemble member as the truth and verifying the ensemble mean of the remaining members against that. The predictability estimates fall between 35-37 days (taken as forecast lead when the correlation reaches 0.5) and are rather insensitive to the initial MJO phase. The model shows slightly higher skill when the initial conditions contain strong MJO events compared to weak events, although the difference in skill is evident only from lead 1 to 20. Similar to other models, the RMM-index-based skill arises mostly from the circulation components of the index. The skill of the convective component of the index drops to 0.5 by day 20 as opposed to day 30 for circulation fields. The propagation of the MJO anomalies over the Maritime Continent does not appear problematic in the GEOS-5 hindcasts implying that the Maritime Continent predictability barrier may not be a major concern in this model. Finally, the MJO prediction skill in this version of GEOS-5 is superior to that of the current seasonal prediction system at the GMAO; this could be partly attributed to a slightly better representation of the MJO in the free running version of this model and partly to the improved atmospheric initialization from MERRA-2.
Over the past several decades, a tremendous effort has been made to improve model performance in the simulation of the climate system. The cold or warm sea surface temperature (SST) bias in the tropics is still a problem common to most coupled ocean atmosphere general circulation models (CGCMs). The precipitation biases in CGCMs are also accompanied by SST and surface wind biases. The deficiencies and biases over the equatorial oceans through their influence on the Walker circulation likely contribute the precipitation biases over land surfaces. In this study, we introduce an approach in the CGCM modeling to correct model biases. This approach utilizes the history of the model's short-term forecasting errors and their seasonal dependence to modify model's tendency term and to minimize its climate drift. The study shows that such an approach removes most of model climate biases. A number of other aspects of the model simulation (e.g. extratropical transient activities) are also improved considerably due to the imposed pre-processed initial 3-hour model drift corrections. Because many regional biases in the GEOS-5 CGCM are common amongst other current models, our approaches and findings are applicable to these other models as well.
A new version of the modeling and analysis system used to produce sub-seasonal to seasonal forecasts has just been released by the NASA Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution (approximately 12 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a free-running coupled simulation with the new system and results from a series of retrospective seasonal forecasts. Results from retrospective forecasts show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. The precipitation forecast skill is comparable to previous S2S systems, and the only trade off is an increased double ITCZ, which is expected as we go to higher atmospheric resolution.
In MERRA-2, observed precipitation is inserted in place of model-generated precipitation at the land surface. The use of observed precipitation was originally developed for MERRA-Land(a land-only replay of MERRA with model-generated precipitation replaced with observations).Previously shown that the land hydrology in MERRA-2 and MERRA-Land is better than MERRA. We test whether the improved land surface hydrology in MERRA-2 leads to the expected improvements in the land surface energy fluxes and 2 m air temperatures (T2m).