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.

448 records · Page 25

Leaf traits and canopy structure together explain canopy functional diversity: an airborne remote sensing approach

Plant functional diversity is strongly connected to photosynthetic carbon assimilation in terrestrial ecosystems. However, many of the plant functional traits that regulate photosynthetic capacity, including foliar nitrogen concentration and leaf mass per area, vary significantly between and within plant functional types and vertically through forest canopies, resulting in considerable landscape-scale heterogeneity in three dimensions. Hyperspectral imagery has been used extensively to quantify functional traits across a range of ecosystems but is generally limited to providing information for top-of-canopy leaves only. On the other hand, lidar data can be used to retrieve the vertical structure of forest canopies. Because these data are rarely collected at the same time, there are unanswered questions about the effect of forest structure on the 3-D spatial patterns of functional traits across ecosystems. In the United States, the National Ecological Observatory Network's Airborne Observation Platform (NEON AOP) provides an opportunity to address this structure-function relationship by collecting lidar and hyperspectral data together across a variety of ecoregions. With a fusion of hyperspectral and lidar data from the NEON AOP and field-collected foliar trait data, we assessed the impacts of forest structure on spatial patterns of N. In addition, we examine the influence of abiotic gradients and management regimes on top of canopy %N and total canopy N (i.e. the total amount of N (g/m 2 ) within a forest canopy) at a NEON site consisting of a mosaic of open longleaf pine and dense broadleaf deciduous forests. Here, our resulting maps suggest that in contrast with top-of-canopy values, total canopy N variation is dampened across this landscape resulting in relatively homogeneous spatial patterns. At the same time, we found that leaf functional diversity and canopy structural diversity showed distinct dendritic patterns related to the spatial distribution of plant functional types.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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.↗

Prior biosphere model impact on global terrestrial CO2 fluxes estimated from OCO-2 retrievals

This study assesses the impact of different state of the art global biospheric CO2 flux models, when applied as prior information, on inverse model “top-down” estimates of terrestrial CO2 fluxes obtained when assimilating Orbiting Carbon Observatory 2 (OCO-2) observations. This is done with a series of observing system simulation experiments (OSSEs) using synthetic CO2 column-average dry air mole fraction (XCO2) retrievals sampled at the OCO-2 satellite spatiotemporal frequency. The OSSEs utilized a 4-D variational (4D-Var) assimilation system with the GEOS-Chem global chemical transport model (CTM) to estimate CO2 net ecosystem exchange (NEE) fluxes using synthetic OCO-2 observations. The impact of biosphere models in inverse model estimates of NEE is quantified by conducting OSSEs using the NASA-CASA, CASA-GFED, SiB-4, and LPJ models as prior estimates and using NEE from the multi-model ensemble mean of the Multiscale Synthesis and Terrestrial Model Intercomparison Project as the “truth”. Results show that the assimilation of simulated XCO2 retrievals at OCO-2 observing modes over land results in posterior NEE estimates which generally reproduce “true” NEE globally and over terrestrial TransCom-3 regions that are well-sampled. However, we find larger spread among posterior NEE estimates, when using different prior NEE fluxes, in regions and seasons that have limited OCO-2 observational coverage and a large range in “bottom-up” NEE fluxes. Seasonally averaged posterior NEE estimates had standard deviations (SD) of ∼10 % to ∼50 % of the multi-model-mean NEE for different TransCom-3 land regions with significant NEE fluxes (regions/seasons with a NEE flux ≥0.5 PgC/yr). On a global average, the seasonally averaged residual impact of the prior model NEE assumption on the posterior NEE spread is ∼10 %–20 % of the posterior NEE mean. Additional OCO-2 OSSE simulations demonstrate that posterior NEE estimates are also sensitive to the assumed prior NEE flux uncertainty statistics, with spread in posterior NEE estimates similar to those when using variable prior model NEE fluxes. In fact, the sensitivity of posterior NEE estimates to prior error statistics was larger than prior flux values in some regions/times in the tropics and Southern Hemisphere where sufficient OCO-2 data were available and large differences between the prior and truth were evident. Overall, even with the availability of spatiotemporally dense OCO-2 data, noticeable residual differences (up to ∼20 %–30 % globally and 50 % regionally) in posterior NEE flux estimates remain that were caused by the choice of prior model flux values and the specification of prior flux uncertainties.

Sajeev Philip↗

