Search NASA⌕ Search

SEARCH · Search NASA

Results for “forecast skill”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Convective Weather Forecast Quality Metrics for Air Traffic Management Decision-Making

Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.

Chatterji, Gano B.↗
NASA NTRS2021

Description of the NASA GEOS Composition Forecast Modeling System GEOS-CF v1.0

The Goddard Earth Observing System composition forecast (GEOS-CF) system is a high-resolution (0.25 degree) global constituent prediction system from NASA’s Global Modeling and Assimilation Office (GMAO). GEOS-CF offers a new tool for atmospheric chemistry research, with the goal to supplement NASA’s broad range of space-based and in-situ observation sand to support flight campaign planning, support of satellite observations, and air quality research. GEOS-CF expands on the GEOS weather and aerosol modeling system by introducing the GEOS-Chem chemistry module to provide analyses and 5-day forecasts of atmospheric constituents including ozone (O3), carbon monoxide (CO), nitrogen dioxide (NO2), and fine particulate matter (PM2.5). The chemistry module integrated in GEOS-CF is identical to the offline GEOS-Chem model and readily benefits from the innovations provided by the GEOS-Chem community.Evaluation of GEOS-CF against satellite, ozone sonde and surface observations show realistic simulated concentrations of O3, NO2, and CO, with normalized mean biases of -0.1 to -0.3, normalized root mean square errors (NRMSE) between 0.1-0.4, and correlations between 0.3-0.8. Comparisons against surface observations highlight the successful representation of air pollutants under a variety of meteorological conditions, yet also highlight current limitations, such as an over prediction of summertime ozone over the Southeast United States. GEOS-CFv1.0 generally overestimates aerosols by 20-50% due to known issues in GEOS-Chem v12.0.1 that have been addressed in later versions.The 5-day hourly forecasts have skill scores comparable to the analysis. Model skills can be improved significantly by applying a bias-correction to the surface model output using a machine-learning approach.

GEOS-CF↗
NASA NTRS2021

Drought Relief and Reversal Over North America from 1500 to 2016

Season-to-season persistence of soil moisture drought varies across North America. Such inter-seasonal autocorrelation can have modest skill in forecasting future conditions several months in advance. Because robust instrumental observations of precipitation span less than 100 years, the temporal stability of the relationship between seasonal moisture anomalies is uncertain. The North American Seasonal Precipitation Atlas (NASPA) is a gridded network of separately reconstructed cool (December-April) and warm (May-July) season precipitation series and offers new insights on the intra-annual changes in drought for up to 2,000 years. Here, the NASPA precipitation reconstructions are rescaled to represent the long-term soil moisture balance during the cool season and three-month long atmospheric moisture during the warm season. These rescaled seasonal reconstructions are then used to quantify the frequency, magnitude, and spatial extent of cool season drought that was relieved or reversed during the following summer months. The adjusted seasonal reconstructions reproduce the general patterns of large-scale drought amelioration and termination in the instrumental record during the 20th century and are used to estimate relief and reversals for the most skillfully reconstructed past 500-years. Sub-continental to continental scale reversals of cool season drought in the following warm season have been rare, but the reconstructions display periods prior to the instrumental data of increased reversal probabilities for the Mid-Atlantic region and the US Southwest. Drought relief at the continental scale may arise in part from macroscale ocean-atmosphere processes, while the smaller scale regional reversals may reflect land-surface feedbacks and stochastic variability.

Drought↗
DOE OSTI2023

Optimizing Seasonal‐To‐Decadal Analog Forecasts With a Learned Spatially‐Weighted Mask

Abstract Seasonal‐to‐decadal climate prediction is crucial for decision‐making in a number of industries, but forecasts on these timescales have limited skill. Here, we develop a data‐driven method for selecting optimal analogs for seasonal‐to‐decadal analog forecasting. Using an interpretable neural network, we learn a spatially‐weighted mask that quantifies how important each grid point is for determining whether two climate states will evolve similarly. We show that analogs selected using this weighted mask provide more skillful forecasts than analogs that are selected using traditional spatially‐uniform methods. This method is tested on two prediction problems using the Max Planck Institute for Meteorology Grand Ensemble: multi‐year prediction of North Atlantic sea surface temperatures, and seasonal prediction of El Niño Southern Oscillation. This work demonstrates a methodical approach to selecting analogs that may be useful for improving seasonal‐to‐decadal forecasts and understanding their sources of skill.

Geology↗
NASA NTRS1974

Short-range forecasts with the GISS model of the global atmosphere

