Development of a composite drought indicator for operational drought monitoring in the MENA region

This paper presents the composite drought indicator (CDI) that Jordanian, Lebanese, Moroccan, and Tunisian government agencies now produce monthly to support operational drought management decision making, and it describes their iterative co-development processes. The CDI is primarily intended to monitor agricultural and ecological drought on a seasonal time scale. It uses remote sensing and modelled data inputs, and it reflects anomalies in precipitation, vegetation, soil moisture, and evapotranspiration. Following quantitative and qualitative validation assessments, engagements with policymakers, and consideration of agencies' technical and institutional capabilities and constraints, we made changes to CDI input data, modelling procedures, and integration to tailor the system for each national context. We summarize validation results, drought modelling challenges and how we overcame them through CDI improvements, and we describe the monthly CDI production process and outputs. Finally, we synthesize procedural and technical aspects of CDI development and reflect on the constraints we faced as well as trade-offs made to optimize the CDI for operational monitoring to support policy decision-making-including aspects of salience, credibility, and legitimacy-within each national context.

Applications of Dynamic Land Surface Information for Passive Microwave Precipitation Retrieval

Accurate, physically based precipitation retrieval over global land surfaces is an important goal of the NASA/JAXA Global Precipitation Measurement Mission (GPM). This is a difficult problem for the passive microwave constellation, as the signal over radiometrically warm land surfaces in the microwave frequencies means that the measurements used are indirect and typically require inferring some type of relationship between an observed scattering signal and precipitation at the surface. GPM, with collocated radiometer and dual-frequency radar, is an excellent tool for tackling this problem and improving global retrievals. In the years following the launch of the GPM Core Observatory satellite, physically based passive microwave retrieval of precipitation over land continues to be challenging. Validation efforts suggest that the operational GPM passive microwave algorithm, the Goddard profiling algorithm (GPROF), tends to overestimate precipitation at the low (<5 mm h^-1) end of the distribution over land. In this work, retrieval sensitivities to dynamic surface conditions are explored through enhancement of the algorithm with dynamic, retrieved information from a GPM-derived optimal estimation scheme. The retrieved parameters describing surface and background characteristics replace current static or ancillary GPROF information including emissivity, water vapor, and snow cover. Results show that adding this information decreases probability of false detection by 50% and, most importantly, the enhancements with retrieved parameters move the retrieval away from dependence on ancillary datasets and lead to improved physical consistency.

Earth observations and integrative models in support of food and water security

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries.

The instantaneous retrieval of precipitation over land by temporal variation at 19 GHz

The primary signal used in all current passive microwave precipitation retrieval algorithms over land is the depression of the instantaneous brightness temperature (TB) caused by ice scattering. This study presents a new methodology to retrieve instantaneous precipitation rate over land by using TB temporal variation (ΔTB) at 19 GHz, which primarily reflects the surface emissivity variation due to the precipitation impact. As a proof-of-concept, we exploit observations from five polar-orbiting satellites over the Southern Great Plains (SGP) of the United States. Results show that ΔTB at 19 GHz correlate well with the instantaneous precipitation rate. Further analysis shows that ΔTB at 19 GHz is better correlated with the precipitation rate when multiple satellite observations are used due to the much shorter re-visit time for a certain location. The retrieved instantaneous precipitation rate over SGP from ΔTB at 19 GHz reasonably agrees with the surface radar observations, with the correlation, the root mean square error and the bias being 0.49, 2.39 mm/hr and 6.54%, respectively. Future work seeks to combine the ice scattering signal at high frequencies and this surface emissivity variation signal at low frequencies to achieve an optimal retrieval performance.

Impacts of aerosol-monsoon interaction on rainfall and circulation over Northern India and the Himalaya Foothills

