Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 1,333 records · Page 74

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Second ARM Aerosol Chemical Speciation Monitor Users’ Meeting Report

The aerosol chemical speciation monitor (ACSM) was developed to adapt the technology of the aerosol mass spectrometer (AMS) to routine, long-term, standalone monitoring. The calculation of particulate mass concentration from ACSM data requires the measurement of the response of the instrument to aerosol of specific size and composition as well as assumptions about the instrument response based on laboratory measurements and field experience acquired over more than two decades of operation of AMS and a decade of operation of the ACSM. Three parameters in the concentration calculations that are particularly important are the NO3 response factor (RFNO3), relative ionization efficiency (RIE), and the collection efficiency (CE). The values of RFNO3 and RIEs are determined from calibration, however the jump calibration method previously used in calibration can result in errors in the RIE for sulfate. This has been corrected by implementing a continuous calibration method. The default collection efficiency is 0.5. This has been shown to result in mass loadings that do not agree with mass determinations from other instruments because of effects of composition on the vaporization of the particles. The previous work of investigators addressing this issue is discussed. After preliminary work on ACSM and scanning mobility particle sizer (SMPS) data from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) observatory collected in late 2016 and 2017 produced a parameterization of composition dependent collection efficiency very different from the results of previous studies, SMPS data were examined and we determined that there was significant mass that the instrument did not capture because the particles with diameters larger than 465 nm are not counted by this instrument. Data for the ultra-high-sensitivity aerosol spectrometer (UHSAS), SMPS, and ACSM are available for nearly all of 2019. The data from UHSAS and SMPS collected in 2019 were compared. We found that the UHSAS data has particle counts and total volumes significantly less than measured by the condensation particle counter (CPC) and SMPS. The SMPS data were extended by fitting the average volume distribution with a log normal curve and using this relationship to estimate a mass value for the SMPS over extended diameter range. The extended SMPS mass values result in a CDCE parameterization that is in better agreement with the results of other investigators, but is still different from other formulations. The working group recommends that the ACSM data be processed with a collective efficiency (CE)=1, that this be documented clearly in the metadata, and the use of the default CE of 0.5 or a formulation of composition-dependent collection efficiencies (CDCE) chosen by the user should be implemented based on the ammonium nitrate mass fraction. This is clearly necessary for the wintertime SGP ACSM data because of the high nitrate concentrations.

47 OTHER INSTRUMENTATION↗

VIC-Global Parameter Dataset Sensitivity with the Variable Infiltration Capacity Model: Evaluating the importance of dynamic land surface parameters when using the VIC-Global parameter dataset

Accurate prediction of runoff is essential to water resources management, flood risk assessment, and ecosystem protection. However, many hydrological models still have relatively substantial limitations when representing the influence of land use and land cover (LULC) on runoff generation and routing. Changes in LULC, such as deforestation, urban expansion, agricultural intensification, and wetland loss, have been shown to alter the water balance at the land surface through fundamental hydrologic processes (e.g., interception, infiltration, evapotranspiration, and soil storage). However, it remains an open question what the exact magnitude and timing of these impacts are for the spatial and temporal scales commonly used in engineering applications. In this analysis we focus on one aspect of recent LULC change for assessing human impacts, which is urbanization. Specifically we seek to determine the impacts of urbanization on the magnitude and timing of surface runoff and baseflow in HUC-12 basins in Clark County, Nevada which has experienced rapid urbanization. We use the Variable Infiltration Capacity (VIC) hydrology model with a widely used off-the-shelf dataset of land surface parameters, VIC-Global, both of which have been commonly used in the past for water and energy balance modeling for large scale hydrologic studies. We examine two scenarios where the first scenario removes all urbanized land cover and parameterizes those areas of the basins as barren or open shrubland. The second scenario tests the opposite case where all areas of the basins are classified as urban regardless of their present classification. The results from the VIC model show there is a low sensitivity for daily surface runoff between scenarios. The daily baseflow values indicate similar low sensitivity to the classification change during specific periods, but then have substantial differences during other period when large precipitation events are occurring. This is likely due to the assumed parameter values for the urban land cover classification made by the VIC-Global dataset. Using a static land cover parameterization is reasonable for large domain hydrology models that are being used for near-term planning horizons (<30 years). However, longer planning horizons where feedbacks between the atmosphere and land surface are important, especially in transient climate situations, considerations for how to update land surface parameters should be incorporated.

