Search NASA⌕ Search

Engineering topics

Andrea Molod

Publications and source records attributed to Andrea Molod.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Tropical Cyclones in the GEOS-S2S-2 Subseasonal Forecasts

This paper analyzes the climatology, prediction skill, and predictability of tropical cyclones (TCs) in NASA’s Global Earth Observing System Subseasonal to Seasonal (GEOS-S2S) forecast system version 2. GEOS reasonably simulates the number and spatial distribution of TCs compared to observations except in the Atlantic where the model simulates too few TCs due to low genesis rates in the Caribbean Sea and Gulf of Mexico. The environmental conditions, diagnosed through a genesis potential index, do not clearly explain model biases in the genesis rates, especially in the Atlantic. At the storm-scale, GEOS reforecasts replicate several key aspects of the thermodynamic and dynamic structure of observed TCs, such as a warm core and the secondary circulation. The model, however, fails to simulate an off-center eyewall when evaluating vertical velocity, precipitation and moisture. The analysis of prediction skill of TC genesis and occurrence shows that GEOS has comparable skill to other global models in WMO S2S archive and that its skill could be further improved by increasing the ensemble size. After calibration, GEOS forecasts are skillful in the Western North Pacific and Southern Indian Ocean up to 20 days in advance. A model-based predictability analysis demonstrates the importance of the Madden-Julian Oscillation (MJO) as a source of predictability of TC occurrence beyond the 14 day lead-time. Forecasts initialized under strong MJO conditions show evidence of predictability beyond week 3. However, due to model biases in the forecast distribution there are notable gaps between MJO-related prediction skill and predictability which require further study.

Jorge L. Garcia-Franco↗
NASA NTRS2023

Seasonal Forecasting Skill for the High Mountain Asia Region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorological relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks-to-months is an area of active research and development. NASA’s 15 Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1-3 month lead times in the HMA region, including a portion of the Indian Subcontinent, during the retrospective forecast period, 1981-2016. To assess forecast skill, we evaluate 2-m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and 20 Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are 25 evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system’s ability to forecast HMA hydrometeorology.

GEOS↗
NASA NTRS2022

Local Air-Sea Interactions at Ocean Mesoscale in Western Boundary Currents

We present results from a new, global, high-resolution (∼3-km for ocean and ∼6-km for atmosphere) realistic earth system simulation. This simulation allows us to examine aspects of small-scale air-sea interaction beyond what previous studies have reported. Our study focuses on recurring intermittent wind events in the Gulf Stream region. These events induce local air-sea heat fluxes above Sea Surface Temperature (SST) anomalies with horizontal scales smaller than 500km. In particular, strong latent heat bursts above warm SST anomalies are observed during these wind events. We show that such wind events are associated with a secondary circulation that acts to fuel the latent heat bursts by transferring dry air and momentum down to the surface. The intensity of this secondary circulation is related to the strength of small-scale SST fronts that border SST anomalies. The study of such phenomena requires high-resolution in both the atmospheric and oceanic components of the model.

Ehud Strobach↗
NASA NTRS2022

Earth System Model Parameter Adjustment Using a Green's Functions Approach

We demonstrate the practicality and effectiveness of using a Green's functions estimation approach for adjusting uncertain parameters in an Earth system model (ESM). This estimation approach has previously been applied to an intermediate-complexity climate model and to individual ESM components, e.g., ocean, sea ice, or carbon cycle components. Here, the Green's functions approach is applied to a state-of-the-art ESM that comprises a global atmosphere/land configuration of the Goddard Earth Observing System (GEOS) coupled to an ocean and sea ice configuration of the Massachusetts Institute of Technology general circulation model (MITgcm). Horizontal grid spacing is approximately 110 km for GEOS and 37–110 km for MITgcm. In addition to the reference GEOS-MITgcm simulation, we carried out a series of model sensitivity experiments, in which 20 uncertain parameters are perturbed. These “control” parameters can be used to adjust sea ice, microphysics, turbulence, radiation, and surface schemes in the coupled simulation. We defined eight observational targets: sea ice fraction, net surface shortwave radiation, downward longwave radiation, near-surface temperature, sea surface temperature, sea surface salinity, and ocean temperature and salinity at 300 m. We applied the Green's functions approach to optimize the values of the 20 control parameters so as to minimize a weighted least-squares distance between the model and the eight observational targets. The new experiment with the optimized parameters resulted in a total cost reduction of 9 % relative to a simulation that had already been adjusted using other methods. The optimized experiment attained a balanced cost reduction over most of the observational targets. We also report on results from a set of sensitivity experiments that are not used in the final optimized simulation but helped explore options and guided the optimization process. These experiments include an assessment of sensitivity to the number of control parameters and to the selection of observational targets and weights in the cost function. Based on these sensitivity experiments, we selected a specific definition for the cost function. The sensitivity experiments also revealed a decreasing overall cost as the number of control variables was increased. In summary, we recommend using the Green's functions estimation approach as an additional fine-tuning step in the model development process. The method is not a replacement for modelers' experience in choosing and adjusting sensitive model parameters. Instead, it is an additional practical and effective tool for carrying out final adjustments of uncertain ESM parameters.

Green's Function↗
NASA NTRS2021

Impacts of the Eruption of Mount Pinatubo on Surface Temperatures and Precipitation Forecasts With the NASA GEOS Subseasonal-to-Seasonal System

A contemporary seasonal forecasting system is used to study the impacts of a volcanic sulfate injection into the stratosphere on the seasonal forecasts for surface temperatures, the El Niño Southern Oscillation (ENSO), and precipitation. The focus is a case study of the June 1991 eruption of Mt. Pinatubo, Philippines and the period from July 1991 to February 1992. Version 2 of the Goddard Earth Observing System (GEOS) subseasonal-to-seasonal (S2S) forecasting system is used in this study. GEOS-S2S includes the GOddard Chemistry, Aerosols, Radiation and Transport (GOCART) aerosol module, which allows to prognostically simulate aerosol distributions. GOCART is coupled to the radiation and cloud modules to include the impact of the eruption on forecasted radiation and precipitation. The coupled GEOS-S2S system was initialized in May 1991 with fields based on observations to produce ten-member 9-month forecasts with and without the volcanic sulfur injection. The results of these ensemble experiments demonstrate that including Mt. Pinatubo in seasonal forecasts would improve the forecasts of the 1991–1992 global mean temperature and precipitation but worsen the forecast of ENSO by strengthening of El Niño beyond what showed in observations. Most significant changes in the forecasts of temperatures and precipitation are limited to the tropics. The only land area where the inclusion of Pinatubo significantly lowered the forecasted precipitation is tropical Africa.

Mt. Pinatubo↗
NASA NTRS2021

Tendency Bias Correction in the GEOS AGCM

In this report, we summarize the results of global and regional TBC applied to the GEOS AGCM(the same model used to generate MERRA-2 though run at lower resolution, referred to as the M2_AGCM) employing the analysis increments (analysis minus forecast) generated by MERRA-2 to calculate the TBC terms. Extending the results of Chang et al. (2019)and Schubert et al. (2019), we examine in more detail the seasonality of the TBC impacts, including a deeper look into the reasons for why TBC appears to be least effective during boreal winter. We also look into the ability of TBC to correct other diagnostic fields such as cloudiness and surface fluxes and, more generally, attempt to provide some insight into why TBC is effective in correcting long term climate biases. Looking beyond TBC, we also present some preliminary results from an extension of the TBC approach in which we include a state-dependent term

Yehui Chang↗
NASA NTRS2021

A Phenomenon-Based Decomposition of Model-Based Estimates of Boreal Winter ENSO Variability

Climate models are now routinely being used to simulate and predict climate variability on time scales ranging from sub-seasonal to seasonal and longer. As such, there are now long histories of such simulations and predictions spanning multiple decades and multiple ensemble members, both of which are crucial for separating climate signal from climate noise. A key focus of such runs has been the El Niño-Southern Oscillation (ENSO), spurred by recent improvements in our ability to predict such events, though questions remain as to how well climate models do beyond simply always predicting the “canonical” atmospheric response to an ENSO event—something simple statistical models already do reasonably well. This is a critical issue that needs addressing, given the importance of event-to-event differences for predicting regional impacts of ENSO teleconnections, and the need to justify the expense of running sophisticated climate models. Unfortunately, current diagnostic tools are not well suited for quantifying the different sources of variability associated with specific phenomena such as ENSO. More generally, while much effort has focused on addressing model bias, less has been done to address errors in second moment statistics—an issue whose importance is gaining increased attention particularly as we build climate prediction systems capable of taking advantage of forecasts of opportunity—a capability that requires reliable estimates of forecast uncertainty. In this report, we outline a phenomenon-based statistical decomposition of climate variance(in essence a detailed variance budget)that is specifically tailored to address the above questions by separating the variability (both the signal and noise) into that tied to the long-term average impact of a particular phenomenon(the composite mean) and the event-to-event(E2E) variability about the composite mean. In addition, we provide related decompositions of the correlations that allow us to quantify how much of the agreement with observations (the skill) comes from the composite mean as opposed to from the E2Evariability. As an example, we present the results of such a decomposition for ENSO based on simulations with the GEOS atmospheric general circulation model (AGCM), with a focus on the monthly mean impacts over North America during boreal winter(December –March). Here we take advantage of existing GEOS AGCM simulations that were produced as companion simulations to MERRA-2for the period 1980-2016. Comparisons are made throughout with MERRA-2.

ENSO↗
NASA NTRS2020

Windows of Opportunity for Skillful Forecasts Subseasonal to Seasonal and Beyond

There is high demand and a growing expectation for predictions of environmental conditions that go beyond 0–14-day weather forecasts with outlooks extending to one or more seasons and beyond. This is driven by the needs of the energy, water management, and agriculture sectors, to name a few. There is an increasing realization that, unlike weather forecasts, prediction skill on longer time scales can leverage specific climate phenomena or conditions for a predictable signal above the weather noise. Currently, it is understood that these conditions are intermittent in time and have spatially heterogeneous impacts on skill, hence providing strategic windows of opportunity for skillful forecasts. Research points to such windows of opportunity, including El Niño or La Niña events, active periods of the Madden–Julian oscillation, disruptions of the stratospheric polar vortex, when certain large-scale atmospheric regimes are in place, or when persistent anomalies occur in the ocean or land surface. Gains could be obtained by increasingly developing prediction tools and metrics that strategically target these specific windows of opportunity. Across the globe, reevaluating forecasts in this manner could find value in forecasts previously discarded as not skillful. Users’ expectations for prediction skill could be more adequately met, as they are better aware of when and where to expect skill and if the prediction is actionable. Given that there is still untapped potential, in terms of process understanding and prediction methodologies, it is safe to expect that in the future forecast opportunities will expand. Process research and the development of innovative methodologies will aid such progress.

Annarita Mariotti↗
NASA NTRS2019

Annual Cycle of Planetary Boundary Layer Heights estimated from NOAA Profiler Network Data

The Planetary Boundary Layer (PBL) height over the Unites States Great and Central Plains during 1992-2012 is examined here using a combination of Wind Profiler-derived (WP) PBL height estimates and reanalysis fields from the Modern Era Retrospective Reanalysis Version 2 (MERRA-2). The combined analysis allows process study of the reasons behind the monthly mean behavior of the observed PBL heights under clear-sky conditions as well as the PBL height variability. WP PBL height monthly mean annual cycles were grouped into general categories of behavior, each analyzed using MERRA-2 fields of sensible and latent heat flux, surface temperature, net radiation and soil moisture. In the ’canonical’ category the latent heat plays little role in the determination of the monthly mean PBL height, and it follows the annual cycle of the net radiation. In the other categories, precipitation and latent heat flux had more influence in setting the annual cycle. An analysis of variance revealed that the role of latent heat in determining the PBL height variations is large (explaining up to 40% of PBL height variability) even in the ’canonical’ category for which latent heat played no role in setting the monthly mean. In other categories the latent heat explained up to 80% of the PBL height variations. The amount of that influence is shown to be related to the variability of column soil moisture.

MERRA-2↗
NASA NTRS2018

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗
NASA NTRS2017

Direct Estimation of the Global Distribution of Vertical Velocity Within Cirrus Clouds

Cirrus clouds determine the radiative balance of the upper troposphere and the transport of water vapor across the tropopause.The representation of vertical wind velocity, W, in atmospheric models constitutes the largest source of uncertainty in the calculation of the cirrus formation rate. Using global atmospheric simulations with a spatial resolution of 7 kilometers we obtain for the first time a direct estimate of the distribution of W at the scale relevant for cirrus formation, validated against long-term observations at two different ground sites. The standard deviation in W, sigma (sub W), varies widely over the globe with the highest values resulting from orographic uplift and convection, and the lowest occurring in the Arctic. Globally about 90 of the simulated sigma (sub W) values are below 0.1 meters per second and about one in 10 (sup 4) cloud formation events occur in environments with sigma (sub W) greater than 0.8 meters per second. Combining our estimate with reanalysis products and an advanced cloud formation scheme results in lower homogeneous ice nucleation frequency than previously reported, and a decreasing average ice crystal concentration with decreasing temperature. These features are in agreement with observations and suggest that the correct parameterization of sigma (sub W) is critical to simulate realistic cirrus properties.

Radiative Balance↗