Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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

Techno-Economic Analysis and Life Cycle Assessment for Guiding Bioenergy Pathway Development

The Center for Bioenergy Innovation (CBI) is a bioenergy research center recently launched by the U.S. Department of Energy Office of Science's Biological and Environmental Research program. CBI's objective is to enable innovations across the bioenergy supply chain through development of poplar and switchgrass phenotypes, microorganisms for consolidated bio-processing (CBP), in which biomass fermentation and conversion is accomplished in a single unit operation, and lignin valorization processes. In addition to basic science, CBI's activities involve the application of techno-economic analysis (TEA) and life cycle assessment (LCA) to the bioenergy pathways under study. TEA and LCA is used to quantify the impacts that research and development breakthroughs have on the economic viability and sustainability of a bioenergy pathway, and to identify drivers for reducing costs and improving sustainability through targeted research efforts. Combining TEA and LCA allows for understanding the economic and sustainability benefits and trade-offs of these integrated bioenergy systems. The focus of this presentation is a parameterized, integrated TEA-LCA model, which will be discussed along with an application to bioenergy pathways currently under study. The model combines data and information developed by CBI researchers into a superstructure of potential bioenergy pathways and their life cycles. Model parameters include feedstock species and phenotype, farming, harvesting, storage and transportation practices, and the microbes used for biomass conversion and lignin valorization, among others. The parameterized TEA-LCA model allows for the rapid quantification of pathway economic and sustainability metrics such as minimum fuel selling price, payback period, fossil fuel consumption and net greenhouse gas emissions. To demonstrate the model, economic and sustainability metrics for a bioenergy pathway utilizing CBP will be quantified and compared to the same metrics for a conventional bioenergy pathway. Results of the comparison will be used to identify and prioritize high-impact research areas within the CBP pathway. fermentation and conversion is accomplished in a single unit operation, and lignin valorization processes. In addition to basic science, CBI's activities involve the application of techno-economic analysis (TEA) and life cycle assessment (LCA) to the bioenergy pathways under study. TEA and LCA is used to quantify the impacts that research and development breakthroughs have on the economic viability and sustainability of a bioenergy pathway, and to identify drivers for reducing costs and improving sustainability through targeted research efforts. Combining TEA and LCA allows for understanding the economic and sustainability benefits and trade-offs of these integrated bioenergy systems. The focus of this presentation is a parameterized, integrated TEA-LCA model, which will be discussed along with an application to bioenergy pathways currently under study. The model combines data and information developed by CBI researchers into a superstructure of potential bioenergy pathways and their life cycles. Model parameters include feedstock species and phenotype, farming, harvesting, storage and transportation practices, and the microbes used for biomass conversion and lignin valorization, among others. The parameterized TEA-LCA model allows for the rapid quantification of pathway economic and sustainability metrics such as minimum fuel selling price, payback period, fossil fuel consumption and net greenhouse gas emissions. To demonstrate the model, economic and sustainability metrics for a bioenergy pathway utilizing CBP will be quantified and compared to the same metrics for a conventional bioenergy pathway. Results of the comparison will be used to identify and prioritize high-impact research areas within the CBP pathway.

09 BIOMASS FUELS↗

The role of polar regions in global climate, and a new parameterization of global heat transport

The effects of the transport of heat between polar regions and other latitudes on climate sensitivity and stability are examined within the framework of simple energy balance models. New heat transport parameterizations adjust radiative equilibrium distributions of temperature with latitude on the basis of Hadley cells and baroclinically unstable eddies; including the effects of static stability changes with latitude eliminates the possible error in estimating the pole-equator temperature difference. It is found that climate sensitivity and stability for the new transport parameterizations can differ from other models, and is capable of simulating the sensitivity required by existing climate data.

Lindzen, R. S.↗

Modeling Interfaces to Support Low-Level Waste Disposal System Performance Assessments - 20366