42 ENGINEERING↗

WFIP3

The Wind Forecasting Improvement Project 3 (WFIP-3) is the first offshore-based wind resource characterization project within the WFIP construct, funded by the U.S. Department of Energy. WFIP-3 will provide a unique field study that will deliver the comprehensive suite of data needed to inform a series of modeling efforts that will develop and evaluate parameterization schemes suited to offshore environments and improved industry-targeted applications. The field study has two goals: (1) detailed sampling of the vertical structure of the Marine Atmospheric Boundary Layer (MABL) at key observational areas, creating a rich dataset that will be used to refine and validate parameterization schemes, and (2) wide-area sampling of the MABL to create a multi-scale array of observations informing and guiding models of resource characterization. We will deploy a multi-platform array of measurements that span the MABL and create a multi-scale observational array stretching south from Marth’s Vineyard across the wind energy areas.

17 WIND ENERGY↗

Infinite quantum signal processing

Quantum signal processing (QSP) represents a real scalar polynomial of degree d using a product of unitary matrices of size 2 × 2 , parameterized by ( d + 1 ) real numbers called the phase factors. This innovative representation of polynomials has a wide range of applications in quantum computation. When the polynomial of interest is obtained by truncating an infinite polynomial series, a natural question is whether the phase factors have a well defined limit as the degree d → ∞ . While the phase factors are generally not unique, we find that there exists a consistent choice of parameterization so that the limit is well defined in the ℓ 1 space. This generalization of QSP, called the infinite quantum signal processing, can be used to represent a large class of non-polynomial functions. Our analysis reveals a surprising connection between the regularity of the target function and the decay properties of the phase factors. Our analysis also inspires a very simple and efficient algorithm to approximately compute the phase factors in the ℓ 1 space. The algorithm uses only double precision arithmetic operations, and provably converges when the ℓ 1 norm of the Chebyshev coefficients of the target function is upper bounded by a constant that is independent of d . This is also the first numerically stable algorithm for finding phase factors with provable performance guarantees in the limit d → ∞ .

Dong, Yulong [Department of Mathematics, Universit↗

Comparative Study and Validation of Photovoltaic Model Formulations for the IBPSA Modelica Library based on Rooftop Measurement Data

Domain-overarching system models are crucial to investigate sector coupling concepts. Specifically, the coupling of building and electrical energy systems becomes crucial to integrate renewable energy sources such as photovoltaic power systems (PV). For such interdisciplinary simulation models, Modelica is a suitable language. However, most open-source Modelica libraries are either domain-specific or lack simple-to-parameterize PV models. We close this gap by developing a PV model for the IBPSA Modelica library. The model comprises two I-V-characteristic models and three mounting-dependent approaches to calculate the cell temperature. The I-V-characteristic models follow a single- and two-diodes approach. This study uses measurement data from a rooftop PV system in Berlin, Germany, for validation. The focus lies on comparing the implemented single- and two-diodes approach. Results prove that both models accurately calculate the modules’ DC power output and cell temperature. The two-diodes approach slightly outperforms the single-diode one at the expense of a higher parameterization effort.

Maier, Laura↗

Trait-Based Representation of Biological Nitrification: Model Development, Testing, and Predicted Community Composition

Trait-based microbial models show clear promise as tools to represent the diversity and activity of microorganisms across ecosystem gradients. These models parameterize specific traits that determine the relative fitness of an “organism” in a given environment, and represent the complexity of biological systems across temporal and spatial scales. In this study we introduce a microbial community trait-based modeling framework (MicroTrait) focused on nitrification (MicroTrait-N) that represents the ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA) and nitrite-oxidizing bacteria (NOB) using traits related to enzyme kinetics and physiological properties.We used this model to predict nitrifier diversity, ammonia (NH3) oxidation rates, and nitrous oxide (N2O) production across pH, temperature, and substrate gradients. Predicted nitrifier diversity was predominantly determined by temperature and substrate availability, the latter was strongly influenced by pH. The model predicted that transient N2O production rates are maximized by a decoupling of the AOB and NOB communities, resulting in an accumulation and detoxification of nitrite to N2O by AOB. However, cumulative N2O production (over 6 month simulations) is maximized in a system where the relationship between AOB and NOB is maintained.When the reactions uncouple, the AOB become unstable and biomass declines rapidly, resulting in decreased NH3 oxidation and N2O production. We evaluated this model against site level chemical datasets from the interior of Alaska and accurately simulated NH3 oxidation rates and the relative ratio of AOA:AOB biomass. The predicted community structure and activity indicate (a) parameterization of a small number of traits may be sufficient to broadly characterize nitrifying community structure and (b) changing decadal trends in climate and edaphic conditions could impact nitrification rates in ways that are not captured by extant biogeochemical models.

