SEARCH · Engineering Papers

Results for “Covariance Generation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

GeoNEX-SUBSETS: Cutouts of geostationary satellite data over long-term monitoring sites

Satellite remote sensing data are important tool for extrapolating the knowledge obtained at the Fluxnet towers. We introduce our GEO-NEX geostationary satellite subset products, which make it easy to compare between ground observation and Geo-NEX products. The MODIS subset has been frequently used for validating ecosystem models developed at flux towers and upscaling the observed flux data to regional scales. However, MODIS on the polar orbiting satellites can observe target regions only once a day, while the Fluxnet eddy-covariance data are compiled as sub-hourly datasets. As a results, summarizing the sub-hourly flux data into daily statistics is necessary for the comparison between MODIS and Flux data. Here, the new generation geostationary satellite sensors (GOES-16/17 ABI and Himawari-8/9 AHI) has the high-frequent observation feature (10 minutes) in addition to similar spectral band and spatial resolution with MODIS. The high frequent observation allows us to compare the flux diurnal cycle with geostationary satellite sensor data. Some studies have already shown the effective utilization of time series of geostationary satellite data for ecosystem modeling. We are producing NEX Level-1G products, which is the gridded Top-of-Atmosphere reflectance and brightness temperature data of geostationary satellite sensors. We cut out the NEX Level-1G data using the same file format with the MODIS subset except for the projection. The other data products (e.g., surface reflectance, land surface temperature, vegetation indices and climate data) will be also added upon their availability. We selected the ground observation sites from Fluxnet, Phenocam, and AERONET networks. The NEX subset data will be provided through NASA NEX data portal.

GeoNEX↗

NASA NTRS2002

SNC Meteorites and Martian Reservoirs

Jones first suggested that the inverse covariation of initial epsilon (Nd-143) and Sr-87/Sr-86 of the shergottites could be explained by interaction between mantle-derived magmas with another isotopic reservoir(s) (i.e., assimilation or contamination). In that model, magmas were generated in a source region that was isotopically similar to the Nakhla source and the second reservoir(s) was presumed to be crust. The text also permitted the second reservoir to be another type of mantle, but I can confirm that a second mantle reservoir was never seriously considered by that author. Other features of this model were that (i) it occurred at a particular time, 180 m.y. ago, and (ii) the interacting reservoirs had been separated at approximately 4.5 b.y. In a later paper Jones explored this mixing model more quantitatively and concluded that magmas from a Nakhla-like source region at 180 m.y. would fall on or near an isotopic Nd-Sr-Pb hyperplane defined by the shergottites. This criterion was a necessary prerequisite for the parent magma(s) of the shergottites to have initially been Nakhla-like isotopically. At this juncture, it is perhaps worthwhile to note that this mixing model was not presented to explain geochemical variations but as a justification for a 180 m.y. crystallization age for the shergottites and a 1.3 b.y. crystallization age for the nakhlites. In the mid-1980's crystallization ages estimated for Nakhla ranged from approximately 1.3 b.y to 4.5 b.y. Similarly, preferred crystallization ages for the shergottites ranged from 360 m.y., to 1.3 b.y., to 4.5 b.y. In all these models, the 180 m.y. event seen in the shergottites was deemed to be metamorphic. The fit between the Nakhla-like source region and the shergottite hyperplane was a validation both of the 1.3 b.y. igneous age of Nakhla and the 180 m.y. igneous age of the shergottites.

Jones, J. H.↗

NASA NTRS2013

The Aqua-Planet Experiment (APE): CONTROL SST Simulation

Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE.Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies.The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed.The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behavior and investigate convergence of the aqua-planet climate with increasing resolution.

AGCM↗

NASA NTRSDate not supplied

Characterizing Spatiotemporal Uncertainty in Interpolated Meteorological Data