In low-level waste (LLW) disposal sites, interfaces between cementitious materials, used as waste forms and/or engineered barriers, and the surrounding soil or backfill material are often encountered. Reactions across these interfaces may lead to chemical and structural alteration of the cementitious and backfill materials that may ultimately affect long-term performance. For example, the ingress of carbon dioxide from soil gas into a waste tank concrete shell at the Hanford Site or Savannah River Site (SRS) or a low-level waste disposal vault at SRS leads to carbonation of cement hydration products (e.g. portlandite, calcium-silicate hydrate or C-S-H, and ettringite). The result of carbonation is a decrease in pH in the cement paste portion of the tank wall or vault concrete structure that could ultimately lead to de-passivation of the embedded structural steel and, eventually, to cracking. Cracking can lead to increased ingress of water into the structure, corresponding increased release of constituents of concern, and to subsequent increased transport of these constituents to the surrounding environment. The mobility of trace constituents (e.g., radionuclides of concern) in waste forms may increase in response to changes in pH, pore structure, and mineralogical gradients. Depletion of portlandite results in subsequent decalcification of C-S-H that leads to changes in the cement strength and may lead to structural failure. At the boundary between waste forms and concrete vaults, the migration of sulfate ions from the salt waste form (e.g., the saltstone waste form used for LLW disposal at the Savannah River Site in Saltstone Disposal Units, or SDUs) into the barrier or vault concrete has been predicted to result in expansive mineral phase formation (or 'sulfate attack'). Expansive mineral formation could result in cracking and potential loss of structural integrity in the vault concrete. In the performance assessment (PA) for the SRS SDUs, prediction of both the carbonation and sulfate attack fronts are critical to assessing the long-term performance of these LLW disposal vaults. A general purpose geochemical reactive transport model has been developed in LeachXS/ORCHESTRA to evaluate typical interfaces for LLW disposal environments, including SDUs at SRS. The model may be used to predict the long-term performance of interfaces between cementitious materials and saltstone or backfill materials with respect to primary phases by considering mass transport, cement chemistry, geochemical speciation, and multi-ionic diffusion over interfaces between different materials. The present study compares the performance and interactions between vault concrete and saltstone waste form for a representative SDU scenario. The materials were characterized using pH-dependence leaching test (EPA 1313) and semi-dynamic transport tests (EPA 1315) with test results used to parameterize the interface model. The development of a credible inter model provides an important basis for supporting future PAs where carbonation and sulfate attack are important aging and degradation mechanisms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Diagnosis of the Marine Low Cloud Simulation in the NCAR Community Earth System Model (CESM) and the NCEP Global Forecast System (GFS)-Modular Ocean Model v4 (MOM4) coupled model

We present a diagnostic analysis of the marine low cloud climatology simulated by two state-of-the-art coupled atmosphere-ocean models: the NCAR Community Earth System Model (CESM) and the NCEP Global Forecasting System (GFS). In both models, the shallow convection and boundary layer turbulence parameterizations have been recently updated: both models now use a mass-flux scheme for the parameterization of shallow convection, and a turbulence parameterization capable of handling Stratocumulus (Sc)-topped Planetary Boundary Layers (PBLs). For shallow convection, both models employ a convective trigger function based on the concept of convective inhibition and both include explicit convective overshooting/penetrative entrainment formulation. For Sc-topped PBL, both models treat explicitly turbulence mixing and cloud-top entrainment driven by cloud-top radiative cooling. Our focus is on the climatological transition from Sc to shallow Cumulus (Cu)-topped PBL in the subtropical eastern oceans. We show that in the CESM the coastal Sc-topped PBLs in the subtropical Eastern Pacific are well-simulated but the climatological transition from Sc to shallow Cu is too abrupt and happens too close to the coast. By contrast, in the GFS coupled simulation the coastal Sc amount and PBL depth are severely underestimated while the transition from Sc to shallow Cu is ³delayed² and offshore Sc cover is too extensive in the subtropical Eastern Pacific. We discuss the possible connections between such differences in the simulations and differences in the parameterizations of shallow convection and boundary layer turbulence in the two models.

Xiao, Heng↗

Constraints on the Profiles of Total Water PDF in AGCMs from AIRS and a High-Resolution Model

Atmospheric general circulation model (AGCM) cloud parameterizations generally include an assumption about the subgrid-scale probability distribution function (PDF) of total water and its vertical profile. In the present study, the Atmospheric Infrared Sounder (AIRS) monthly-mean cloud amount and relative humidity fields are used to compute a proxy for the second moment of an AGCM total water PDF called the RH01 diagnostic, which is the AIRS mean relative humidity for cloud fractions of 0.1 or less. The dependence of the second moment on horizontal grid resolution is analyzed using results from a high-resolution global model simulation.The AIRS-derived RH01 diagnostic is generally larger near the surface than aloft, indicating a narrower PDF near the surface, and varies with the type of underlying surface. High-resolution model results show that the vertical structure of profiles of the AGCM PDF second moment is unchanged as the grid resolution changes from 200 to 100 to 50 km, and that the second-moment profiles shift toward higher values with decreasing grid spacing.Several Goddard Earth Observing System, version 5 (GEOS-5), AGCM simulations were performed with several choices for the profile of the PDF second moment. The resulting cloud and relative humidity fields were shown to be quite sensitive to the prescribed profile, and the use of a profile based on the AIRS-derived proxy results in improvements relative to observational estimates. The AIRS-guided total water PDF profiles, including their dependence on underlying surface type and on horizontal resolution, have been implemented in the version of the GEOS-5 AGCM used for publicly released simulations.

Molod, Andrea↗

Experimental and Computational Interrogation of Fast SCR Mechanism and Active Sites on H-Form SSZ-13

Experiment and density functional theory (DFT) models are combined to develop a unified, quantitative model of the mechanism and kinetics of fast selective catalytic reduction (SCR) of NO/NO2 mixtures over H-SSZ-13 zeolite. Rates, rate orders, and apparent activation energies collected under differential conditions reveal two distinct kinetic regimes. First-principles thermodynamics simulations are used to determine the relative coverages of free Brønsted sites, chemisorbed NH4+ and physisorbed NH3 as a function of reaction conditions. First-principles metadynamics calculations show that all three sites can contribute to the rate-limiting N-N bond forming step in fast SCR. The results are used to parameterize a kinetic model that encompasses the full range of reaction conditions and recovers observed rate orders and apparent activation energies. Observed kinetic regimes are related to changes in most-abundant surface intermediates. Financial support was provided by the National Science Foundation GAOLI program under award number 1258690-CBET. We thank the Center for Research Computing at Notre

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental evidence for dynamic friction on rock fractures from frequency-dependent nonlinear hysteresis and harmonic generation: Low-Frequency Torsional Dynamic Friction

Frictional properties affect the propagation of high-amplitude seismic waves across rock fractures and faults. Laboratory evidence suggests that these properties can be measured in active seismic surveys, potentially offering a route to characterizing friction in situ. Here, we present experimental results from a subresonance torsional modulus and attenuation apparatus that utilizes micron-scale sinusoidal oscillations to probe the nonlinear stress-strain relation at a range of strain amplitudes and rates. Nonlinear effects are further quantified using harmonic distortion; however, time series data best illuminate underlying physical processes. The low-frequency stress-strain hysteretic loops show stiffening at the sinusoid's static ends, but stiffening is reduced above a threshold frequency. This shape is determined by harmonic generation in the strain; the stress signal has no harmonics, confirming that the fractured sample is the source of the nonlinearity. These qualitative observations suggest the presence of rate-dependent friction and are consistent between fractures in three different rock types. We propose that static friction at the low strain rate part of the cycle, when given sufficient “healing” time at low oscillation frequencies, causes this stiffening cusp shape in the hysteresis loop. While rate-and-state friction is commonly used to represent dynamic friction, it cannot capture static friction or negative slip velocities. So we implement another dynamic friction model, based on the work of Dahl, which describes this process and produces similar results. Since the two models have a similar form, parameterizations of field data could constraint fault model inputs, such as specific location velocity strengthening or weakening properties.

58 GEOSCIENCES↗

Can mesoscale models capture the effect from cluster wakes offshore?

Long wakes from offshore wind turbine clusters can extend tens of kilometers downstream, affecting the wind resource of a large area. Given the ability of mesoscale numerical weather prediction models to capture important atmospheric phenomena and mechanisms relevant to wake evolution, they are often used to simulate wakes behind large wind turbine clusters and their impact over a wider region. Yet, uncertainty persists regarding the accuracy of representing cluster wakes via mesoscale models and their wind turbine parameterizations. Here, we evaluate the accuracy of the Fitch wind farm parameterization in the Weather Research and Forecasting model in capturing cluster-wake effects using two different options to represent turbulent mixing in the planetary boundary layer. To this end, we compare operational data from an offshore wind farm in the North Sea that is fully or partially waked by an upstream array against high-resolution mesoscale simulations. In general, we find that mesoscale models accurately represent the effect of cluster wakes on front-row turbines of a downstream wind farm. However, the same models may not accurately capture cluster-wake effects on an entire downstream wind farm, due to misrepresenting internal-wake effects.

17 WIND ENERGY↗

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗

Theoretical and methodological challenges in hierarchical Bayesian inference for model-form uncertainty