59 BASIC BIOLOGICAL SCIENCES↗

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

trait-based modeling↗

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

Does Marine Surface Tension Have Global Biogeography? Addition for the OCEANFILMS Package

We apply principles of Gibbs phase plane chemistry across the entire ocean-atmosphere interface to investigate aerosol generation and geophysical transfer issues. Marine surface tension differences comprise a tangential pressure field controlling trace gas fluxes, primary organic inputs, and sea spray salt injections, in addition to heat and momentum fluxes. Mapping follows from the organic microlayer composition, now represented in ocean system models. Organic functional variations drive the microforcing, leading to (1) reduced turbulence and (by extension) laminar gas-energy diffusion; plus (2) altered bubble film mass emission into the boundary layer. Interfacial chemical behaviors are, therefore, closely reviewed as the background. We focus on phase transitions among two dimensional “solid, liquid, and gaseous” states serving as elasticity indicators. From the pool of dissolved organic carbon (DOC) only proteins and lipids appear to occupy significant atmospheric interfacial areas. The literature suggests albumin and stearic acid as the best proxies, and we distribute them through ecodynamic simulation. Consensus bulk distributions are obtained to control their adsorptive equilibria. We devise parameterizations for both the planar free energy and equation of state, relating excess coverage to the surface pressure and its modulus. Constant settings for the molecular surrogates are drawn from laboratory study and successfully reproduce surfactant solid-to-gas occurrence in compression experiments. Since DOC functionality measurements are rare, we group them into super-ecological province tables to verify aqueous concentration estimates. Outputs are then fed into a coverage, tension, elasticity code. The resulting two dimensional pressure contours cross a critical range for the regulation of precursor piston velocity, bubble breakage, and primary aerosol sources plus ripple damping. Concepts extend the water-air adsorption theory currently embodied in our OCEANFILMS aerosol emissions package, and the two approaches could be inserted into Earth System Models together. Uncertainties in the logic include kinetic and thermochemical factors operating at multiple scales.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Influence of CAM5 Aerosol Configuration on Simulated Tropical Cyclones in the North Atlantic

This study examines the influence of prescribed and prognostic aerosol model configurations on the formation of tropical cyclones (TCs) in the North Atlantic Ocean in Community Atmosphere Model version 5 (CAM5). The impact of aerosol parameterization is examined by investigating storm track density, genesis density, potential intensity, and genesis potential index. This work shows that both CAM5 configurations simulate reduced storm frequency when compared to observations and that differences in TC climatology between the model configurations can be explained by differences in the large-scale environment. The analysis shows that simulation with the prognostic aerosol parameterization scheme reasonably captures the observed interannual variability in tropical cyclones and aerosols (i.e., dust) in the North Atlantic, while simulation with the prescribed configuration (climatology) is less favorable. The correlation between dust and TCs in observations (i.e., reanalysis and satellite datasets) is shown to be negative, and this relationship was also found for the prognostic aerosol configuration despite an overall decrease in the frequency of TCs. This indicates that, to accurately replicate certain aspects of TC interannual variability, the aerosol configuration within CAM5 needs to account for the appropriate dust variability.

Huff, J. Jacob A. (ORCID:0000000329143817)↗

Application of an Online-Coupled Regional Climate Model, WRF-CAM5, over East Asia for Examination of Ice Nucleation Schemes: Part I. Comprehensive Model Evaluation and Trend Analysis for 2006 and 2011