Interpolated meteorological data invariably contain errors. These errors have structure in time and space, particularly autocorrelation, which can cause the effects of errors to compound when model outputs are aggregated temporally or spatially. One way to account for this uncertainty is with a probabilistic model from which samples can be drawn that are coherent with respect to underlying spatial and temporal covariance structure. This work describes a probabilistic method for spatial interpolation of point-wise meteorological time series. Observational data from weather stations are generally sparse in space and dense in time (but sometimes missing). The method works by projecting time series onto orthogonal basis vectors and spatially interpolating each resulting component independently. Under suitable assumptions, and data transformations to better satisfy those assumptions, Gaussian process regression provides a complete description of the joint predictive distribution over a Gaussian random field. Spatiotemporally coherent realizations are generated as the sum of conditional (spatial) simulations of each orthogonal (temporal) component. Data-derived and generic orthogonal bases are considered. In addition to spatial interpolation, imputation of missing observational data is examined. The method is applied using near-surface air temperature over the Western United States and validated by comparing theoretical versus actual coverage of predictive distributions and analyzing the degree to which spatial and temporal covariance structure is reproduced. Computational considerations, relating to conditional simulation of random fields, are also addressed.

Conor T Doherty↗

NASA NTRS1992

An assessment of gravity model improvements using TOPEX/Poseidon TDRSS observations

The contribution of TOPEX/Poseidon (T/P) TDRSS data to geopotential model recovery is assessed. Simulated TDRSS one-way and Bilateration Ranging Transponder System (BRTS) observations have been generated and orbitally reduced to form normal equations for geopotential parameters. These normals have been combined with those of the latest prelaunch T/P gravity model solution using data from over 30 satellites. A study of the resulting solution error covariance shows that TDRSS can make important contributions to geopotential recovery, especially for improving T/P specific effects like those arising from orbital resonance. It is argued that future effort is desirable both to establish TDRSS orbit determination limits in a reference frame compatible with that used for the precise laser/DORIS orbits, and the reduction of these TDRSS data for geopotential recovery.

Putney, B. H.↗

NASA NTRS2009

Efficient Method for Analyzing Hyperspectral Remote Sensing Data

The current and next generation satellite infrared sensors such as AIRS, IASI, CrIS, and CLARREO all have one thing in common: large number of spectral channels. In order to fully utilize the information content of these sensors, a large number of Radiative Transfer (RT) calculations through the inhomogeneous atmosphere are needed. It is also computationally intensive to invert atmospheric and surface properties using all the channels. Usually, only subsets of channels are used to perform physical inversions for atmospheric profiles. We will explore ways to speed up RT calculations and inversions and test the algorithms using data from the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument. We will describe a retrieval algorithm using a Principal Component-based Radiative Transfer Model (PCRTM) for generating atmospheric temperature/moisture/ozone profiles and surface properties. The retrieval algorithm performs both the radiative transfer calculations and inversions in the Principal Component (PC) domain. The inversion algorithm is based on a non-linear Levenberg- Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the ultraspectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels. We will also use data collected during the Joint Airborne IASI Validation Experiment (JAIVEx) field campaign to validate the algorithm and IASI retrievals.

Liu, Xu↗

NASA NTRS2017

Estimating Evaporative Fraction From Readily Obtainable Variables in Mangrove Forests of the Everglades, U.S.A.

A remote-sensing-based model to estimate evaporative fraction (EF) the ratio of latent heat (LE; energy equivalent of evapotranspiration -ET-) to total available energy from easily obtainable remotely-sensed and meteorological parameters is presented. This research specifically addresses the shortcomings of existing ET retrieval methods such as calibration requirements of extensive accurate in situ micro-meteorological and flux tower observations, or of a large set of coarse-resolution or model-derived input datasets. The trapezoid model is capable of generating spatially varying EF maps from standard products such as land surface temperature [T(sub s)] normalized difference vegetation index (NDVI)and daily maximum air temperature [T(sub a)]. The 2009 model results were validated at an eddy-covariance tower (Fluxnet ID: US-Skr) in the Everglades using T(sub s) and NDVI products from Landsat as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors. Results indicate that the model accuracy is within the range of instrument uncertainty, and is dependent on the spatial resolution and selection of end-members (i.e. wet/dry edge). The most accurate results were achieved with the T(sub s) from Landsat relative to the T(sub s) from the MODIS flown on the Terra and Aqua platforms due to the fine spatial resolution of Landsat (30 m). The bias, mean absolute percentage error and root mean square percentage error were as low as 2.9% (3.0%), 9.8% (13.3%), and 12.1% (16.1%) for Landsat-based (MODIS-based) EF estimates, respectively. Overall, this methodology shows promise for bridging the gap between temporally limited ET estimates at Landsat scales and more complex and difficult to constrain global ET remote-sensing models.

