NASA-Unified Weather Research And Forecasting Research Articles

High-resolution mapping of nitrogen oxide emissions in large US cities from TROPOMI retrievals of tropospheric nitrogen dioxide columns

Abstract. Satellite-derived spatiotemporal patterns of nitrogen oxide (NOx) emissions can improve accuracy of emission inventories to better support air quality and climate research and policy studies. In this study, we develop a new method by coupling the chemical transport Model-Independent SATellite-derived Emission estimation Algorithm for Mixed-sources (MISATEAM) with a divergence method to map high-resolution NOx emissions across US cities using TROPOspheric Monitoring Instrument (TROPOMI) tropospheric nitrogen dioxide (NO2) retrievals. The accuracy of the coupled method is validated through application to synthetic NO2 observations from the NASA-Unified Weather Research and Forecasting (NU-WRF) model, with a horizontal spatial resolution of 4 km × 4 km for 33 large and mid-size US cities. Validation reveals excellent agreement between inferred and NU-WRF-provided emission magnitudes (R= 0.99, normalized mean bias, NMB = −0.01) and a consistent spatial pattern when comparing emissions for individual grid cells (R=0.88±0.06). We then develop a TROPOMI-based database reporting annual emissions for 39 US cities at a horizontal spatial resolution of 0.05° × 0.05° from 2018 to 2021. This database demonstrates a strong correlation (R= 0.90) with the National Emission Inventory (NEI) but reveals some bias (NMB = −0.24). There are noticeable differences in the spatial patterns of emissions in some cities. Our analysis suggests that uncertainties in TROPOMI-based emissions and potential misallocation of emissions and/or missing sources in bottom-up emission inventories both contribute to these differences.

A new method for inferring city emissions and lifetimes of nitrogen oxides from high-resolution nitrogen dioxide observations: a model study

Abstract. We present a new method to infer emissions and lifetimes of nitrogen oxides (NOx) based on tropospheric nitrogen dioxide (NO2) observations together with reanalysis wind fields for cities located in polluted backgrounds. Since the accuracy of the method is difficult to assess due to lack of “true values” that can be used as a benchmark, we apply the method to synthetic NO2 observations derived from the NASA-Unified Weather Research and Forecasting (NU-WRF) model at a high horizontal spatial resolution of 4 km × 4 km for cities over the continental United States. We compare the inferred emissions and lifetimes with the values given by the NU-WRF model to evaluate the method. The method is applicable to 26 US cities. The derived results are generally in good agreement with the values given by the model, with the relative differences of 2 % ± 17 % (mean ± standard deviation) and 15 % ± 25 % for lifetimes and emissions, respectively. Our investigation suggests that the use of wind data prior to the satellite overpass time improves the performance of the method. The correlation coefficients between inferred and NU-WRF lifetimes increase from 0.56 to 0.79 and for emissions increase from 0.88 to 0.96 when comparing results based on wind fields sampled simultaneously with satellite observations and averaged over 9 h data prior to satellite observations, respectively. We estimate that uncertainties in NOx lifetime and emissions arising from the method are approximately 15 % and 20 %, respectively, for typical (US) cities. The total uncertainties reach up to 43 % (lifetimes) and 45 % (emissions) by considering the additional uncertainties associated with satellite NO2 observations and wind data. We expect this new method to be applicable to NO2 observations from the TROPOspheric Monitoring Instrument (TROPOMI) and geostationary satellites, such as Geostationary Environment Monitoring Spectrometer (GEMS) or the Tropospheric Emissions: Monitoring Pollution (TEMPO) instrument, to estimate urban NOx emissions and lifetimes globally.

Understanding the Impacts of Land Surface and PBL Observations on the Terrestrial and Atmospheric Legs of Land–Atmosphere Coupling

