Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

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↗
NASA NTRS2012

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↗
NASA NTRS2014

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↗
DOE OSTI2022

Representation of Leaf-to-Canopy Radiative Transfer Processes Improves Simulation of Far-Red Solar-Induced Chlorophyll Fluorescence in the Community Land Model Version 5

Recent advances in satellite observations of solar-induced chlorophyll fluorescence (SIF) provide a new opportunity to evaluate and constrain the simulation of terrestrial gross primary productivity (GPP). Accurate representation of the processes driving SIF emission and the radiative transfer of SIF to remote sensing sensors is an essential prerequisite for the evaluation and data assimilation. Recently, SIF simulations have been incorporated into several land surface models, but the scaling of SIF from leaf-level to canopy level is usually not well-represented. In this work, we incorporate the simulation of far-red SIF observed at nadir into the Community Land Model version 5 (CLM5). An efficient and accurate method based on escape probability is developed to scale SIF from leaf-level to top-of-canopy while taking clumping and the radiative transfer processes into account. SIF simulated by CLM5 and a canopy-level model agreed well at sites except one in needle leaf forest (R2>0.91, root-mean-square error < 0.19W m -2 sr -1 um -1 ), and captured the day-to-day variation of tower-measured SIF at temperate forest sites (R 2 >0.68). At the global scale, simulated SIF generally captured the spatial and seasonal pat37 terns of satellite-observed SIF (R2 > 0.76 except for tropical forest). Factors including the fluorescence emission model, clumping, bidirectional effect, and canopy properties (leaf optical properties and leaf area index) had considerable impacts on SIF simulation, and the discrepancies between simulated and observed SIF varied with plant functional type. By improving the representation of radiative transfer for SIF simulation, our model allows better comparisons between simulated and observed SIF towards constraining and evaluating GPP simulations.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2023

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↗
NASA NTRS1980

Iron and titanium distribution on the moon from orbital gamma ray spectrometry with implications for crustal evolutionary models

A set of Fe and Ti maps and regional values are obtained from the Apollo 15 and 16 orbital gamma ray data by energy band analysis. High-Ti basalts predominate the early and late stages of mare volcanism with high-Fe basaltic volcanism in the interim. The first evidence of a high-Ti-KREEP basalt association is found in the Aristarchus region. A N-S asymmetry for Fe and Ti in the east limb and farside highlands complicates the E-W asymmetry for Th but substantiates crustal inhomogeneity. The observed crustal inhomogeneity adds an additional objection to the primitive source model for crustal evolution. The high-Ti-KREEP basalt association and the general trend of decreasing mare basalt Ti with time lend support to the cumulate source model; however, this model cannot account for young, high-Ti maria. The dynamic assimilation model better accounts for chemical variations observed on the moon.

Davis, P. A., Jr.↗
NASA NTRS2006

Detection of Ice Polar Stratospheric Clouds from Assimilation of Atmospheric Infrared Sounder Data

A novel technique is presented for detection of ice polar stratospheric clouds (PSCs) that form at extremely low temperatures in the lower polar stratosphere during winter. Temperature is a major factor in determining abundance of PSCs, which in turn provide surfaces for heterogeneous chemical reactions leading to ozone loss and radiative cooling. The technique infers the presence of ice PSCs using radiances from the Atmospheric Infrared Sounder (AIRS) in the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. Brightness temperatures are computed from short-term GEOS-5 forecasts for several hundred AIRS channels, using a radiation transfer module. The differences between collocated AIRS observations and these computed values are the observed-minus-forecast (O-F) residuals in the assimilation system. Because the radiation model assumes clear-sky conditions, we hypothesize that these O-F residuals contain quantitative information about PSCs. This is confirmed using sparse data from the Polar Ozone and Aerosol Measurement (POAM) III occultation instrument. The analysis focuses on 0-F residuals for the 6.79pm AIRS moisture channel. At coincident locations, when POAM III detects ice clouds, the AIRS O-F residuals for this channel are lower than -2K. When no ice PSCs are evident in POAM III data, the AIRS 0-F residuals are larger. Given this relationship, the high spatial density of AIRS data is used to construct maps of regions where 0-F residuals are lower than -2K, as a proxy for ice PSCs. The spatial scales and spatio-temporal variations of these PSCs in the Antarctic and Arctic are discussed on the basis of these maps.

Stajner, Ivanka↗
NASA NTRS1994

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗
NASA NTRS1987