trapezoid↗

NASA NTRS2018

A Global Spatially Contiguous Solar-Induced Fluorescence (CSIF) Dataset Using Neural Networks

Satellite-retrieved solar-induced chlorophyll fluorescence (SIF) has shown great potential to monitor the photosynthetic activity of terrestrial ecosystems. However, several issues, including low spatial and temporal resolution of the gridded datasets and high uncertainty of the individual retrievals, limit the applications of SIF. In addition, inconsistency in measurement footprints also hinders the direct comparison between gross primary production (GPP) from eddy covariance (EC) flux towers and satellite-retrieved SIF. In this study, by training a neural network (NN) with surface reflectance from the MODerate-resolution Imaging Spectroradiometer (MODIS) and SIF from Orbiting Carbon Observatory-2 (OCO-2), we generated two global spatially contiguous SIF (CSIF) datasets at moderate spatiotemporal (0.05° 4-day) resolutions during the MODIS era, one for clear-sky conditions (2000–2017) and the other one in all-sky conditions (2000–2016). The clear-sky instantaneous CSIF (CSIF(sub clear-inst)) shows high accuracy against the clear-sky OCO-2 SIF and little bias across biome types. The all-sky daily average CSIF (CSIF(sub all-daily)) dataset exhibits strong spatial, seasonal and interannual dynamics that are consistent with daily SIF from OCO-2 and the Global Ozone Monitoring Experiment-2 (GOME-2). An increasing trend (0.39 %) of annual average CSIFall-daily is also found, confirming the greening of Earth in most regions. Since the difference between satellite-observed SIF and CSIF is mostly caused by the environmental down-regulation on SIF(sub yield), the ratio between OCO-2 SIF and CSIF(sub clear-inst) can be an effective indicator of drought stress that is more sensitive than the normalized difference vegetation index and enhanced vegetation index. By comparing CSIF(sub all-daily) with GPP estimates from 40 EC flux towers across the globe, we find a large cross-site variation (c.v. = 0.36) of the GPP–SIF relationship with the highest regression slopes for evergreen needleleaf forest. However, the cross-biome variation is relatively limited (c.v. = 0.15). These two contiguous SIF datasets and the derived GPP–SIF relationship enable a better understanding of the spatial and temporal variations of the GPP across biomes and climate.

Zhang, Yao↗

NASA NTRS1992

Maximum likelihood techniques applied to quasi-elastic light scattering

There is a necessity of having an automatic procedure for reliable estimation of the quality of the measurement of particle size from QELS (Quasi-Elastic Light Scattering). Getting the measurement itself, before any error estimates can be made, is a problem because it is obtained by a very indirect measurement of a signal derived from the motion of particles in the system and requires the solution of an inverse problem. The eigenvalue structure of the transform that generates the signal is such that an arbitrarily small amount of noise can obliterate parts of any practical inversion spectrum. This project uses the Maximum Likelihood Estimation (MLE) as a framework to generate a theory and a functioning set of software to oversee the measurement process and extract the particle size information, while at the same time providing error estimates for those measurements. The theory involved verifying a correct form of the covariance matrix for the noise on the measurement and then estimating particle size parameters using a modified histogram approach.

Edwards, Robert V.↗

NASA NTRS1994

On the Lighthill relationship and sound generation from isotropic turbulence

In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.

Zhou, YE↗

NASA NTRS1997

Relative Navigation In Elliptical Orbits Using An Iterative Nonlinear Filter

