Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational data assimilation”

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.

At least 307 records · Page 17

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

Diagnostic Studies with GLA Fields

Assessments of the NASA Goddard Earth Observing System-1 Data Assimilation System(GEOS-1 DAS), regarding heating rates, energetics, and angular momentum quantities were made. These diagnostics can be viewed as measures of climate variability. Comparisons with the NOAA/NCEP reanalysis system of momentum and energetics diagnostics are included. Water vapor and angular momentum are diagnosed in many models, including those of NASA, as part of the Atmospheric Model Intercomparison Project. 'Me GEOS-I and NOAA/NCEP global atmospheric angular momentum values are coherent on time scales down to about three days. Furthermore, they agree with the series of Earth angular momentum, as measured by tiny fluctuations in the rotation rate of the Earth, as variations in the length of day. The torques that effect such changes in atmospheric and Earth momentum are dominated by the influence of particular mountain systems, including the Rockies, Himalayas, and Andes, upon mountain torques on time scales shorter than about two weeks. Other project areas included collaboration with Goddard Space Flight Center to examine the impact of mountainous areas and the treatments of parameterizations on diagnoses of the atmosphere. Relevant preprints are included herein.

Salstein, David A.↗

The Use of MERRA-2 Near Surface Meteorology to Understand the Behavior of Planetary Boundary Layer Heights Derived from Wind Profiler Data over the US Great Plains

The atmospheric general circulation model (GCM) that underlies the MERRA-2 reanalysis includes a suite of physical parameterizations that describe the processes that occur in the planetary boundary layer (PBL). The data assimilation system assures that the atmospheric state variables used as input to these parameterizations are constrained to the best fit to all of the available observations. Many studies, however, have shown that the GCM-based estimates of MERRA-2 PBL heights are biased high, and so are not reliable for application related to constituent transport or the carbon cycle. A new 20-year record of PBL heights was derived from Wind Profiler (WP) backscatter data measured at a wide network of stations throughout the US Great Plains and has been validated against independent estimates. The behavior of these PBL heights shows geographical and temporal variations that are difficult to attribute to particular physical processes without additional information that are not part of the observational record. In the present study, we use information on physical processes from MERRA-2 to understand the behavior of the WP derived PBL heights. The behavior of the annual cycle of both MERRA-2 and WP PBL heights shows three classes of behavior: (i) canonical, where the annual cycle follows the annual cycle of the sun, (ii) delayed, where the PBL height reaches its annual maximum after the annual maximum of the solar insolation, and (iii) double maxima, where the PBL height begins to rise with the solar insolation but falls sometimes during the summer and then rises again. Although the magnitude of these types of variations is described by the WP PBL record, the explanation for these behaviors and the relationship to local precipitation, temperature, hydrology and sensible and latent heat fluxes is articulated using information from MERRA-2.

MERRA↗

SST Improvements in the NASA GMAO Weather Analysis and Prediction System

It is well known that the sea surface temperature (SST) exhibits a diurnal cycle. Since January 2017, the NASA GMAO's Weather Analysis & Prediction System (https://gmao.gsfc.nasa.gov/weather_prediction/) includes a model for this diurnal cycle, a cool-skin layer and assimilates for skin SST. Data assimilation for the skin SST is carried out along with the 3-D state of the atmosphere, primarily using infrared radiance measurements, that are sensitive to the cool-skin and diurnal warming. In the context of Numerical Weather Prediction (NWP), this implies that the skin SST from our system is tightly coupled to the overlaying atmosphere and constrained by all the available atmospheric and surface-sensitive observations (in situ and satellite brightness temperatures). Comparisons of the SST diurnal cycle with other SST retrievals (SEVIRI (Spinning Enhanced Visible & InfraRed Imager) and AMSR-2 (Advanced Microwave Scanning Radiometer-2)) has shown mostly similar diurnal variation, however, the former decays at a faster rate. This poster presents our SST plus atmosphere assimilation methodology, and recent improvements to the skin SST model to address the known issues. Global, hourly estimates of the cool-skin and diurnal warm layer fields would be of interest to the SST retrieval community.

Akella, Santha↗

Global analysis and forecasts of carbon monoxide on Mars

