Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

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 361 records · Page 20

Rapid Optimization of Total Variation with Applications in Imaging, Additive Manufacturing, and Qualification

Total Variation optimization penalizes the gradient of a control variable or state. While this work focuses on image processing in particular, it has also found applications in inverse problems and topology optimization. In image processing, the goal is to maintain faithfulness to the original image while denoising and/or deblurring. Additionally, bilevel optimization over the spatially varying regularization weights can illuminate interfaces such as damage regions and other anomalies. We will address two fundamental challenges with TV-optimization: (i) the typical slow convergence of existing TV-optimization methods, and (ii) the selection of spatially varying TV parameters to promote interface detection. Additionally, we will apply such techniques to image data collected in additive manufacturing. In said context, stochasticity in build events induces flaws in the manufactured piece, compromising the integrity of said part. There is a critical need for in-situ monitoring to spot anomalies once they form, and in this setting we apply our total variation and hyperparameter solvers. We will develop a customized algorithm based on for extreme-scale TV-optimization that achieves super-linear or quadratic-convergence, a critical property for real-time, image-by-image analysis. A worst-case outcome is a preprocessing step that enhances image quality in-situ, specifically for out-of-focus and noisy images.

36 MATERIALS SCIENCE↗

Effects of SRPPF ARS Effluent Variations on Magnesium Oxysulfate Solidification

A magnesium oxysulfate (MOS) grout formulation has been identified to solidify the liquid effluent from the Savannah River Plutonium Processing Facility (SRPPF) Aqueous Recovery System (ARS). The formulation uses reactive, light burnt MgO (28 wt%), anhydrous MgSO 4 (10 wt%), and dead burnt MgO (62 wt%), resulting in a grout with good mixability, acceptable density, acceptable setting times, and a leachate pH around 9.4, within the assumed Waste Isolation Pilot Plant (WIPP) brine pH range. A series of tests were then designed to observe how this formula reacted to variations in the liquid effluent. Liquid effluent temperature, the caustic pH level, concentration of NaNO 3 , and the presence of neutralization products, trace metals, and/or sodium sulfate (all potential residuals found in the effluent from upstream processing), were evaluated to see how each impacts the grout mixing time, peak temperature, setting time, presence of bleed water, leachate pH, and density. The standard mix has an average mixing time of 18 minutes, peak temperature of 100.5 °C, and set time of 60 minutes. Increasing the temperature of the liquid effluent increased the reaction rate of the mix and reduced mixing and setting times. Between 30 and 40°C the change was relatively small, raising no more than 6 °C compared to the standard mix maximum temperature. Increasing the pH of the liquid effluent had no significant impact on the mix. The addition of neutralization products and trace metals had an overall impact of decreasing the reactivity of the mix. The addition of more ions in the liquid effluent, such as an increase in the concentration of NaNO 3 and the addition of Na 2 SO 4 , generally decreased the reactivity of the mix. None of the variations to the liquid effluent hindered solidification as indicated by the lack of bleed water found on all the samples. The leachate pH values for all mixes tested did not significantly vary from the expected pH of 9.4 and none of the density values fell below 1.8 g/cm 3 . Overall, changes to the liquid effluent were found to have a minimal impact on the formulation suggesting the grout's ability to properly form despite increased temperatures and the presence of residuals typically found in the liquid effluent from the ARS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

VRN3P: Variational Recurrent Neural Network Based Net-Load Prediction under High Solar Penetration

This is the final technical report for the SETO-funded VRN3P project (PNNL# 76914). The goal of this project, led by Pacific Northwest National Laboratory (PNNL), in collaboration with Lawrence Livermore National Laboratory (LLNL) and Portland General Electric (PGE), was to develop and validate a deep variational recurrent neural network-based net-load prediction (VRN3P) framework for probabilistic time-series forecasting of day-ahead net-load under high solar penetration scenarios. The project team reports successful design of a novel probabilistic net-load forecasting architecture, comprising of a variational autoencoder and a recurrent neural network, which demonstrates 30% improvement in forecast performance, 60% improvement in training time, and consumes 44% less memory, when compared with conventional baseline models. The team tested the VRN3P model performance on GridLAB-D test-cases representing varying BTM solar penetration levels of 20%, 30%, and 50%, with integrated time-series net-load profiles provided by the utility partner (PGE). The VRN3P model demonstrate <2% hourly MAPE (averaged over the year) for day- ahead net-load forecast on the test scenario with 20% BTM solar. Transfer learning extension of the VRN3P model has demonstrated 8.33× speed-up in training, while still achieving acceptable forecast performance of 2.24% hourly MAPE on the 30% BTM solar penetration test-scenario. A preliminary version of the VRN3P GridAPPS-D™has been developed, along with a web-based interactive user-interface (named ‘Forte’) which has made available on GitHub for public use.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Variational Quantum Linear Solver

Previously proposed quantum algorithms for solving linear systems of equations cannot be implemented in the near term due to the re quired circuit depth. Here, we propose a hybrid quantum-classical algorithm, called Variational Quantum Linear Solver (VQLS), for solving linear systems on near-term quantum computers. VQLS seeks to variationally prepare |x$\rangle$ such that A|x$\rangle$ ∝ |b$\rangle$. We derive an operationally meaningful termination condition for VQLS that allows one to guarantee that a desired solution precision ϵ is achieved. Specifically, we prove that C $&GreaterSlantEqual;$ ϵ 2 /κ 2 , where C is the VQLS cost function and κ is the condition number of A. We present efficient quantum circuits to estimate C, while providing evidence for the classical hardness of its estimation. Using Rigetti’s quantum computer, we success fully implement VQLS up to a problem size of 1024 × 1024. Finally, we numerically solve nontrivial problems of size up to 2 50 × 2 50 . For the specific examples that we consider, we heuristically find that the time complexity of VQLS scales efficiently in ϵ, κ, and the system size N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarking Adaptive Variational Quantum Eigensolvers

By design, the variational quantum eigensolver (VQE) strives to recover the lowest-energy eigenvalue of a given Hamiltonian by preparing quantum states guided by the variational principle. In practice, the prepared quantum state is indirectly assessed by the value of the associated energy. Novel adaptive derivative-assembled pseudo-trotter (ADAPT) ansatz approaches and recent formal advances now establish a clear connection between the theory of quantum chemistry and the quantum state ansatz used to solve the electronic structure problem. Here we benchmark the accuracy of VQE and ADAPT-VQE to calculate the electronic ground states and potential energy curves for a few selected diatomic molecules, namely H 2 , NaH, and KH. Using numerical simulation, we find both methods provide good estimates of the energy and ground state, but only ADAPT-VQE proves to be robust to particularities in optimization methods. Another relevant finding is that gradient-based optimization is overall more economical and delivers superior performance than analogous simulations carried out with gradient-free optimizers. The results also identify small errors in the prepared state fidelity which show an increasing trend with molecular size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of the Aeroelastic Dynamics of Lightweight Flexible Variations of the SNL-NRT Turbine

Current trends show that wind turbines are growing in size to meet a rising demand for renewable energy generation, and their upscaled rotors have inherently become more flexible to maintain a proportionally lighter design. This is because larger rotors must be less massive relative to their diameter to minimize the levelized cost of energy (LCOE), which means that blades that are notably less stiff are produced as a result. These structural changes to blades are often reflected in their compromised aeroelastic stability and amplified deformation during operation, which has the potential to decrease the blade’s expected lifetime and the performance of the machine overall. Variations in blade flexibility are also known to influence vortex-wake structures downstream of the turbine, causing patterns of velocity deficit to evolve in ways that affect the performance of other turbines in the farm. This research explores how the increased flexibility of modern utility-scale wind turbine blades influences rotor aeroelastic behavior and interactions with farm flow. High-fidelity simulations of Sandia National Laboratories’ (SNL) National Rotor Testbed (NRT) wind turbine are presented. Flexible variations of the NRT baseline blade are simulated in a variety of realistic operational conditions typically expected at the SNL’s SWiFT facility in Lubbock, Texas. Solutions are then compared to investigate how specific changes to the structural properties of the NRT baseline blade’s design and construction can influence its aeroelastic response at the rotor and the evolution of the turbine’s wake.

Farrell, Alayna (ORCID:000000023555720X)↗

Evaluation of the Four-Dimensional Ensemble-Variational Hybrid Data Assimilation with Self-Consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts

54 ENVIRONMENTAL SCIENCES↗

Genomic Analysis of the Natural Variation of Fatty Acid Composition in Seed Oils of Camelina sativa

Camelina sativa is an oilseed crop that has shown strong promise as a biofuel feedstock. The profile of fatty acids greatly influences the oil quality; however, genetic mechanisms that determine the natural variation of fatty acid composition in camelina are not fully understood. A genome wide association study (GWAS) was performed to uncover genetic loci that may contribute to the contents of major fatty acids such as oleic and linolenic acids in camelina seed. Two approaches were taken to improve the GWAS efficiency. First, growing a diversity panel of 212 accessions in four locations and two nitrogen fertilization conditions revealed great variation in fatty acid contents in seeds. Second, using an improved reference genome, abundant markers, including 203,320 single nucleotide polymorphisms (SNPs) and 99,067 insertions/deletions (indels), were developed, which refined the population structure of the diversity panel. GWAS resulted in 118 genetic markers across 31 trait/treatment conditions. Closely linked markers were determined based on linkage decay and by comparing secondarily associated markers when highly associated ones were removed. Candidate genes were examined by comparing the pangenomes of 12 high-quality reference genomes. This study provides new resources to understand seed lipid metabolism and improve camelina oils through molecular breeding.

Life Sciences & Biomedicine - Other Topics↗

Dimensionality Reduction with Variational Encoders Based on Subsystem Purification

Efficient methods for encoding and compression are likely to pave the way toward the problem of efficient trainability on higher-dimensional Hilbert spaces, overcoming issues of barren plateaus. Here, we propose an alternative approach to variational autoencoders to reduce the dimensionality of states represented in higher dimensional Hilbert spaces. To this end, we build a variational algorithm-based autoencoder circuit that takes as input a dataset and optimizes the parameters of a Parameterized Quantum Circuit (PQC) ansatz to produce an output state that can be represented as a tensor product of two subsystems by minimizing $Tr(ρ^2)$. The output of this circuit is passed through a series of controlled swap gates and measurements to output a state with half the number of qubits while retaining the features of the starting state in the same spirit as any dimension-reduction technique used in classical algorithms. The output obtained is used for supervised learning to guarantee the working of the encoding procedure thus developed. We make use of the Bars and Stripes (BAS) dataset for an 8 × 8 grid to create efficient encoding states and report a classification accuracy of 95% on the same. Thus, the demonstrated example provides proof for the working of the method in reducing states represented in large Hilbert spaces while maintaining the features required for any further machine learning algorithm that follows.

97 MATHEMATICS AND COMPUTING↗

Crystal Chemistry and Structural Variations for Zircon Samples from Various Localities

This study investigates the variations of structural parameters and chemistry of a partially metamict and seven detrital zircon samples from different localities using single-crystal X-ray diffraction, synchrotron high-resolution powder X-ray diffraction, and electron-probe micro-analysis techniques. The unit-cell parameters for the eight zircon samples vary linearly with increasing unit-cell volume, V. A zircon sample from the Canadian Arctic Islands has the smallest unit-cell parameters, bond distances, ideal stoichiometric composition, unaffected by α-radiation damage, and is chemically pure. A zircon sample from Jemaa, Nigeria has the largest unit-cell parameters because of the effect of α-radiation doses received over a long time (2384 Ma). All the samples show good correlations between Zr and Si apfu (atom per formula unit) versus unit-cell volume, V. The α-radiation doses in the samples are lower than ~3.5 × 1015 α-decay events/mg. Substitutions of other cations at the Zr and Si sites control the variations of the structural parameters. Relatively large unit-cell parameters and bond distances occur because the Zr site accommodates other cations that have larger ionic radii than the Zr atom. Geological age increases the radiation doses in zircon and it is related to V.

58 GEOSCIENCES↗

Comprehensive Genome-Wide Natural Variation and Expression Analysis of Tubby-like Proteins Gene Family in Brachypodium distachyon

The Tubby-like proteins (TLPs) gene family is a group of transcription factors found in both animals and plants. In this study, we identified twelve B. distachyon TLPs, divided into six groups based on conserved domains and evolutionary relationships. We predicted cis-regulatory elements involved in light, hormone, and biotic and abiotic stresses. The expression patterns in response to light and hormones revealed that BdTLP3, 4, 7, and 14 are involved in light responses, and BdTLP1 is involved in ABA responses. Furthermore, BdTLP2, 7, 9, and 13 are expressed throughout vegetative and reproductive stages, whereas BdTLP1, 3, 5, and 14 are expressed at germinating grains and early vegetative development, and BdTLP4, 6, 8, and 10 are expressed at the early reproduction stage. The natural variation in the eleven most diverged B. distachyon lines revealed high conservation levels of BdTLP1-6 to high variation in BdTLP7-14 proteins. Based on diversifying selection, we identified amino acids in BdTLP1, 3, 8, and 13, potentially substantially affecting protein functions. This analysis provided valuable information for further functional studies to understand the regulation, pathways involved, and mechanism of BdTLPs.

Plant Sciences↗

Efficient Application of the Factorized form of the Unitary Coupled-Cluster Ansatz for the Variational Quantum Eigensolver Algorithm by Using Linear Combination of Unitaries

The variational quantum eigensolver is one of the most promising algorithms for near-term quantum computers. It has the potential to solve quantum chemistry problems involving strongly correlated electrons with relatively low-depth circuits, which are otherwise difficult to solve on classical computers. The variational eigenstate is constructed from a number of factorized unitary coupled-cluster terms applied onto an initial (single-reference) state. Current algorithms for applying one of these operators to a quantum state require a number of operations that scale exponentially with the rank of the operator. We exploit a hidden SU(2) symmetry to allow us to employ the linear combination of unitaries approach, Our Prepare subroutine uses n+2 ancilla qubits for a rank-n operator. Our Select (U) scheme uses O(n) Cnot gates. This results in a full algorithm that scales like the cube of the rank of the operator n 3 , a significant reduction in complexity for rank five or higher operators. This approach, when combined with other algorithms for lower-rank operators (when compared to the standard implementation), will make the factorized form of the unitary coupled-cluster approach much more efficient to implement on all types of quantum computers.

97 MATHEMATICS AND COMPUTING↗

Global Characteristics of the IBEX Ribbon and Its Temporal Variation: Implications for Relationships with Solar Cycle Variability and the Interstellar Magnetic Field

The Interstellar Boundary Explorer (IBEX) has revealed the existence of the energetic neutral atom (ENA) ribbon, a narrow and nearly circular band of enhanced flux that enables probing plasma and magnetic field conditions outside the heliopause (HP). In this study, we investigate the global structure and temporal variation of the IBEX ribbon using 14 yr of IBEX-Hi ENA observations. We derive geometric parameters of the ribbon: center location, angular radius, eccentricity, and major axis orientation. We find that the ribbon exhibits systematic variations tied to the solar cycle. The ribbon center shows a solar cycle–modulated latitudinal motion at low energies. The angular radius increases with ENA energy and displays a pronounced enlargement between 2013 and 2018, coinciding with changes in high-latitude fast solar wind from polar coronal holes. The eccentricity remains consistently around 0.3 with some temporal evolution, while the major axis rotates gradually over time, becoming increasingly perpendicular to the ecliptic meridional plane. A notable hemispheric asymmetry has emerged over time: the ribbon in the northern hemisphere remains remarkably stable in position, whereas the southern ribbon exhibits temporal migration. By combining these results with solar wind observations, global heliosphere, and ribbon modeling, we suggest that this asymmetry is explained by a coupled effect of magnetic field draping and the latitude-dependent timing of fast solar wind from polar coronal holes. These results establish the most detailed long-term geometric characterization of the IBEX ribbon to date and provide observational support for models of secondary ENA production and magnetic field draping outside the HP.

Noh, Sung Jun [New Mexico Consortium, Los Alamos, ↗

CV4Quantum: Reducing the Sampling Overhead in Probabilistic Error Cancellation Using Control Variates

Quasiprobabilistic decompositions (QPDs) play a key role in maximizing the utility of near-term quantum hardware. For example, Probabilistic Error Cancellation (PEC) (an error mitigation technique) and circuit cutting (which enables large quantum computations to be performed on quantum hardware with a limited number of qubits) both involve QPDs. Computations based on QPDs typically incur large sampling overheads that grow exponentially with the number of error-terms mitigated or number of circuit-cuts employed, limiting their practical feasibility. In this work, we adapt the control variates variance reduction technique from the statistics literature in order to reduce the sampling overhead in QPD-based computations. We demonstrate our method using simulation experiments that mimic a realistic PEC scenario. In our experiments, we observed a more than 50% reduction in the number of samples needed to achieve a given precision, in more than 50% of the PEC-based estimations performed in the study when using our approach. We discuss how future research on constructing good control variates can lead to even stronger sampling overhead reduction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Seasonal variations in fire conditions are important drivers in the trend of aerosol optical properties over the south-eastern Atlantic

Abstract. From June to October, southern Africa produces one-third of the global biomass burning (BB) emissions by widespread fires. BB aerosols are transported westward over the south-eastern Atlantic with the mid-tropospheric winds, resulting in significant radiative effects. Ascension Island (ASI) is located midway between Africa and South America. From June 2016 to October 2017, a 17-month in situ observation campaign on ASI found a low single-scattering albedo (SSA) as well as a high mass absorption cross-section of black carbon (MACBC), demonstrating the strong absorbing marine boundary layer in the south-eastern Atlantic. Here we investigate the monthly variations of critical optical properties of BB aerosols, i.e. SSA and MACBC, during the BB seasons and the driving factors behind these variations. Both SSA and MACBC increase from June to August and decrease in September and October. The average SSA during the BB seasons is 0.81 at 529 nm wavelength, with the highest mean ∼ 0.85 in October and the lowest ∼ 0.78 in August. The absorption enhancement (Eabs) derived from the MACBC shows similar trends with SSA, with the average during the whole of the BB seasons at ∼ 1.96 and ∼ 2.07 in 2016 and 2017, respectively. As the Eabs is higher than the ∼ 1.5 commonly adopted value by climate models, this result suggests the marine boundary layer in the south-eastern Atlantic is more absorbing than model simulations. We find the enhanced ratio of BC to CO (ΔBC/ΔCO, equal to BC/ΔCO as the BC background concentration is considered to be 0) is well correlated with SSA and MACBC, providing a simple way to estimate the aerosol optical characteristics in the south-eastern Atlantic. The exponential function we proposed can approximate SSA and MACBC with BC/ΔCO, and when BC/ΔCO is small it can capture the rapid growth of SSA as BC/ΔCO decreases. BC/ΔCO is influenced by combustion conditions and aerosol scavenging. From the analysis of the location of BB, the primary source fuel, the water content in the fuel, combined with the mean cloud cover and precipitation in the transport areas of the BB plume, we conclude that the increase in BC/ΔCO from June to August is likely to be caused by burning becoming more flaming. The reduction in the water content of fuels may be responsible for the change in the burning conditions from June to August. The decrease in BC/ΔCO in September and October may be caused by two factors, one being a lower proportion of flaming conditions, possibly associated with a decrease in mean surface wind speed in the burning area, and the other being an increase in precipitation in the BB transport pathway, leading to enhanced aerosol scavenging, which ultimately results in an increase in SSA and MACBC.

54 ENVIRONMENTAL SCIENCES↗

The sources and diurnal variations of submicron aerosols in a coastal–rural environment near Houston, US

Aerosol properties were characterized at a rural site southwest of Houston from May to September 2022 during the intensive operation periods (IOPs) of the Tracking Aerosol Convection Interactions ExpeRiment (TRACER). Backward trajectory analysis reveals three major air mass types: marine air mass from the Gulf, urban air mass influenced by urban emissions, and regional air mass. Marine aerosols typically show a bimodal size distribution and have the lowest particle number and mass concentrations of PM 1 (particulate matter with an aerodynamic diameter of less than 1 µm), while aerosols from air masses strongly influenced by urban emissions exhibit the highest concentrations. Organic aerosol (OA) accounts for more than 50 % of PM 1 for urban and regional air masses, whereas sulfate is comparable to OA in marine air masses. Positive matrix factorization (PMF) analysis of aerosol mass spectra identifies 6 OA factors: hydrocarbon-like OA (HOA), OA from the oxidation of monoterpenes (MT-SOA), OA from the reactive uptake of isoprene epoxydiols by acidic sulfate particles (isoprene-SOA), oxygenated OA arising from shipping emissions (shipping-OOA), and two oxygenated OA factors with high O : C ratios (OOA1 and OOA2). OOA2 has the highest O : C ratio and exhibits elevated mass concentration in the afternoon. Similar diurnal variation of highly oxidized OA factors was commonly observed in the Houston area during previous studies and attributed to the SOA formation by photochemistry and mixing from aloft. Here, using air mass backward trajectories and a 1-D box model, we show the diurnal trend of OOA2 mass concentration is instead driven by changes in air mass arriving at the rural site. The air mass changes are likely caused by the shift between land breezes and sea/bay breezes. Within the same air mass type (e.g., either urban or marine air mass), OOA2 mass concentration is largely independent of wind direction and shows essentially no diurnal variation, suggesting OOA2 is related to aged OA with minimal influence by local emissions. This study helps identify the major sources of OA in the Houston region and highlights the impacts of both atmospheric chemistry and meteorology on aerosol properties in the coastal–rural environment.

54 ENVIRONMENTAL SCIENCES↗

DeepPhenoMem V1.0: deep learning modelling of canopy greenness dynamics accounting for multi-variate meteorological memory effects on vegetation phenology

Abstract. Vegetation phenology plays a key role in controlling the seasonality of ecosystem processes that modulate carbon, water and energy fluxes between the biosphere and atmosphere. Accurate modelling of vegetation phenology in the interplay of Earth's surface and the atmosphere is thus crucial to understand how the coupled system will respond to and shape climatic changes. Phenology is controlled by meteorological conditions at different timescales: on the one hand, changes in key meteorological variables (temperature, water, radiation) can have immediate effects on the vegetation development; on the other hand, phenological changes can be driven by past environmental conditions, known as memory effects. However, the processes governing meteorological memory effects on phenology are not completely understood, resulting in their limited performance of vegetation phenology represented in land surface models. A deep learning model, specifically a long short-term memory network (LSTM), has the potential to capture and model the meteorological memory effects on vegetation phenology. Here, we apply the LSTM to model the vegetation phenology using meteorological drivers and high-temporal-resolution canopy greenness observations through digital repeat photography by the PhenoCam network. We compare a multiple linear regression model, a no-memory-effect LSTM model and a full-memory-effect LSTM model to predict the whole seasonal greenness trajectory and the corresponding phenological transition dates across 50 sites and 317 site years during 2009–2018, covering deciduous broadleaf forests, evergreen needleleaf forests and grasslands. Results show that the deep learning model outperforms the multiple linear regression model, and the full-memory-effect LSTM model performs better than the no-memory-effect model for all three plant function types (median R2 of 0.878, 0.957 and 0.955 for broadleaf forests, evergreen needleleaf forests and grasslands). We also find that the full-memory-effect LSTM model is capable of predicting the seasonal dynamic variations of canopy greenness and reproducing trends in shifting phenological transition dates. We also performed a sensitivity analysis of the full-memory-effect LSTM model to assess its plausibility, revealing its coherence with established knowledge of vegetation phenology sensitivity to meteorological conditions, particularly changes in temperature. Our study highlights that (1) multi-variate meteorological memory effects play a crucial role in vegetation phenology, and (2) deep learning opens up new avenues for improving the representation of vegetation phenological processes in land surface models via a hybrid modelling approach.

Geology↗

3D Continuous Forcing Dataset from 3D Constrained Variational Analysis at SGP

The continuous 3D large-scale forcing (VARANAL3D) data set derived from 3D constrained variational analysis (3DCVA) extends the conventional constrained variational analysis method by incorporating multiple sub-columns within the analysis domain. This advancement introduces spatial variability into the large-scale forcing fields, thereby enriching the data set’s applicability. The VARANAL3D data set spans from 2004 to 2018 and covers a region of 5˚×4.5˚ domain around the ARM SGP site. The analysis domain is divided into 10×9 sub-columns with 0.5˚ resolution. The 3D large-scale forcing data provides necessary variables to drive and evaluate single-column models (SCM), cloud-resolving models (CRM) ,and large-eddy simulations (LES), as well as information for testing model sensitivity to spatial variability of the large-scale forcing data, facilitating more rigorous testing and refinement of physical processes in SCM/CRM/LES.

54 ENVIRONMENTAL SCIENCES↗