The two step filter is applied to process intersatellite radar measurements to determine the motion of one satellite relative to another in close elliptical orbits. This filter breaks a nonlinear estimation problem into two state vectors. The 'first step' state is chosen so as to have a linear measurement equation. This is nonlinearly related to the 'second step' state which describes the dynamics. Two different forms are used. In one, the first step state is the second step state vector augmented by the measurement equation. In the other, the first step and second step state vectors are of equal dimension. The two step filter is compared against an iterated extended Kalman filter and a Kalman filter using a change of variables. Analytical differences between the two step estimator and these conventional filters are highlighted. Special concerns for initializing the first step state covariance matrix and handling the possibility of numerically rank deficient covariance matrices are addressed. Numerical simulations are performed which show that the Two Step estimator produces a lower estimation bias under two circumstances; large apriori initial error; and small dimension observation vectors which require a longer arc of measurements to generate observability of the state.

Garrison, James L.↗

NASA NTRS1988

Airborne radar technology for windshear detection

The objectives and accomplishments of the two-and-a-half year effort to describe how returns from on-board Doppler radar are to be used to detect the presence of a wind shear are reported. The problem is modeled as one of first passage in terms of state variables, the state estimates are generated by a bank of extended Kalman filters working in parallel, and the decision strategy involves the use of a voting algorithm for a series of likelihood ratio tests. The performance issue for filtering is addressed in terms of error-covariance reduction and filter divergence, and the performance issue for detection is addressed in terms of using a probability measure transformation to derive theoretical expressions for the error probabilities of a false alarm and a miss.

Hibey, Joseph L.↗

NASA NTRS2017

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

NASA NTRS1995

On-line estimation of error covariance parameters for atmospheric data assimilation

A simple scheme is presented for on-line estimation of covariance parameters in statistical data assimilation systems. The scheme is based on a maximum-likelihood approach in which estimates are produced on the basis of a single batch of simultaneous observations. Simple-sample covariance estimation is reasonable as long as the number of available observations exceeds the number of tunable parameters by two or three orders of magnitude. Not much is known at present about model error associated with actual forecast systems. Our scheme can be used to estimate some important statistical model error parameters such as regionally averaged variances or characteristic correlation length scales. The advantage of the single-sample approach is that it does not rely on any assumptions about the temporal behavior of the covariance parameters: time-dependent parameter estimates can be continuously adjusted on the basis of current observations. This is of practical importance since it is likely to be the case that both model error and observation error strongly depend on the actual state of the atmosphere. The single-sample estimation scheme can be incorporated into any four-dimensional statistical data assimilation system that involves explicit calculation of forecast error covariances, including optimal interpolation (OI) and the simplified Kalman filter (SKF). The computational cost of the scheme is high but not prohibitive; on-line estimation of one or two covariance parameters in each analysis box of an operational bozed-OI system is currently feasible. A number of numerical experiments performed with an adaptive SKF and an adaptive version of OI, using a linear two-dimensional shallow-water model and artificially generated model error are described. The performance of the nonadaptive versions of these methods turns out to depend rather strongly on correct specification of model error parameters. These parameters are estimated under a variety of conditions, including uniformly distributed model error and time-dependent model error statistics.

Dee, Dick P.↗

NASA NTRS1995

Latitudinal distribution of black carbon soot in the upper troposphere and lower stratosphere

Black carbon soot from the upper troposphere and lower stratosphere has been systematically collected at latitudes from 90 deg N to 45 deg S. The measured latitudinal distribution of this soot at 10 to 11 km altitude is found to covary with commercial air traffic fuel use, suggesting that aircraft fuel combustion at altitude is the principal source. In addition, at latitudes where the commercial air traffic is high, measured black carbon soot values are high even at 20 km altitude, suggesting that aircraft-generated soot injected just above the tropopause may be transported to higher altitudes. During the volcanically influenced period in which these samples were collected, the number abundances, total mass, and calculated total surface area of black carbon soot are 2-3 orders of magnitude lower than similar measures of sulfuric acid aerosol. During volcanically quiescent periods, the calculated total surface area of black carbon soot aerosol is of the same order of magnitude as that of the background sulfuric acid aerosol. It appears from this comparison that black carbon soot is only capable of influencing lower stratosphere or upper troposphere chemistry during periods when the aerosol budget is not dominated by volcanic activity. It remains to determine the extent to which black carbon soot particles act as nuclei for sulfuric acid aerosol formation. However, mass balance calculations suggest that aircraft soot injected at altitude does not represent a significant source of condensation nuclei for sulfuric acid aerosols.

