Engineering PapersSearch

SEARCH · Engineering Papers

Results for “data assimilation numerical modeling ocean circulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

34 records · Page 2

Evaluating the Contribution of NASA Remotely-Sensed Data Sets on a Convection-Allowing Forecast Model

The Short-term Prediction Research and Transition (SPoRT) Center is a collaborative partnership between NASA and operational forecasting partners, including a number of National Weather Service forecast offices. SPoRT provides real-time NASA products and capabilities to help its partners address specific operational forecast challenges. One challenge that forecasters face is using guidance from local and regional deterministic numerical models configured at convection-allowing resolution to help assess a variety of mesoscale/convective-scale phenomena such as sea-breezes, local wind circulations, and mesoscale convective weather potential on a given day. While guidance from convection-allowing models has proven valuable in many circumstances, the potential exists for model improvements by incorporating more representative land-water surface datasets, and by assimilating retrieved temperature and moisture profiles from hyper-spectral sounders. In order to help increase the accuracy of deterministic convection-allowing models, SPoRT produces real-time, 4-km CONUS forecasts using a configuration of the Weather Research and Forecasting (WRF) model (hereafter SPoRT-WRF) that includes unique NASA products and capabilities including 4-km resolution soil initialization data from the Land Information System (LIS), 2-km resolution SPoRT SST composites over oceans and large water bodies, high-resolution real-time Green Vegetation Fraction (GVF) composites derived from the Moderate-resolution Imaging Spectroradiometer (MODIS) instrument, and retrieved temperature and moisture profiles from the Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI). NCAR's Model Evaluation Tools (MET) verification package is used to generate statistics of model performance compared to in situ observations and rainfall analyses for three months during the summer of 2012 (June-August). Detailed analyses of specific severe weather outbreaks during the summer will be presented to assess the potential added-value of the SPoRT datasets and data assimilation methodology compared to a WRF configuration without the unique datasets and data assimilation.

Zavodsky, Bradley T.

Evaluation of Arctic Sea Ice Thickness Simulated by Arctic Ocean Model Intercomparison Project Models

Six Arctic Ocean Model Intercomparison Project model simulations are compared with estimates of sea ice thickness derived from pan-Arctic satellite freeboard measurements (2004-2008); airborne electromagnetic measurements (2001-2009); ice draft data from moored instruments in Fram Strait, the Greenland Sea, and the Beaufort Sea (1992-2008) and from submarines (1975-2000); and drill hole data from the Arctic basin, Laptev, and East Siberian marginal seas (1982-1986) and coastal stations (1998-2009). Despite an assessment of six models that differ in numerical methods, resolution, domain, forcing, and boundary conditions, the models generally overestimate the thickness of measured ice thinner than approximately 2 mand underestimate the thickness of ice measured thicker than about approximately 2m. In the regions of flat immobile landfast ice (shallow Siberian Seas with depths less than 25-30 m), the models generally overestimate both the total observed sea ice thickness and rates of September and October ice growth from observations by more than 4 times and more than one standard deviation, respectively. The models do not reproduce conditions of fast ice formation and growth. Instead, the modeled fast ice is replaced with pack ice which drifts, generating ridges of increasing ice thickness, in addition to thermodynamic ice growth. Considering all observational data sets, the better correlations and smaller differences from observations are from the Estimating the Circulation and Climate of the Ocean, Phase II and Pan-Arctic Ice Ocean Modeling and Assimilation System models.

sea ice

Calibration of Sea Ice Motion from QuikSCAT with those from SSM/I and Buoy

QuikSCAT backscatter and DMSP SSM/I radiance data are used to derive sea ice motion for both the Arctic and Antarctic region using wavelet analysis method. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing satellite data from Synthetic Aperture Radar (SAR). Sea ice motion of the Arctic for the period from October 1999 to March 2000 derived from QuikSCAT and SSM/I data agrees well with that derived from ocean buoys quantitatively. Thus the ice tracking results from QuikSCAT and SSM/I are complement to each other, Then, three sea-ice drift daily results from QuikSCAT, SSM/I, and buoy data can be merged to generate composite maps with more complete coverage of sea ice motion than those from single data source. A series of composite sea ice motion maps for December 1999 show that the major circulation patterns of sea ice motion are changing and shifting significantly within every four days and they are dominated by wind forcing. Sea-ice drift in the summer can not be derived from NSCAT and SSM/I data. In later summer of 1999 (in September), however, QuikSCAT data can provide good sea ice motion information in the Arctic. QuiksCAT can also provide at least partial sea ice motion information until June 15 in early summer 1999. For the Antarctic, case study shows that sea ice motion derived from QuikSCAT data is predominantly forced by and is consistent with wind field derived from QuikSCAT around the polar region. These calibrated/validated results indicate that QuikSCAT, SSM/I, and buoy merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.