Estimation of Mesoscale Atmospheric Latent Heating Profiles from TRMM Rain Statistics Utilizing a Simple One-Dimensional Model

In this study, a model is developed to estimate mesoscale-resolution atmospheric latent heating (ALH) profiles. It utilizes rain statistics deduced from Tropical Rainfall Measuring Mission (TRMM) data, and cloud vertical velocity profiles and regional surface thermodynamic climatologies derived from other available data sources. From several rain events observed over tropical ocean and land, ALH profiles retrieved by this model in convective rain regions reveal strong warming throughout most of the troposphere, while in stratiform rain regions they usually show slight cooling below the freezing level and significant warming above. The mesoscale-average, or total, ALH profiles reveal a dominant stratiform character, because stratiform rain areas are usually much larger than convective rain areas. Sensitivity tests of the model show that total ALH at a given tropospheric level varies by less than +/- 10 % when convective and stratiform rain rates and mesoscale fractional rain areas are perturbed individually by +/- 15 %. This is also found when the non-uniform convective vertical velocity profiles are replaced by one that is uniform. Larger variability of the total ALH profiles arises when climatological ocean- and land-surface temperatures (water vapor mixing ratios) are independently perturbed by +/- 1.0 K (+/- 5%) and +/- 5.0 K (+/- 15%), respectively. At a given tropospheric level, such perturbations can cause a +/- 25% variation of total ALH over ocean, and a factor-of-two sensitivity over land. This sensitivity is reduced substantially if perturbations of surface thermodynamic variables do not change surface relative humidity, or are not extended throughout the entire model evaporation layer. The ALH profiles retrieved in this study agree qualitatively with tropical total diabatic heating profiles deduced in earlier studies. Also, from January and July 1999 ALH-profile climatologies generated separately with TRMM Microwave Imager and Precipitation Radar rain statistics, it is shown that ALH profiles can be retrieved utilizing diverse satellite-derived rain products that offer convective and stratiform discrimination. Therefore, the ALH retrieval model developed in this study can be used to make regional estimates of total diabatic heating profiles in the future Global Precipitation Measurement mission, and to assimilate these profiles into numerical weather forecast and climate models.

Iacovazzi, Robert A., Jr.↗

Estimation of Mesoscale Atmospheric Latent Heating Profiles from TRMM Rain Statistics Utilizing a Simple One-Dimensional Model

In this study, a model is developed to estimate mesoscale-resolution atmospheric latent heating (ALH) profiles. It utilizes rain statistics deduced from Tropical Rainfall Measuring Mission (TRMM) data, and cloud vertical velocity profiles and regional surface thermodynamic climatologies derived from other available data sources. From several rain events observed over tropical ocean and land, ALH profiles retrieved by this model in convective rain regions reveal strong warming throughout most of the troposphere, while in stratiform rain regions they usually show slight cooling below the freezing level and significant warming above. The mesoscale-average, or total, ALH profiles reveal a dominant stratiform character, because stratiform rain areas are usually much larger than convective rain areas. Sensitivity tests of the model show that total ALH at a given tropospheric level varies by less than +/- 10 % when convective and stratiform rain rates and mesoscale fractional rain areas are perturbed individually by 1 15 %. This is also found when the non-uniform convective vertical velocity profiles are replaced by one that is uniform. Larger variability of the total ALH profiles arises when climatological ocean- and land-surface temperatures (water vapor mixing ratios) are independently perturbed by +/- 1.0 K (+/- 5 %) and +/- 5.0 K (+/- 15 %), respectively. At a given tropospheric level, such perturbations can cause a +/- 25 % variation of total ALH over ocean, and a factor-of-two sensitivity over land. This sensitivity is reduced substantially if perturbations of surface thermodynamic variables do not change surface relative humidity, or are not extended throughout the entire model evaporation layer. The ALH profiles retrieved in this study agree qualitatively with tropical total diabatic heating profiles deduced in earlier studies. Also, from January and July 1999 ALH-profile climatologies generated separately with TRMM Microwave Imager and Precipitation Radar rain statistics, it is shown that ALH profiles can be retrieved utilizing diverse satellite-derived rain products that offer convective and stratiform discrimination. Therefore, the ALH retrieval model developed in this study can be used to make regional estimates of total diabatic heating profiles in the future Global Precipitation Measurement mission, and to assimilate these profiles into numerical weather forecast and climate models.