Variational four-dimensional analysis using quasi-geostrophic constraints

A variational four-dimensional analysis technique using quasi-geostrophic models as constraints is examined using gridded fields as data. The analysis method uses a standard iterative nonlinear minimization technique to find the solution to the constraining forecast model which best fits the data as measured by a predefined functional. The minimization algorithm uses the derivative of the functional with respect to each of the initial condition values. This derivative vector is found by inserting the weighted differences between the model solution and the inserted data into a backwards integrating adjoint model. The four-dimensional analysis system was examined by applying it to fields created from a primitive equations model forecast and to fields created from satellite retrievals. The results show that the technique has several interesting characteristics not found in more traditional four-dimensional assimilation techniques. These features include a close fit of the model solution to the observations throughout the analysis interval and an insensitivity to the frequency of data insertion or the amount of data. The four-dimensional analysis technique is very versatile and can be extended to more complex problems with little theoretical difficulty.

Derber, John C.↗
NASA NTRS2006

A Test of Sensitivity to Convective Transport in a Global Atmospheric CO2 Simulation

Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.

Bian, H.↗
NASA NTRS2004

Alternative Approaches to Land Initialization for Seasonal Precipitation and Temperature Forecasts

The seasonal prediction system of the NASA Global Modeling and Assimilation Office is used to generate ensembles of summer forecasts utilizing realistic soil moisture initialization. To derive the realistic land states, we drive offline the system's land model with realistic meteorological forcing over the period 1979-1993 (in cooperation with the Global Land Data Assimilation System project at GSFC) and then extract the state variables' values on the chosen forecast start dates. A parallel series of forecast ensembles is performed with a random (though climatologically consistent) set of land initial conditions; by comparing the two sets of ensembles, we can isolate the impact of land initialization on forecast skill from that of the imposed SSTs. The base initialization experiment is supplemented with several forecast ensembles that use alternative initialization techniques. One ensemble addresses the impact of minimizing climate drift in the system through the scaling of the initial conditions, and another is designed to isolate the importance of the precipitation signal from that of all other signals in the antecedent offline forcing. A third ensemble includes a more realistic initialization of the atmosphere along with the land initialization. The impact of each variation on forecast skill is quantified.

Koster, Randal↗
NASA NTRS2018

Studying Diurnal Variations of Aerosols with NASA MERRA-2 Reanalysis Data

Aerosols play an important role in atmospheric dynamics, climate variations, and Earth's energy cycle by altering the radiation balance in the atmosphere through interaction with clouds, providing fertilizer for forests and canopy, and as a supply of iron to the ocean over long time periods. Studies suggest that much of the feedback between dust aerosols and dynamics is associated with diurnal and synoptic scale variability. However, the lack of sub-daily resolution of aerosols from satellite observations makes it difficult to study the diurnal characteristics, especially over tropical and subtropical regions. Investigation of this topic utilizes over 37 years of simulated global aerosol products from NASA atmospheric reanalysis, in the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) data set, available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). MERRA-2 covers the period 1980-present, and is continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using data from MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated using the MERRA-2 aerosol model, which interacts directly with radiation parameterization, and is radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Hourly, monthly, and monthly diurnal data are available at spatial resolution of 0.5o x 0.625o (latitude x longitude). By using MERRA-2 hourly and monthly diurnal products, different aerosol diurnal variabilities are observed over North America, Africa, Asia, and Australia, that may be due to different meteorological conditions and aerosol sources. The presentation will also provide an overview of MERRA-2 data services at GES DISC, such as how to find and download data, and how to quickly visualize and analyze data online with Giovanni.

aerosol↗
NASA NTRS2016

Maps Suggest Transport and Source Processes of PM2.5 at 1 km x 1 km for the Whole San Joaquin Valley, Winter 2011 (Generalizations from DISCOVER-AQ)