U.S. Eastern Continental Shelf Carbon Cycling (USECoS): Modeling, Data Assimilation, and Analysis

Although the oceans play a major role in the uptake of fossil fuel CO2 from the atmosphere, there is much debate about the contribution from continental shelves, since many key shelf fluxes are not yet well quantified: the exchange of carbon across the land-ocean and shelf-slope interfaces, air-sea exchange of CO2, burial, and biological processes including productivity. Our goal is to quantify these carbon fluxes along the eastern U.S. coast using models quantitatively verified by comparison to observations, and to establish a framework for predicting how these fluxes may be modified as a result of climate and land use change. Our research questions build on those addressed with previous NASA funding for the USECoS (U.S. Eastern Continental Shelf Carbon Cycling) project. We have developed a coupled biogeochemical ocean circulation model configured for this study region and have extensively evaluated this model with both in situ and remotely-sensed data. Results indicate that to further reduce uncertainties in the shelf component of the global carbon cycle, future efforts must be directed towards 1) increasing the resolution of the physical model via nesting and 2) making refinements to the biogeochemical model and quantitatively evaluating these via the assimilation of biogeochemical data (in situ and remotely-sensed). These model improvements are essential for better understanding and reducing estimates of uncertainties in current and future carbon transformations and cycling in continental shelf systems. Our approach and science questions are particularly germane to the carbon cycle science goals of the NASA Earth Science Research Program as well as the U.S. Climate Change Research Program and the North American Carbon Program. Our interdisciplinary research team consists of scientists who have expertise in the physics and biogeochemistry of the U.S. eastern continental shelf, remote-sensing data analysis and data assimilative numerical models.

Mannino, Antonio

Tropical Rainfall Measuring Mission (TRMM) and the Future of Rainfall Estimation from Space

Tropical rainfall is important in the hydrological cycle and to the lives and welfare of humans. Three-fourths of the energy that drives the atmospheric wind circulation comes from the latent heat released by tropical precipitation. Recognizing the importance of rain in the tropics, NASA for the U.S.A. and NASDA for Japan have partnered in the design, construction and flight of a satellite mission to measure tropical rainfall and calculate the associated latent heat release. The Tropical Rainfall Measuring Mission (TRMM) satellite was launched on November 27, 1997, and data from all the instruments first became available approximately 30 days after launch. Since then, much progress has been made in the calibration of the sensors, the improvement of the rainfall algorithms and applications of these results to areas such as Data Assimilation and model initialization. TRMM has reduced the uncertainty of climatological rainfall in tropics by over a factor of two, therefore establishing a standard for comparison with previous data sets and climatologies. It has documented the diurnal variation of precipitation over the oceans, showing a distinct early morning peak and this satellite mission has shown the utility of precipitation information for the improvement of numerical weather forecasts and climate modeling. This paper discusses some promising applications using TRMM data and introduces a measurement concept being discussed by NASA/NASDA and ESA for the future of rainfall estimation from space.

Kakar, Ramesh

Arctic Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of DMSP SSM/I (Special Sensor Microwave/Imager) 85 GHz and 37 GHz radiance data, SMMR (Scanning Multichannel Microwave Radiometer) 37 GHz, and NSCAT (NASA Scatterometer) 13.9 GHZ data can be used to obtain daily sea ice drift information for both the northern and southern polar regions. The derived maps of sea ice drift provide both improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars (SAR). Examples of derived ice-drift maps in the Arctic illustrate large-scale circulation reversals within a period of a couple weeks. Comparisons with ice displacements derived from buoys show good quantitative agreement. NASA Scatterometer (NSCAT) 13.9 GHZ data have been also used for wavelet analysis to derive sea-ice drift. First, the 40' incidence-angle, sigma-zero (surface roughness) daily map of whole Arctic region with 25 km of pixel size from satellite's 600 km swath has been constructed. Then, the similar wavelet transform procedure to SSM/I data can be applied. Various scales of wavelet transform and threshold have been tested. By overlaying , neighbor filtering, and block-averaging the results of multiscale wavelet transforms, the final sea ice drift vectors are much smooth and representative to the sea ice motion. This wavelet analysis procedure is robust and can make a major contribution to the understanding of ice motion over large areas at relatively high temporal resolutions. The results of wavelet analysis of SSM/I and NSCAT images and buoy data can be merged by some data fusion techniques and will help to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of surface evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote geographic sources of surface evaporation for precipitation, based on the implementation of three-dimensional constituent tracers of regional water vapor sources (termed water vapor tracers, WVT) in a general circulation model. The major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In the WVT approach, each tracer is associated with an evaporative source region for a prognostic three-dimensional variable that represents a partial amount of the total atmospheric water vapor. The physical processes that act on a WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be predicted within the model simulation, and can be validated against the model's prognostic water vapor. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional sources, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In most North American continental regions, the local source of precipitation is correlated with total precipitation. There is a general positive correlation between local evaporation and local precipitation, but it can be weaker because large evaporation can occur when precipitation is inhibited. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.

