Search NASASearch

Engineering topics

Andrea M Molod

Publications and source records attributed to Andrea M Molod.

At least 19 records

The Simulation and Subseasonal Forecasting of Hydrological Variables: Insights from a Simple Water Balance Model

Past work has shown that a land surface model’s (LSM’s) implicit (not explicitly coded) relationships between soil moisture and both evapotranspiration (ET) and runoff largely determine the LSM’s hydrological behavior. Here we estimate the relationships that appear to be operating in the real world and compare them to those of the LSM component of a state-of-the-art Earth system model (ESM). The two sets of relationships are determined by calibrating them within a simple water balance model (WBM): once using stream gauge observations from small, unregulated rivers over the eastern half of the U.S., and once using the runoffs generated by the LSM as part of a state-of-the-art atmospheric reanalysis. Hydrological simulations and subseasonal hydrological forecasts performed with the two calibrated versions of the WBM provide two key results. First, the version calibrated to the LSM-generated runoffs does successfully reproduce, to first order, the hydrological behavior of the full LSM within its ESM environment. Second, of the two WBM versions, the one calibrated to the observations reproduces more accurately a broad collection of fully independent streamflow observations as well as a similarly broad collection of in-situ soil moisture measurements. Taken together, the two results suggest that the observations-calibrated ET and runoff efficiency functions do successfully represent, at least to some degree, soil moisture controls over hydrological variability in Nature and can serve as potentially useful targets for further LSM development.

Water Balance Model

Asymmetry in Subseasonal Surface Air Temperature Forecast Error with Respect to Soil Moisture Initialization

Soil moisture (W) helps control evapotranspiration (ET), and ET variations can in turn have a distinct impact on 2-m air temperature (T2M), given that increases in evaporative cooling encourage reduced temperatures. Soil moisture is accordingly linked to T2M, and realistic soil moisture initialization has, in previous studies, been shown to improve the skill of subseasonal T2M forecasts. The relationship between soil moisture and evapotranspiration, however, is distinctly nonlinear, with ET tending to increase with soil moisture in drier conditions and to be insensitive to soil moisture variations in wetter conditions. Here, through an extensive analysis of subseasonal forecasts produced with a state-of-the-art seasonal forecast system, this nonlinearity is shown to imprint itself on T2M forecast error in the conterminous United States in two unique ways: (i) the T2M forecast bias (relative to independent observations) induced by a negative precipitation bias tends to be larger for dry initializations, and (ii) on average, the unbiased root-mean-square error (ubRMSE) tends to be larger for dry initializations. Such findings can aid in the identification of forecasts of opportunity; taken a step further, they suggest a pathway for improving bias correction and uncertainty estimation in subseasonal T2M forecasts by conditioning each on initial soil moisture state.

Air Temperature

Improving the Representation of Subtropical Boundary Layer Clouds in the NASA GEOS Model with the Eddy-Diffusivity/Mass-Flux Parameterization

A systematic underestimation of subtropical planetary boundary layer (PBL) stratocu-mulus clouds by the GEOS model has been significantly improved by a new Eddy-Diffusivity/Mass-Flux (EDMF) parameterization. The EDMF parameterization representsthe subgrid-scale transport in the dry and moist parts of the PBL in a unified mannerand it combines an adjusted eddy-diffusivity PBL scheme from GEOS with a stochasticmulti-plume mass-flux model. The new EDMF version of the GEOS model is first com-pared against the CONTROL version in a single-column model (SCM) framework for twobenchmark cases representing subtropical stratocumulus and shallow cumulus clouds, andvalidated against large eddy simulations. Global simulations are performed and comparedagainst observations and reanalysis data. The results show that the EDMF version of theGEOS model produces more realistic subtropical PBL clouds. Noteworthy, the EDMF im-provements first detected in the SCM framework translate into similar improvements of theglobal GEOS model.

Kay Suselj

Subseasonal-to-Seasonal Hindcast Skill Assessment of Ridging Events Related to Drought Over the Western United States

Persistent atmospheric ridging events centered near the western United States are associated with widespread precipitation deficits and meteorological drought. Due to the relatively low skill of- dynamical models in forecasting precipitation on subseasonal‐to‐seasonal (S2S) time scales across this region, forecasts of ridging are explored in this study as a potential bridge for early warning drought prediction. To assess skill, we evaluate deterministic and probabilistic S2S hindcasts (out to 6 week lead time) of ridging events and geopotential height anomalies over the western United States in five ensemble hindcast systems. Prediction skill for ridging events is shown to be highly variable across models. For some models, longer‐time‐averaged patterns of geopotential height anomalies across the first 6 weeks can be skillfully simulated, when evaluated probabilistically against climatology. The most skillful models show modest skill in forecasting above normal ridging occurrences at lead times of Weeks 3–4 and Weeks 5–6,with some sensitivity to the specific ridge location and method for determining skill. Using the European Centre for Medium‐Range Weather Forecasts (ECMWF) model as a case study, longer lead time forecast busts are shown to often occur under extended periods of highly active La Nina‐like tropical convection. Model errors at simulating these tropical convection features may have a disproportionately large impact and degrade downstream forecasts over the western United States. Despite the documented model shortcomings, our results highlight an opportunity to target skillful ridging forecasts from dynamical models for improving early warning drought forecasting at S2S lead times.