We present interpreted data analysis using MAIAC (Multiangle implementation of Atmospheric Correction) retrievals and appropriate RAPid Update Cycle (RAP) meteorology to map respirable aerosol (PM2.5) for the period January and February, 2011. The San Joaquin Valley is one of the unhealthiest regions in the USA for PM2.5 and related morbidity. The methodology evaluated can be used for the entire moderate-resolution imaging spectrometer (MODIS, VIIRS) data record. Other difficult areas of the West: Riverside, CA, Salt Lake City, UT, and Doa Ana County, NM share similar difficulties and solutions. The maps of boundary layer depth for 1116 hr local time from RAP allows us to interpret aerosol optical thickness as a concentration of particles in a nearly well-mixed box capped by clean air. That mixing is demonstrated by DISCOVER-AQ data and afternoon samples from the airborne measurements, P3B (on-board) and B200 (HSRL2 lidar). This data and the PM2.5 gathered at the deployment sites allowed us to estimate and then evaluate consistency and daily variation of the AOT to PM2.5 relationship. Mixed-effects modeling allowed a refinement of that relation from day to day; RAP mixed layers explained the success of previous mixed-effects modeling. Compositional, size-distribution, and MODIS angle-of-regard effects seem to describe the need for residual daily correction beyond ML depth.We report on an extension method to the entire San Joaquin Valley for all days with MODIS imagery using the permanent PM2.5 stations, evaluated for representativeness. Resulting map movies show distinct sources, particularly Interstate-5 (at approx. 1km x 1km resolution) and the broader Bakersfield area. Accompanying winds suggest transport effects and variable pathways of pollution cleanout. Such estimates should allow morbiditymortality studies. They should be also useful for actual model assimilations, where composition and sources are uncertain. We conclude with a description of new work to extend these insights to similar regions, e.g. interior valleys of California, the Po Valley, the Mediterranean litoral, and the Ganges Plain.This work show generalizable use of remote sensing, a major goal of DISCOVER-AQ, Deriving Information on Surface Conditions from COlumn and VERtically Resolved Observations Relevant to Air Quality.

Chatfield, R.↗
DOE OSTI2022

Beyond carbon flux partitioning: Carbon allocation and nonstructural carbon dynamics inferred from continuous fluxes

Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Highly Siderophile Element Abundances in Martian Meteorites

Critical evaluation of new and literature data for highly siderophile elements (HSE) in Martian (SNC) meteorites allows several first order conclusions to be drawn. (i) Re concentrations in SNC meteorites are nearly constant (within a factor of two) and do not correlate with rock type. Exceptions to this rule are Chassigny and Dar al Gani (DaG) 476, both of which are inferred to have experienced terrestrial Re contamination. (ii) Fractionations between Rh and Pd are small. Excluding Shergotty, the Rh/Pd ratio of the SNC suite is 0.22\pm0.05. (iii) Os and Ir contents vary by about four orders of magnitude; and positive correlations with MgO, Cr, and Ni suggest that these variations are not controlled by sulfide fractionation. A possible exception is the orthopyroxenite ALH84001, whose HSE's (including Ni, which is compatible in opx) are very low. (iv) Zagami, Shergotty, and Nakhla have nearly identical HSE signatures. Shergotty and Zagami have experienced assimilation-fractional crystallization (AFC) and have "crustal" Sr and Nd isotopic signatures. Conversely, the Nakhla parent was a small degree partial melt of a depleted mantle that interacted little with the Martian crust. These observations suggest that "evolved" HSE signatures can be produced by either fractional crystallization or small degrees of partial melting. (v) Chassigny and other mafic SNC's have HSE signatures that are very distinct from those of Nakhla-Zagami-Shergotty. The HSE elemental ratios of mafic SNC's approach chondritic, implying that the Martian mantle has nearly chondritic relative abundances of the HSE's. (vi) This chondritic HSE signature is observed in SNC's of various ages, suggesting that this is an ancient feature that has not evolved over time. (vii) No correlation is observed between HSE's and signatures of crustal contamination (e.g., Sr isotopes), indicating that the HSE signatures of the SNC suite are not derived from the crust. (vii) The Ru/Pd for the SNC suite ratio is about 0.5 chondritic. We attribute this depletion of Ru to the presence of spinel in the SNC source regions. (viii) Assignment of this depletion of Ru to spinel seems to also imply that spinel is not the cause of Os and Ir fractionations, requiring some other mafic phase to account for the observed variations. -----.

Jones, J. H.↗
NASA NTRS1979

The origin of KREEP