Model Simulations of CO2 Transport Using Assimilated Meteorological Fields

The numerical simulation of CO2 transport (and other tracers such as CO, CH4, and biomass burning tracers) in the atmosphere is required to determine the fate of anthropogenic source gases. Estimation of the CO2 exchange between the ocean surface, the terrestrial biosphere, and the atmosphere is of first-order importance to understanding the global carbon cycle and the processes that are most crucial in determining the atmospheric CO2 concentration. Forward transport simulations have been conducted using two-dimensional, time-dependent grids of average surface fluxes (from TRANSCOM) and three-dimensional wind data from a prototype data assimilation system (FV-DAS) run by the Goddard Data Assimilation Office. The objective is to better understand the contribution of meteorological variability to changes in CO2 and other constituents, By accurately accounting for meteorological variability, through use of assimilated winds, we hope to better characterize the distribution of surface sources and sinks (and chemistry where applicable). With assimilated meteorology such chemistry/transport runs provide the basic framework to analyze existing (and proposed) measurement data on a point-by-point, real-time basis. We compare with measured CO2 concentration gradients on a daily, seasonal, regional, and interhemispheric basis to examine the consistency of sources, sinks, and transport formulation. We will also examine the inter-annual variability of atmospheric CO2 due to atmospheric circulation changes using longer runs with assimilated winds.

Kaea, S. R.

Oceanic Observations

For many years, merchant ships and the naval fleets of various countries have been the major source of data over and in the open ocean. Oceanographic research experiments and process studies in the field have also contributed to the climatological data bases for the global ocean, but, for the most part, these have been limited in duration and extent. However, over the last 10 years under the auspices of the World Climate Research Program and the International Geosphere Biosphere Program the role of the oceans in global and climate change has taken on increased significance. This has created a need for a considerably improved understanding of the seasonal, interannual, decadal and longer time-scale variability of the physical and biogeochemical attributes of the global ocean. As a result, over the past 10 years several major international field programs have been implemented and have had a tremendous impact on the number of in situ observations obtained for the global ocean. The Tropical Ocean Global Atmosphere (TOGA) program, the World Ocean Circulation Experiment (WOCE), and the Joint Global Ocean Flux Study (JGOFS) were designed with observational, modelling, and process study components aimed at analyzing different aspects of the ocean's role in the coupled climate system. In parallel with the field programs, continuous space-based observations of sea surface temperature, sea surface topography, and sea surface winds spanning nearly a decade or longer have become a reality. During this same time period, numerical ocean models and computational power have advanced to the point where the oceanographic observations, both in situ and remotely sensed, can be assimilated into numerical ocean models in order to provide a four-dimensional (x-y-z-t) depiction of the evolving state of the global ocean.

Busalacchi, Antonio J.

Water Vapor Tracers as Diagnostics of the Regional Hydrologic Cycle