Spatial and temporal variations in the Martian carbon monoxide (CO) cycle have been investigated through combining Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) retrievals of carbon monoxide mixing ratio and Mars Climate Sounder (MCS) temperature profiles with a Martian global circulation model (GCM) to produce the first global reanalysis of the Martian CO cycle. The reanalysis reduces the root mean square error between the forecast and CRISM CO retrievals by a factor of 2–4, dependent on the time of year. Forecasts initiated from the reanalysis show an improved match to standalone CRISM CO retrievals from northern winter solstice to northern summer solstice, indicating the benefit of data assimilation in simulating the CO abundance. The northern summer solstice CO minimum between 10°S – 50°S in the CRISM CO retrievals is found to be caused by a suppression of CO-enriched air leaking from the Hellas and Argyre basins, and is also strongly influenced by the amount of carbon dioxide condensing at the time. The GCM is the first reported to simulate the local CO minimum however it is weaker in strength, as CO-enriched air is also released northward in the area of Argyre basin in the GCM because of a shifted boundary between the southern polar vortex and mid-latitudes. The reanalysis, as a result of the assimilation of MCS temperature profiles, indicates that the polar vortex boundary is northward of Argyre basin and hence no northerly transport of CO-enriched air should be present over this location. Differences in local CO abundance between the reanalysis and GCM are also evident in spatio-temporal regions where no nearby CRISM CO retrievals are available as a result of alterations in local circulation and the sublimation/condensation of carbon dioxide through the combined assimilation of CRISM CO retrievals alongside MCS temperature profiles, illustrating how constraints can be imposed indirectly on the CO cycle through the powerful technique of data assimilation.

James A. Holmes↗

Regional-Scale Modeling at NASA Goddard Space Flight Center

Over the past decade, the Goddard Mesoscale Modeling and Dynamics Group has used a popular regional scale model, MM5, to study precipitation processes. Our group is making contributions to the MM5 by incorporating the following physical and numerical packages: improved Goddard cloud processes, a land processes model (Parameterization for Land-Atmosphere-Cloud Exchange - PLACE), efficient but sophisticated radiative processes, conservation of hydrometeor mass (water budget), four-dimensional data assimilation for rainfall, and better computational methods for trace gas transport. At NASA Goddard, the MM5 has been used to study: (1) the impact of initial conditions, assimilation of satellite-derived rainfall, and cumulus parameterizations on rapidly intensifying oceanic cyclones, hurricanes and typhoons, (2) the dynamic and thermodynamic processes associated with the development of narrow cold frontal rainbands, (3) regional climate and water cycles, (4) the impact of vertical transport by clouds and lightning on trace gas distributiodproduction associated with South and North American mesoscale convective systems, (5) the development of a westerly wind burst (WWB) that occurred during the TOGA COARE and the diurnal variation of precipitation in the tropics, (6) a Florida sea breeze convective event and a Mid-US flood event using a sophisticated land surface model, (7) the influence of soil heterogeneity on land surface energy balance in the southwest GCIP region, (8) explicit simulations (with 1.33 to 4 km horizontal resolution) of hurricanes Bob (1991) and Bonnie (1998), (9) a heavy precipitation event over Taiwan, and (10) to make real time forecasts for a major NASA field program. In this paper, the modifications and simulated cases will be described and discussed.

Tao, W.-K.↗

Evaluating Nighttime CALIOP 0.532 micron Aerosol Optical Depth and Extinction Coefficient Retrievals

NASA Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) Version 3.01 5-km nighttime 0.532 micron aerosol optical depth (AOD) datasets from 2007 are screened, averaged and evaluated at 1 deg X 1 deg resolution versus corresponding/co-incident 0.550 micron AOD derived using the US Navy Aerosol Analysis and Prediction System (NAAPS), featuring two-dimensional variational assimilation of quality-assured NASA Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging Spectroradiometer (MISR) AOD. In the absence of sunlight, since passive radiometric AOD retrievals rely overwhelmingly on scattered radiances, the model represents one of the few practical global estimates available from which to attempt such a validation. Daytime comparisons, though, provide useful context. Regional-mean CALIOP vertical profiles of night/day 0.532 micron extinction coefficient are compared with 0.523/0.532 micron ground-based lidar measurements to investigate representativeness and diurnal variability. In this analysis, mean nighttime CALIOP AOD are mostly lower than daytime (0.121 vs. 0.126 for all aggregated data points, and 0.099 vs. 0.102 when averaged globally per normalised 1 deg. X 1 deg. bin), though the relationship is reversed over land and coastal regions when the data are averaged per normalised bin (0.134/0.108 vs. 0140/0.112, respectively). Offsets assessed within single bins alone approach +/- 20 %. CALIOP AOD, both day and night, are higher than NAAPS over land (0.137 vs. 0.124) and equal over water (0.082 vs. 0.083) when averaged globally per normalised bin. However, for all data points inclusive, NAAPS exceeds CALIOP over land, coast and ocean, both day and night. Again, differences assessed within single bins approach 50% in extreme cases. Correlation between CALIOP and NAAPS AOD is comparable during both day and night. Higher correlation is found nearest the equator, both as a function of sample size and relative signal magnitudes inherent at these latitudes. Root mean square deviation between CALIOP and NAAPS varies between 0.1 and 0.3 globally during both day/night. Averaging of CALIOP along-track AOD data points within a single NAAPS grid bin improves correlation and RMSD, though day/night and land/ocean biases persist and are believed systematic. Vertical profiles of extinction coefficient derived in the Caribbean compare well with ground-based lidar observations, though potentially anomalous selection of a priori lidar ratios for CALIOP retrievals is likely inducing some discrepancies. Mean effective aerosol layer top heights are stable between day and night, indicating consistent layer-identification diurnally, which is noteworthy considering the potential limiting effects of ambient solar noise during day.