AbstractAccurately representing land–atmosphere (LA) interactions and coupling in NWP systems remains a challenge. New observations, incorporated into models via assimilation or calibration, hold the promise of improved forecast skill, but erroneous model coupling can hinder the benefits of such activities. To better understand model representation of coupled interactions and feedbacks, this study demonstrates a novel framework for coupled calibration of the single column model (SCM) capability of the NASA Unified Weather Research and Forecasting (NU-WRF) system coupled to NASA’s Land Information System (LIS). The local land–atmosphere coupling (LoCo) process chain paradigm is used to assess the processes and connections revealed by calibration experiments. Two summer case studies in the U.S. Southern Great Plains are simulated in which LSM parameters are calibrated to diurnal observations of LoCo process chain components including 2-m temperature, 2-m humidity, surface fluxes (Bowen ratio), and PBL height. Results show a wide range of soil moisture and hydraulic parameter solutions depending on which LA variable (i.e., observation) is used for calibration, highlighting that improvement in either soil hydraulic parameters or initial soil moisture when not in tandem with the other can provide undesirable results. Overall, this work demonstrates that a process chain calibration approach can be used to assess LA connections, feedbacks, strengths, and deficiencies in coupled models, as well as quantify the potential impact of new sources of observations of land–PBL variables on coupled prediction.

Object-Based Comparison of Data-Driven and Physics-Driven Satellite Estimates of Extreme Rainfall.

The Global Precipitation Measurement (GPM) constellation of spaceborne sensors provides a variety of direct and indirect measurements of precipitation processes. Such observations can be employed to derive spatially and temporally consistent gridded precipitation estimates either via data-driven retrieval algorithms or by assimilation into physically based numerical weather models. We compare the data-driven Integrated Multisatellite Retrievals for GPM (IMERG) and the assimilation-enabled NASA-Unified Weather Research and Forecasting (NU-WRF) model against Stage IV reference precipitation for four major extreme rainfall events in the southeastern United States using an object-based analysis framework that decomposes gridded precipitation fields into storm objects. As an alternative to conventional "grid-by-grid analysis," the object-based approach provides a promising way to diagnose spatial properties of storms, trace them through space and time, and connect their accuracy to storm types and input data sources. The evolution of two tropical cyclones are generally captured by IMERG and NU-WRF, while the less organized spatial patterns of two mesoscale convective systems pose challenges for both. NU-WRF rain rates are generally more accurate, while IMERG better captures storm location and shape. Both show higher skill in detecting large, intense storms compared to smaller, weaker storms. IMERG's accuracy depends on the input microwave and infrared data sources; NU-WRF does not appear to exhibit this dependence. Findings highlight that an object-oriented view can provide deeper insights into satellite precipitation performance and that the satellite precipitation community should further explore the potential for "hybrid" data-driven and physics-driven estimates in order to make optimal usage of satellite observations.

Assimilation of Satellite-Derived Soil Moisture for Improved Forecasts of the Great Plains Low-Level Jet

AbstractIn the context of forecasting societally impactful Great Plains low-level jets (GPLLJs), the potential added value of satellite soil moisture (SM) data assimilation (DA) is high. GPLLJs are both sensitive to regional soil moisture gradients and frequent drivers of severe weather, including mesoscale convective systems. An untested hypothesis is that SM DA is more effective in forecasts of weakly synoptically forced, or uncoupled GPLLJs, than in forecasts of cyclone-induced coupled GPLLJs. Using the NASA Unified Weather Research and Forecasting (NU-WRF) Model, 75 GPLLJs are simulated at 9-km resolution both with and without NASA Soil Moisture Active Passive SM DA. Differences in modeled SM, surface sensible (SH) and latent heat (LH) fluxes, 2-m temperature (T2), 2-m humidity (Q2), PBL height (PBLH), and 850-hPa wind speed (W850) are quantified for individual jets and jet-type event subsets over the south-central Great Plains, as well as separately for each GPLLJ sector (entrance, core, and exit). At the GPLLJ core, DA-related changes of up to 5.4 kg m−2in SM can result in T2, Q2, LH, SH, PBLH, and W850 differences of 0.68°C, 0.71 g kg−2, 59.9 W m−2, 52.4 W m−2, 240 m, and 4 m s−1, respectively. W850 differences focus along the jet axis and tend to increase from south to north. Jet-type differences are most evident at the GPLLJ exit where DA increases and decreases W850 in uncoupled and coupled GPLLJs, respectively. Data assimilation marginally reduces negative wind speed bias for all jets, but the correction is greater for uncoupled GPLLJs, as hypothesized.