S2S

Differences in Tropical High Clouds Among Reanalyses: Origins and Radiative Impacts

We examine differences among reanalysis highcloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifthgeneration European Centre for Medium-range Weather Forecasts – ECMWF – reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanalysis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Climate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanalyses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of longwave, short-wave, and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and short-wave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radiative heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERAInterim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and feed back to convective behaviour through their relationships with thermodynamic stability.

High Clouds

Using a Simple Water Balance Framework to Quantify the Impact of Soil Moisture Initialization on Subseasonal Evapotranspiration and Air Temperature Forecasts

Past studies have shown that accurate soil moisture initialization can contribute significant skill to near-surface air temperature (T2M) forecasts at subseasonal leads. The mechanisms by which soil moisture contributes such skill are examined here with a simple water balance-based model that captures the essence of soil moisture behavior in a state-of-the-art subseasonal-to-seasonal (S2S) forecasting system. The simple model successfully transforms initial soil moisture contents into average “forecasted” ET values at 16-30 day lead that agree well, during summer, with the values forecasted by the full NASA GEOS S2S system, indicating that soil moisture initialization dominates over forecasted meteorology in determining ET fluxes at subseasonal leads. When the simple model’s ET anomalies are interpreted in terms of T2M anomalies, a similar conclusion is reached for T2M: soil moisture initialization explains much (about 50% in the eastern half of the continental US) of the T2M anomalies produced by the full GEOS S2S system at 16-30 day lead, and the T2M forecasts produced by the simple model capture about half of the skill attained by the full system. The simple model’s framework is particularly conducive to an analysis of uncertainty in forecasts. Drier soils are generally found to induce larger uncertainty in ET (and thus T2M) forecasts, a result linked to the functional form relating ET to soil moisture in the simple model and verified by an analysis of the ensemble spreads within the forecasts produced by the full GEOS S2S system

Randal D Koster

Representation of Tropical Storms by the Modern-Era Retrospective Analysis for Research and Applications Version 2

This study examines the veracity of the tropical cyclone (TC) statistics estimated from the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) global atmospheric reanalysis, focusing on the climatological-mean genesis regions, tracks and their lifetime maximum intensity, as well as the interannual and intraseasonal variations in TC activity. The results are validated against the International Best Track Archive for Climate Stewardship (IBTrACS) data for 1980-2016. MERRA-2 represents the spatial distribution of the TC genesis location and the tracks realistically well over all main development regions (MDRs), but with notable discrepancies including too weak TC activity over the Arabian Sea and the South China Sea and too strong activity over the North Indian Ocean. Over the western North Pacific and the North Atlantic, MERRA-2 underestimates recurving TCs steered by background westerlies in the mid-latitudes and thereby exaggerates northwestward moving TCs, resulting in an overestimate of the landfall probability in East Asia and North America. Excessive development of TCs over the Bay of Bengal also tends to exaggerate the landfall probability in India. In spite of the discrepancies in the annual TC number, the seasonal variation of TC genesis is realistic in MERRA-2. MERRA-2 also captures the TC intensity relationship between the minimum pressure and the maximum surface wind speed at the mature stage, although the maximum intensity is weaker than in the observations. While MERRA-2 tends to reproduce the interannual variations of the observed TC number and the power dissipation index (PDI), the level of accuracy depends on the region. MERRA-2 describes the changes in the TC genesis region and tracks realistically according to the different phases of El Nino and the Southern Oscillation (ENSO) and the Madden-Julian Oscillation (MJO), although it is less realistic over the North Indian Ocean.

International Best Track Archive for Climate Stewa

Convective Entrainment Rates Estimated from Aura CO and CloudSat/CALIPSO Observations and Comparison with GEOS-5

Entrainment rate (λ) in convective parameterizations remains a sensitive parameter with much uncertainty in model simulations. This study estimates λ using carbon monoxide (CO) measurements jointly from the Microwave Limb Sounder (MLS) and Tropospheric Emission Spectrometer (TES) onboard the Aura satellite, associated with deep convective cases identified by CloudSat and CALIPSO observations. CO is treated as a conserved quantity over convective transport time scales and a simple entraining-plume model is used to derive entrainment rates. The relationships of the observational estimates of λ as a function of convective height, environmental relative humidity and convective available potential energy (CAPE) derived from Atmospheric Infrared Sounder (AIRS) data are compared with those from Goddard Earth Observing System Model (GEOS-5) simulations. Bulk statistics of λ show that the values of λ are predominately below 20% km(exp-1) for deep convection and the occurrence frequency of any λ decreases with increasing λ. Composite λ values are generally lower in the tropics compared to northern mid-latitudes in both observations and the GEOS-5 model. A decrease of λ with increasing convective height is found in both observations and model simulations. We also find that λ tends to decrease with increasing CAPE in the observation-based λ’s and plume-based GEOS-5 λ’s, although the model given λ’s have a non-monotonic relation with CAPE. The observed λ’s have a weak relation with lower-to-mid tropospheric RH, while both GEOS-5 plume-based and given λ increases with increasing RH.