aerosol optical depth↗

Satellite-Based Spatiotemporal Trends in PM2.5 Concentrations: China 2004-2013

Three decades of rapid economic development is causing severe and widespread PM2.5(particulate matter (is) less than 2.5 ) pollution in China. However, research on the health impacts of PM2.5 exposure has been hindered by limited historical PM2.5 concentration data. We estimated ambient PM2.5 concentrations from 2004 to 2013 in China at 0.1 deg resolution using the most recent satellite data and evaluated model performance with available ground observations. We developed a two-stage spatial statistical model using the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 aerosol optical depth (AOD) and assimilated meteorology, land use data, and PM2.5 concentrations from China's recently established ground monitoring network. An inverse variance weighting (IVW) approach was developed to combine MODIS Dark Target and Deep Blue AOD to optimize data coverage. We evaluated model predicted PM2.5 concentrations from 2004 to early 2014 using ground observations. The overall model cross-validation R(sup 2) and relative prediction error were 0.79 and 35.6%, respectively. Validation beyond the model year (2013) indicated that it accurately predicted PM(sub 2.5) concentrations with little bias at the monthly (R(sup 2) = 0.73), regression slope = 0.91) and seasonal (R(sup 2) = 0.79), regression slope = 0.92) levels. Seasonal variations revealed that winter was the most polluted season and that summer was the cleanest season. Analysis of predicted PM2.5 levels showed a mean annual increase of 1.97 micro-g/cu cm between 2004 and 2007 and a decrease of 0.46 micro-g/cu cm between 2008 and 2013. Our satellite-driven model can provide reliable historical PM2.5 estimates in China at a resolution comparable to those used in epidemiologic studies on the health effects of long-term PM2.5 exposure in North America. This data source can potentially advance research on PM2.5 health effects in China.

satellite remote sensing↗

Global 3-D Ionospheric Electron Density Reanalysis Based on Multisource Data Assimilation

We report preliminary results of a global 3-D ionospheric electron density reanalysis demonstration study during 2002-2011 based on multisource data assimilation. The monthly global ionospheric electron density reanalysis has been done by assimilating the quiet days ionospheric data into a data assimilation model constructed using the International Reference Ionosphere (IRI) 2007 model and a Kalman filter technique. These data include global navigation satellite system (GNSS) observations of ionospheric total electron content (TEC) from ground-based stations, ionospheric radio occultations by CHAMP, GRACE, COSMIC, SAC-C, Metop-A, and the TerraSAR-X satellites, and Jason-1 and 2 altimeter TEC measurements. The output of the reanalysis are 3-D gridded ionospheric electron densities with temporal and spatial resolutions of 1 h in universal time, 5deg in latitude, 10deg in longitude, and approx.30 km in altitude. The climatological features of the reanalysis results, such as solar activity dependence, seasonal variations, and the global morphology of the ionosphere, agree well with those in the empirical models and observations. The global electron content derived from the international GNSS service global ionospheric maps, the observed electron density profiles from the Poker Flat Incoherent Scatter Radar during 2007-2010, and foF2 observed by the global ionosonde network during 2002-2011 are used to validate the reanalysis method. All comparisons show that the reanalysis have smaller deviations and biases than the IRI-2007 predictions. Especially after April 2006 when the six COSMIC satellites were launched, the reanalysis shows significant improvement over the IRI predictions. The obvious overestimation of the low-latitude ionospheric F region densities by the IRI model during the 23/24 solar minimum is corrected well by the reanalysis. The potential application and improvements of the reanalysis are also discussed.

