Results for “Global ocean state estimation” — Page 6

NASA NTRS2003

Merging Ocean Color Data From Multiple Missions

Oceanic phytoplankton may play an important role in the cycling of carbon on the Earth, through the uptake of carbon dioxide in the process of photosynthesis. Although they are ubiquitous in the global oceans, their abundances and dynamics are difficult to estimate, primarily due to the vast spatial extent of the oceans and the short time scales over which their abundances can change. Consequently, the effects of oceanic phytoplankton on biogeochemical cycling, climate change, and fisheries are not well known. In response to the potential importance of phytoplankton in the global carbon cycle and the lack of comprehensive data, NASA and the international community have established high priority satellite missions designed to acquire and produce high quality ocean color data (Table 6.1). Ten of the missions are routine global observational missions: the Ocean Color and Temperature Sensor (OCTS), the Polarization and Directionality of the Earth's Reflectances sensor (POLDER), Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectrometer-AM (MODIS-AM), Medium Resolution Imaging Spectrometer (MERIS), Global Imager (GLI), MODIS-PM, Super-GLI (S-GLI), and the Visible/Infrared Imager and Radiometer Suite (VIIRS) on the NPOESS Preparatory Project (NPP) and the National Polar-orbiting Operational Environmental Satellite System (NPOESS). In addition, there are several other missions capable of providing ocean color data on smaller scales. Most of these missions contain the spectral band complement considered necessary to derive oceanic chlorophyll concentrations and other related parameters. Many contain additional bands that can provide important ancillary information about the optical and biological state of the oceans.

Gregg, Watson W.↗

NASA NTRSDate not supplied

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↗

NASA NTRSDate not supplied

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↗

NASA NTRS2011

The NASA Soil Moisture Active Passive (SMAP) Mission - Algorithm and Cal/Val Activities and Synergies with SMOS and Other L-Band Missions

NASA's Soil Moisture Active Passive (SMAP) mission, planned for launch in late 2014, has as its key measurement objective the frequent, global mapping of near-surface soil moisture and its freeze-thaw state. SMAP soil moisture and freeze/thaw measurements at 10 km and 3 km resolutions respectively, would enable significantly improved estimates of water, energy and carbon transfers between the land and atmosphere. Soil moisture control of these fluxes is a key factor in the performance of atmospheric models used for weather forecasts and climate projections Soil moisture measurements are also of great importance in assessing floods and for monitoring drought. In addition, observations of soil moisture and freeze/thaw timing over the boreal latitudes can help reduce uncertainties in quantifying the global carbon balance. The SMAP measurement concept utilizes an L-band radar and radiometer sharing a rotating 6-meter mesh reflector antenna. The SMAP radiometer and radar flight hardware and ground processing designs are incorporating approaches to identify and mitigate potential terrestrial radio frequency interference (RFI). The radar and radiometer instruments are planned to operate in a 680 km polar orbit, viewing the surface at a constant 40-degree incidence angle with a 1000-km swath width, providing 3-day global coverage. Data from the instruments would yield global maps of soil moisture and freeze/thaw state to be provided at 10 km and 3 km resolutions respectively, every two to three days. Plans are to provide also a radiometer-only soil moisture product at 40-km spatial resolution. This product and the underlying brightness temperatures have characteristics similar to those provided by the Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are unique opportunities for common data product development and continuity between the two missions. SMAP also has commonalities with other satellite missions having L-band radiometer and/or radar sensors applicable to soil moisture measurement, such as Aquarius, SAO COM, and ALOS-2. The algorithms and data products for SMAP are being developed in the SMAP Science Data System (SDS) Testbed. The algorithms are developed and evaluated in the SDS Testbed using simulated SMAP observations as well as observational data from current airborne and spaceborne L-band sensors including SMOS. The SMAP project is developing a Calibration and Validation (Cal/Val) Plan that is designed to support algorithm development (pre-launch) and data product validation (post-launch). A key component of the Cal/Val Plan is the identification, characterization, and instrumentation of sites that can be used to calibrate and validate the sensor data (Level I) and derived geophysical products (Level 2 and higher). In this presentation we report on the development status of the SMAP data product algorithms, and the planning and implementation of the SMAP Cal/Val program. Several components of the SMAP algorithm development and Cal/Val plans have commonality with those of SMOS, and for this reason there are shared activities and resources that can be utilized between the missions, including in situ networks, ancillary data sets, and long-term monitoring sites.