Evaluation of NU-WRF model performance on air quality simulation under various model resolutions – an investigation within the framework of MICS-Asia Phase III

Abstract. Horizontal grid resolution has a profound effect on model performances on meteorology and air quality simulations. In contribution to MICS-Asia Phase III, one of whose goals was to identify and reduce model uncertainty in air quality prediction, this study examined the impact of grid resolution on meteorology and air quality simulation over East Asia, focusing on the North China Plain (NCP) region. The NASA Unified Weather Research and Forecasting (NU-WRF) model has been applied with the horizontal resolutions at 45, 15, and 5 km. The results revealed that, in comparison with ground observations, no single resolution can yield the best model performance for all variables across all stations. From a regional average perspective (i.e., across all monitoring sites), air temperature modeling was not sensitive to the grid resolution but wind and precipitation simulation showed the opposite. NU-WRF with the 5 km grid simulated the wind speed best, while the 45 km grid yielded the most realistic precipitation as compared to the site observations. For air quality simulations, finer resolution generally led to better comparisons with observations for O3, CO, NOx, and PM2.5. However, the improvement of model performance on air quality was not linear with the resolution increase. The accuracy of modeled surface O3 of the 15 km grid was greatly improved over the one from the 45 km grid. A further increase in grid resolution to 5 km, however, showed diminished impact on model performance improvement on O3 prediction. In addition, a 5 km resolution grid showed large advantage in better capturing the frequency of high-pollution occurrences. This was important for the assessment of noncompliance with ambient air quality standards, which was key to air quality planning and management. Balancing the modeling accuracy and resource limitation, a 15 km grid resolution was suggested for future MICS-Asia air quality modeling activity if the research region remained unchanged. This investigation also found a large overestimate of ground-level O3 and an underestimate of surface NOx and CO, likely due to missing emissions of NOx and CO.

Evaluation of the CABLEv2.3.4 Land Surface Model Coupled to NU‐WRFv3.9.1.1 in Simulating Temperature and Precipitation Means and Extremes Over CORDEX AustralAsia Within a WRF Physics Ensemble

AbstractThe Community Atmosphere Biosphere Land Exchange (CABLE) model is a third‐generation land surface model (LSM). CABLE is commonly used as a stand‐alone LSM, coupled to the Australian Community Climate and Earth Systems Simulator global climate model and coupled to the Weather Research and Forecasting (WRF) model for regional applications. Here, we evaluate an updated version of CABLE within a WRF physics ensemble over the COordinated Regional Downscaling EXperiment (CORDEX) AustralAsia domain. The ensemble consists of different cumulus, radiation and planetary boundary layer (PBL) schemes. Simulations are carried out within the NASA Unified WRF modeling framework, NU‐WRF. Our analysis did not identify one configuration that consistently performed the best for all diagnostics and regions. Of the cumulus parameterizations the Grell‐Freitas cumulus scheme consistently overpredicted precipitation, while the new Tiedtke scheme was the best in simulating the timing of precipitation events. For the radiation schemes, the RRTMG radiation scheme had a general warm bias. For the PBL schemes, the YSU scheme had a warm bias, and the MYJ PBL scheme a cool bias. Results are strongly dependent on the region of interest, with the northern tropics and southwest Western Australia being more sensitive to the choice of physics options compared to southeastern Australia which showed less overall variation and overall better performance across the ensemble. Comparisons with simulations using the Unified Noah LSM showed that CABLE in NU‐WRF has a more realistic simulation of evapotranspiration when compared to GLEAM estimates.

Expanding the Goddard CSH Algorithm for GPM: New Extratropical Retrievals.