Iacovazzi, Robert A., Jr.↗

The Re-Analysis of Ozone Profile Data from a 41-Year Series of SBUV Instruments

In this study we present the validation of ozone profiles from a number of Solar Back Scattered Ultra Violet (SBUV) and SBUV/2 instruments that were recently reprocessed using an updated (Version 8.6) algorithm. The SBUV dataset provides the longest available record of global ozone profiles, spanning a 41-year period from 1970 to 2011 (except a 5-year gap in the 1970s) and includes ozone profile records obtained from the Nimbus-4 BUV and Nimbus-7 SBUV instruments, and a series of SBUV(/2) instruments launched on NOAA operational satellites (NOAA 09, 11, 14, 16, 17, 18, 19). Although modifications in instrument design were made in the evolution from the BUV instrument to the modern SBUV(/2) model, the basic principles of the measurement technique and retrieval algorithm remain the same. The long term SBUV data record allows us to create a consistent, calibrated dataset of ozone profiles that can be used for climate studies and trend analyses. In particular, we focus on estimating the various sources of error in the SBUV profile ozone retrievals using independent observations and analysis of the algorithm itself. For the first time we include in the metadata a quantitative estimate of the smoothing error, defined as the error due to profile variability that the SBUV observing system cannot inherently measure. The magnitude of the smoothing error varies with altitude, latitude, season and solar zenith angle. Between 10 and 1 hPa the smoothing errors for the SBUV monthly zonal mean retrievals are of the order of 1 %, but start to increase above and below this layer. The largest smoothing errors, as large as 15-20%, were detected in in the troposphere. The SBUV averaging kernels, provided with the ozone profiles in version 8.6, help to eliminate the smoothing effect when comparing the SBUV profiles with high vertical resolution measurements, and make it convenient to use the SBUV ozone profiles for data assimilation and model validation purposes. The smoothing error can also be minimized by combining layers of data, and we will discuss recommendations for this approach as well. The SBUV ozone profiles have been intensively validated against satellite profile measurements obtained from the Microwave Limb Sounders (MLS) (on board the UARS and AURA satellites), Stratospheric Aerosol and Gas Experiment (SAGE) and Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Also, we compare coincident and collocated SBUV ozone retrievals with observations made by ground-based instruments, such as microwave spectrometers, lidars, Umkehr instruments and balloon-borne ozonosondes. Finally, we compare the SBUV ozone profiles with output from the NASA GSFC GEOS-CCM model. In the stratosphere between 25 and 1 hPa the mean biases and standard deviations are within 5% for monthly mean ozone profiles. Above and below this layer the vertical resolution of the SBUV algorithm decreases and the effects of vertical smoothing should be taken into account. Though the SBUV algorithm has a coarser vertical resolution in the lower stratosphere and troposphere, it is capable of precisely estimating the integrated ozone column between the surface and 25 hPa. The time series of the tropospheric - lower stratospheric ozone column derived from SBUV agrees within 5% with the corresponding values observed by an ensemble of ozone sonde stations in North Hemisphere. Drift of the ozone time series obtained from each SBUV(/2) instrument relative to ground based and satellite measurements are evaluated and some features of individual SBUV(l2) instruments are discussed. In addition to evaluating individual instruments against independent observations, we also focus on the instrument to instrument consistency in the series. Overall, Version 8.6 ozone profiles obtained from two different SBUV(l2) instruments compare within a couple of percent during overlap periods and are consistently varying in time, with some exceptions. Some of the noted discrepancies might bssociated with ozone diurnal variations, since the difference in the local time of the observations for a pair of SBUV(l2) instruments could be several hours. Other issues include the potential short-term drift in measurements as the instrument orbit drifts, and measurements are obtained at high solar zenith angles (>85 ). Based on the results of the validation, a consistent, calibrated dataset of SBUV ozone profiles has been created based on internal calibration only.

Kramarova, Natalya↗

Volatiles in High-K Lunar Basalts