Njoku, Eni↗

NASA NTRS2020

Development and Evaluation of the Raindrop Size Distribution Parameters for the NASA Global Precipitation Measurement Mission Ground Validation Program

The National Aeronautics and Space Administration Global Precipitation Measurement (GPM) mission ground validation program uses dual-polarization radar moments to estimate raindrop size distribution (DSD) parameters, the mass-weighted mean drop diameter Dmass, and normalized intercept parameter NW, to validate the GPM Core Observatory–derived DSD parameters. The disdrometer-based Dmass and NW are derived through empirical relationships between Dmass and differential reflectivity ZDR, and between NW, reflectivity ZH, and Dmass. This study employs large datasets collected from two-dimensional video disdrometers (2DVD) during six different field studies to derive the requisite empirical relationships. The uncertainty of the derived Dmass(ZDR) relationship is evaluated through comparisons of 2DVD-calculated and ZDR-estimated Dmass, where ZDR is calculated directly from 2DVD observations. Similarly, the uncertainty of the NW(ZH, Dmass) relationship is evaluated through 2DVD-calculated and Dmass and ZH-estimated NW, where Dmass and ZH are directly calculated from 2DVD observations. This study also presents the sensitivity of Dmass(ZDR) relationships to climate regime and to disdrometer type after developing three additional Dmass(ZDR) relationships from second-generation Particle Size Velocity (PARSIVEL2) disdrometer (P2) observations collected in the Pacific Northwest, in Iowa, and at Kwajalein Atoll in the tropical Pacific Ocean. The application of P2-derived Dmass(ZDR) relationship based on precipitation in the northwestern United States to P2 observations collected over the tropical ocean resulted in the highest error among comparisons of the three datasets.

Ali Tokay↗

NASA NTRS2010

Comparing Global Atmospheric CO2 Flux and Transport Models with Remote Sensing (and Other) Observations

We report recent progress derived from comparison of global CO2 flux and transport models with new remote sensing and other sources of CO2 data including those from satellite. The overall objective of this activity is to improve the process models that represent our understanding of the workings of the atmospheric carbon cycle. Model estimates of CO2 surface flux and atmospheric transport processes are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, to provide the basic framework for carbon data assimilation, and ultimately for future projections of carbon-climate interactions. Models can also be used to test consistency within and between CO2 data sets under varying geophysical states. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 2000 through 2009. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at 1x1 degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-3), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to remote sensing observations from TCCON, GOSAT, and AIRS as well as relevant in situ observations. Examples of the influence of key process representations are shown from both forward and inverse model comparisons. We find that the model can resolve much of the synoptic, seasonal, and interannual variability in the observations, although reasons for persistent discrepancies in northern hemisphere vegetation uptake are examined. At this time, we do not find any serious shortcomings in the model transport representation, but this is still the subject of close scrutiny. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in CO2 fluxes and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗

NASA NTRS2011

TRMM-Based Lightning Climatology