The Goddard Convective-Stratiform Heating (CSH) algorithm has been used to retrieve latent heating (LH) associated with clouds and cloud systems in support of the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) mission. The CSH algorithm required the use of a cloud-resolving model (CRM) to simulate LH profiles to build look-up tables (LUTs). However, the current LUTs in the CSH algorithm are not suitable for retrieving LH profiles at high latitudes or winter conditions that are needed for GPM. The NASA Unified-Weather Research and Forecasting (NU-WRF) model is used to simulate three eastern continental US (CONUS) synoptic winter and three western coastal/offshore events. The relationship between LH structures (or profiles) and other precipitation properties (radar reflectivity, freezing level height, echo-top height, maximum radar reflectivity height and surface precipitation rate) is examined, and a new classification system is adopted with varying ranges for each of these precipitation properties to create LUTs representing high latitude/winter conditions. The performance of the new LUTs is examined using a self-consistency check for one CONUS and one West Coast offshore event by comparing LH profiles retrieved from the LUTs using model-simulated precipitation properties with those originally simulated by the model. The results of the self-consistency check validate the new classification and LUTs. High latitude retrievals from the new LUTs are merged with those from the CSH algorithm to retrieve LH profiles over the GPM domain using precipitation properties retrieved from the GPM combined algorithm.

NU-WRF Aerosol Transport Simulation over West Africa: Effects of Biomass Burning on Smoke Aerosol Distribution

AbstractSeries of aerosol transport hindcasts for West Africa were conducted using the Weather Research and Forecasting (WRF) Model coupled to chemistry within the NASA-Unified WRF (NU-WRF) framework. The transport of biomass-burning aerosols in April and December 2009 was investigated over two types of simulation domains. One-month simulations with 9-km grid spacing for April or December 2009 covered most of North and West Africa and were evaluated by comparison with measurements of the total-column aerosol optical depth, Ångström exponent, and horizontal wind components at various pressure levels. The horizontal wind components at 700 hPa were identified as key factors in determining the transport patterns of biomass-burning aerosols from sub-Saharan West Africa to the Sahel. The vertical accumulation of biomass-burning aerosols close to 700 hPa was demonstrated in 1-day simulations with 1-km horizontal grid spacing. A new simple parameterization for the effects of heat release by biomass burning was designed for this resolution and tested together with the conventional parameterization based on fixed smoke injection heights. The aerosol vertical profiles were somewhat sensitive to the selection of parameterization, except for cases with the assumption of excessive heating by biomass burning. The new parameterization works reasonably well and offers flexibility to relate smoke transport to biomass-burning plume rise that can be correlated with the satellite fire radiative power measurements, which is advantageous relative to the conventional parameterization.

Microphysics and Radiation Effect of Dust on Saharan Air Layer: An HS3 Case Study

A Saharan air layer (SAL) event associated with a nondeveloping African easterly wave (AEW) over the main development region of the eastern Atlantic was sampled by the NASA Global Hawk aircraft on 24–25 August 2013 during the NASA Hurricane and Severe Storm Sentinel (HS3) campaign and was simulated with the NASA Unified Weather Research and Forecasting (NU-WRF) Model. Airborne, ground-based, and spaceborne measurements were used to evaluate the model performance. The microphysical and radiative effects of dust and other aerosols on the SAL structure and environment were investigated with the factor-separation method. The results indicate that relative to a simulation without dust–radiative and microphysical impacts, Saharan dust and other aerosols heated the SAL air mainly through shortwave heating by the direct aerosol–radiation (AR) effect, resulting in a warmer (up to 0.6 K) and drier (up to 5% RH reduction) SAL and maintaining the strong temperature inversion at the base of the SAL in the presence of predominant longwave cooling. Radiative heating of the dust accentuated a vertical circulation within the dust layer, in which air rose (sank) in the northern (southern) portions of the dust layer. Furthermore, above and to the south of the dust layer, both the microphysical and radiative impacts of dust tended to counter the vertical motions associated with the Hadley circulation, causing a small weakening and southward shift of convection in the intertropical convergence zone (ITCZ) and reduced anvil cloud to the north. Changes in moisture and cloud/precipitation hydrometeors were largely driven by the dust-induced changes in vertical motion. Dust strengthened the African easterly jet by up to ~1 m s−1 at the southern edge of the jet, primarily through the AR effect, and produced modest increases in vertical wind shear within and in the vicinity of the dust layer. These modulations of the SAL and AEW environment clearly contributed to the nondevelopment of this AEW.

Evaluation of Cloud Microphysical Schemes for a Warm Frontal Snowband during the GPM Cold Season Precipitation Experiment (GCPEx)