The boreal summer of 2008 was unusual for the Indian monsoon, featuring exceptional heavy loading of dust aerosols over the Arabian Sea and northern-central India, near normal all-India rainfall, but excessive heavy rain, causing disastrous flooding in the Northern Indian Himalaya Foothills (NIHF) regions, accompanied by persistent drought conditions in central and southern India. Using NASA Unified-physics Weather Research Forecast (NUWRF) model with fully interactive aerosol physics and dynamics, we carried out three sets of 7-day ensemble model forecast experiments: 1) control with no aerosol, 2) aerosol radiative effect only and 3) aerosol radiative and aerosol-cloud-microphysics effects, to study the impacts of aerosol-monsoon interactions on monsoon variability over the NIHF during the summer of 2008. Results show that aerosol-radiation interaction (ARI), i.e., dust aerosol transport, and dynamical feedback processes induced by aerosol-radiative heating, plays a key role in altering the large-scale monsoon circulation system, reflected by an increased north-south tropospheric temperature gradient, a northward shift of heavy monsoon rainfall, advancing the monsoon onset by 1-5 days over the HF, consistent with the EHP hypothesis (Lau et al. 2006). Additionally, we found that dust aerosols, via the semi-direct effect, increase atmospheric stability, and cause the dissipation of a developing monsoon onset cyclone over northeastern India/northern Bay of Bengal. Eventually, in a matter of several days, ARI transforms the developing monsoon cyclone into meso-scale convective cells along the HF slopes. Aerosol-Cloud-microphysics Interaction (ACI) further enhances the ARI effect in invigorating the deep convection cells and speeding up the transformation processes. Results indicate that even in short-term (up to weekly) numerical forecasting of monsoon circulation and rainfall, effects of aerosol-monsoon interaction can be substantial and cannot be ignored.

Tradeoff between cost and accuracy in large-scale surface water dynamic modeling

Recent efforts have led to the development of the local inertia formulation (INER) for an accurate but still cost-efficient representation of surface water dynamics, compared to the widely used kinematic wave equation (KINE). In this study, both formulations are evaluated over the Amazon basin in terms of computational costs and accuracy in simulating streamflows and water levels through synthetic experiments and comparisons against ground-based observations. Varying time steps are considered as part of the evaluation and INER at 60-second time step is adopted as the reference for synthetic experiments. Five hybrid (HYBR) realizations are performed based on maps representing the spatial distribution of the two formulations that physically represent river reach flow dynamics within the domain. Maps have fractions of KINE varying from 35.6% to 82.8%. KINE runs show clear deterioration along the Amazon river and main tributaries, with maximum RMSE values for streamflow and water level reaching 7827m3.s-1 and 1379cm near the basin's outlet. However, KINE is at least 25% more efficient than INER with low model sensitivity to longer time steps. A significant improvement is achieved with HYBR, resulting in maximum RMSE values of 3.9-292m3.s-1 for streamflows and 1.1-28.5cm for water levels, and cost reduction of 6-16%, depending on the map used. Optimal results using HYBR are obtained when the local inertia formulation is used in about one third of the Amazon basin, reducing computational costs in simulations while preserving accuracy. However, that threshold may vary when applied to different regions, according to their hydrodynamics and geomorphological characteristics.

Operational Hydrological Forecasting during the IPHEx-IOP Campaign - Meet the Challenge