Gridded climatologies of total lightning flash rates seen by the spaceborne Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) have been updated. OTD collected data from May 1995 to March 2000. LIS data (equatorward of about 38 deg) has been added for 1998-2010. Flash counts from each instrument are scaled by the best available estimates of detection efficiency. The long LIS record makes the merged climatology most robust in the tropics and subtropics, while the high latitude data is entirely from OTD. The mean global flash rate from the merged climatology is 46 flashes per second. The peak annual flash rate at 0.5 deg scale is 160 fl/square km/yr in eastern Congo. The peak monthly average flash rate at 2.5 scale is 18 fl/square km/mo, from early April to early May in the Brahmaputra Valley of far eastern India. Lightning decreases in this region during the monsoon season, but increases further north and west. A monthly average peak from early August to early September in northern Pakistan also exceeds any monthly averages from Africa, despite central Africa having the greatest yearly average. Most continental regions away from the equator have an annual cycle with lightning flash rates peaking in late spring or summer. The main exceptions are India and southeast Asia, with springtime peaks in April and May. For landmasses near the equator, flash rates peak near the equinoxes. For many oceanic regions, the peak flash rates occur in autumn. This is particularly noticeable for the Mediterranean and North Atlantic. Landmasses have a strong diurnal cycle of lightning, with flash rates generally peaking between 3-5 pm local solar time. The central United States flash rates peak later, in late evening or early night. Flash rates peak after midnight in northern Argentina. These regions are known for large, intense, long-lived mesoscale convective systems.

Cecil, Daniel J.↗

NASA NTRS2013

North American Tropical Cyclone Landfall and SST: A Statistical Model Study

A statistical-stochastic model of the complete life cycle of North Atlantic (NA) tropical cyclones (TCs) is used to examine the relationship between climate and landfall rates along the North American Atlantic and Gulf Coasts. The model draws on archived data of TCs throughout the North Atlantic to estimate landfall rates at high geographic resolution as a function of the ENSO state and one of two different measures of sea surface temperature (SST): 1) SST averaged over the NA subtropics and the hurricane season and 2) this SST relative to the seasonal global subtropical mean SST (termed relSST). Here, the authors focus on SST by holding ENSO to a neutral state. Jackknife uncertainty tests are employed to test the significance of SST and relSST landfall relationships. There are more TC and major hurricane landfalls overall in warm years than cold, using either SST or relSST, primarily due to a basinwide increase in the number of storms. The signal along the coast, however, is complex. Some regions have large and significant sensitivity (e.g., an approximate doubling of annual major hurricane landfall probability on Texas from -2 to +2 standard deviations in relSST), while other regions have no significant sensitivity (e.g., the U.S. mid-Atlantic and Northeast coasts). This geographic structure is due to both shifts in the regions of primary TC genesis and shifts in TC propagation.

El Nino↗

NASA NTRS2015

The Observed State of the Water Cycle in the Early Twenty-First Century

This study quantifies mean annual and monthly fluxes of Earth's water cycle over continents and ocean basins during the first decade of the millennium. To the extent possible, the flux estimates are based on satellite measurements first and data-integrating models second. A careful accounting of uncertainty in the estimates is included. It is applied within a routine that enforces multiple water and energy budget constraints simultaneously in a variational framework in order to produce objectively determined optimized flux estimates. In the majority of cases, the observed annual surface and atmospheric water budgets over the continents and oceans close with much less than 10% residual. Observed residuals and optimized uncertainty estimates are considerably larger for monthly surface and atmospheric water budget closure, often nearing or exceeding 20% in North America, Eurasia, Australia and neighboring islands, and the Arctic and South Atlantic Oceans. The residuals in South America and Africa tend to be smaller, possibly because cold land processes are negligible. Fluxes were poorly observed over the Arctic Ocean, certain seas, Antarctica, and the Australasian and Indonesian islands, leading to reliance on atmospheric analysis estimates. Many of the satellite systems that contributed data have been or will soon be lost or replaced. Models that integrate ground-based and remote observations will be critical for ameliorating gaps and discontinuities in the data records caused by these transitions. Continued development of such models is essential for maximizing the value of the observations. Next-generation observing systems are the best hope for significantly improving global water budget accounting.

Numerical analysis/modeling↗

NASA NTRS2009

Developing a Model-based Capability to Analyze Requirements for the Climate Observing System