Detailed observations from the Global Precipitation Measurement (GPM) mission Cold Season Precipitation Experiment (GCPEx) of an intense warm frontal band on 18 February 2012 were used to evaluate several bulk microphysical parameterizations within the NASA-Unified Weather Research and Forecasting (NU-WRF) Model. These included the Predicted Particle Properties (P3), Morrison (MORR), Stony Brook University (SBU), and Goddard four-class ice (4ICE) microphysics schemes. All schemes were able to predict the snowband, but the simulated intensities varied because of various assumptions in these schemes. The saturation adjustment scheme within MORR promoted excessive amounts of cloud water evaporational cooling in the warm sector, which contributed to a decrease in midlevel instability approaching the frontal band and thus a weaker band. In contrast, the explicit calculation of cloud water condensation/evaporation in the P3 scheme produced limited amounts of evaporational cooling, which allowed for greater midlevel instability to support band development. The P3 and SBU schemes produced moderate rime/graupel mass within the band that was confirmed by observations, while the MORR and 4ICE schemes drastically underpredicted the graupel mass. The high-density, fast-falling rimed particles in P3 underwent weak sublimation and melting, which helped promote a stronger horizontal temperature gradient and greater low-level instability along the frontal band compared to the other schemes. Overall, the schemes that use specified thresholds for converting between the predefined ice-phase categories of cloud ice, snow, and graupel had the most unrepresentative hydrometeor types. These results highlight the advantage of predicting ice particle properties and explicitly calculating cloud water condensation/evaporation in the P3 scheme.

Impact of Assimilated Precipitation-Sensitive Radiances on the NU-WRF Simulation of the West African Monsoon

AbstractThis work assimilates multisensor precipitation-sensitive microwave radiance observations into a storm-scale NASA Unified Weather Research and Forecasting (NU-WRF) Model simulation of the West African monsoon. The analysis consists of a full description of the atmospheric states and a realistic cloud and precipitation distribution that is consistent with the observed dynamic and physical features. The analysis shows an improved representation of monsoon precipitation and its interaction with dynamics over West Africa. Most significantly, assimilation of precipitation-affected microwave radiance has a positive impact on the distribution of precipitation intensity and also modulates the propagation of cloud precipitation systems associated with the African easterly jet. Using an ensemble-based assimilation technique that allows state-dependent forecast error covariance among dynamical and microphysical variables, this work shows that the assimilation of precipitation-sensitive microwave radiances over the West African monsoon rainband enables initialization of storms. These storms show the characteristics of continental tropical convection that enhance the connection between tropical waves and organized convection systems.

Biogenic isoprene emissions driven by regional weather predictions using different initialization methods: case studies during the SEAC<sup>4</sup>RS and DISCOVER-AQ airborne campaigns

Abstract. Land and atmospheric initial conditions of the Weather Research and Forecasting (WRF) model are often interpolated from a different model output. We perform case studies during NASA's SEAC4RS and DISCOVER-AQ Houston airborne campaigns, demonstrating that using land initial conditions directly downscaled from a coarser resolution dataset led to significant positive biases in the coupled NASA-Unified WRF (NUWRF, version 7) surface and near-surface air temperature and planetary boundary layer height (PBLH) around the Missouri Ozarks and Houston, Texas, as well as poorly partitioned latent and sensible heat fluxes. Replacing land initial conditions with the output from a long-term offline Land Information System (LIS) simulation can effectively reduce the positive biases in NUWRF surface air temperature by ∼ 2 °C. We also show that the LIS land initialization can modify surface air temperature errors almost 10 times as effectively as applying a different atmospheric initialization method. The LIS-NUWRF-based isoprene emission calculations by the Model of Emissions of Gases and Aerosols from Nature (MEGAN, version 2.1) are at least 20 % lower than those computed using the coarser resolution data-initialized NUWRF run, and are closer to aircraft-observation-derived emissions. Higher resolution MEGAN calculations are prone to amplified discrepancies with aircraft-observation-derived emissions on small scales. This is possibly a result of some limitations of MEGAN's parameterization and uncertainty in its inputs on small scales, as well as the representation error and the neglect of horizontal transport in deriving emissions from aircraft data. This study emphasizes the importance of proper land initialization to the coupled atmospheric weather modeling and the follow-on emission modeling. We anticipate it to also be critical to accurately representing other processes included in air quality modeling and chemical data assimilation. Having more confidence in the weather inputs is also beneficial for determining and quantifying the other sources of uncertainties (e.g., parameterization, other input data) of the models that they drive.