Online-coupled climate and chemistry models are necessary to realistically represent the interactions between climate variables and chemical species and accurately simulate aerosol direct and indirect effects on cloud, precipitation, and radiation. In this Part I of a two-part paper, simulations from the Weather Research and Forecasting model coupled with the physics package of Community Atmosphere Model (WRF-CAM5) are conducted with the default heterogeneous ice nucleation parameterization over East Asia for two full years: 2006 and 2011. A comprehensive model evaluation is performed using satellite and surface observations. The model shows an overall acceptable performance for major meteorological variables at the surface and in the boundary layer, as well as column variables (e.g., precipitation, cloud fraction, precipitating water vapor, downward longwave and shortwave radiation). Moderate to large biases exist for cloud condensation nuclei over oceanic areas, cloud variables (e.g., cloud droplet number concentration, cloud liquid and ice water paths, cloud optical depth, longwave and shortwave cloud forcing). These biases indicate a need to improve the model treatments for cloud processes, especially cloud droplets and ice nucleation, as well as to reduce uncertainty in the satellite retrievals. The model simulates well the column abundances of chemical species except for column SO 2 but relatively poor for surface concentrations of several species such as CO, NO 2 , SO 2 , PM2.5, and PM10. Several reasons could contribute to the underestimation of major chemical species in East Asia including underestimations of anthropogenic emissions and natural dust emissions, uncertainties in the spatial and vertical distributions of the anthropogenic emissions, as well as biases in meteorological, radiative, and cloud predictions. Despite moderate to large biases in the chemical predictions, the model performance is generally consistent with or even better than that reported for East Asia with only a few exceptions. The model generally reproduces the observed seasonal variations and the difference between 2006 and 2011 for most variables or chemical species. Overall, these results demonstrate promising skills of WRF-CAM5 for long-term simulations at a regional scale and suggest several areas of potential improvements.

54 ENVIRONMENTAL SCIENCES↗

Modeling Sea Ice Effects for Wave Energy Resource Assessments

Wave-generated power has potential as a valuable coastal resource, but the wave climate needs to be mapped for feasibility before wave energy converters are installed. Numerical models are used for wave resource assessments to quantify the amount of available power and its seasonality. Alaska is the U.S. state with the longest coastline and has extensive wave resources, but it is affected by seasonal sea ice that dampens the wave energy and the full extent of this dampening is unknown. To accurately characterize the wave resource in regions that experience seasonal sea ice, coastal wave models must account for these effects. The aim of this study is to determine how the dampening effects of sea ice change wave energy resource assessments in the nearshore. Here, we show that by combining high-resolution sea ice imagery with a sea ice/wave dampening parameterization in an unstructured grid, the Simulating Waves Nearshore (SWAN) model improves wave height predictions and demonstrates the extent to which wave power decreases when sea ice is present. The sea ice parametrization decreases the bias and root mean square errors of wave height comparisons with two wave buoys and predicts a decrease in the wave power of up to 100 kW/m in areas around Prince William Sound, Alaska. The magnitude of the improvement of the model/buoy comparison depends on the coefficients used to parameterize the wave–ice interaction.

16 TIDAL AND WAVE POWER↗

The Impact of Model Resolution on the Southern Hemisphere in Community Climate System Model Version 4 Idealized Climate Simulations

Model resolution plays a large role in accurately simulating the Southern Hemisphere circulation in both the ocean and atmosphere. Resolving the mesoscale field is important as it has been shown to have a significant impact on the large-scale climate in eddy-rich regions, which are regions of large CO2 absorption. The presence of ocean and atmospheric mesoscale features can affect sea surface temperatures, the strength and location of storm tracks, and many other air-sea processes. Additionally, with an improvement in resolution, the eddy kinetic energy in the ocean can be expected to change considerably. The significance model resolution has on the Southern Hemisphere is examined using the Community Climate System Model, Version 4, eddy-parameterizing and eddy-resolving simulations. The CO2 concentrations and ozone levels are specified independently to better understand how the mesoscale field responds to extreme changes in external forcing and the resulting climate impacts. Overall, in the eddy-parameterizing simulations, the ozone forcing is found to be more important than the changes in CO2 concentrations. However, in the case of the eddy-resolving simulations, the CO2 concentrations are found to be more dominant, especially in eddy-rich regions. These results demonstrate the need for an increase in model resolution for climate prediction.