Models are foundational for estimating states of the earth's climate system, both as tools to extrapolate information in time and space, and as observation 'operators' used to relate what is analyzed and predicted to what is observed. Expanding the simulation approach further, observing system simulation experiments (OSSEs) are designed to mimic the complete process of analyzing the climate state by replacing real observations with entirely simulated ones determined from a model-based depiction of nature. OSSEs provide a framework to 'fly' simulated satellite instruments through a synthetic atmosphere and investigate the trade-spaces of measurements for various satellite configurations and sampling strategies, and assess their measurement impact on modeling and forecasting capabilities. Such a tool is a crucial but as yet unfulfilled need for future mission selection and design. The components of a state-of-the-art OSSE system are being assembled at the Global Modeling and Assimilation Office (GMAO, Code 610.1) at NASA/GSFC, leveraging on the GMAO's existing modeling and data assimilation infrastructure for numerical weather prediction (NWP). The OSSE framework is based on the GMAO's Goddard Earth Observing System atmospheric general circulation model, version 5 (GEOS-5) and the Gridpoint Statistical Interpolation (GSI) observational analysis scheme, combined with the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model developed by the Atmospheric Chemistry and Dynamics Branch (Code 613.3). This system is an evolving, key component of Goddard's planned development of an Integrated Earth System Analysis (IESA) capability, which will bring together into a single, fully interactive system Goddard's modeling and assimilation efforts in atmosphere, ocean and chemistry and aerosols to provide a comprehensive analysis and prediction system for weather and climate In addition to providing a state-of-the-art capability for assimilating current observation types, GEOS-5, and the future IESA, provide the capability to identify the need for, and assess the potential impact of, future observing systems under consideration for improving weather and climate prediction.

Gelaro, Ronald↗

NASA NTRS2020

Toward Coupled Data Assimilation in NASA’s GEOS: Developments in the Ocean Context

The Global Modeling & Assimilation Office (GMAO) at NASA GSFC produces analyses and predictions of the Earth system using various configurations of the Goddard Earth Observing System (GEOS) model and assimilation system. The current sub-seasonal-to-seasonal prediction system (GEOS-S2S) is based on a coupled atmosphere-ocean-land-ice configuration of GEOS which includes the Modular Ocean Model version 5 (MOM5) run at approximately 50-km resolution and a de-coupled OI-based ocean analysis that uses an initialization of MOM5 forced by the MERRA-2 reanalysis. GMAO will soon implement an updated GEOS-S2S system that will run at 25-km resolution and adopt aspects of the hybrid four-dimensional ensemble-variational (H4DEnVar) system already running in the production-version atmospheric analysis system, including a Local Ensemble Transform Kalman Filter (LETKF) to provide initial conditions for the oceanic state. This presentation will focus on developments to sustain the GMAO's systems on longer time horizons, where more radical transformations will be required to adapt to advanced computing environments, higher resolution and more diverse model components, and new observations for the Earth system. Results will describe progress toward a version of the GEOS coupled system that will be based around the Joint Effort for Data assimilation Integration (JEDI) framework being developed within Joint Center for Satellite Data Assimilation (JCSDA) and include an updated ocean model, MOM6. Discussion will focus specifically on the use of a Unified Forward Operator (UFO) for simulating observations and the Object Oriented Prediction System (OOPS) for providing the state estimate. These features are being developed as a multi-agency effort under the auspices of the JCSDA and are being adopted in the GMAO for all its applications of coupled data assimilation including S2S, numerical weather prediction, and reanalysis.

Mahajan, Rahul↗

NASA NTRS2011

Global Electric Circuit Implications of Combined Aircraft Storm Electric Current Measurements and Satellite-Based Diurnal Lightning Statistics