A nine-layer, primitive equation (PE) model of the global atmosphere developed at the Goddard Institute for Space Studies (GISS) has been used to generate six 48-hr forecasts during December 1972 and January 1973. Operational analyses north of 18 N and experimental global analyses made available by the National Meteorological Center (NMC), NOAA, were used as the initial conditions; the operational analyses were used to verify the forecasts at 12-hr intervals over the Northern Hemisphere north of 22 N. The combined analyses were used to verify the forecasts in the global domain. RMS errors of the sea-level pressure, 1000-mb heights, and vector geostrophic winds, and 500-mb heights and vector geostrophic winds indicate that the GISS forecasts have skill comparable to those made by operational PE models.

Druyan, L. M.↗
NASA NTRS1989

Dynamical predictability in a simple general circulation model - Average error growth

Average predictability and error growth in a simple realistic two-level general circulation model (GCM) were investigated using a series of Monte Carlo experiments for fixed external forcing (perpetual winter in the Northern Hemisphere). It was found that, for realistic initial errors, the dependence of the limit of dynamic predictability on total wavenumber was similar to that found for the ECMWF model for the 1980/1981 winter conditions, with the lowest wavenumbers showing significant skill for forecast ranges of more than 1 month. On the other hand, for very small amplitude errors distributed according to the climate spectrum, the total error growth was superexponential, reaching a maximum growth rate (2-day doubling time) in about 1 week. A simple empirical model of error variance, which involved two broad wavenumber bands and incorporating a 3/2 power saturation term, was found to provide an excellent fit to the GCM error growth behavior.

Schubert, Siegfried D.↗
NASA NTRS1990

Simulations of winter blocking episodes using observed sea surface temperatures

The formation of major Northern Hemisphere blockings during January 1979 is studied in numerical simulations of the global atmosphere with a sea surface temperature (SST) field updated with observations during the model integration. Both the standard 4 x 5 degrees latitude-longitude and high-resolution 2 x 2.5 degrees versions of the Goddard Laboratory for Atmospheres general circulation models (GCMs) are employed. The SST field is provided by a blend analysis of in situ and satellite-retrieved data. The simulations by the high-resolution GCM with a realistic SST field of the ocean surface are shown to be capable of producing two successive realistic major blockings in the Pacific and Atlantic through a one-month period. The skill in predicting the blocking formation as obtained in the simulations is due to improved skill in forecasting of ultralog waves. Although the results are encouraging, only one initial state is involved in this series of simulations; therefore, the results should not be generalized yet at this stage of investigation.

Kung, Ernest C.↗
NASA NTRS2009

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗
NASA NTRS2018

Local Versus Remote Contributions of Soil Moisture to Near-Surface Temperature Variability

Soil moisture variations have a straightforward impact on overlying air temperatures, wetter soils can induce higher evaporative cooling of the soil and thus, locally, cooler temperatures overall. Not known, however, is the degree to which soil moisture variations can affect remote air temperatures through their impact on the atmospheric circulation. In this talk we describe a two-pronged analysis that addresses this question. In the first segment, an extensive ensemble of NASA/GSFC GEOS-5 atmospheric model simulations is analyzed statistically to isolate and quantify the contributions of various soil moisture states, both local and remote, to the variability of air temperature at a given local site. In the second segment, the relevance of the derived statistical relationships is evaluated by applying them to observations-based data. Results from the second segment suggest that the GEOS-5-based relationships do, at least to first order, hold in nature and thus may provide some skill to forecasts of air temperature at subseasonal time scales, at least in certain regions.

Koster, R.↗
NASA NTRS2021

The 2019 Southern Hemisphere Stratospheric Polar Vortex Weakening and Its Impacts

This study offers an overview of the low-frequency (i.e., monthly to seasonal) evolution, dynamics, predictability, and surface impacts of a rare Southern Hemisphere (SH) stratospheric warming that occurred in austral spring 2019. Between late August to mid- September 2019, the stratospheric circumpolar westerly jet weakened rapidly, and Antarctic stratospheric temperatures rose dramatically. The deceleration of the vortex at 10 hPa was as drastic as that of the first ever observed major sudden stratospheric warming in the SH during 2002, while the mean Antarctic warming over the course of spring 2019 broke the previous record of 2002 by ~50% in the mid-stratosphere. This event was preceded by a poleward shift of the SH polar night jet in the uppermost stratosphere in early winter, which was then followed by record-strong planetary wave-one activity propagating upward from the troposphere in August that acted to dramatically weaken the polar vortex throughout the depth of the stratosphere. The weakened vortex winds and elevated temperatures moved downward to the surface from mid-October to December, promoting a record strong swing of the Southern Annular Mode (SAM) to its negative phase. This record-negative SAM appeared to be a primary driver of the extreme hot and dry conditions over subtropical eastern Australia that accompanied the severe wildfires that occurred in late spring 2019. State-of-the-art dynamical seasonal forecast systems skillfully predicted the significant vortex weakening of spring 2019 and subsequent development of negative SAM from as early as late July.