Daher, Houraa (ORCID:0000000200177346)↗

Assessment of canopy chlorophyll content retrieval in maize and soybean: Implications of hysteresis on the development of generic algorithms

Here, canopy chlorophyll content (Chl) closely relates to plant photosynthetic capacity, nitrogen status and productivity. The goal of this study is to develop remote sensing techniques for accurate estimation of canopy Chl during the entire growing season without re-parameterization of algorithms for two contrasting crop species, maize and soybean. These two crops represent different biochemical mechanisms of photosynthesis, leaf structure and canopy architecture. The relationships between canopy Chl and reflectance, collected at close range and resampled to bands of the Multi Spectral Instrument (MSI) aboard Sentinel-2, were analyzed in samples taken across the entirety of the growing seasons in three irrigated and rainfed sites located in eastern Nebraska between 2001 and 2005. Crop phenology was a factor strongly influencing the reflectance of both maize and soybean. Substantial hysteresis of the reflectance vs. canopy Chl relationship existed between the vegetative and reproductive stages. The effect of the hysteresis on vegetation indices (VI), applied for canopy Chl estimation, depended on the bands used and their formulation. The hysteresis greatly affected the accuracy of canopy Chl estimation by widely-used VIs with near infrared (NIR) and red reflectance (e.g., normalized difference vegetation index (NDVI), enhanced vegetation index (EVI) and simple ratio (SR)). VIs that use red edge and NIR bands (e.g., red edge chlorophyll index (CIred edge), red edge NDVI and the MERIS terrestrial chlorophyll index (MTCI)) were minimally affected by crop phenology (i.e., they exhibited little hysteresis) and were able to accurately estimate canopy Chl in two crops without algorithm re-parameterization and, thus, were found to be the best candidates for generic algorithms to estimate crop Chl using the surface reflectance products of MSI Sentinel-2.

59 BASIC BIOLOGICAL SCIENCES↗

Connecting CO intensity mapping to molecular gas and star formation in the epoch of galaxy assembly

Intensity mapping, which images a single spectral line from unresolved galaxies across cosmological volumes, is a promising technique for probing the early universe. Here we present predictions for the intensity map and power spectrum of the CO(1–0) line from galaxies at z∼2.4–2.8, based on a parameterized model for the galaxy–halo connection, and demonstrate the extent to which properties of high-redshift galaxies can be directly inferred from such observations. We find that our fiducial prediction should be detectable by a realistic experiment. Motivated by significant modeling uncertainties, we demonstrate the effect on the power spectrum of varying each parameter in our model. Using simulated observations, we infer constraints on our model parameter space with an MCMC procedure, and show corresponding constraints on the L{sub IR}–L{sub CO} relation and the CO luminosity function. These constraints would be complementary to current high-redshift galaxy observations, which can detect the brightest galaxies but not complete samples from the faint end of the luminosity function. By probing these populations in aggregate, CO intensity mapping could be a valuable tool for probing molecular gas and its relation to star formation in high-redshift galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

PROBING THE EXPANSION HISTORY OF THE UNIVERSE BY MODEL-INDEPENDENT RECONSTRUCTION FROM SUPERNOVAE AND GAMMA-RAY BURST MEASUREMENTS

To probe the late evolution history of the universe, we adopt two kinds of optimal basis systems. One of them is constructed by performing the principle component analysis, and the other is built by taking the multidimensional scaling approach. Cosmological observables such as the luminosity distance can be decomposed into these basis systems. These basis systems are optimized for different kinds of cosmological models that are based on different physical assumptions, even for a mixture model of them. Therefore, the so-called feature space that is projected from the basis systems is cosmological model independent, and it provides a parameterization for studying and reconstructing the Hubble expansion rate from the supernova luminosity distance and even gamma-ray burst (GRB) data with self-calibration. The circular problem when using GRBs as cosmological candles is naturally eliminated in this procedure. By using the Levenberg–Marquardt technique and the Markov Chain Monte Carlo method, we perform an observational constraint on this kind of parameterization. The data we used include the “joint light-curve analysis” data set that consists of 740 Type Ia supernovae and 109 long GRBs with the well-known Amati relation.

79 ASTRONOMY AND ASTROPHYSICS↗