Using rotating vane electric field mills and Gerdien capacitors, we measured the electric field profile and conductivity during 850 overflights of thunderstorms and electrified shower clouds (ESCs) spanning regions including the Southeastern United States, the Western Atlantic Ocean, the Gulf of Mexico, Central America and adjacent oceans, Central Brazil, and the South Pacific. The overflights include storms over land and ocean, and with positive and negative fields above the storms. Over three-quarters (78%) of the land storms had detectable lightning, while less than half (43%) of the oceanic storms had lightning. Integrating our electric field and conductivity data, we determined total conduction currents and flash rates for each overpass. With knowledge of the storm location (land or ocean) and type (with or without lightning), we determine the mean currents by location and type. The mean current for ocean thunderstorms is 1.7 A while the mean current for land thunderstorms is 1.0 A. The mean current for ocean ESCs 0.41 A and the mean current for land ESCs is 0.13 A. We did not find any significant regional or latitudinal based patterns in our total conduction currents. By combining the aircraft derived storm currents and flash rates with diurnal flash rate statistics derived from the Lightning Imaging Sensor (LIS) and Optical Transient Detector (OTD) low Earth orbiting satellites, we reproduce the diurnal variation in the global electric circuit (i.e., the Carnegie curve) to within 4% for all but two short periods of time. The agreement with the Carnegie curve was obtained without any tuning or adjustment of the satellite or aircraft data. Given our data and assumptions, mean contributions to the global electric circuit are 1.1 kA (land) and 0.7 kA (ocean) from thunderstorms, and 0.22 kA (ocean) and 0.04 (land) from ESCs, resulting in a mean total conduction current estimate for the global electric circuit of 2.0 kA. Mean storm counts are 1100 for land thunderstorms, 530 for ocean ESCs, 390 for ocean thunderstorms, and 330 for land ESCs.

Mach, Douglas M.↗

NASA NTRS2011

SMOS/SMAP Synergy for SMAP Level 2 Soil Moisture Algorithm Evaluation

Soil Moisture Active Passive (SMAP) satellite has been proposed to provide global measurements of soil moisture and land freeze/thaw state at 10 km and 3 km resolutions, respectively. SMAP would also provide a radiometer-only soil moisture product at 40-km spatial resolution. This product and the supporting brightness temperature observations are common to both SMAP and European Space Agency's Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are opportunities for synergies between the two missions. These include exploiting the data for calibration and validation and establishing longer term L-band brightness temperature and derived soil moisture products. In this investigation we will be using SMOS brightness temperature, ancillary data, and soil moisture products to develop and evaluate a candidate SMAP L2 passive soil moisture retrieval algorithm. This work will begin with evaluations based on the SMOS product grids and ancillary data sets and transition to those that will be used by SMAP. An important step in this analysis is reprocessing the multiple incidence angle observations provided by SMOS to a global brightness temperature product that simulates the constant 40 degree incidence angle observations that SMAP will provide. The reprocessed brightness temperature data provide a basis for evaluating different SMAP algorithm alternatives. Several algorithms are being considered for the SMAP radiometer-only soil moisture retrieval. In this first phase, we utilized only the Single Channel Algorithm (SCA), which is based on the radiative transfer equation and uses the channel that is most sensitive to soil moisture (H-pol). Brightness temperature is corrected sequentially for the effects of temperature, vegetation, roughness (dynamic ancillary data sets) and soil texture (static ancillary data set). European Centre for Medium-Range Weather Forecasts (ECMWF) estimates of soil temperature for the top layer (as provided as part of the SMOS ancillary data) were used to correct for surface temperature effects and to derive microwave emissivity. ECMWF data were also used for precipitation forecasts, presence of snow, and frozen ground. Vegetation options are described below. One year of soil moisture observations from a set of four watersheds in the U.S. were used to evaluate four different retrieval methodologies: (1) SMOS soil moisture estimates (version 400), (2) SeA soil moisture estimates using the SMOS/SMAP data with SMOS estimated vegetation optical depth, which is part of the SMOS level 2 product, (3) SeA soil moisture estimates using the SMOS/SMAP data and the MODIS-based vegetation climatology data, and (4) SeA soil moisture estimates using the SMOS/SMAP data and actual MODIS observations. The use of SMOS real-world global microwave observations and the analyses described here will help in the development and selection of different land surface parameters and ancillary observations needed for the SMAP soil moisture algorithms. These investigations will greatly improve the quality and reliability of this SMAP product at launch.

