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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

NASA Team 2 Sea Ice Concentration Algorithm Retrieval Uncertainty

Satellite microwave radiometers are widely used to estimate sea ice cover properties (concentration, extent, and area) through the use of sea ice concentration (IC) algorithms. Rare are the algorithms providing associated IC uncertainty estimates. Algorithm uncertainty estimates are needed to assess accurately global and regional trends in IC (and thus extent and area), and to improve sea ice predictions on seasonal to interannual timescales using data assimilation approaches. This paper presents a method to provide relative IC uncertainty estimates using the enhanced NASA Team (NT2) IC algorithm. The proposed approach takes advantage of the NT2 calculations and solely relies on the brightness temperatures (TBs) used as input. NT2 IC and its associated relative uncertainty are obtained for both the Northern and Southern Hemispheres using the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) TB. NT2 IC relative uncertainties estimated on a footprint-by-footprint swath-by-swath basis were averaged daily over each 12.5-km grid cell of the polar stereographic grid. For both hemispheres and throughout the year, the NT2 relative uncertainty is less than 5%. In the Southern Hemisphere, it is low in the interior ice pack, and it increases in the marginal ice zone up to 5%. In the Northern Hemisphere, areas with high uncertainties are also found in the high IC area of the Central Arctic. Retrieval uncertainties are greater in areas corresponding to NT2 ice types associated with deep snow and new ice. Seasonal variations in uncertainty show larger values in summer as a result of melt conditions and greater atmospheric contributions. Our analysis also includes an evaluation of the NT2 algorithm sensitivity to AMSR-E sensor noise. There is a 60% probability that the IC does not change (to within the computed retrieval precision of 1%) due to sensor noise, and the cumulated probability shows that there is a 90% chance that the IC varies by less than +/-3%. We also examined the daily IC variability, which is dominated by sea ice drift and ice formation/melt. Daily IC variability is the highest, year round, in the MIZ (often up to 20%, locally 30%). The temporal and spatial distributions of the retrieval uncertainties and the daily IC variability is expected to be useful for algorithm intercomparisons, climate trend assessments, and possibly IC assimilation in models.

Sea Ice↗
NASA NTRS2000

Global Geophysical Fluids Center of IERS

The Global Geophysical Fluids Center (GGFC) and its seven Special Bureaus (SB, for Atmosphere, Oceans, Tides, Hydrology, Mantle, Core and Gravity/Geocenter) were establishes by the International Earth Rotation Service in 1998, to support global geodynamic research. Mass transports in the geophysical fluids of the Earth system will cause observable geodynamic effects on a broad time scale.These include (1) variations in the solid Earth's rotation (in length-of-day and polar motion/nutation) via the conservation of angular momentum and effected by torques at the fluid-solid Earth interface; (2) changes in the global gravitational field according to Newton's gravitational law; and (3) motion in the center of mass of the solid Earth relative to that of the whole Earth ("geocenter") via the conservation of linear momentum. These minute signals have become observable by space geodetic techniques, primarily VLBI, SLR, GPS, and DORIS, and new exciting data will be available by space gravity, altimetry, SAR, and magnetic missions. In this sense the precise space geodetic techniques have become effective means of remote sensing of global mass transports. The GGFC and its SBs have the responsibility of supporting, facilitating, and providing services to the worldwide research community in the related research areas. We compute, analyze, compare, archive, and disseminate the time series of the angular momenta and the related torques, gravitational coefficients, and geocenter shift for all geophysical fluids, based on global observational data, and/or products from state-of-the-art models some of which assimilate such data. The computed quantities, algorithm and data formats are standardized. This paper reviews our activities, reports the status, and looks forward into the future.

Chao, Benjamin F.↗
NASA NTRS2007

Salinity Remote Sensing and the Study of the Global Water Cycle