Antarctic Oscillation↗
NASA NTRS2014

Prediction of the Arctic Oscillation in Boreal Winter by Dynamical Seasonal Forecasting Systems

This study assesses the skill of boreal winter Arctic Oscillation (AO) predictions with state-of-the-art dynamical ensemble prediction systems (EPSs): GloSea4, CFSv2, GEOS-5, CanCM3, CanCM4, and CM2.1. Long-term reforecasts with the EPSs are used to evaluate how well they represent the AO and to assess the skill of both deterministic and probabilistic forecasts of the AO. The reforecasts reproduce the observed changes in the large-scale patterns of the Northern Hemispheric surface temperature, upper level wind, and precipitation associated with the different phases of the AO. The results demonstrate that most EPSs improve upon persistence skill scores for lead times up to 2 months in boreal winter, suggesting some potential for skillful prediction of the AO and its associated climate anomalies at seasonal time scales. It is also found that the skill of AO forecasts during the recent period (1997-2010) is higher than that of the earlier period (1983-1996).

Arctic Oscillation↗
NASA NTRS2009

An Anomaly Correlation Skill Score for the Evaluation of the Performance of Hyperspectral Infrared Sounders

With the availability of very accurate forecasts, the metric of accuracy alone for the evaluation of the performance of a retrieval system can produce misleading results. A useful characterization of the quality of a retrieval system and its potential to contribute to an improved weather forecast is its skill, which we define as the ability to make retrievals of geophysical parameters which are closer to the truth than the six hour forecast, when the truth differs significantly from the forecast. We illustrate retrieval skill using one day of AMSU and AIRS data with three different retrieval algorithms, which result in retrievals for more than 90% of the potential retrievals under clear and cloudy conditions. Two of the three algorithms have better than 1 K rms "RAOB quality" accuracy on the troposphere, but only one has skill between 900 and 100 mb. AIRS was launched on the EOS Aqua spacecraft in May 2002 into a 705 km polar sun-synchronous orbit with accurately maintained 1:30 PM ascending node. Essentially uninterrupted data are freely available since September 2002.

retrieval↗
DOE OSTI2026

Subseasonal Forecasting and MJO Teleconnections in Machine Learning Weather Prediction Models

Abstract In recent years, machine‐learning (ML) models trained on reanalysis data have rivaled physics‐based forecast models in terms of performance skill for global weather forecasting. With increased rollout stability, the question of how these models perform for subseasonal to seasonal (S2S, week 3–8) forecasting has emerged. In this study we run a large set of subseasonal hindcasts over 2004–2023 to evaluate two ML weather forecast models at the S2S time scale, SFNO‐HENS (Nvidia, fully ML) and NeuralGCM (Google Research, hybrid). Corresponding hindcasts from the European Centre for Medium‐Range Weather Forecasts (ECMWF) are used as a baseline for comparison to a physics‐based model. Because our focus is on predicting moisture transport over the Western United States between October and March, we evaluate the models' prediction skill for the Madden‐Julian Oscillation (MJO) and its associated teleconnections in the North Pacific. We find that both ML models are competitive with the ECWMF model, with comparable skill in predicting the North Pacific large‐scale circulation and the MJO at week 3 and beyond. Even though overall the mid‐latitude subseasonal prediction skill remains low, the ML models exhibit interesting behavior such as a realistic propagation of the MJO across the Maritime Continent and realistic teleconnections. A SFNO‐HENS sensitivity experiment with altered initial conditions in the tropics demonstrates the stability of the model, and it illustrates the capability of ML models to represent important physical processes of the atmosphere at the S2S time scale. Plain Language Summary Predicting weather patterns and precipitation a few weeks in advance (subseasonal time scale) is of great interest for stakeholders such as water managers in the Southwest United States (US), where arid conditions prevail. Subseasonal forecasts from traditional weather forecast models exhibit low skill in the region, limiting their applicability. Here we examine whether the recent breakthrough in weather forecasting made with machine learning/artificial intelligence models can translate to improved subseasonal forecasts. Recently‐developed machine learning models exhibit comparable skill to a state‐of‐the‐art physics‐based model for predicting weather patterns in the North Pacific/North America region, and associated moisture transport. The same applies to their skill in predicting the tropical pattern, the Madden‐Julian Oscillation, and its important remote perturbations over the midlatitude East Pacific and Southwest US. Additionally, a perturbation experiment carried out with one of the machine learning models illustrates their ability to not only predict the evolution of atmospheric fields, but also to learn and represent physical processes such as tropics‐extratropics Rossby wave propagation. Key Points Two machine learning weather forecast models exhibit state‐of‐the‐art prediction skill at the subseasonal time scale in the Pacific sector The models equal ECWMF in terms of Madden‐Julian oscillation (MJO) prediction skill, and they accurately predict the MJO propagation and associated teleconnections The two machine‐learning models represent key physical processes for subseasonal prediction, despite being trained for weather forecasting