Numerous studies suggest that local feedback of evaporation on precipitation, or recycling, is a significant source of water for precipitation. Quantitative results on the exact amount of recycling have been difficult to obtain in view of the inherent limitations of diagnostic recycling calculations. The current study describes a calculation of the amount of local and remote sources of water for precipitation, based on the implementation of passive constituent tracers of water vapor (termed water vapor tracers, WVT) in a general circulation model. In this case, the major limitation on the accuracy of the recycling estimates is the veracity of the numerically simulated hydrological cycle, though we note that this approach can also be implemented within the context of a data assimilation system. In this approach, each WVT is associated with an evaporative source region, and tracks the water until it precipitates from the atmosphere. By assuming that the regional water is well mixed with water from other sources, the physical processes that act on the WVT are determined in proportion to those that act on the model's prognostic water vapor. In this way, the local and remote sources of water for precipitation can be computed within the model simulation, and can be validated against the model's prognostic water vapor. Furthermore, estimates of precipitation recycling can be compared with bulk diagnostic approaches. As a demonstration of the method, the regional hydrologic cycles for North America and India are evaluated for six summers (June, July and August) of model simulation. More than 50% of the precipitation in the Midwestern United States came from continental regional tracers, and the local source was the largest of the regional tracers (14%). The Gulf of Mexico and Atlantic 2 regions contributed 18% of the water for Midwestern precipitation, but further analysis suggests that the greater region of the Tropical Atlantic Ocean may also contribute significantly. In general, most North American land regions showed a positive correlation between evaporation and recycling ratio (except the Southeast United States) and negative correlations of recycling ratio with precipitation and moisture transport (except the Southwestern United States). The Midwestern local source is positively correlated with local evaporation, but it is not correlated with water vapor transport. This is contrary to bulk diagnostic estimates of precipitation recycling. In India, the local source of precipitation is a small percentage of the precipitation owing to the dominance of the atmospheric transport of oceanic water. The southern Indian Ocean provides a key source of water for both the Indian continent and the Sahelian region.

Bosilovich, Michael G.

Assimilation of Altimeter Data into a Quasigeostrophic Model of the Gulf Stream System: Assimilation Results - Part 2

The improvement in the climatological behavior of a numerical model as a consequence of the assimilation of surface data is investigated. The model used for this study is a quasigeostrophic (QG) model of the Gulf Stream region. The data that have been assimilated are maps of sea surface height that have been obtained as the superposition of sea surface height variability deduced from the Geosat altimeter measurements and a mean field constructed from historical hydrographic data. The method used for assimilating the data is the nudging technique. Nudging has been implemented in such a way as to achieve a high degree of convergence of the surface model fields toward the observations. Comparisons of the assimilation results with available in situ observations show a significant improvement in the degree of realism of the climatological model behavior, with respect to the model in which no data are assimilated. The remaining discrepancies in the model mean circulation seem to be mainly associated with deficiencies in the mean component of the surface data that are assimilated. On the other hand, the possibility of building into the model more realistic eddy characteristics through the assimilation of the surface eddy field proves very successful in driving components of the mean model circulation that are in relatively good agreement with the available observations. Comparisons with current meter time series during a time period partially overlapping the Geosat mission show that the model is able to 'correctly' extrapolate the instantaneous surface eddy signals to depths of approximately 1500 m. The correlation coefficient between current meter and model time series varies from values close to 0.7 in the top 1500 m to values as low as 0.1-0.2 in the deep ocean.

Capotondi, Antonietta

The Aquarius Mission: Sea Surface Salinity from Space

Aquarius is a new satellite mission concept to study the impact of the global water cycle on the ocean, including the response of the ocean to buoyancy forcing and the subsequent feedback of the ocean on the climate. The measurement objective of Aquarius is sea surface salinity, which reflects the concentration of freshwater at the ocean surface. Salinity affects the dielectric constant of sea water and, consequently, the radiometric emission of the sea surface to space. Rudimentary space observations with an L-band radiometer were first made from Skylab in the mid-70s and numerous aircraft missions of increasing quality and improved technology have been conducted since then. Technology is now available to carry out a global mission, which includes both an accurate L band (1.413 Ghz) radiometer and radar system in space and a global array of in situ observations for calibration and validation, in order to address key NASA Earth Science Enterprise questions about the global cycling of water and the response of the ocean circulation to climate change. The key scientific objectives of Aquarius examine the cycling of water at the ocean's surface, the response of the ocean circulation to buoyancy forcing, and the impact of buoyancy forcing on the ocean's thermal feedback to the climate. Global surface salinity will also improve our ability to model the surface solubility chemistry needed to estimate the air-sea exchange of CO2. In order to meet these science objectives, the NASA Salinity Sea Ice Working Group over the past three years has concluded that the mission measurement goals should be better than 0.2 practical salinity units (psu) accuracy, 100 km resolution, and weekly to revisits. The Aquarius mission proposes to meet these measurement requirements through a real aperture dual-polarized L band radiometer and radar system. This system can achieve the less than 0.1 K radiometric temperature measurement accuracy that is required. A 3 m antenna at approx. 600km altitude in a sun-synchronous orbit and 300 km swath can provide the desired 100 km resolution global coverage every week. Within this decade, it may be possible to combine satellite sea surface salinity measurements with ongoing satellite observations of temperature, surface height, air-sea fluxes; vertical profiles of temperature and salinity from the Argo program; and modern ocean/atmosphere modeling and data assimilation tools, in order to finally address the complex influence of buoyancy on the ocean circulation and climate.