The SMOS and AquariusISAC-D satellite missions will begin a new era to map the global sea surface salinity (SSS) field and its variability from space within the next twothree years. They will provide critical data needed to study the interactions between the ocean circulation, global water cycle and climate. Key scientific issues to address are (1) mapping large expanses of the ocean where conventional SSS data do not yet exist, (2) understanding the seasonal and interannual SSS variations and the link to precipitation, evaporation and sea-ice patterns, (3) links between SSS and variations in the oceanic overturning circulation, (4) air-sea coupling processes in the tropics that influence El Nino, and (4) closing the marine freshwater budget. There is a growing body of oceanographic evidence in the form of salinity trends that portend significant changes in the hydrologic cycle. Over the past several decades, highlatitude oceans have become fresher while the subtropical oceans have become saltier. This change is slowly spreading into the subsurface ocean layers and may be affecting the strength of the ocean's therrnohaline overturning circulation. Salinity is directly linked to the ocean dynamics through the density distribution, and provides an important signature of the global water cycle. The distribution and variation of oceanic salinity is therefore attracting increasing scientific attention due to the relationship to the global water cycle and its influence on circulation, mixing, and climate processes. The oceans dominate the water cycle by providing 86% of global surface evaporation (E) and receiving 78% of global precipitation (P). Regional differences in E-P, land runoff, and the melting or freezing of ice affect the salinity of surface water. Direct observations of E-P over the ocean have large uncertainty, with discrepancies between the various state-of-the-art precipitation analyses of a factor of two or more in many regions. Quantifying the climatic influence of the oceanic water cycle requires more accurately resolving the net air-sea water flux. Measuring global SSS trends on seasonal to interannual timescales by satellite is fundamental to this problem because the SSS trends represent detectable time-integrated signals of the variable marine hydrological cycle. Satellite measurements, coupled with an array of in situ observations, will provide global synoptic SSS fields for the first time history. These data will provide a strong constraint on climate models and data assimilation efforts, which must properly represent the freshwater budget in terms of E-P, ocean advection and surface layer mixing in order to accurately simulate the true ocean state. The SSS fields will allow us to quantify the covariability between the SSS and the strong seasonal E-P cycle in the tropics and high latitudes. Field measurement campaigns to exploit satellite and in situ measurements to close the seasonal E-P cycle over an ocean region are being considered. Lastly the satellite systems will monitor and trace the large long-lived SSS anomalies from year to year that have the potential to influence El Nino and the large scale ocean circulation.

Lagerloef, G. S. E.↗
NASA NTRS1996

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗
NASA NTRS2010

Using Leaf Chlorophyll to Parameterize Light-Use-Efficiency Within a Thermal-Based Carbon, Water and Energy Exchange Model