Ryan E Stanfield

The Land Ice Representation in the NASA Goddard Seasonal Forecasting System

Models used in seasonal forecasting systems have aimed for a level of sophistication comparable to Earth system models by incorporating complex physical processes, including those relevant to the polar regions. Here, we examine the polar surface climate in version 3 of the NASA Goddard Earth Observing System Sub-seasonal to Seasonal prediction system (GEOS S2S v3). The model is composed of the ⅟₂̊ resolution Jason4.0version of the GEOS AGCM and the MOM5 ocean model at ⅟₄̊ resolution. The GEOSS2S model incorporates interactive aerosols and two-moment cloud microphysics. As compared to version 2, the model incorporates improved radiative transfer and are designed diurnal cycle atmosphere-ocean interface layer. This study aims to understand how the version 3 model compares with version 2, and with contemporary Earth system models.

Richard I Cullather

Inferred Sea Level Prediction in the NASA GMAO Seasonal Forecasting System

Reliable predictions of sea level anomalies on seasonal timescales with lead times of 1 to 9 months may have relevance to stakeholders – for example, in the advance deployment of resources for coastal flood mitigation. Routine prediction and analysis may also highlight physical processes associated with sea level change and modeling capabilities on seasonal and other timescales. These forecasts may represent interannual changes in the seasonal slope of the ocean surface, teleconnection effects such as the El Niño/Southern Oscillation phenomenon, and variations in seasonal hydrology including precipitation and coastal runoff. Coupled atmosphere/ocean models are routinely used in the seasonal prediction of temperature anomalies, precipitation anomalies, sea ice cover, and climate indices such as the Niño3.4 predictions under the North American Multi-Model Ensemble (NMME) protocol. Within the limits of their configuration, these complex Earth-system models have a potential for depicting regional changes in oceanic column properties, including the sea surface height. Seasonal prediction models generally have no representation of long-term mass contributions from melting land ice, or changes in vertical land motion; their output may be more specifically characterized as predictions of the ocean dynamic sea level. In practice however, the sea surface height prognostic variable is substantially compromised by the forecast model response to initial conditions. Imbalances between the initial, observed hydrologic cycle and the forecast model state produce abrupt adjustments in the model sea surface height. As a result, most seasonal prediction systems employ a constraint on the globally-averaged sea surface height that is applied at each time step. This essentially renders the prognostic sea surface height variable as unserviceable. Several approaches have previously been used to retrieve sea level information from seasonal forecasts beyond the use of the sea surface height variable. Here, we extend a method of relating other prognostic values, including ocean circulation and climate indices, to observed sea level variations. We use the merged altimetry record of the NASA MEaSUREs Gridded Sea Surface Height Anomalies data set and monthly revised local reference gauge observations from the National Oceanography Centre Permanent Service for Mean Sea Level (PSMSL) to evaluate derived prognostic variables from the NASA Global Modeling and Assimilation Office subseasonal-to-seasonal system version 2.1 (GMAO S2S v2.1). We focus on results for the midlatitudes with particular emphasis on US gauge locations. As shown in previous studies, prognostic ENSO-related indices in boreal winter are well correlated with gauge observations for the US west coast, but also for other locations in the southeastern US. Other forecast climate indices such as the North Atlantic Oscillation have relations to sea level that are limited both seasonally and spatially. As expected, surface atmospheric pressure (e.g., inverse barometer effect) is found to be particularly well correlated with observed sea level. We provide a characterization of forecast skill for seasonal sea level with this method.

Richard I Cullather

Improving Boundary Layer Data Assimilation Using Observation Data from Multiple Observing Systems in the NASA GEOS System

The Planetary Boundary Layer (PBL) is a complex interface that mediates energy and moisture exchanges between the Earth’s surface and atmosphere. Accurate simulation and observation of PBL characteristics, such as PBL height and thermodynamic structure, have proven to be a challenge. In our latest efforts, we have focused on improving PBL thermodynamic structure using data from multiple observing systems in the Goddard Earth Observing System (GEOS), developed by the NASA Global Modeling and Assimilation Office (GMAO). We present strategies and results from assimilating PBL height data derived from radiosondes, GNSS radio occultation, space-based lidar (CALIPSO, CATS, IceSat-2), ground-based lidar (MPLNET), and radar wind profilers, including a novel global PBL height analysis dataset. We also discuss the impacts of better representing capping inversions by using PBL height data together with other observations in GEOS through adjustments to the background error covariance. Long-term statistics of the impact of assimilating and utilizing PBL height data in GEOS are presented. In addition, we explore an innovative approach to assimilate GNSS-RO refractivity data in the lower troposphere.

Eun-Gyeong Yang