This report describes challenges associated with the hierarchical Bayesian approach to inform model-form uncertainty (MFU) representations, which are parameterized modifications to a mathematical models’ governing equations to express uncertainty in form of the equations. To inform model-form uncertainties, hierarchical Bayesian inference is often employed. Here, the MFU parameters are distributed parametrically, and the hyperparameters of the parametric distribution are informed through Bayesian inference, with the aim of determining the MFU parameter distribution that best agrees with calibration data. In practice, however, we have found the hierarchical Bayesian approach falls short of this aim. We discuss theoretical and methodological challenges of the approach, and we present several numerical demonstrations of these challenges. To conclude, we suggest promising alternative approaches for future investigation.

97 MATHEMATICS AND COMPUTING↗

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗

An Extended Objective Evaluation of the 29-km Eta Model for Weather Support to the United States Space Program

This report describes the Applied Meteorology Unit's objective verification of the National Centers for Environmental Prediction 29-km eta model during separate warm and cool season periods from May 1996 through January 1998. The verification of surface and upper-air point forecasts was performed at three selected stations important for 45th Weather Squadron, Spaceflight Meteorology Group, and National Weather Service, Melbourne operational weather concerns. The statistical evaluation identified model biases that may result from inadequate parameterization of physical processes. Since model biases are relatively small compared to the random error component, most of the total model error results from day-to-day variability in the forecasts and/or observations. To some extent, these nonsystematic errors reflect the variability in point observations that sample spatial and temporal scales of atmospheric phenomena that cannot be resolved by the model. On average, Meso-Eta point forecasts provide useful guidance for predicting the evolution of the larger scale environment. A more substantial challenge facing model users in real time is the discrimination of nonsystematic errors that tend to inflate the total forecast error. It is important that model users maintain awareness of ongoing model changes. Such changes are likely to modify the basic error characteristics, particularly near the surface.

Nutter, Paul↗

Performance of the Goddard Multiscale Modeling Framework with Goddard Ice Microphysical Schemes

The multiscale modeling framework (MMF), which replaces traditional cloud parameterizations with cloud-resolving models (CRMs) within a host atmospheric general circulation model (GCM), has become a new approach for climate modeling. The embedded CRMs make it possible to apply CRM-based cloud microphysics directly within a GCM. However, most such schemes have never been tested in a global environment for long-term climate simulation. The benefits of using an MMF to evaluate rigorously and improve microphysics schemes are here demonstrated. Four one-moment microphysical schemes are implemented into the Goddard MMF and their results validated against three CloudSat/CALIPSO cloud ice products and other satellite data. The new four-class (cloud ice, snow, graupel, and frozen drops/hail) ice scheme produces a better overall spatial distribution of cloud ice amount, total cloud fractions, net radiation, and total cloud radiative forcing than earlier three-class ice schemes, with biases within the observational uncertainties. Sensitivity experiments are conducted to examine the impact of recently upgraded microphysical processes on global hydrometeor distributions. Five processes dominate the global distributions of cloud ice and snow amount in long-term simulations: (1) allowing for ice supersaturation in the saturation adjustment, (2) three additional correction terms in the depositional growth of cloud ice to snow, (3) accounting for cloud ice fall speeds, (4) limiting cloud ice particle size, and (5) new size-mapping schemes for snow and graupel. Despite the cloud microphysics improvements, systematic errors associated with subgrid processes, cyclic lateral boundaries in the embedded CRMs, and momentum transport remain and will require future improvement.

snow↗

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (MEME-CCS)

This effort employs and adapts a rigorous multiscale electricity modeling platform at the National Renewable Energy Laboratory (NREL) to evaluate carbon capture and sequestration (CCS) and negative emissions technologies (NET) from the ARPA-E FLECCS program. NREL's modeling platform includes the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model, which projects future electricity generation mixes at sub-state-level resolution that are downscaled to the unit-level to enable hourly, zonal or nodal electricity production cost modeling in the PLEXOS model. The ReEDS-PLEXOS modeling suite is well-established for examining electric sector futures with high renewable energy penetrations, energy storage, electrification, and distributed generation. The key advancement proposed herein utilizes collaboration with CCS experts at the University of Wyoming (U.WY) and the FLECCS technology development teams to create innovative methods for representing CCS and NET in the ReEDS and PLEXOS models. Expanded technology options and new operational parameterizations are integrated into these models to allow an unprecedented combination of scope and resolution for exploring the future of CCS and NET. The resulting capabilities permit wide-ranging scenario analysis to assess CCS and NET deployment under alternative scenarios of CO2 prices and competitiveness of flexible CCS and NET technologies. We demonstrate sample deployment and operational outcomes to show how these models are being used to assess the future potential for FLECCS technologies and their impacts on the U.S. electricity system. These products will help an emerging CCS/NET industry in the United States by providing economics-driven guidance to technology developers while informing policy and investment decisions in the public and private sectors.