Chlorine is an unusual isotopic system, being essentially unfractionated ((delta)Cl-37 approximately 0 per mille ) between bulk terrestrial samples and chondritic meteorites and yet showing large variations in lunar (approximately -4 to +81 per mille), martian, and vestan (HED) samples. Among lunar samples, the volatile-bearing mineral apatite (Ca5(PO4)3[F,Cl,OH]) has been studied for volatiles in K-, REE-, and P (KREEP), very high potassium (VHK), low-Ti and high-Ti basalts, as well as samples from the lunar highlands. These studies revealed a positive correlation between in-situ (delta)Cl-37 measurements and bulk incompatible trace elements (ITEs) and ratios. Such trends were interpreted to originate from Cl isotopic fractionation during the degassing of metal chlorides during or shortly after the differentiation of the Moon via a magma ocean. In this study, we investigate the volatile inventories of a group of samples for which new-era volatile data have yet to be reported - the high-K (greater than 2000 ppm bulk K2O), high-Ti, trace element-rich mare basalts. We used isotope imaging on the Cameca NanoSIMS 50L at JSC to obtain the Cl isotopic composition [((Cl-37/(35)Clsample/C-37l/(35)Clstandard)-1)×1000, to get a value in per thousand (per mille)] which ranges from approximately -2.7 +/- 2 per mille to +16.1 +/- 2 per mille (2sigma), as well as volatile abundances (F & Cl) of apatite in samples 10017, 10024 & 10049. Simply following prior models, as lunar rocks with high bulk-rock abundances of ITEs we might expect the high-K, high-Ti basalts to contain apatite characterized by heavily fractionated (delta)Cl-37 values, i.e., Cl obtained from mixing between unfractionated mantle Cl (approximately 0 per mille) and the urKREEP reservoir (possibly fractionated to greater than +25 per mille.). However, the data obtained for the studied samples do not conform to either the early degassing or mixing models. Existing petrogentic models for the origin of the high-K, high-Ti basalts do not include urKREEP assimilation into their LMO cumulate sources. Therefore, Cl in these basalts either originated from source region heterogeneity or through assimilation or metasomatism by volatile and incompatible trace element rich materials. The new data presented here could provide evidence for the existence of region(s) in the lunar interior that are ITE-enriched and contain Cl that does not share isotopic affinities with lunar urKREEP, possibly representing the composition of the purported 'neuKREEP'.

Barnes, Jessica J.↗

Satellite Gravimetry Applied to Drought Monitoring

Near-surface wetness conditions change rapidly with the weather, which limits their usefulness as drought indicators. Deeper stores of water, including root-zone soil wetness and groundwater, portend longer-term weather trends and climate variations, thus they are well suited for quantifying droughts. However, the existing in situ networks for monitoring these variables suffer from significant discontinuities (short records and spatial undersampling), as well as the inherent human and mechanical errors associated with the soil moisture and groundwater observation. Remote sensing is a promising alternative, but standard remote sensors, which measure various wavelengths of light emitted or reflected from Earth's surface and atmosphere, can only directly detect wetness conditions within the first few centimeters of the land s surface. Such sensors include the Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E) C-band passive microwave measurement system on the National Aeronautic and Space Administration's (NASA) Aqua satellite, and the combined active and passive L-band microwave system currently under development for NASA's planned Soil Moisture Active Passive (SMAP) satellite mission. These instruments are sensitive to water as deep as the top 2 cm and 5 cm of the soil column, respectively, with the specific depth depending on vegetation cover. Thermal infrared (TIR) imaging has been used to infer water stored in the full root zone, with limitations: auxiliary information including soil grain size is required, the TIR temperature versus soil water content curve becomes flat as wetness increases, and dense vegetation and cloud cover impede measurement. Numerical models of land surface hydrology are another potential solution, but the quality of output from such models is limited by errors in the input data and tradeoffs between model realism and computational efficiency. This chapter is divided into eight sections, the next of which describes the theory behind satellite gravimetry. Following that is a summary of the GRACE mission and how hydrological information is gleaned from its gravity products. The fourth section provides examples of hydrological science enabled by GRACE. The fifth and sixth sections list the challenging aspects of GRACE derived hydrology data and how they are being overcome, including the use of data assimilation. The seventh section describes recent progress in applying GRACE for drought monitoring, including the development of new soil moisture and drought indicator products, and that is followed by a discussion of future prospects in satellite gravimetry based drought monitoring.

Rodell, Matthew↗

Offshore wind energy forecasting sensitivity to sea surface temperature input in the Mid-Atlantic