Blake, David F.↗

NASA NTRS1992

The processing of the Viking Orbiter range data and its contribution to Mars gravity solutions

The processing of Doppler data has been the primary method for deriving models of the Mars gravity field. Since the Mariner 9 and Viking spacecraft were placed in orbit about Mars, many models from degree and order 6 to degree and order 50 have been developed. However, during the Viking mission, some 26,000 range measurements to the two Viking Orbiters were also obtained. These data have not previously been used in the derivation of Mars gravity models. A portion of these range data have been processed simultaneously with the Doppler data. Normal equations were generated for both sets of data and were used to create two solutions complete to degree and order 30: a nominal solution including both the range and the Doppler data (MGM-R100), and another solution including only the Doppler data (MGM-R101). Tests with the covariances of these solutions, as well as with orbit overlap tests indicate that the interplanetary range data can be used to improve the modeling of the Mars gravity field.

Lemoine, Frank G.↗

NASA NTRS1992

Correlation, Cost Risk, and Geometry

The geometric viewpoint identifies the choice of a correlation matrix for the simulation of cost risk with the pairwise choice of data vectors corresponding to the parameters used to obtain cost risk. The correlation coefficient is the cosine of the angle between the data vectors after translation to an origin at the mean and normalization for magnitude. Thus correlation is equivalent to expressing the data in terms of a non orthogonal basis. To understand the many resulting phenomena requires the use of the tensor concept of raising the index to transform the measured and observed covariant components into contravariant components before vector addition can be applied. The geometric viewpoint also demonstrates that correlation and covariance are geometric properties, as opposed to purely statistical properties, of the variates. Thus, variates from different distributions may be correlated, as desired, after selection from independent distributions. By determining the principal components of the correlation matrix, variates with the desired mean, magnitude, and correlation can be generated through linear transforms which include the eigenvalues and the eigenvectors of the correlation matrix. The conversion of the data to a non orthogonal basis uses a compound linear transformation which distorts or stretches the data space. Hence, the correlated data does not have the same properties as the uncorrelated data used to generate it. This phenomena is responsible for seemingly strange observations such as the fact that the marginal distributions of the correlated data can be quite different from the distributions used to generate the data. The joint effect of statistical distributions and correlation remains a fertile area for further research. In terms of application to cost estimating, the geometric approach demonstrates that the estimator must have data and must understand that data in order to properly choose the correlation matrix appropriate for a given estimate. There is a general feeling by employers and managers that the field of cost requires little technical or mathematical background. Contrary to that opinion, this paper demonstrates that a background in mathematics equivalent to that needed for typical engineering and scientific disciplines at the masters or doctorate level is appropriate within the field of cost risk.

Dean, Edwin B.↗

NASA NTRS2001

DAO's Next Generation Physical-Space/Finite-Volume Data Assimilation System: Formulation and Initial Evaluation

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will describe the general system formulation, the adaptive quality control and general aspects of the error covariance modeling. The NASA-NCAR GCM is a completely new model which replaces the CEOs GCM used in the previous GEOS-1/2/3 Data Assimilation systems. A particular configuration of adaptive Statistical Quality Control and the Physical-space Statistical Analysis System (PSAS) are currently implemented in DAO's operational Data Assimilation System. However, the unique finite-volume formulation of the NASA-NCAR GCM, combined with the generality of the observation-space formulation of PSAS, provides for a very simple and accurate model-analysis interface. The system assimilates a variety of conventional and satellite observations. In particular, TOVS Level 1B radiances are assimilated using a 1-D variational scheme, both in clear sky and cloudy conditions. Computationally, the fvDAS runs approximately 10 times faster than the operational GEOS-Terra system. We will show that the next-generation fvDAS has much improved observation-minus-6hr forecast (O-F) statistics, as well as 5-day forecast skills. Top of the atmosphere radiation fields are in closer agreement with CERES measurements, with realistic precipitation and moisture fields. We will also show that the finite-volume formulation of the fvDAS produce assimilated fields which are more suitable for driving constituent transport models.

daSilva, A.↗