Chlorophylls absorb photosynthetically active radiation and thus function as vital pigments for photosynthesis, which makes leaf chlorophyll content (C(sub ab) useful for monitoring vegetation productivity and an important indicator of the overall plant physiological condition. This study investigates the utility of integrating remotely sensed estimates of C(sub ab) into a thermal-based Two-Source Energy Balance (TSEB) model that estimates land-surface CO2 and energy fluxes using an analytical, light-use-efficiency (LUE) based model of canopy resistance. The LUE model component computes canopy-scale carbon assimilation and transpiration fluxes and incorporates LUE modifications from a nominal (species-dependent) value (LUE(sub n)) in response to short term variations in environmental conditions, However LUE(sub n) may need adjustment on a daily timescale to accommodate changes in plant phenology, physiological condition and nutrient status. Day to day variations in LUE(sub n) were assessed for a heterogeneous corn crop field in Maryland, U,S.A. through model calibration with eddy covariance CO2 flux tower observations. The optimized daily LUE(sub n) values were then compared to estimates of C(sub ab) integrated from gridded maps of chlorophyll content weighted over the tower flux source area. The time continuous maps of daily C(sub ab) over the study field were generated by focusing in-situ measurements with retrievals generated with an integrated radiative transfer modeling tool (accurate to within +/-10%) using at-sensor radiances in green, red and near-infrared wavelengths acquired with an aircraft imaging system. The resultant daily changes in C(sub ab) within the tower flux source area generally correlated well with corresponding changes in daily calibrated LUE(sub n) derived from the tower flux data, and hourly water, energy and carbon flux estimation accuracies from TSEB were significantly improved when using C(sub ab) for delineating spatio-temporal variations in LUE(sub n). The results demonstrate the synergy between thermal infrared and shortwave reflective wavebands in producing valuable remote sensing data for operational monitoring of carbon and water fluxes.

Houlborg, Rasmus↗
NASA NTRS2020

Stratospheric Impact on the Northern Hemisphere Winter and Spring Ozone Interannual Variability in the Troposphere

In this study we use ozone and stratospheric ozone tracer simulations from the high-resolution (0.5°x0.5°) Goddard Earth Observing System, Version 5 (GEOS-5) in a replay mode to study the impact of stratospheric ozone on tropospheric ozone interannual variability (IAV). We use these simulations in conjunction with ozone sonde measurements from 1990 to 2016 during the winter and spring seasons. The simulations include a stratospheric ozone tracer (StratO3) to aid 15in the evaluation of the impact of stratospheric ozone IAV on the IAV of tropospheric ozone at different altitudes and location. The model is in good agreement with the observed interannual variation of tropospheric ozone, except for the post-Pinatubo period(1992-1994)over the region of North America. Ozone sonde data show a negative ozone anomaly in 1992-1994 following the Pinatubo eruption, with recovery thereafter. The simulated anomaly is only half the magnitude of that observed. Our analysis suggests that the simulated Stratosphere-troposphere exchange (STE)flux deduced from the analysis 20might be too strong over the North American (50°N-70°N) region after the Mt. Pinatubo eruption in the early 1990s,maskingthe impact of lower stratospheric ozone concentration on tropospheric ozone. European ozone sonde measurements show a similar but weaker ozone depletion after the Mt. Pinatubo eruption, which is fully reproduced by the model. Analysis based on the stratospheric ozone tracer (StratO3) identifies differences in strength and vertical extent of stratospheric ozone impact on the tropospheric ozone interannual variation (IAV)between North America and Europe. Over NorthAmerica25stations, the StratO3IAV has a significant impact on tropospheric ozone from the upper to lower troposphere and explains about 60% and 66% of the simulated ozone IAV at 400 hPa, and ~11% and 34% at 700 hPa in winter and spring respectively. Over European stations, the influence is limited to the middle to upper troposphere, and becomes much smaller at 700 hPa. The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2)assimilated fields exhibit strong longitudinal variations over northern hemisphere (NH) mid-high latitudes, with lower tropopause height and 30lower geopotential height over North America than Europe. These variations associated with the relevant variations in the location of tropospheric jet flows are responsible for the longitudinal differences in the stratospheric ozone impact, with stronger effects over North America than over Europe

tropospheric ozone IAV, Model, ozonesonde, stratos↗
NASA NTRS1999

Systematic Differences between Satellite-Based Presipitation Climatologies over the Tropical Oceans

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. It is clear however that there still remain significant uncertainties with satellite precipitation retrievals which limit their usefulness for many purposes. Systematic differences i'A tropical precipitation estimates have been brought to light in comparison activities such as the GPCP Algorithm Intercomparison Project and more recent Wetnet Precipitation Intercomparison Project 3. These uncertainties are assuming more importance because of the demands for validation associated with global climate modeling and data assimilation methodologies. The objective of the present study is to determine the physical basis for systematic differences in spatial structure of tropical precipitation as portrayed by several different satellite-based data sets. The study is limited to oceanic regions only and deals primarily with aspects of spatial variability. We are specifically interested in why MSU channel 1 and GPI precipitation differences are so striking over the Eastern Pacific ITCZ and why they both differ from other microwave emission-based precipitation estimates from SSM/I and a scattering-based deep convective ice index from MSU channel 2. Our results to date have shown that MSU channel I precipitation estimates are biased high over the Eastern Pacific ITCZ because of two factors: (1) the hypersensitivity of this frequency to cloud water in contrast to falling rain drops, and (2) unaccounted for scattering effects by precipitation-size ice which depresses the signal of the liquid water emission. Likewise, cold cloud top climatologies such as the GPI show an excess (a deficit) in estimated rainfall over the E. Pacific ITCZ (Warm Pool region). We show that these algorithms need to account for regionally varying heights (or temperatures) at which tropical convection detrains to form cirrus shields. A second objective we pursue is to identify variations in the macroscale cloud physical and thermodynamic properties of precipitation regimes" and relate these differences to tropical dynamical mechanisms of tropical heat and moisture balance. Finally, we interpret the algorithm differences and their associations with tropical dynamics in terms of WCRP GPCP goals for constructing precipitation climatologies.