Koblinsky, Chester

Assimilation of Soil Moisture Observations Over Land Improves Analysis and Prediction of Tropical Cyclone Idai

Soil moisture conditions can impact the circulation and structure of a tropical cyclone (TC) when part or all of the circulation is over land. Dry land surface conditions may lead to faster dissipation of a TC over land, whereas very wet conditions may lead to a prolonged maintenance of its intensity. While this relationship is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tbs significantly improves modeled land surface states. Here we evaluate: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP Tb observations. We find that in the analysis with SMAP assimilation, the TC has a better-defined, more aligned vertical structure over land relative to the control run; moreover, the analyzed TC size, as measured by the wind speed radius, better matches the observed TC size. We further find significant reductions in the forecast intensity error and the forecast along-track error, measured against observations. The largest error reductions occur at lead times of 36 to 72 hours, suggesting that the land with its longer memory gains in importance as a source of predictability at this timescale. An investigation of the underlying mechanisms leading to the skill improvements from SMAP data assimilation revealed that the assimilation of SMAP leads to wetter soil moisture conditions and an increased latent heat flux in the SMAP analysis, which results in a TC with higher column-integrated total moisture content and total energy compared to the control analysis.

Jana Kolassa

Hawaii Ocean Mixing Experiment: Program Summary

It is becoming apparent that insufficient mixing occurs in the pelagic ocean to maintain the large scale thermohaline circulation. Observed mixing rates fall a factor of ten short of classical indices such as Munk's "Abyssal Recipe." The growing suspicion is that most of the mixing in the sea occurs near topography. Exciting recent observations by Polzin et al., among others, fuel this speculation. If topographic mixing is indeed important, it must be acknowledged that its geographic distribution, both laterally and vertically, is presently unknown. The vertical distribution of mixing plays a critical role in the Stommel Arons model of the ocean interior circulation. In recent numerical studies, Samelson demonstrates the extreme sensitivity of flow in the abyssal ocean to the spatial distribution of mixing. We propose to study the topographic mixing problem through an integrated program of modeling and observation. We focus on tidally forced mixing as the global energetics of this process have received (and are receiving) considerable study. Also, the well defined frequency of the forcing and the unique geometry of tidal scattering serve to focus the experiment design. The Hawaiian Ridge is selected as a study site. Strong interaction between the barotropic tide and the Ridge is known to take place. The goals of the Hawaiian Ocean Mixing Experiment (HOME) are to quantify the rate of tidal energy loss to mixing at the Ridge and to identify the mechanisms by which energy is lost and mixing generated. We are challenged to develop a sufficiently comprehensive picture that results can be generalized from Hawaii to the global ocean. To achieve these goals, investigators from five institutions have designed HOME, a program of historic data analysis, modeling and field observation. The Analysis and Modeling efforts support the design of the field experiments. As the program progresses, a global model of the barotropic (depth independent) tide, and two models of the baroclinic (depth varying) tide, all validated with near-Ridge data, will be applied, to reveal the mechanisms of tidal energy conversion along the Ridge, and allow spatial and temporal integration of the rate of conversion. Field experiments include a survey to identify "hot spots" of enhanced mixing and barotropic to baroclinic conversion, a Nearfield study identifying the dominant mechanisms responsible for topographic mixing, and a Farfield program which quantifies the barotropic energy flux convergence at the Ridge and the flux divergence associated with low mode baroclinic waves radiation. The difference is a measure of the tidal power available for mixing at the Ridge. Field work is planned from years 2000 through 2002, with analysis and modeling efforts extending through early 2006. If successful, HOME will yield an understanding of the dominant topographic mixing processes applicable throughout the global ocean. It will advance understanding of two central problems in ocean science, the maintenance of the abyssal stratification, and the dissipation of the tides. HOME data will be used to improve the parameterization of dissipation in models which presently assimilate TOPEX-POSEIDON observations. The improved understanding of the dynamics and spatial distribution of mixing processes will benefit future long-term programs such as CLIVAR.

Ray, Richard D.

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa