SEARCH · Engineering Papers
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.
Co-Author Voon Hui Lai, Australian National University, will be attending IUGG and presenting the work. The session of interested is "S05 - Advances in Earthquake and Explosion Monitoring Using Distributed Acoustic Sensing"
For sustainable biomass production of Miscanthus × giganteus (hereafter miscanthus), understanding the impact of stand age and nitrogen (N) fertilization on biomass yield is crucial. This study investigated the effects of varying N fertilization rates (0, 56, 112, and 168 kg N ha-1) on yield components (tiller height, density, and weight) and their correlations with end-of-season biomass yield in miscanthus. We also explored end-of-season biomass yield prediction using in-season traits (canopy height, leaf area index (LAI), and leaf chlorophyll content (LCC)). The study was conducted at two sites in Illinois: a previously unfertilized 10-year-old miscanthus research stand at Urbana and a 16-year-old commercial stand at Pesotum with a history of annual 56N application. Results from 2018-2021 in Urbana and 2020-2021 in Pesotum showed increased biomass yields with N fertilization, varying by rate, year, and location. Biomass yield in Pesotum peaked at 56N, while in Urbana, it increased significantly at 112 kg N ha-1. Biomass yield was strongly correlated with tiller height and weight measured at Urbana across N rates. Morphological traits measured every 2-3 weeks during the 2020 and 2021 growing seasons showed that canopy height was the strongest single predictor of miscanthus biomass yield, followed by LCC. Mid-August to September measurements of these traits were the best predictors of biomass yield. Multiple regressions involving the canopy height and LCC further improved yield predictions. We conclude that while N enhances biomass yields of aging miscanthus, the optimum rate depends on the site, environmental conditions, and management.
The Land Surface data include bihemispherical and directional-hemispherical reflectances (albedo), hemispherical directional and bidirectional reflectance factors (BRF), BRF model parameters, leaf-area index (LAI), fraction of photosynthetically active radiation (FPAR), and normalized difference vegetation index (NDVI) on a 1.1 km grid. The land surface data include hemispherical directional reflectance factor, bihemispherical reflectance (i.e., albedo), bidirectional reflectance factor, directional hemispherical reflectance, BRF model parameters, FPAR, and terrain-referenced view and illumination angles. [Location=GLOBAL LAND] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1.1 km - 17.6 km; Longitude_Resolution=1.1 km - 17.6 km; Horizontal_Resolution_Range=1 km - < 10 km or approximately .01 degree - < .09 degree; Temporal_Resolution=about 15 orbits/day].
The Land Surface data include bihemispherical and directional-hemispherical reflectances (albedo), hemispherical directional and bidirectional reflectance factors (BRF), BRF model parameters, leaf-area index (LAI), fraction of photosynthetically active radiation (FPAR), and normalized difference vegetation index (NDVI) on a 1.1 km grid. The land surface data include hemispherical directional reflectance factor, bihemispherical reflectance (i.e., albedo), bidirectional reflectance factor, directional hemispherical reflectance, BRF model parameters, FPAR, and terrain-referenced view and illumination angles. [Location=GLOBAL LAND] [Temporal_Coverage: Start_Date=2000-02-24; Stop_Date=] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Latitude_Resolution=1.1 km - 17.6 km; Longitude_Resolution=1.1 km - 17.6 km; Horizontal_Resolution_Range=1 km - < 10 km or approximately .01 degree - < .09 degree; Temporal_Resolution=about 15 orbits/day; Temporal_Resolution_Range=Daily - < Weekly, Daily - < Weekly].
NorthWest Research Associates (NWRA) has developed an Inverse Model for inverting aircraft wake vortex data. The objective of the inverse modeling is to obtain estimates of the vortex circulation decay and crosswind vertical profiles, using time history measurements of the lateral and vertical position of aircraft vortices. The Inverse Model performs iterative forward model runs using estimates of vortex parameters, vertical crosswind profiles, and vortex circulation as a function of wake age. Iterations are performed until a user-defined criterion is satisfied. Outputs from an Inverse Model run are the best estimates of the time history of the vortex circulation derived from the observed data, the vertical crosswind profile, and several vortex parameters. The forward model, named SHRAPA, used in this inverse modeling is a modified version of the Shear-APA model, and it is described in Section 2 of this document. Details of the Inverse Model are presented in Section 3. The Inverse Model was applied to lidar-observed vortex data at three airports: FAA acquired data from San Francisco International Airport (SFO) and Denver International Airport (DEN), and NASA acquired data from Memphis International Airport (MEM). The results are compared with observed data. This Inverse Model validation is documented in Section 4. A summary is given in Section 5. A user's guide for the inverse wake vortex model is presented in a separate NorthWest Research Associates technical report (Lai and Delisi, 2007a).
NorthWest Research Associates (NWRA) has developed an inverse model for inverting landing aircraft vortex data. The data used for the inversion are the time evolution of the lateral transport position and vertical position of both the port and starboard vortices. The inverse model performs iterative forward model runs using various estimates of vortex parameters, vertical crosswind profiles, and vortex circulation as a function of wake age. Forward model predictions of lateral transport and altitude are then compared with the observed data. Differences between the data and model predictions guide the choice of vortex parameter values, crosswind profile and circulation evolution in the next iteration. Iterations are performed until a user-defined criterion is satisfied. Currently, the inverse model is set to stop when the improvement in the rms deviation between the data and model predictions is less than 1 percent for two consecutive iterations. The forward model used in this inverse model is a modified version of the Shear-APA model. A detailed description of this forward model, the inverse model, and its validation are presented in a different report (Lai, Mellman, Robins, and Delisi, 2007). This document is a User's Guide for the Wake Vortex Inverse Model. Section 2 presents an overview of the inverse model program. Execution of the inverse model is described in Section 3. When executing the inverse model, a user is requested to provide the name of an input file which contains the inverse model parameters, the various datasets, and directories needed for the inversion. A detailed description of the list of parameters in the inversion input file is presented in Section 4. A user has an option to save the inversion results of each lidar track in a mat-file (a condensed data file in Matlab format). These saved mat-files can be used for post-inversion analysis. A description of the contents of the saved files is given in Section 5. An example of an inversion input file, with preferred parameters values, is given in Appendix A. An example of the plot generated at a normal completion of the inversion is shown in Appendix B.
In this paper, we explore an alternative approach for estimating LAI for remotely sensed data from AVIRIS based on estimates of canopy liquid water.
Many secondary deciduous forests of eastern North America are approaching a transition in which maturing early successional tree species are declining, resulting in an uncertain future for these forests’ century-long carbon (C) sink. We initiated the Forest Accelerated Succession Experiment (FASET) at the University of Michigan Biological Station to examine the patterns and mechanisms underlying forest C cycling following the stem-girdling induced mortality of >6,700 early successional Populus spp. (aspen) and Betula papyrifera (paper birch). Since 2008, meteorological flux tower-based C cycling observations from the 33-ha treatment forest have been paired with those from a nearby unmanipulated forest. Following over a decade of observations, herein we revisit our core hypothesis: that net ecosystem production (NEP) would increase following the transition to mid-late successional species dominance due to increased canopy structural complexity. Supporting our hypothesis, NEP was stable, briefly declined, and then increased relative to the control in the decade following disturbance; however, increasing NEP was not associated with rising structural complexity but rather the rapid 1-year recovery of total LAI as mid-late successional Acer, Quercus, and Pinus assumed canopy dominance. The transition to mid-late successional species dominance improved carbon-use efficiency (CUE = NEP/gross primary production) as ecosystem respiration declined. Similar soil respiration rates in control and treatment forests, along with species differences in leaf physiology and the rising relative growth rates of mid-late successional species in the treatment forest, suggest changes in aboveground plant respiration and growth, respectively, were primarily responsible for increases in NEP. We conclude that deciduous forests transitioning from early to middle succession are capable of sustained or increased NEP, even when experiencing extensive tree mortality. This adds to mounting evidence that aging deciduous forests in the region will function as C sinks for decades to come.
Vegetation phenology has received increasing attention in climate change research. Near-surface sensing using digital repeat photography has proven to be useful for ecosystem-scale monitoring of vegetation phenology. However, our understanding of the link between phenological metrics derived from digital repeat photography and the phenology of forest canopy photosynthesis is still incomplete, especially for evergreen plant species. Using 49 site-years of digital images from the PhenoCam Network from eight evergreen needleleaf forest (ENF) and six deciduous broadleaf forest (DBF) sites in North America, we explored the potential of the green chromatic (GCC) and red chromatic coordinates (RCC) in tracking forest canopy photosynthesis by comparing camera-derived start- and end-of-growing season (SOS and EOS, respectively) with corresponding estimates derived from eddy covariance-derived daily gross primary productivity (GPP). We found that for DBF sites, both GCC and RCC performed comparable in capturing SOS and EOS. However, similar to earlier studies, GCC had limited potential in capturing GPP-based SOS or EOS for ENF sites. In contrast, we found RCC was a better predictor of both GPP-based SOS and EOS for ENF sites. Environmental and ecological explanations were both provided that phenological transitions derived from RCC were highly correlated with spring and autumn meteorological conditions, as well as having overall higher correlations with phenology based on LAI, a critical variable for describing canopy development. Our results demonstrate that RCC is an underappreciated metric for tracking vegetation phenology, especially for ENF sites where GCC failed to provide reliable estimates for GPP-based SOS or EOS. Our results improve confidence in using digital repeat photography to characterize the phenology of canopy photosynthesis across forest types.
Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.
This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.
We report many studies have shown that solar-induced fluorescence (SIF) has a good potential to predict gross primary production (GPP) of vegetation. What we measured by remote sensing or near-surface platforms is top-of-canopy (TOC) SIF (SIF toc ), which is not necessarily equal or proportional to total emitted SIF (SIF tot ) from the entire canopy due to the (re)absorption and scattering effects. However, photosynthesis, the process that plants use to fix carbon from the atmosphere, occurs at the entire vertical canopy. Here, by using the recollision theory, we calculated SIF tot at 760 nm from the measured SIF toc , hyperspectral reflectance (R), canopy interception (i) and leaf albedo (ω l ). Among them, both SIF toc and R can be obtained from concurrent TOC spectral measurements; i and ω l in the near-infrared region can be estimated from the open access datasets with a good accuracy. Our result confirms that the measured SIF toc only accounts for a small fraction of SIF tot : the above-the-canopy sensor can only “see” on average 22.9% of SIF tot at Harvard Forest. SIF tot has the following advantages over SIF toc in estimating GPP: (1) SIF tot improves the diurnal estimate of canopy GPP, especially capable to capture the midday depression of photosynthesis which may cause the large discrepancies between SIF toc and GPP on a diurnal basis, (2) SIF tot produces a stronger correlation with GPP from plants with complex canopy structure or under sky conditions with more diffuse irradiance, and (3) the SIF tot -GPP relationship shows a stronger resilience to environmental stresses. The fluorescence escape ratio (f esc ), the ratio between SIF toc multiplied by π and SIF tot , is mostly determined by the sun-canopy-sensor geometry and leaf inclination distribution. The effect of LAI and the leaf chlorophyll concentration on f esc is marginal at the 760 nm wavelength. Our results suggest that converting SIF toc into SIF tot provides a better tool to understand and estimate GPP.
Abstract Increasing temperature and water cycle changes due to warming climate may increase the frequency and intensity of wildfires. Fire emission projections are useful for informing strategies for adaptation and mitigation of fire impacts on societies and ecosystems. Here, we construct a neural network (NN) model explained by the Shapley Additive explanation to predict fire PM 2.5 emissions change and understand their drivers over the contiguous US (CONUS) in the mid‐21st century under a high greenhouse gas emissions scenario (SSP5‐8.5). Using future meteorology and leaf area index (LAI) simulated by eight global climate models from the Coupled Model Intercomparison Project Phase 6, future population density, and present‐day land use and land cover (LULC) as input to the NN model, the total fire PM 2.5 emissions over CONUS are projected to increase by 4%–75% (model spread). Among different regions, fire emissions in the western US are projected to increase more significantly in June‐July‐August than in other seasons and regions, with the median ratios of future to present‐day fire emissions ranging from 1.67 to 2.86. The increases in fire emissions are mainly driven by increasing normalized temperature (23%–29%) and decreasing soil moisture (2%–10%) in the future. When future LULC change is considered, the projected fire emissions further increase by 58%–83% over the western US compared to projections without LULC change because of future increases in vegetation fraction. The results highlight the important role of warmer temperature, decreasing soil moisture, and LULC change in increasing fire emissions in the future.
Here, we investigate the spin and orbital correlations of a superexchange model with spin $\textit{S}$ = 1 and orbital $\textit{L}$ = 1 relevant for $5d^4$ transition-metal Mott insulators, using exact diagonalization and density matrix renormalization group (DMRG). For spin-orbit coupling λ = 0, the orbitals are in an entangled state that is decoupled from the spins. We find two phases with increasing λ: (I) the $\textit{S}$2 phase with two peaks in the structure factor for λ ≤ λ$_{c1}$ ≈ 0.34$\textit{J}$ where $\textit{J}$ is the ferromagnetic exchange; and, (II) the $\textit{S}$1 phase for λ$_{c1}$ < λ ≤ λ$_{c2}$ ≈ 1.2$\textit{J}$ with emergent antiferromagnetic correlations. Both $\textit{S}$1 and $\textit{S}$2 phases are shown to exhibit power-law correlations, indicative of a gapless spectrum. Upon increasing λ > λ$_{c2}$ leads to a product state of local spin-orbital singlets that exhibit exponential decay of correlations, indicative of a gapped phase. We obtain insights into the phases from the well-known Uimin-Lai-Sutherland model in an external field that provides an approximate description of our model within mean-field theory.
Using a matrix product state algorithm with infinite boundary conditions, we compute high-resolution dynamic spin and quadrupolar structure factors in the thermodynamic limit to explore the low-energy excitations of isotropic bilinear-biquadratic spin-1 chains. Haldane mapped the spin-1 Heisenberg antiferromagnet to a continuum field theory, the nonlinear sigma model (NLσM). In this work, we find that the NLσM fails to capture the influence of the biquadratic term and provides only an unsatisfactory description of the Haldane phase physics. But several features in the Haldane phase can be explained by noninteracting multimagnon states. The physics at the Uimin-Lai-Sutherland point is characterized by multisoliton continua. Moving into the extended critical phase, we find that these excitation continua contract, which we explain using a field-theoretic description. New excitations emerge at higher energies and, in the vicinity of the purely biquadratic point, they show simple cosine dispersions. Using block fidelities, we identify them as elementary one-particle excitations and relate them to the integrable Temperley-Lieb chain.
We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field h. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor S(q,ω) as a function of h. At h=0, the envelope of the lowest-energy excitations goes soft at two points, q 1 =2π/3 and q 2 =4π/3, dubbed the phase A. With increasing field, the spectral peaks at each of the gapless points bifurcate, making in total four soft modes, and combine to form a new set of excitations that soften at a single point q=π at h c1 ≈0.94. Beyond h c1 the system enters another gapless B phase until the transition at h c2 =4 to the fully polarized phase. We compare the ULS model results with those for the Affleck-Kennedy-Lieb-Tasaki model as a representative of the gapped Haldane phase. We explain the mechanism of the gapless to gapless transition in the ULS model using its conserved charges and a spinon band picture. We also discuss the universality of central charges of the BLBQ family of models subjected to a magnetic field.
The availability of global high-resolution land cover maps provides promising a priori knowledge for characterizing subpixel heterogeneity and improving predictions of directional reflectance of coarse-resolution pixels. Due to mutual shadowing and sheltering effects between the adjacent forest and cropland patches, the spectral nonlinear mixing of patchy ecotones is significant, especially when the sun illuminates the ecotone from the forest side with high solar zenith angle. The spectral linear mixture (SLM) approach leads to overestimation of the bidirectional reflectance factor (BRF) in the red band in the principal plane (PP), with a maximum absolute error (MAE) of 0.0063 and a maximum relative error (MRE) of 52.5%, and to underestimation in the near-infrared band in PP with an MAE of 0.0940 and an MRE of 14.5%. In a scenario with randomly distributed boundary orientations, the overestimation of SLM increases with the degree of fragmentation and the view zenith angle. We propose a Radiative Transfer model for patchy ECotones (RTEC). which improves R² from 0.61 to 0.94 in the red band of Landsat-8 directional reflectance at the validation site. The RTEC model provides an efficient and analytical approach for directional reflectance predictions over heterogeneous patchy landscapes at coarse resolution and will be used for biophysical parameter retrievals [e.g., the leaf area index (LAI)] in future applications.
Here, we present a de novo, extensible, modular, research-oriented, single-player reimplementation of the "AlphaZero" class of algorithms originally described in [1]. We also present a novel task called the Zoning Game, designed as a toy problem to test the applicability of various AlphaZero-based paradigms to infrastructure planning problems, along with several other games. We implement specializations of AlphaZero to play several games. We also include code that serves as a starting point to explore language generation with AlphaZero, as well as documentation and other related materials. [1] Silver, D., Schrittwieser, J., Simonyan, K., Antonoglou, I., Huang, A., Guez, A., Hubert, T., Baker, L., Lai, M., Bolton, A., Chen, Y., Lillicrap, T., Hui, F., Sifre, L., van den Driessche, G., Graepel, T., & Hassabis, D. (2017). Mastering the game of Go without human knowledge. Nature, 550(7676), 354–359. https://doi.org/10.1038/nature24270