Bindlish, Rajat↗

NASA NTRS2022

Detecting the Greatest Changes in Global Satellite-Based Precipitation Observations

In recent years, the analysis of abrupt and non-abrupt changes in precipitation has received much attention due to the importance of climate change-related issues (e.g., extreme climate events). In this study, we used a novel segmentation algorithm, DBEST (Detecting Breakpoints and Estimating Segments in Trend), to analyze the greatest changes in precipitation using a monthly pixel-based satellite precipitation dataset (TRMM 3B43) at three different scales: (i) global, (ii) continental, and (iii) climate zone, during the 1998–2019 period. We found significant breakpoints, 14.1%, both in the form of abrupt and non-abrupt changes, in the global scale precipitation at the 0.05 significance level. Most of the abrupt changes were observed near the Equator in the Pacific Ocean and Asian continent, relative to the rest of the globe. Most detected breakpoints occurred during the 1998–1999 and 2009–2011 periods on the global scale. The average precipitation change for the detected breakpoint was ±100 mm, with some regions reaching ±3000 mm. For instance, most portions of northern Africa and Asia experienced major changes of approximately +100 mm. In contrast, most of the South Pacific and South Atlantic Ocean experienced changes of −100 mm during the studied period. Our findings indicated that the larger areas of Africa (23.9%), Asia (22.9%), and Australia (15.4%) experienced significant precipitation breakpoints compared to North America (11.6%), South America (9.3%), Europe (8.3%), and Oceania (9.6%). Furthermore, we found that the majority of detected significant breakpoints occurred in the arid (31.6%) and polar (24.1%) climate zones, while the least significant breakpoints were found for snow-covered (11.5%), equatorial (7.5%), and warm temperate (7.7%) climate zones. Positive breakpoints’ temporal coverage in the arid (54.0%) and equatorial (51.9%) climates were more than those in other climates zones. Here, the findings indicated that large areas of Africa and Asia experienced significant changes in precipitation (−250 to +250 mm). Compared to the average state (trend during a specific period), the greatest changes in precipitation were more abrupt and unpredictable, which might impose a severe threat to the ecology, environment, and natural resources.

Breakpoint↗

NASA NTRS2023

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust↗

NASA NTRS1999

GPS Altimetry

The advent of satellite altimetry has greatly improved our ability to observe global ocean circulation. However, the swath of a single, nadir-viewing satellite altimeter is only a few km and the track spacing is several hundred km to resolve the two-dimensional structure of ocean eddies. Our goal is to increase spatial and temporal coverage by monitoring Global Positioning System (GPS) signals reflected from the ocean. A constellation of spacecraft would each carry a GPS receiver capable of recording 8 reflections simultaneously. The reflections are well distributed in azimuth and elevation and can be tracked continuously while the satellite is in view, and another is then acquired, as illustrated below. The diagram depicts a new approach at altimetry measurements where ocean surface reflected GPS signals are simultaneously tracked and processed in a GPS flight receiver in space. The reflected GPS signals from the ocean must be compared precisely with the direct GPS signals in order to infer the characteristics of the ocean from the combined data set. Understanding the features and accuracy of GPS altimetry measurement is crucial to establishing its suitability for oceanography. Preliminary work has enabled us to theoretically model the signal output of the correlator for a variety of system parameters such as wind speed (sea roughness), receiver height, incidence angle, receiver range and Doppler filter bandwidth and antenna gain. Expected signal-to-noise ratio has been estimated from which we have inferred, to a first approximation, the basic receiver gain requirements for a space-based altimeter and the expected range raw error. In 1998, work on a different task led to the extraction of the first reflected GPS signal observed from a spaceborne receiver during the 1995 Space Transportation System-68 (STS-68) Shuttle Radar Laboratory-2 (SRL-2) high resolution synthetic aperture radar mission. Good comparisons with our signal models have been obtained. Having established that only modest signal-to-noise ratios are obtainable unless very high gain antennas are used, we expect that single measurements might not provide an estimate of sea state parameters as accurate as that obtainable with traditional remote sensing instruments, if costs are limited. Therefore spatial and temporal averaging of many measurements is required. Since the receiving satellite tracks do not repeat, measurements in a given area will be used to refine the local solution as a function of time, and will define the spatial resolution. In order to make this instrument viable in space, we need to detect and process many scattered signals. In 1998 we made a first step at understanding what configurations of antenna gain and orientations capture the largest number of viable signals to be used in the subsequent spatial/temporal averaging process. Additional information is contained in the original.