KREEP is a lunar material having very high concentrations of incompatible elements; its name is an acronym for the incompatibles K, rare-earth elements (REE), and P. Although a few pristine (endogenously igneous) KREEPy samples were returned from the Apollo 15 and 17 sites, most KREEPy samples are polymict breccias. Most models of KREEP petrogenesis have been based on partial melting of a variety of sources. Such models fail to explain the veritable absence of variations in incompatible element patterns over the sampled portion of the moon. We have defined a KREEP component based on the average composition of Apollo 14 breccias having extremely high concentrations of incompatible elements. Normalization of accurate incompatible data for KREEPy samples from the Apollo 12, 14, 15, 16, and 17 sites to this component virtually always shows no resolvable fractionation (e.g., <10% variation in the La/Lu ratio), whereas partial melting models typically produce larger fractionations (±20–25% in La/Lu) from a factor of 2 difference in degree of partial melting. Required is a single major source that could provide KREEP to widely separated locations on the nearside of the moon. The anorthositic crust of the moon is commonly attributed to the flotation of plagioclase on a deep, moon-wide magma ocean. Fractional crystallization of this magma ocean would have produced large enrichments of incompatibles in a residual liquid. No other plausible major source of incompatibles has been proposed. We borrow the German prefix ur—meaning primeval and designate this residual liquid ‘urKREEP.’ We propose that all KREEPy rocks originated by dilution of urKREEP with crustal or mantle materials during assimilation, or zone-refining (pristine samples), or impact-induced brecciation (breccias and melt rocks). The formation of urKREEP cannot be dated precisely. Correction of breccia Rb-Sr model ages for Rb loss or gain during the early intense bombardments yields ages that cluster in the range 4.4–4.5 Gy. This implies that crystallization of the magma ocean was essentially complete at this time and is in general agreement with U-Pb evidence indicating crustal formation at 4.4 Gy. Assuming that the moon had the composition of an H-group chondrite depleted in Fe-Ni and FeS and that half the incompatibles fractionated into materials other than urKREEP, the thickness of a moon-wide urKREEP layer was <2 km. Thorium concentrations determined by gamma ray spectroscopy indicate that about 4% of the incompatibles in an H chondritic moon are now in the outermost kilometer.

Paul H. Warren↗
NASA NTRS1999

Monitoring Global Geophysical Fluids by Space Geodesy

Since its establishment on 1/1/1998 by the International Earth Rotation Service, the Coordinating Center for Monitoring Global Geophysical Fluids (MGGF) and its seven Special Bureaus have engaged in an effort to support and facilitate the understanding of the geophysical fluids in global geodynamics research. Mass transports in the atmosphere-hydrosphere-solid Earth-core system (the "global geophysical fluids") will cause the following geodynamic effects on a broad time scale: (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, with ever increasing precision/accuracy and temporal/spatial resolution. Each of the seven Special Bureaus within MGGF is responsible for calculations related to a specific Earth component or aspect -- Atmosphere, Ocean, Hydrology, Ocean Tides, Mantle, Core, and Gravity/Geocenter. Angular momenta and torques, gravitational coefficients, and geocenter shift will be computed for geophysical fluids based on global observational data, and from state-of-the-art models, some of which assimilate such data. The computed quantities, algorithm and data formats are standardized. The results are archived and made available to the scientific research community. This paper reports the status of the MGGF activities and current results.

Chao, Benjamin F.↗
NASA NTRS2000

Climatology of the Finite-Volume CCM and Sensitivity Studies

A new model has been developed in the Data Assimilation Office at NASA Goddard Space Flight Center. It is intended for data assimilation and general circulation modeling, including representations of radiatively active trace gases such as stratospheric ozone. The model is based on the flux-form semi-Lagrangian' dynamical core, developed by Lin and Rood. In collaboration with the National Center for Atmospheric Research, the physical parametrization package from the Community Climate model, Version 3, has been included. This paper will discuss the climatology of a multi-year integration of the model, with observed boundary conditions (sea-surface temperature, sea ice, and incoming solar radiation). The standard model version has a resolution of 2.5 degrees longitude, 2 degrees latitude and has 55 levels between the surface and 80km. The mean climate and interannual variability of this model will be discussed, with emphasis on the middle atmosphere. The sensitivity to horizontal resolution will be examined, as will the impact of increasing and decreasing the vertical resolution. Changes which differ significantly from the interannual variations of the control run will be highlighted. Impacts of using long-term mean boundary conditions or values which vary from year to year will be discussed. Changes in the tropospheric circulation induced by lowering the upper boundary of the model will be discussed. The emphasis of the paper is on the climatology of the middle atmosphere and its impacts on the troposphere, with particular attention given to the climatologically important region near the tropopause.

Pawson, Steven↗