Peings, Yannick↗
NASA NTRS2013

Prediction of the Arctic Oscillation in Boreal Winter by Dynamical Seasonal Forecasting Systems

This study assesses the prediction skill of the boreal winter Arctic Oscillation (AO) in the state-of-the-art dynamical ensemble prediction systems (EPSs): the UKMO GloSea4, the NCEP CFSv2, and the NASA GEOS-5. Long-term reforecasts made with the EPSs are used to evaluate representations of the AO, and to examine skill scores for the deterministic and probabilistic forecast of the AO index. The reforecasts reproduce the observed changes in the large-scale patterns of the Northern Hemispheric surface temperature, upper-level wind, and precipitation according to the AO phase. Results demonstrate that all EPSs have better prediction skill than the persistence prediction for lead times up to 3-month, suggesting a great potential for skillful prediction of the AO and the associated climate anomalies in seasonal time scale. It is also found that the deterministic and probabilistic forecast skill of the AO in the recent period (1997-2010) is higher than that in the earlier period (1983-1996).

Coupled Models↗
DOE OSTI2024

Challenges of COVID-19 Case Forecasting in the US, 2020–2021

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub ( https://covid19forecasthub.org ). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1–4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naïve baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2020

Evaluation of Global Fire Weather Database Reanalysis and Short-Term Forecast Products

Daily Fire Weather Index (FWI) System components calculated from the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), are compared to FWI calculations from a global network of weather stations over 2004–2018, and short-term, experimental (8 d) daily FWI forecasts are evaluated for their skill across the Terrestrial Ecoregions of the World for 2018. FWI components from MERRA-2 were, in general, biased low compared to station data, but this reflects a mix of coherent low and high biases of different magnitudes. Biases in different MERRA-2 FWI components were related to different biases in weather input variables for different regions, but temperature and relative humidity biases were the most important overall. FWI forecasts had high skill for 1–2 d lead times for most of the world. For longer lead times, forecast skill decreased most quickly at high latitudes and was most closely related to decreasing skill of relative humidity forecasts. These results provide a baseline for the evaluation and use of fire weather products calculated from global analysis and forecast fields.

Global Fire Weather Database↗
NASA NTRS2021

Developing a Hydrological Monitoring and Sub-Seasonal to Seasonal Forecasting System for South and Southeast Asian River Basins

South and Southeast Asia is subject to significant hydrometeorological extremes, including drought. Under rising temperatures, growing populations, and an apparent weakening of the South Asian monsoon in recent decades, concerns regarding drought and its potential impacts on water and food security are on the rise. Reliable sub-seasonal to seasonal (S2S) hydrological forecasts could, in principle, help governments and international organizations to better assess risk and act in the face of an oncoming drought. Here, we leverage recent improvements in S2S meteorological forecasts and the growing power of Earth Observations to provide more accurate monitoring of hydrological states for forecast initialization. Information from both sources is merged in a South and Southeast Asia sub-seasonal to seasonal hydrological forecasting system (SAHFS-S2S), developed collaboratively with the NASA SERVIR program and end-users across the region. This system applies the Noah-MultiParameterization (NoahMP) Land Surface Model (LSM) in the NASA Land Information System (LIS), driven by downscaled meteorological fields from the Global Data Assimilation System (GDAS) and Climate Hazards InfraRed Precipitation products (CHIRP and CHIRPS) to optimize initial conditions. The NASA Goddard Earth Observing System Model - sub-seasonal to seasonal (GEOS-S2S) forecasts, downscaled using the National Center for Atmospheric Research (NCAR) General Analog Regression Downscaling (GARD) tool and quantile mapping, are then applied to drive 5-km resolution hydrological forecasts to a 9-month forecast time horizon. Results show that the skillful predictions of root zone soil moisture can be made one to two months in advance for forecasts initialized in rainy seasons and up to 8 months when initialized in dry seasons. The memory of accurate initial conditions can positively contribute to forecast skills throughout the entire 9-month prediction period in areas with limited precipitation. This SAHFS-S2S has been operationalized at the International Centre for Integrated Mountain Development (ICIMOD) to support drought monitoring and warning needs in the region.

Yifan Zhou↗