Xinan Yue↗

Comparing Two Approaches for Assessing Observation Impact

Langland and Baker introduced an approach to assess the impact of observations on the forecasts. In that approach, a state-space aspect of the forecast is defined and a procedure is derived ultimately relating changes in the aspect with changes in the observing system. Some features of the state-space approach are to be noted: the typical choice of forecast aspect is rather subjective and leads to incomplete assessment of the observing system, it requires availability of a verification state that is in practice correlated with the forecast, and it involves the adjoint operator of the entire data assimilation system and is thus constrained by the validity of this operator. This article revisits the topic of observation impacts from the perspective of estimation theory. An observation-space metric is used to allow inferring observation impact on the forecasts without the limitations just mentioned. Using differences of observation-minus-forecast residuals obtained from consecutive forecasts leads to the following advantages: (i) it suggests a rather natural choice of forecast aspect that directly links to the data assimilation procedure, (ii) it avoids introducing undesirable correlations in the forecast aspect since verification is done against the observations, and (iii) it does not involve linearization and use of adjoints. The observation-space approach has the additional advantage of being nearly cost free and very simple to implement. In its simplest form it reduces to evaluating the statistics of observationminus- background and observation-minus-analysis residuals with traditional methods. Illustrations comparing the approaches are given using the NASA Goddard Earth Observing System.

Variational Analysis↗

Real Time Land-Surface Hydrologic Modeling Over Continental US

The land surface component of the hydrological cycle is fundamental to the overall functioning of the atmospheric and climate processes. Spatially and temporally variable rainfall and available energy, combined with land surface heterogeneity cause complex variations in all processes related to surface hydrology. The characterization of the spatial and temporal variability of water and energy cycles are critical to improve our understanding of land surface-atmosphere interaction and the impact of land surface processes on climate extremes. Because the accurate knowledge of these processes and their variability is important for climate predictions, most Numerical Weather Prediction (NWP) centers have incorporated land surface schemes in their models. However, errors in the NWP forcing accumulate in the surface and energy stores, leading to incorrect surface water and energy partitioning and related processes. This has motivated the NWP to impose ad hoc corrections to the land surface states to prevent this drift. A proposed methodology is to develop Land Data Assimilation schemes (LDAS), which are uncoupled models forced with observations, and not affected by NWP forcing biases. The proposed research is being implemented as a real time operation using an existing Surface Vegetation Atmosphere Transfer Scheme (SVATS) model at a 40 km degree resolution across the United States to evaluate these critical science questions. The model will be forced with real time output from numerical prediction models, satellite data, and radar precipitation measurements. Model parameters will be derived from the existing GIS vegetation and soil coverages. The model results will be aggregated to various scales to assess water and energy balances and these will be validated with various in-situ observations.

Houser, Paul R.↗

An Assessment of Actual and Potential Building Climate Zone Change and Variability From the Last 30 Years Through 2100 Using NASA's MERRA and CMIP5 Simulations

Background: In the US, residential and commercial building infrastructure combined consumes about 40% of total energy usage and emits about 39% of total CO2 emission (DOE/EIA "Annual Energy Outlook 2013"). Building codes, as used by local and state enforcement entities are typically tied to the dominant climate within an enforcement jurisdiction classified according to various climate zones. These climate zones are based upon a 30-year average of local surface observations and are developed by DOE and ASHRAE. Establishing the current variability and potential changes to future building climate zones is very important for increasing the energy efficiency of buildings and reducing energy costs and emissions in the future. Objectives: This paper demonstrates the usefulness of using NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA) atmospheric data assimilation to derive the DOE/ASHRAE building climate zone maps and then using MERRA to define the last 30 years of variability in climate zones for the Continental US. An atmospheric assimilation is a global atmospheric model optimized to satellite, atmospheric and surface in situ measurements. Using MERRA as a baseline, we then evaluate the latest Climate Model Inter-comparison Project (CMIP) climate model Version 5 runs to assess potential variability in future climate zones under various assumptions. Methods: We derive DOE/ASHRAE building climate zones using surface and temperature data products from MERRA. We assess these zones using the uncertainties derived by comparison to surface measurements. Using statistical tests, we evaluate variability of the climate zones in time and assess areas in the continental US for statistically significant trends by region. CMIP 5 produced a data base of over two dozen detailed climate model runs under various greenhouse gas forcing assumptions. We evaluate the variation in building climate zones for 3 different decades using an ensemble and quartile statistics to provide an assessment of potential building climate zone changes relative to the uncertainties demonstrated using MERRA. Findings and Conclusions: These results show that there is a statistically significant increase in the area covered by warmer climate zones and a tendency for a reduction of area in colder climate zones in some limited regions. The CMIP analysis shows that models vary from relatively little building climate zone change for the least sensitive and conservation assumptions to a warming of at most 3 zones for certain areas, particularly the north central US by the end of the 21st century.

Stackhouse, Paul W., Jr.↗

High-Frequency Planetary Waves in the Polar Middle Atmosphere as seen in a data Assimilation System

This study examines the winter southern hemisphere vortex of 1998 using four times daily output from a data assimilation system to focus on the polar 2-day, wave number 2 component of the 4-day wave. The data assimilation system products are from a test version of the finite volume data assimilation system (fvDAS) being developed at Goddard Space Flight Center (GSFC) and include an ozone assimilation system. Results show that the polar 2-day wave dominates during July 1998 at 70 degrees. The period of the quasi 2-day wave is somewhat shorter than 2 days (about 1.7 days) during July 1998 with an average perturbation temperature amplitude for the month of over 2.5 K. The 2-day wave propagates more slowly than the zonal mean zonal wind, consistent with Rossby wave theory, and has EP flux divergence regions associated with regions of negative horizontal potential vorticity gradients, as expected from linear instability theory. Results for the assimilation-produced ozone mixing ratio show that the 2-day wave represents a major source of ozone variation in this region. The ozone wave in the assimilation system is in good agreement with the wave seen in the POAM (Polar Ozone and Aerosol Measurement) ozone observations for the same time period. Some differences with linear instability theory are noted as well as spectral peaks in the ozone field, not seen in the temperature field, that may be a consequence of advection.

Coy, L.↗

How Well are Recent Climate Variability Signals Resolved by Satellite Radiative Flux Estimates?

One notable aspect of Earth s climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. Understanding the character of radiative flux estimates and relating them to variations in other energy fluxes and climate state variables is key to improving our understanding of climate. In this work we will evaluate several recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project FD radiative flux profiles are available from rnid-1983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and HlRS operational soundings profiles. Full and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NASA/Goddard Earth Observing System assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. Several estimates of downward LW flux at the surface inferred from microwave data are also examined. Since these products have been evaluated with Baseline Surface Radiation Network data over land we focus over ocean regions and use the DOE/NOAA/NASA Shipboard Ocean Atmospheric Radiation (SOAR) surface flux measurements to characterize performance of these data sets under both clear and cloudy conditions. Some aspects of performance are stratified according to SST and vertical motion regimes. Comparisons to the TRMM/CERES SRB data in 1998 are also interpreted. These radiative fluxes are then analyzed to determine how surface (and TOA) radiative exchanges respond to interannual signals of ENS0 warm and cold events. Our analysis includes regional changes as well as integrated signals over land, ocean and various latitude bands. Changes in water vapor and cloud forcing signatures are prominent on interannual time scales. Prominent signals are also found in the SW fluxes for the Pinatubo volcanic event. These systematic changes in fluxes are related to changes in large-scale circulations and energy transport in the atmosphere and ocean. Some estimates of signal-to-noise and reliability are discussed to place our results in context.

Robertson, Franklin R.↗

Characterization of Turbulent Latent and Sensible Heat Flux Exchange Between the Atmosphere and Ocean in MERRA

Turbulent fluxes of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth's energy and water balance. Characterizing both the spatiotemporal variability and the fidelity of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. This study examines the veracity of the recently completed NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) product with respect to its representation of the surface turbulent heat fluxes. A validation of MERRA turbulent heat fluxes and near-surface bulk variables at local, high-resolution space and time scales is achieved by making comparisons to a large suite of direct observations. Both in situ and satellite-observed gridded surface heat flux estimates are employed to investigate the spatial and temporal variability of the surface fluxes with respect to their annual mean climatologies, their seasonal covariability of near-surface bulk parameters, and their representation of extremes. The impact of data assimilation on the near-surface parameters is assessed through evaluation of incremental analysis update tendencies produced by the assimilation procedure. It is found that MERRA turbulent surface heat fluxes are relatively accurate for typical conditions but have systematically weak vertical gradients in moisture and temperature and have a weaker covariability between the near-surface gradients and wind speed than found in observations. This results in an underestimate of the surface latent and sensible heat fluxes over the western boundary current and storm track regions. The assimilation of observations mostly acts to bring MERRA closer to observational products by increasing moisture and temperature near the surface and decreasing the near-surface wind speeds. The major patterns of spatial and temporal variability of the turbulent heat fluxes produced by MERRA compare favorably to observationally based estimates. However, MERRA is distinct in terms of amplitude. These results suggest that MERRA is likely to be a valuable resource for a number of research applications though, as with all turbulent flux estimates, systematic issues should be taken into account

Robert, J. Brent↗

Characterization of Turbulent Latent and Sensible Heat Flux Exchange Between the Atmosphere and Ocean in MERRA

Turbulent fluxes of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth s energy and water balance. Characterizing both the spatiotemporal variability and the fidelity of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. This study examines the veracity of the recently completed NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) product with respect to its representation of the surface turbulent heat fluxes. A validation of MERRA turbulent heat fluxes and near-surface bulk variables at local, high-resolution space and time scales is achieved by making comparisons to a large suite of direct observations. Both in situ and satellite-observed gridded surface heat flux estimates are employed to investigate the spatial and temporal variability of the surface fluxes with respect to their annual mean climatologies, their seasonal covariability of near-surface bulk parameters, and their representation of extremes. The impact of data assimilation on the near-surface parameters is assessed through evaluation of incremental analysis update tendencies produced by the assimilation procedure. It is found that MERRA turbulent surface heat fluxes are relatively accurate for typical conditions but have systematically weak vertical gradients in moisture and temperature and have a weaker covariability between the near-surface gradients and wind speed than found in observations. This results in an underestimate of the surface latent and sensible heat fluxes over the western boundary current and storm track regions. The assimilation of observations mostly acts to bring MERRA closer to observational products by increasing moisture and temperature near the surface and decreasing the near-surface wind speeds. The major patterns of spatial and temporal variability of the turbulent heat fluxes produced by MERRA compare favorably to observationally based estimates. However, MERRA is distinct in terms of amplitude. These results suggest that MERRA is likely to be a valuable resource for a number of research applications though, as with all turbulent flux estimates, systematic issues should be taken into account.

Roberts, J. Brent↗

Comparison of Land Skin Temperature from a Land Model, Remote Sensing, and In-situ Measurement

Land skin temperature (Ts) is an important parameter in the energy exchange between the land surface and atmosphere. Here hourly Ts from the Community Land Model Version 4.0, MODIS satellite observations, and in-situ observations in 2003 were compared. Compared with the in-situ observations over four semi-arid stations, both MODIS and modeled Ts show negative biases, but MODIS shows an overall better performance. Global distribution of differences between MODIS and modeled Ts shows diurnal, seasonal, and spatial variations. Over sparsely vegetated areas, the model Ts is generally lower than the MODIS observed Ts during the daytime, while the situation is opposite at nighttime. The revision of roughness length for heat and the constraint of minimum friction velocity from Zeng et al. [2012] bring the modeled Ts closer to MODIS during the day, and have little effect on Ts at night. Five factors contributing to the Ts differences between the model and MODIS are identified, including the difficulty in properly accounting for cloud cover information at the appropriate temporal and spatial resolutions, and uncertainties in surface energy balance computation, atmospheric forcing data, surface emissivity, and MODIS Ts data. These findings have implications for the cross-evaluation of modeled and remotely sensed Ts, as well as the data assimilation of Ts observations into Earth system models.

Modeling↗

The Ionospheric Three-Dimensional Electron Density Variations Induced by the 21 August 2017 Total Solar Eclipse by Using Global Ionospheric Specification

Global Ionospheric Specification (GIS) is based on the Gauss–Markov Kalman filter to assimilate the slant total electron content (TEC) observed from ground-based GPS receivers and space-based radio occultation instrumentations in order to reconstruct three-dimensional (3D) ionospheric electron density structure, and it can remotely sense and monitor the weather condition in space. In this study, five minutes of high temporal resolution GIS is implemented in order to reconstruct the 3D electron density structure on the 21 August 2017 total solar eclipse and analyze the variations induced by the moon’s shadow. To obtain more information of the ionosphere, from the extend 2200 GPS stations on the continental United States, are added for assimilation. The results show the ionosphere peak height (hmF2) uplift was 30–50 km altitude in latitude 25–40°N, and that the electron density depletion at higher altitudes (400 km) has a more noticeable time delay than at low altitudes (200 km), especially in low-latitude regions.

Ionosphere↗