Abstract. As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25∘ × 0.25∘ resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054∘ × 0.054∘ resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02∘ × 0.02∘ resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 m) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (< 30 h each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10−1 or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

17 WIND ENERGY↗

An Integrated Approach to Weather Radar Calibration and Monitoring Using Ground Clutter and Satellite Comparisons

Abstract The stability and accuracy of weather radar reflectivity calibration are imperative for quantitative applications, such as rainfall estimation, severe weather monitoring and nowcasting, and assimilation in numerical weather prediction models. Various radar calibration and monitoring techniques have been developed, but only recently have integrated approaches been proposed, that is, using different calibration techniques in combination. In this paper the following three techniques are used: 1) ground clutter monitoring, 2) comparisons with spaceborne radars, and 3) the self-consistency of polarimetric variables. These techniques are applied to a C-band polarimetric radar (CPOL) located in the Australian tropics since 1998. The ground clutter monitoring technique is applied to each radar volumetric scan and provides a means to reliably detect changes in calibration, relative to a baseline. It is remarkably stable to within a standard deviation of 0.1 dB. To obtain an absolute calibration value, CPOL observations are compared to spaceborne radars on board TRMM and GPM using a volume-matching technique. Using an iterative procedure and stable calibration periods identified by the ground echoes technique, we improve the accuracy of this technique to about 1 dB. Finally, we review the self-consistency technique and constrain its assumptions using results from the hybrid TRMM–GPM and ground echo technique. Small changes in the self-consistency parameterization can lead to 5 dB of variation in the reflectivity calibration. We find that the drop-shape model of Brandes et al. with a standard deviation of the canting angle of 12° best matches our dataset.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the 7-km GEOS-5 Nature Run

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

GEOS-5↗

A Systems Biology Approach to Energy Flow in H 2 Producing Microbial Communities: Microbial Mats

This report describes the results of work conducted at NASA Ames Research Center, in collaboration with investigators at Lawrence Livermore National Laboratory, during the time period 12/15/2014 to 12/15/2017. The work was an extension of a collaboration with LLNL that began in 2008. Work this reporting period focused primarily on tasks related to microbial mats from both the Elkhorn Slough and Baja California field sites, including the analysis and publication of data produced prior to 12/15/2014. During this reporting period, in addition to the microbial mat work, we also conducted a series of experiments with a particular algal-bacterial association (C. sorokiniana and A. brasilense). A number of aspects of this well-characterized algal-bacterial association enabled us not only to further our understanding of interactions between these particular organisms, but to apply that knowledge and techniques to our collection of microorganisms isolated from microbial mats. Microbial mat experiments focused both on intact mat communities and on cultures of organisms isolated from the mats (the Microbial Mat Living Library). In a series of co-culture experiments, we documented both positive and negative interactions amongst heterotrophic and autotrophic isolates of microbial mat microorganisms. One particularly important and novel mat cyanobacterium (ESFC-1) was studied in detail. We developed techniques to pursue an investigation of functional redundancy (the idea that closely related organisms perform identical roles in a community) in cyanobacterial autotrophs in the Elkhorn Slough microbial mat. In those experiments, we documented numerous instances of increases in growth rate between co-cultures of two organisms (relative to their growth rate in isolation) performing what would seem to be identical functional roles. Simplified microbial mat communities consisting of organisms from the Microbial Mat Living Library were constructed and studied in laboratory experiments. The “constructed mats” were grown in standard laboratory glassware and plasticware, but cultures were also combined with artificial materials in an effort to re-create the three-dimensional structure of naturally occurring microbial mats. These constructed mats re-created a number of functions attributable to natural microbial mats (N-fixation, H2-production) and simplified their study using molecular ecological techniques due to the smaller number of microbial “species” present. Cyanobacteria which migrated in response to vertical variations in light intensity (irradiance) obtained characteristic positions in constructed mats that were similar to their positions in natural mats. A detailed investigation of nitrogen cycling in both the Baja California and Elkhorn Slough microbial mats was initiated during this reporting period. The nitrogen cycle in the microbial mats has not been studied in depth, despite the fact that nitrogen usually limits growth in marine environments. The organisms conducting nitrogen cycling in both microbial mats were studied using the genomes available from earlier work on these mats, as well as with a targeted approach using molecular ecological methods developed here. Quantitative measurements of rates of nitrogen transformations were made using stable isotopically labeled nitrogen. The processes quantified and examined using molecular ecological methods included: nitrogen fixation, ammonification, ammonium and nitrate assimilation, denitrification, anaerobic oxidation of ammonium, and dissimilatory nitrate reduction to ammonium. In a series of experiments designed to examine a potential application for microbial mats, we demonstrated that both Elkhorn Slough and Baja California microbial mats were very effective in removing large quantities of both nitrate and ammonium from waste streams containing these compounds at levels where they would be considered to be pollutants. Considerable effort was expended during this reporting period to analyze, and prepare for publication, metagenomes and metatranscriptomes collected from two experimental manipulations of microbial mats. The first, a diel study, was conducted in November 2011, while the second, an investigation of the effects of sulfate on community composition, was conducted during this reporting period. For the earlier study, reference-based and reference-free methods were used to assess microbial ecology and genetic partitioning in these complex microbial systems, and that work is available as a pre-print. Four metagenomes from the same experiment were published. In the second, longer (nearly 200 day) manipulation, rRNA iTag libraries for cDNA and DNA were analyzed to determine the effects of sulfate concentration on microbial mat community composition over the course of the experiment. Work on the algal-bacterial association (Chlorella sorokiniana and Azospirillum brasilense) focused initially on demonstrating the transfer of carbon and nitrogen between the partners, using stable isotopic labeling and the LLNL NanoSIMS. Subsequently, the potential role of diazotrophy (nitrogen fixation) by the bacterial partner in supporting the growth of the association was assessed in a series of experiments. It appears to be the case that, although the bacterial partner demonstrably increases the fitness of the algae under nitrogen limitation, and cells of C. sorokiniana receive nitrogen originally fixed from atmospheric nitrogen by A. brasilense, the amount of nitrogen obtained by diazotrophy is too small to explain the increase in fitness.