Lichten, Stephen M.↗

NASA NTRS1989

Gravitational model improvement at the Goddard Space Flight Center

Major new computations of terrestrial gravitational field models were performed by the Geodynamics Branch of Goddard Space Flight Center (GSFC). This development has incorporated the present state of the art results in satellite geodesy and have relied upon a more consistent set of reference constants than was heretofore utilized in GSFC's GEM models. The solutions are complete in spherical harmonic coefficients out to degree 50 for the gravity field parameters. These models include adjustment for a subset of 66 ocean tidal coefficients for the long wavelength components of 12 major ocean tides. This tidal adjustment was made in the presence of 550 other fixed ocean tidal terms representing 32 major and minor ocean tides and the Wahr frequency dependent solid earth tidal model. In addition 5-day averaged values for Earth rotation and polar motion were derived for the time period of 1980 onward. Two types of models were computed. These are satellite only models relying exclusively on tracking data and combination models which have incorporated satellite altimetry and surface gravity data. The satellite observational data base consists of over 1100 orbital arcs of data on 31 satellites. A large percentage of these observations were provided by third generation laser stations (less than 5 cm). A calibration of the model accuracy of the GEM-T2 satellite only solution indicated that it was a significant improvement over previous models based solely upon tracking data. The rms geoid error for this field is 110 cm to degree and order 36. This is a major advancement over GEM-T1 whose errors were estimated to be 160 cm. An error propagation using the covariances of the GEM-T2 model for the TOPEX radial orbit component indicates that the rms radial errors are expected to be 12 cm. The combination solution, PGS-3337, is a preliminary effort leading to the development of GEM-T3. PGS-3337 has incorporated global sets of surface gravity data and the Seasat altimetry to produce a model complete to (50,50). A solution for the dynamic ocean topography to degree and order 10 was included as part of this adjustment.

Marsh, J. G.↗

NASA NTRSDate not supplied

Exploring Radiation, Water and Energy Closure in Key Deep Tropical Regions Over the Last 3 Decades

The closure of global energy and water budgets at various time scales depends heavily upon the on the processes of water and energy exchange in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance in key tropical regions is observed of water budget terms between observations and atmospheric reanalysis systems, implying that processes are not simulated well. For instance, periods of imbalance in the water budget may be at least partly attributed to precipitation estimates due to the transition between cloud and precipitation systems in which ice cloud affects microwave radar estimates (Lietmann-Niimi et al. 2024, submitted). A recent paper by Hsaio et al (2023, pre-published version) found that the transition process is characterized by wind shear and longwave cloud radiative feedback signatures as well. This is just one example how these processes also affect the radiative fluxes and thus the energy budget at the TOA. In this presentation, we assess the net radiative budget and components in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES EBAF/SYN1deg Ed4.2, ISSCP FH, etc.) and describing the temporal variability of energy balance terms across the deep tropics for the period from 1988 to near present in six key oceanic regions. These six key tropical oceanic regions are West Indian Ocean, the south Indian Ocean east, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean. In these regions we assess the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear to assess the uncertainties of energy closure in these regions. These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. These differences also affect closure terms and inferences on systematic uncertainties within the observational estimates. We conclude with recommendations for continuing work

Radiation budget↗