Evaluating hourly rainfall characteristics over the U.S. Great Plains in dynamically downscaled climate model simulations using NASA‐Unified WRF

AbstractAccurate simulation of extreme precipitation events remains a challenge in climate models. This study utilizes hourly precipitation data from ground stations and satellite instruments to evaluate rainfall characteristics simulated by the NASA‐Unified Weather Research and Forecasting (NU‐WRF) regional climate model at horizontal resolutions of 4, 12, and 24 km over the Great Plains of the United States. We also examined the sensitivity of the simulated precipitation to different spectral nudging approaches and the cumulus parameterizations. The rainfall characteristics in the observations and simulations were defined as an hourly diurnal cycle of precipitation and a joint probability distribution function (JPDF) between duration and peak intensity of precipitation events over the Great Plains in summer. We calculated a JPDF for each data set and the overlapping area between observed and simulated JPDFs to measure the similarity between the two JPDFs. Comparison of the diurnal precipitation cycles between observations and simulations does not reveal the added value of high‐resolution simulations. However, the performance of NU‐WRF simulations measured by the JPDF metric strongly depends on horizontal resolution. The simulation with the highest resolution of 4 km shows the best agreement with the observations in simulating duration and intensity of wet spells. Spectral nudging does not affect the JPDF significantly. The effect of cumulus parameterizations on the JPDFs is considerable but smaller than that of horizontal resolution. The simulations with lower resolutions of 12 and 24 km show reasonable agreement but only with the high‐resolution observational data that are aggregated into coarse resolution and spatially averaged.

Sensitivity of CONUS Summer Rainfall to the Selection of Cumulus Parameterization Schemes in NU-WRF Seasonal Simulations

Abstract This study investigates the sensitivity of daily rainfall rates in regional seasonal simulations over the contiguous United States (CONUS) to different cumulus parameterization schemes. Daily rainfall fields were simulated at 24-km resolution using the NASA-Unified Weather Research and Forecasting (NU-WRF) Model for June–August 2000. Four cumulus parameterization schemes and two options for shallow cumulus components in a specific scheme were tested. The spread in the domain-mean rainfall rates across the parameterization schemes was generally consistent between the entire CONUS and most subregions. The selection of the shallow cumulus component in a specific scheme had more impact than that of the four cumulus parameterization schemes. Regional variability in the performance of each scheme was assessed by calculating optimally weighted ensembles that minimize full root-mean-square errors against reference datasets. The spatial pattern of the seasonally averaged rainfall was insensitive to the selection of cumulus parameterization over mountainous regions because of the topographical pattern constraint, so that the simulation errors were mostly attributed to the overall bias there. In contrast, the spatial patterns over the Great Plains regions as well as the temporal variation over most parts of the CONUS were relatively sensitive to cumulus parameterization selection. Overall, adopting a single simulation result was preferable to generating a better ensemble for the seasonally averaged daily rainfall simulation, as long as their overall biases had the same positive or negative sign. However, an ensemble of multiple simulation results was more effective in reducing errors in the case of also considering temporal variation.

Development of high-resolution dynamic dust source function - A case study with a strong dust storm in a regional model

A high-resolution dynamic dust source has been developed in the NASA Unified-Weather Research and Forecasting (NU-WRF) model to improve the existing coarse static dust source. In the new dust source map, topographic depression is in 1-km resolution and surface bareness is derived using the Normalized Difference Vegetation Index (NDVI) data from Moderate Resolution Imaging Spectroradiometer (MODIS). The new dust source better resolves the complex topographic distribution over the Western United States where its magnitude is higher than the existing, coarser resolution static source. A case study is conducted with an extreme dust storm that occurred in Phoenix, Arizona in 02–03 UTC July 6, 2011. The NU-WRF model with the new high-resolution dynamic dust source is able to successfully capture the dust storm, which was not achieved with the old source identification. However the case study also reveals several challenges in reproducing the time evolution of the short-lived, extreme dust storm events.

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.

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.

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.

Introducing multisensor satellite radiance-based evaluation for regional Earth System modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.