08 HYDROGEN↗

Thermodynamic constraints on the textural evolution of eucrite, EET 90020

Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.

Eucrites↗

Data from: "Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine"

This archive contains data used to support conclusions drawn in “Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine”, by Reinhardt et al., 2011. Data were collected over one summer season in plots within the Alpine Treeline Warming Experiment (ATWE), before climate manipulations began. The experiment was located on Niwot Ridge, in the Front Range of the Colorado Rocky Mountains. This data package includes five comma-separated-values (.csv) files, five Microsoft Excel (.xlsx) files, one .pdf file, and two types of geospatial files: keyhole markup language (.kml), and ESRI shapefiles (.shp). .csv files can be opened using any simple text-editing software (such as Notepad and TextEdit), R, and Microsoft Excel. .xlsx files can only be opened using Microsoft Excel. The .pdf file can be opened using Adobe Acrobat Reader or any other compatible file viewing software. The .kml file can be opened using Google Earth and Google Maps, and shapefiles can be opened using any software compatible with the file type, such as ESRI’s ArcGIS suite and QGIS.Data archived contain gas exchange and plant physiology measurements, non-structural carbohydrate data, among others. Geospatial files are also provided for additional locational context. The files and their contents in this data package are summarized under "Data Summary" in the included Data User's Guide. All files (excluding geospatial) are available in both Microsoft Excel and in .csv format, and are indicated in the Data Summary list as well.-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Climate change is predicted to cause upward shifts in forest tree distributions, which will require seedling recruitment beyond current forest boundaries. However, predicting the likelihood of successful plant establishment beyond current species’ ranges under changing climate is complicated by the interaction of genetic and environmental controls on seedling establishment. To determine how genetics and climate may interact to affect seedling establishment, we transplanted recently germinated seedlings from high- and low-elevation provenances (HI and LO, respectively) of Pinus flexilis in common gardens arrayed along an elevation and canopy gradient from subalpine forest into the alpine zone and examined differences in physiology and morphology between provenances and among sites. Plant dry mass, projected leaf area and shoot:root ratios were 12–40% greater in LO compared with HI seedlings at each elevation. There were no significant changes in these variables among sites except for decreased dry mass of LO seedlings in the alpine site. Photosynthesis, carbon balance (photosynthesis/respiration) and conductance increased >2× with elevation for both provenances, and were 35–77% greater in LO seedlings compared with HI seedlings. There were no differences in dark-adapted chlorophyll fluorescence (Fv/Fm) among sites or between provenances. Our results suggest that for P. flexilis seedlings, provenances selected for above-ground growth may outperform those selected for stress resistance in the absence of harsh climatic conditions, even well above the species’ range limits in the alpine zone. This indicates that forest genetics may be important to understanding and managing species’ range adjustments due to climate change.

54 ENVIRONMENTAL SCIENCES↗