capacity expansion↗

Large-Eddy Simulation of Transition to Turbulence in Boundary Layers

Large-eddy simulation results for laminar-to-turbulent transition in a spatially developing boundary layer are presented. The disturbances are ingested into a laminar flow through an unsteady suction-and-blowing strip. The filtered, three-dimensional time- dependent Navier-Stokes equations are integrated numerically using spectral, high-order finite-difference, and three-stage low-storage Runge-Kutta methods. The buffer-domain technique is used for the outflow boundary condition. The localized dynamic model used to parameterize the subgrid-scale stresses begins to have a significant impact at the beginning of the nonlinear transition (or intermittency) region. The flow structures commonly found in experiments are also observed in the present simulation; the computed linear instability modes and secondary instability lambda-vortex structures are in agreement with the experiments, and the streak-like-structures and turbulent statistics compare with both the experiments and the theory. The physics captured in the present LES are consistent with the experiments and the full Navier-Stokes simulation (DNS), at a significant fraction of the DNS cost. A comparison of the results obtained with several SGS models shows that the localized model gives accurate results both in a statistical sense and in terms of predicting the dynamics of the energy-carrying eddies, without ad hoc adjustments.

Huai, Xiao-Li↗

Intercomparison of the Cloud Water Phase among Global Climate Models

Mixed-phase clouds (clouds that consist of both cloud droplets and ice crystals) are frequently present in the Earth’s atmosphere and influence the Earth’s energy budget through their radiative properties, which are highly dependent on the cloud water phase. In this study, the phase partitioning of cloud water is compared among six global climate models (GCMs) and with Cloud and Aerosol Lidar with Orthogonal Polarization retrievals. It is found that the GCMs predict vastly different distributions of cloud phase for a given temperature, and none of them are capable of reproducing the spatial distribution or magnitude of the observed phase partitioning. While some GCMs produced liquid water paths comparable to satellite observations, they all failed to preserve sufficient liquid water at mixed-phase cloud temperatures. Our results suggest that validating GCMs using only the vertically integrated water contents could lead to amplified differences in cloud radiative feedback. The sensitivity of the simulated cloud phase in GCMs to the choice of heterogeneous ice nucleation parameterization is also investigated. The response to a change in ice nucleation is quite different for each GCM, and the implementation of the same ice nucleation parameterization in all models does not reduce the spread in simulated phase among GCMs. The results suggest that processes subsequent to ice nucleation are at least as important in determining phase and should be the focus of future studies aimed at understanding and reducing differences among the models.

cloud water phase↗

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗

An Overview of NASA’s Newest Engineering Model, ORDEM 4.0

Since the mid-1990s, one of the most important products produced by the NASA Orbital Debris Program Office (ODPO) has been the Orbital Debris Engineering Model (ORDEM). This series of models distills down our knowledge of the orbital debris environment to compute debris fluxes on satellites in a given orbit. This information can be used by spacecraft and upper stage designers and operators to design missions for better protection against the debris environment. The current version of the model is ORDEM 3.2, but the ODPO is working on the next generation of ORDEM, to be designated ORDEM 4.0. ORDEM 4.0 will include many known features from previous models, such as the ability to input a spacecraft orbit and time and to compute the flux as a function of debris size, impact speed, impact direction, and debris material densities, as well as uncertainty information on the flux. ORDEM 4.0 will update debris populations using the most recent measurements, including radar observations by the Haystack Ultrawideband Satellite Imaging Radar (HUSIR), NASA’s Goldstone radar, data from the new Space Surveillance Network Space Fence, and observations of Geosynchronous Earth Orbits (GEO) using the Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT). The latest in situ impact data from returned hardware surfaces will be used. In addition, ORDEM 4.0 will introduce a parameterized debris shape model based on laboratory hypervelocity impact tests, such as DebriSat. This will allow analysts to implement shape characteristics in their damage equations and more accurately predict impact damage risk by debris of different shapes and orientations. This paper provides an overview of some of the new features forthcoming in ORDEM 4.0 and a status report on its development.

Mark Matney↗