Robertson, Frankin R.↗
NASA NTRS2009

Five-Year (2004-2009)Observations of Upper Tropospheric Water Vapor and Cloud Ice from MLS and Comparisons with GEOS-5 Analyses

This paper gives an overview of August 2004 through July 2009 upper tropospheric (UT) water vapor (H2O) and ice water content (IWC) from the Aura Microwave Limb Sounder (MLS) and comparisons with outputs from the NASA Goddard Earth Observing System Version 5 (GEOS-5) data assimilation system. Both MLS and GEOS-5 show that high values of H2O and IWC at 215 to 147 hPa are associated with areas of deep convection. They exhibit good (within approximately 15%) agreement in IWC at these altitudes, but GEOS-5 H2O is approximately 50% (215 hPa) to approximately 30% (147 hPa) larger than MLS, possibility due to its higher temperatures at these altitudes. GOES-5 produces a weaker intertropical convergence zone than MLS, while a seasonally-migrating band of tropical deep convection is clearly evident in both the MLS and GEOS-5 UT H2O and IWC. MLS and GEOS-5 both show spatial anti-correlation between IWC and H2O at 100 hPa, where less H2O is associated with low temperatures in regions of tropical convection. At 100 hPa, GEOS-5 produces 50% less IWC and 15% less H2O in the tropics, and approximately 20% more H2O in the extra-tropics, than does MLS. Behavior of the 100 hPa H2O, which exhibits a quasi-biennial oscillation, appears consistent with it being controlled by temperature. The seasonal cycle in the vertical transport of tropical mean H2O from approximately 147 hPa to approximately 10 hPa appears much stronger in MLS than in GEOS-5. The UT IWC and H2O interannual variations, from both MLS and GEOS-5, show clear imprints of the El Nino-Southern Oscillation.

Jiang, Jonathan H.↗
NASA NTRS2010

A Comparison of Satellite Based, Modeled Derived Daily Solar Radiation Data with Observed Data for the Continental US

Many applications of simulation models and related decision support tools for agriculture and natural resource management require daily meteorological data as inputs. Availability and quality of such data, however, often constrain research and decision support activities that require use of these tools. Daily solar radiation (SRAD) data are especially problematic because the instruments require electronic integrators, accurate sensors are expensive, and calibration standards are seldom available. The Prediction Of Worldwide Energy Resources (NASA/POWER; power.larc.nasa.gov) project at the NASA Langley Research Center estimates daily solar radiation based on data that are derived from satellite observations of outgoing visible radiances and atmospheric parameters based upon satellite observations and assimilation models. The solar data are available for a global 1 degree x 1 degree coordinate grid. SRAD can also be estimated based on attenuation of extraterrestrial radiation (Q0) using daily temperature and rainfall data to estimate the optical thickness of the atmosphere. This study compares daily solar radiation data from NASA/POWER (SRADNP) with instrument readings from 295 stations (SRADOB), as well as with values that were estimated with the WGENR solar generator. WGENR was used both with daily temperature and precipitation records from the stations reporting solar data and records from the NOAA Cooperative Observer Program (COOP), thus providing two additional sources of solar data, SRADWG and SRADCO. Values of SRADNP for different grid cells consistently showed higher correlations (typically 0.85 to 0.95) with SRADOB data than did SRADWG or SRADCO for sites within the corresponding cells. Mean values of SRADOB, SRADWG and SRADNP for sites within a grid cell usually were within 1 MJm-2d-1 of each other, but NASA/POWER values averaged 1.1 MJm-2d-1 lower than SRADOB. The magnitude of this bias was greater at lower latitudes and during summer months and may be at least partially explained by assumptions in ambient aerosol properties. Overall, the NASA/POWER solar radiation data are a promising resource for regional modeling studies where realistic accounting of historic variation is required.

White, Jeffrey W.↗
NASA NTRSDate not supplied

Satellite-Assisted Particulate Matter (SAPM) for the Models, In situ, and Remote Sensing of Aerosols (MIRA) Working Group

Models, In situ, and Remote sensing of Aerosols (MIRA) is an international working group that provides a forum for collaborations amongst these three atmospheric aerosol communities. MIRA consists of a collection of interdisciplinary and independently funded Topic/Project groups with clear goals and generally characterized by requests for additional scientific data. The Satellite-Assisted Particulate Matter (SAPM) Topic group focuses on the study of fine particulate matter (PM2.5), as it is a major contributor to air pollution and negatively impacts human health. Our group aims to provide intercomparisons of various methods and techniques for estimating surface PM2.5 assisted by satellite remote sensors (passive and active), global aerosol models, and in situ aerosol measurements. The overall motivation of our work is to enhance PM2.5 coverage beyond in situ ground stations, which can be limited. In this presentation, we provide an overview of various PM2.5 estimation approaches by current SAPM team members, as each approach has its own strengths and limitations. One approach uses near-surface aerosol extinction retrievals from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and an assumed value of mass extinction efficiency to derive PM2.5 concentrations (Toth et al 2019). A study using this method found a promising agreement (R=0.60; slope=0.89) with U.S. Environmental Protection Agency (EPA) in situ PM2.5 measurements over the contiguous United States (CONUS) for a 12-year (2007-2018) period (Toth et al 2022). In a different approach combining spaceborne lidar and a model, the Cloud Aerosol Transport System (CATS) lidar and Goddard Earth Observing System (GEOS) model are blended in an ensemble variational assimilation scheme to retrieve aerosol extinction and convert speciated mass concentrations to total PM2.5. A recent study (Matus et al in prep.) found agreement within 2 μg/m3 over the CONUS on the 2016 annual mean when compared to EPA ground-based measurements. Another study explored trends in city aerosol optical depth (AOD) from a spaceborne passive remote sensor, Moderate Resolution Imaging Spectroradiometer (MODIS), and found that these trends agree well with trends in city surface PM2.5, which provides confidence in the use of city AOD for city PM2.5 trend studies (Vohra et al 2021). In a follow-on study, steep and significant trends (~2.5 to 7.8% a-1) in AOD are found for several South Asian cities (e.g., Bangalore and Hyderabad), suggesting rapid growth in PM2.5 (Vohra et al 2022). Also in this presentation, we report results on our current collective SAPM study of India, a country that exhibits high levels of PM2.5 pollution. We show the temporal (i.e., seasonal) and spatial variability of PM2.5 over India using the CALIOP and CATS+model approaches, comparisons between the estimated PM2.5 from the retrieval schemes, and validation with in situ PM2.5 ground-based observations from the Central Pollution Control Board in India. This region provides us with an excellent case study to test our PM2.5 retrievals in heavily polluted scenes and over a complex and varying topography. The SAPM Topic group actively seeks international participants/collaborators in our group, including those working with in situ aerosol measurements (e.g., ACTRIS). We are interested in aerosol datasets in order to either improve our PM2.5 estimates from various approaches (e.g., using in situ mass scattering/absorption coefficient and aerosol hygroscopic properties for various aerosol species) and/or validate the PM2.5 estimates (e.g., using in situ ground-based PM2.5 concentrations).

Travis D Toth↗
NASA NTRS2012

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗
NASA NTRS1996

Quality and Control of Water Vapor Winds

Water vapor imagery from the geostationary satellites such as GOES, Meteosat, and GMS provides synoptic views of dynamical events on a continual basis. Because the imagery represents a non-linear combination of mid- and upper-tropospheric thermodynamic parameters (three-dimensional variations in temperature and humidity), video loops of these image products provide enlightening views of regional flow fields, the movement of tropical and extratropical storm systems, the transfer of moisture between hemispheres and from the tropics to the mid- latitudes, and the dominance of high pressure systems over particular regions of the Earth. Despite the obvious larger scale features, the water vapor imagery contains significant image variability down to the single 8 km GOES pixel. These features can be quantitatively identified and tracked from one time to the next using various image processing techniques. Merrill et al. (1991), Hayden and Schmidt (1992), and Laurent (1993) have documented the operational procedures and capabilities of NOAA and ESOC to produce cloud and water vapor winds. These techniques employ standard correlation and template matching approaches to wind tracking and use qualitative and quantitative procedures to eliminate bad wind vectors from the wind data set. Techniques have also been developed to improve the quality of the operational winds though robust editing procedures (Hayden and Veldon 1991). These quality and control approaches have limitations, are often subjective, and constrain wind variability to be consistent with model derived wind fields. This paper describes research focused on the refinement of objective quality and control parameters for water vapor wind vector data sets. New quality and control measures are developed and employed to provide a more robust wind data set for climate analysis, data assimilation studies, as well as operational weather forecasting. The parameters are applicable to cloud-tracked winds as well with minor modifications. The improvement in winds through use of these new quality and control parameters is measured without the use of rawinsonde or modeled wind field data and compared with other approaches.

Jedlovec, Gary J.↗
NASA NTRS2002

A Comparison of the Age-Spectra from Data Assimilation Models

We use kinematic and diabatic back trajectory calculations, driven by winds from a general circulation model (GCM) and two different data assimilation systems (DAS), to compute the age spectrum at three latitudes in the lower stratosphere. The age-spectra are compared to chemical transport model (CTM) calculations, and the mean ages from all of these studies are compared to observations. The age spectra computed using the GCM winds show a reasonably well-isolated tropics in good agreement with observations; however, the age spectra determined from the DAS differ from the GCM spectra. For the diabatic trajectory calculations, the age spectrum is too broad as a result of too much exchange between the tropics and mid-latitudes. The age spectrum determined using the kinematic trajectory calculation is less broad and lacks an age offset; both of these features are due to excessive vertical dispersion of parcels. The tropical and mid-latitude mean age difference between the diabatically and kinematically determined age-spectra is about one year, the former being older. The CTM calculation of the age spectrum using the DAS winds shows the same dispersive characteristics of the kinematic trajectory calculation. These results suggest that the current DAS products will not give realistic trace gas distributions for long integrations; they also help explain why the mean ages determined in a number of previous DAS driven CTM's are too young compared with observations. Finally, we note trajectory-generated age spectra show significant age anomalies correlated with the seasonal cycles, and these anomalies can be linked to year-to-year variations in the tropical heating rate. These anomalies are suppressed in the CTM spectra suggesting that the CTM transport is too diffusive.

Schoeberl, Mark R.↗
NASA NTRS2022

Evaluation of High Mountain Asia-Land Data Assimilation System (version 1) from 2003 to 2016: 2. The impact of assimilating satellite-based snow cover and freeze/thaw observations into a land surface model

This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domainscale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospherichydrological process within a land surface model for use in HMA.

Yuan Xue↗
NASA NTRS2022

Evaluation of High Mountain Asia-Land Data Assimilation System (version 1) from 2003 to 2016: 2. The impact of assimilating satellite-based snow cover and freeze/thaw observations into a land surface model

This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domain-scale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospheric-hydrological process within a land surface model for use in HMA.

Yuan Xue↗
NASA NTRS2003

A Global Aerosol Model Forecast for the ACE-Asia Field Experiment

We present the results of aerosol forecast during the Aerosol Characterization Experiment (ACE-Asia) field experiment in spring 2001, using the Georgia Tech/Goddard Global Ozone Chemistry Aerosol Radiation and Transport (GOCART) model and the meteorological forecast fields from the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The aerosol model forecast provides direct information on aerosol optical thickness and concentrations, enabling effective flight planning, while feedbacks from measurements constantly evaluate the model, making successful model improvements. We verify the model forecast skill by comparing model predicted total aerosol extinction, dust, sulfate, and SO2 concentrations with those quantities measured by the C-130 aircraft during the ACE-Asia intensive operation period. The GEOS DAS meteorological forecast system shows excellent skills in predicting winds, relative humidity, and temperature for the ACE-Asia experiment area as well as for each individual flight, with skill scores usually above 0.7. The model is also skillful in forecast of pollution aerosols, with most scores above 0.5. The model correctly predicted the dust outbreak events and their trans-Pacific transport, but it constantly missed the high dust concentrations observed in the boundary layer. We attribute this missing dust source to the desertification regions in the Inner Mongolia Province in China, which have developed in recent years but were not included in the model during forecasting. After incorporating the desertification sources, the model is able to reproduce the observed high dust concentrations at low altitudes over the Yellow Sea. Two key elements for a successful aerosol model forecast are correct source locations that determine where the emissions take place, and realistic forecast winds and convection that determine where the aerosols are transported. We demonstrate that our global model can not only account for the large-scale intercontinental transport, but also produce the small-scale spatial and temporal variations that are adequate for aircraft measurements planning.

Chin, Mian↗
NASA NTRS1997

Internal Physical Features of a Land Surface Model Employing a Tangent Linear Model

The Earth's land surface, including its biomass, is an integral part of the Earth's weather and climate system. Land surface heterogeneity, such as the type and amount of vegetative covering., has a profound effect on local weather variability and therefore on regional variations of the global climate. Surface conditions affect local weather and climate through a number of mechanisms. First, they determine the re-distribution of the net radiative energy received at the surface, through the atmosphere, from the sun. A certain fraction of this energy increases the surface ground temperature, another warms the near-surface atmosphere, and the rest evaporates surface water, which in turn creates clouds and causes precipitation. Second, they determine how much rainfall and snowmelt can be stored in the soil and how much instead runs off into waterways. Finally, surface conditions influence the near-surface concentration and distribution of greenhouse gases such as carbon dioxide. The processes through which these mechanisms interact with the atmosphere can be modeled mathematically, to within some degree of uncertainty, on the basis of underlying physical principles. Such a land surface model provides predictive capability for surface variables including ground temperature, surface humidity, and soil moisture and temperature. This information is important for agriculture and industry, as well as for addressing fundamental scientific questions concerning global and local climate change. In this study we apply a methodology known as tangent linear modeling to help us understand more deeply, the behavior of the Mosaic land surface model, a model that has been developed over the past several years at NASA/GSFC. This methodology allows us to examine, directly and quantitatively, the dependence of prediction errors in land surface variables upon different vegetation conditions. The work also highlights the importance of accurate soil moisture information. Although surface variables are predicted imperfectly due to inherent uncertainties in the modeling process, our study suggests how satellite observations can be combined with the model, through land surface data assimilation, to improve their prediction.

Yang, Runhua↗
NASA NTRS2002

A Comparison of the Lower Stratospheric Age-Spectra Derived from a General Circulation Model and Two Data Assimilation Systems

We use kinematic and diabatic back trajectory calculations, driven by winds from a general circulation model (GCM) and two different data assimilation systems (DAS), to compute the age spectrum at three latitudes in the lower stratosphere. The age-spectra are compared to chemical transport model (CTM) calculations, and the mean ages from all of these studies are compared to observations. The age spectra computed using the GCM winds show a reasonably isolated tropics in good agreement with observations; however, the age spectra determined from the DAS differ from the GCM spectra. For the DAS diabatic trajectory calculations there is too much exchange between the tropics and mid-latitudes. The age spectrum is thus too broad and the tropical mean age is too old as a result of mixing older mid latitude air with tropical air. Likewise the mid latitude mean age is too young due to the in mixing of tropical air. The DAS kinematic trajectory calculations show excessive vertical dispersion of parcels in addition to excessive exchange between the tropics and mid latitudes. Because air is moved rapidly to the troposphere from the vertical dispersion, the age spectrum is shifted toward the young side. The excessive vertical and meridional dispersion compensate in the kinematic case giving a reasonable tropical mean age. The CTM calculation of the age spectrum using the DAS winds shows the same vertical and meridional dispersive characteristics of the kinematic trajectory calculation. These results suggest that the current DAS products will not give realistic trace gas distributions for long integrations; they also help explain why the extra tropical mean ages determined in a number of previous DAS driven CTM s are too young compared with observations. Finally, we note trajectory-generated age spectra . show significant age anomalies correlated with the seasonal cycles. These anomalies can be linked to year-to-year variations in the tropical heating rate. The anomalies are suppressed in the CTM spectra suggesting that the CTM transport scheme is too diffusive.

Schoeberl, Mark R.↗