An operational streamflow forecasting testbed was implemented during the Intense Observing Period (IOP) of the Integrated Precipitation and Hydrology Experiment (IPHEx-IOP) in May-June 2014 to characterize flood predictability in complex terrain. Specifically, hydrological forecasts were issued daily for 12 headwater catchments in the Southern Appalachians using the Duke Coupled surface-groundwater Hydrology Model (DCHM) forced by hourly atmospheric fields and QPFs (Quantitative Precipitation Forecasts) produced by the NASA-Unified Weather Research and Forecasting (NU-WRF) model. Previous day hindcasts forced by radar-based QPEs (Quantitative Precipitation Estimates) were used to provide initial conditions for present day forecasts. This manuscript first describes the operational testbed framework and workflow during the IPHEx-IOP including a synthesis of results. Second, various data assimilation approaches are explored a posteriori (post-IOP) to improve operational (flash) flood forecasting. Although all flood events during the IOP were predicted by the IPHEx operational testbed with lead times of up to 6 hours, significant errors of over- and, or under-prediction were identified that could be traced back to the QPFs and subgrid-scale variability of radar QPEs. To improve operational flood prediction, three data-merging strategies were pursued post-IOP: 1) the spatial patterns of QPFs were improved through assimilation of satellite-based microwave radiances into NU-WRF; 2) QPEs were improved by merging raingauge observations with ground-based radar observations using bias-correction methods to produce streamflow hindcasts and associated uncertainty envelope capturing the streamflow observations, and 3) river discharge observations were assimilated into the DCHM to improve streamflow forecasts using the Ensemble Kalman Filter (EnKF), the fixed-lag Ensemble Kalman Smoother (EnKS), and the Asynchronous EnKF (i.e. AEnKF) methods. Both flood hindcasts and forecasts were significantly improved by assimilating discharge observations into the DCHM. Specifically, Nash-Sutcliff Efficiency (NSE) values as high as 0.98, 0.71 and 0.99 at 15-min time-scales were attained for three headwater catchments in the inner mountain region demonstrating that the assimilation of discharge observations at the basin's outlet can reduce the errors and uncertainties in soil moisture at very small scales. Success in operational flood forecasting at lead times of 6, 9, 12 and 15hrs was also achieved through discharge assimilation with NSEs of 0.87, 0.78, 0.72 and 0.51, respectively. Analysis of experiments using various data assimilation system configurations indicates that the optimal assimilation time window depends both on basin properties and storm-specific space-time-structure of rainfall, and therefore adaptive, context-aware, configurations of the data assimilation system are recommended to address the challenges of flood prediction in headwater basins.

Evaluation of NU-WRF rainfall forecasts for IFloodS

The Iowa Flood Studies (IFloodS) campaign was conducted in eastern Iowa as a pre-GPM-launch campaign from 1 May to 15 June 2013. During the campaign period, real time forecasts are conducted utilizing NASA-Unified Weather Research and Forecasting (NU-WRF) model to support the everyday weather briefing. In this study, two sets of the NU-WRF rainfall forecasts are evaluated with Stage IV and Multi-Radar Multi-Sensor (MRMS) Quantitative Precipitation Estimation (QPE), with the objective to understand the impact of Land Surface initialization on the predicted precipitation. NU-WRF is also compared with North American Mesoscale Forecast System (NAM) 12 km forecast. In general, NU-WRF did a good job at capturing individual precipitation events. NU-WRF is also able to replicate a better rainfall spatial distribution compare with NAM. Further sensitivity tests show that the high-resolution makes a positive impact on rainfall forecast. The two sets of NU-WRF simulations produce very close rainfall characteristics. The Land surface initialization do not show significant impact on short-term rainfall forecast, and it is largely due to the soil conditions during the field campaign period.

High-resolution NU-WRF simulations of a deep convective-precipitation system during MC3E: Further improvements and comparisons between Goddard microphysics schemes and observations

The Goddard microphysics was recently improved by adding a fourth ice class (frozen drops/hail). This new 4ICE scheme was developed and tested in the Goddard Cumulus Ensemble (GCE) model for an intense continental squall line and a moderate, less organized continental case. Simulated peak radar reflectivity profiles were improved in intensity and shape for both cases, as were the overall reflectivity probability distributions versus observations. In this study, the new Goddard 4ICE scheme is implemented into the regional-scale NASA Unified-Weather Research and Forecasting (NU-WRF) model, modified and evaluated for the same intense squall line, which occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E). NU-WRF simulated radar reflectivities, total rainfall, propagation, and convective system structures using the 4ICE scheme modified herein agree as well as or significantly better with observations than the original 4ICE and two previous 3ICE (graupel or hail) versions of the Goddard microphysics. With the modified 4ICE, the bin microphysics-based rain evaporation correction improves propagation and in conjunction with eliminating the unrealistic dry collection of ice/snow by hail can replicate the erect, narrow, and intense convective cores. Revisions to the ice supersaturation, ice number concentration formula, and snow size mapping, including a new snow breakup effect, allow the modified 4ICE to produce a stronger, better organized system, more snow, and mimic the strong aggregation signature in the radar distributions. NU-WRF original 4ICE simulated radar reflectivity distributions are consistent with and generally superior to those using the GCE due to the less restrictive domain and lateral boundaries.