Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced order representation”

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 163 records · Page 9

Exploiting Satellite Remote-Sensing Data in Fine Particulate Matter Characterization for Serving the Environmental Public Health Tracking Network (EPHTN): The HELIX-Atlanta Experience and NPOESS Implications

As part of the U.S. National Environmental Public Health Tracking Network (EPHTN), the National Center for Environmental Health (NCEH) at the U.S. Centers for Disease Control and Prevention (CDC) led a project in collaboration with the National Aeronautics and Space Administration (NASA) Marshall Space Flight Center (MSFC) called Health and Environment Linked for Information Exchange (HELIX-Atlanta). Under HELIX-Atlanta, pilot projects were conducted to develop methods to better characterize exposure; link health and environmental datasets; and analyze spatial/temporal relationships. This paper describes and demonstrates different techniques for surfacing daily environmental hazards data of particulate matter with aerodynamic diameter less than or equal to 2.5 micrometers (PM(sub 2.5) for the purpose of integrating respiratory health and environmental data for the CDC's pilot study of HELIX-Atlanta. It describes a methodology for estimating ground-level continuous PM(sub 2.5) concentrations using spatial surfacing techniques and leveraging NASA Moderate Resolution Imaging Spectrometer (MODIS) data to complement the U.S. Environmental Protection Agency (EPA) ground observation data. The study used measurements of ambient PM(sub 2.5) from the EPA database for the year 2003 as well as PM(sub 2.5) estimates derived from NASA's MODIS data. Hazard data have been processed to derive the surrogate exposure PM(sub 2.5) estimates. The paper has shown that merging MODIS remote sensing data with surface observations of PM(sub 2.5), may provide a more complete daily representation of PM(sub 2.5), than either data set alone would allow, and can reduce the errors in the PM(sub 2.5) estimated surfaces. Future work in this area should focus on combining MODIS column measurements with profile information provided by satellites like the National Polar-orbiting Operational Environmental Satellite System (NPOESS). The Visible Infrared Imager/Radiometer Suite (VIIRS) and the Aerosol Polarimeter Sensor (APS) NPOESS sensors will provide first-order information on aerosol particle size and are anticipated to provide information on aerosol products at higher resolution and accuracy than MODIS. Use of the NPOESS remote sensing data should result in more robust remotely sensed data that can be coupled with the methods discussed in this paper to generate surface concentrations of PM(2.5) for linkage with health data in Environmental Public Health Tracking.

Al-Hamdan, Mohammad Z.↗

Manifold Learning: What, How, and Why

Manifold learning (ML), also known as nonlinear dimension reduction, is a set of methods to find the low-dimensional structure of data. Dimension reduction for large, high-dimensional data is not merely a way to reduce the data; the new representations and descriptors obtained by ML reveal the geometric shape of high-dimensional point clouds and allow one to visualize, denoise, and interpret them. This review presents the underlying principles of ML, its representative methods, and their statistical foundations, all from a practicing statistician's perspective. It describes the trade-offs and what theory tells us about the parameter and algorithmic choices we make in order to obtain reliable conclusions.

Mathematics↗

Learning the generating functional for variance reduction in lattice QCD

The generating functional in quantum field theory provides the natural framework for constructing correlation functions as derivatives with respect to source operators. We present a methodology that leverages machine-learned normalizing flows to reduce the variance of arbitrary $N$-point correlation functions of bosonic operators in lattice gauge field theory calculations by encoding a representation of the generating functional. We show that it is possible to systematically approach noiseless estimators of correlation functions in this framework. We demonstrate this methodology with applications to calculations of glueball correlation functions and Wilson loops in Quantum Chromodynamics and Yang-Mills theory. The results show up to three orders of magnitude variance reduction.

Abbott, Ryan [Columbia U.] (ORCID:0000000258778005↗

Single-valued representation of unpolarized and polarized semi-inclusive deep inelastic scattering at next-to-next-to-leading order

We revisit the recently published analytic results for unpolarized and polarized semi-inclusive deep inelastic scattering (SIDIS) at next-to-next-to-leading order (NNLO) in quantum chromodynamics (QCD). These expressions for the hard scattering coefficients contain case distinctions in the kinematic (𝑥,𝑧)-plane, splitting the analytic result into four regions. By reexpressing the coefficient functions in terms of single-valued polylogarithms, we remove these case distinctions and can present a unified result valid across the entire kinematic range of SIDIS. This reduces the length of the overall expressions by 30% to 60%.

Deep inelastic scattering↗

A parameter optimization approach to controller partitioning for integrated flight/propulsion control application

A parameter optimization framework is presented to solve the problem of partitioning a centralized controller into a decentralized hierarchical structure suitable for integrated flight/propulsion control implementation. The controller partitioning problem is briefly discussed and a cost function to be minimized is formulated, such that the resulting 'optimal' partitioned subsystem controllers will closely match the performance (including robustness) properties of the closed-loop system with the centralized controller while maintaining the desired controller partitioning structure. The cost function is written in terms of parameters in a state-space representation of the partitioned sub-controllers. Analytical expressions are obtained for the gradient of this cost function with respect to parameters, and an optimization algorithm is developed using modern computer-aided control design and analysis software. The capabilities of the algorithm are demonstrated by application to partitioned integrated flight/propulsion control design for a modern fighter aircraft in the short approach to landing task. The partitioning optimization is shown to lead to reduced-order subcontrollers that match the closed-loop command tracking and decoupling performance achieved by a high-order centralized controller.

Schmidt, Phillip↗

A parameter optimization approach to controller partitioning for integrated flight/propulsion control application

A parameter optimization framework is presented to solve the problem of partitioning a centralized controller into a decentralized hierarchical structure suitable for integrated flight/propulsion control implementation. The controller partitioning problem is briefly discussed and a cost function to be minimized is formulated, such that the resulting 'optimal' partitioned subsystem controllers will closely match the performance (including robustness) properties of the closed-loop system with the centralized controller while maintaining the desired controller partitioning structure. The cost function is written in terms of parameters in a state-space representation of the partitioned sub-controllers. Analytical expressions are obtained for the gradient of this cost function with respect to parameters, and an optimization algorithm is developed using modern computer-aided control design and analysis software. The capabilities of the algorithm are demonstrated by application to partitioned integrated flight/propulsion control design for a modern fighter aircraft in the short approach to landing task. The partitioning optimization is shown to lead to reduced-order subcontrollers that match the closed-loop command tracking and decoupling performance achieved by a high-order centralized controller.

Schmidt, Phillip H.↗

Exotic Phases of a Higgs-Yukawa Model with Reduced Staggered Fermions

We investigate the phase structure of a four dimensional SO(4) invariant lattice Higgs-Yukawa model comprising four reduced staggered fermions interacting with a real scalar field. The fermions belong to the fundamental representation of the symmetry group while the three scalar field components transform in the self-dual representation of SO(4). We explore the phase diagram and find evidence of a continuous transition between a phase where the fermions are massless to one where the fermions acquire mass. This transition is not associated with symmetry breaking and there is no obvious local order parameter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Temporal logics meet telerobotics

The specificity of telerobotics being the presence of a human operator, decision assistance tools are necessary for the operator, especially in hostile environments. In order to reduce execution hazards due to a degraded ability for quick and efficient recovery of unexpected dangerous situations, it is of importance to have the opportunity, amongst others, to simulate the possible consequences of a plan before its actual execution, in order to detect these problematic situations. Hence the idea of providing the operator with a simulator enabling him to verify the temporal and logical coherence of his plans. Therefore, the power of logical formalisms is used for representation and deduction purposes. Starting from the class of situations that are represented, a STRIPS (the STanford Research Institute Problem Solver)-like formalism and its underlying logic are adapted to the simulation of plans of actions in time. The choice of a temporal logic enables to build a world representation, on which the effects of plans, grouping actions into control structures, will be transcribed by the simulation, resulting in a verdict and information about the plan's coherence.

Rutten, Eric↗

Design of a neutral thermal scattering (NeTS) module for hydrogen in light water

The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Project Title: Adaptive Vertical Grid Enhancement for E3SM (Final Technical Report)

During the first phase of the project (2018/05-2021/05), we proved the hypothesis that, without tuning, high vertical resolution improves the representation of unresolved, parameterized low altitude clouds such as stratocumulus clouds in global climate models. We evaluated the effectiveness of the newly developed computational framework, Framework for Improvement by Vertical Enhancement (FIVE) by implementing it into E3SM. FIVE aims to improve representation of the target cloud, low level cloud for this project, with use of a locally high vertical resolution for computation of selected processes. FIVE was assessed via direct comparison among E3SM, E3SM-FIVE, and satellite data as well as a newly developed ensemble artificial neural network methodology. We demonstrated that E3SM-FIVE remarkably improves marine stratocumulus representation. To advance the FIVE methodology, we took a first step toward the adaptive vertical grid in the FIVE method by modifying the E3SM single column model. For future application of FIVE to the global storm resolving model with a few km horizontal grid resolution, we developed the vertical-meridional two-dimensional (2D) Hadley circulation model, which allows one to run simulations with much higher vertical resolution and significantly longer duration than current three-dimensional global storm resolving model. During the second phase of the project (2021/05-2023/11; period for this technical report), we sought a way to appropriately transform FIVE to a regionalized version of FIVE, which allows one to select locations with/without FIVE. In parallel, we demonstrated that E3SM-FIVE with a 25-km mesh significantly reduced the bias associated with coastal stratocumulus. Furthermore, we showed that the same improvement was able to be obtained with the use of both the Regional Refined Model (RRM) targeted to one of the stratocumulus regions with 25-km resolution and FIVE. We explored a new implementation strategy of FIVE with our 2D Hadley circulation model. In parallel we applied the 2D Hadley circulation model to study the double intertropical convergence zone bias with storm resolving model resolution and aerosol-cloud interactions within the Hadley circulation with near large eddy simulation resolution. Last, we have been investigating Arctic mixed-phase stratocumulus clouds with large eddy simulations in order to improve their representation in global models. During the course of the entire project, 6 journal articles were published, and 3 manuscripts are currently in preparation. We gave a total of 30 presentations at various domestic and international conferences. In addition we gave 6 seminar talks at universities and a national laboratory. E3SM- FIVE is available as a branch of E3SM.

54 ENVIRONMENTAL SCIENCES↗

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

Role of the BIPM in UTC Dissemination to the Real Time User

The generation and dissemination of International Atomic Time (TAI), and Coordinated Universal Time (UTC) are explicitly mentioned in the list of principal tasks of the Bureau International des Poids et Mesures (BIPM), that appears in the Compes Rendus of the the 18e Conference Generales des Poids et Measures, in 1987. These time scales are used as the ultimate reference in the most demanding scientific applications and must, therefore, be of the best metrological quality in terms of reliability, long term stability, and conformity of the scale interval with the second, the unit of time of the International System of Units. To meet these requirements, it is necessary that the readings of the atomic clocks, spread all over the world, that are used as basic timing data for TAI and UTC generation, must be combined in the most efficient way possible. In particular, to take full advantage of the quality of each contributing clock calls for observation of its performance over a sufficiently long time. At present, the computation period treats data in blocks covering two months. TAI and UTC are thus deferred-time scales that cannot be immediately available to real-time users. The BIPM can, nevertheless be of help to real-time users. The predictability of UTC is a fundamental attribute of the scale for institutions responsible for the dissemination of real-time time scales. It allows them to improve their local representations of UTC and, thus, implement a more thorough steering of the time scales diffused in real-time. With a view to improving the predicatbility of UTC, the BIPM examines in detail timing techniques and basic theories in order to propose alternative solutions for timing algorithms. This, coupled with a recent improvement of timing data, makes UTC more stable and thus, more predictable. At a more practical level, effort is being devoted to putting in place automatic procedures for reducing the time needed for data collection and treatment: monthly results are already available ten days earlier than before.

Quinn, T. J.↗

Simulation of the Indian Summer Monsoon Using Comprehensive Atmosphere-land Interactions, in the Absence of Two-way Air-sea Interactions

Community Land Model version 2 (CLM2) as a comprehensive land surface model and a simple land surface model (SLM) were coupled to an atmospheric climate model to investigate the role of land surface processes in the development and the persistence of the South Asian summer monsoon. Two-way air-sea interactions were not considered in order to identify the reproducibility of the monsoon evolution by the comprehensive land model, which includes more realistic vertical soil moisture structures, vegetation and 2-way atmosphere-land interactions at hourly intervals. In the monsoon development phase (May and June). comprehensive land-surface treatment improves the representation of atmospheric circulations and the resulting convergence/divergence through the improvements in differential heating patterns and surface energy fluxes. Coupling with CLM2 also improves the timing and spatial distribution of rainfall maxima, reducing the seasonal rainfall overestimation by approx.60 % (1.8 mm/d for SLM, 0.7 mm/dI for CLM2). As for the interannual variation of the simulated rainfall, correlation coefficients of the Indian seasonal rainfall with observation increased from 0.21 (SLM) to 0.45 (CLM2). However, in the mature monsoon phase (July to September), coupling with the CLM2 does not exhibit a clear improvement. In contrast to the development phase, latent heat flux is underestimated and sensible heat flux and surface temperature over India are markedly overestimated. In addition, the moisture fluxes do not correlate well with lower-level atmospheric convergence, yielding correlation coefficients and root mean square errors worse than those produced by coupling with the SLM. A more realistic representation of the surface temperature and energy fluxes is needed to achieve an improved simulation for the mature monsoon period.

Monsoon Precipitation↗

Evaluation of Leading Modes of Climate Variability in the CMIP Archives

The adequate simulation of internal climate variability is key for our understanding of climate as it underpins efforts to attribute historical events, predict on seasonal and decadal time scales, and isolate the effects of climate change. Here the skill of models in reproducing observed modes of climate variability is assessed, both across and within the CMIP3, CMIP5, and CMIP6 archives, in order to document model capabilities, progress across ensembles, and persisting biases. A focus is given to the well-observed tropical and extratropical modes that exhibit small intrinsic variability relative to model structural uncertainty. These include El Niño–Southern Oscillation (ENSO), the Pacific decadal oscillation (PDO), the North Atlantic Oscillation (NAO), and the northern and southern annular modes (NAM and SAM). Significant improvements are identified in models’ representation of many modes. Canonical biases, which involve both amplitudes and patterns, are generally reduced across model generations. For example, biases in ENSO-related equatorial Pacific sea surface temperature, which extend too far westward, and associated atmospheric teleconnections, which are too weak, are reduced. Stronger tropical expression of the PDO in successive CMIP generations has characterized their improvement, with some CMIP6 models generating patterns that lie within the range of observed estimates. For the NAO, NAM, and SAM, pattern correlations with observations are generally higher than for other modes and slight improvements are identified across successive model generations. Finally, for ENSO and PDO spectra and extratropical modes, changes are small compared to internal variability, precluding definitive statements regarding improvement.

54 ENVIRONMENTAL SCIENCES↗

Multidisciplinary Design Optimization and Cruise Mach Number Study of Truss-Braced Wing Aircraft

The Subsonic Ultra Green Aircraft Research (SUGAR) Phase III was led by Dr. Rakesh K. Kapania and Dr. Joseph A. Schetz at the Multidisciplinary Analysis and Design Center for Advanced Vehicles, Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg VA. The research was performed from December 2014 to December 2015. Three major areas were investigated: Multidisciplinary Design Optimization (MDO) studies of truss braced wing (TBW) and strut braced wing (SBW) vehicles at cruise Mach numbers of 0.7 and 0.8 for a flight mission similar to current market single aisle configurations. The performance and the characteristics of the optimized vehicles were compared to the SUGAR Phase II TBW vehicle. These results were obtained without applying any of the extended transonic aerodynamic and aeroelastic tools that will be discussed later. It was found that the cruise Mach number has a large effect on the best “truss” configuration. At Mach 0.7, an SBW has a better fuel consumption and better take-off gross weight (TOGW). However, at Mach 0.8, the TBW is superior because the jury strut aids in satisfying the flutter constraint; Two-dimensional, steady, transonic aerodynamic analysis of the Boeing Airfoil J (BACJ) airfoil was performed for a range of thickness ratios, Mach numbers and lift coefficients. Reynolds-averaged Navier-Stokes (RANS) equations were solved to obtain the lift-curve slope, wave drag coefficient, the location of the center of pressure and to predict the separation at the trailing edge, which may lead to buffeting. One of the goals was to develop a database of lift-curve slope and the location of center of pressure, which could be used in a transonic aeroelastic analysis. Another goal was to compare the wave drag coefficients to those predicted by Lock’s fourth-power law and also to compare the transonic effects obtained from RANS simulations to those predicted by the Korn equations. A third goal was to develop a buffet boundary, which can be integrated into the MDO framework to prevent the optimized designs from probable buffeting; A state-space transonic aeroelastic analysis tool was developed, which can incorporate the nonlinear transonic effects in the unsteady aerodynamics but is yet computationally cheap when used within the VT MDO framework. The aeroelastic analysis uses Leishman- Beddoes (LB) indicial functions, which generated a state-space representation of the aeroelastic system. The indicial functions allow the incorporation of data for steady lift-curve slope and location of the center of pressure. Thus, the steady transonic effects are included, and the unsteady aerodynamic responses are a linearization about the steady results. The aeroelastic approach discretizes the wing into numerous strips, which results in a large eigenvalue problem as each strip has eight augmented aerodynamic states as per the LB theory. Thus, to reduce the computation expense, a reduced order model (ROM) was developed. The approach was validated using a few examples.

Kapania, Rakesh K.↗

Learning high-dimensional parametric maps via reduced basis adaptive residual networks

We propose a scalable framework for the learning of high-dimensional parametric maps via adaptively constructed residual network (ResNet) maps between reduced bases of the inputs and outputs. When just few training data are available, it is beneficial to have a compact parametrization in order to ameliorate the ill-posedness of the neural network training problem. By linearly restricting high-dimensional maps to informed reduced bases of the inputs, one can compress high-dimensional maps in a constructive way that can be used to detect appropriate basis ranks, equipped with rigorous error estimates. A scalable neural network learning framework is thus to learn the nonlinear compressed reduced basis mapping. Unlike the reduced basis construction, however, neural network constructions are not guaranteed to reduce errors by adding representation power, making it difficult to achieve good practical performance. Inspired by recent approximation theory that connects ResNets to sequential minimizing flows, we present an adaptive ResNet construction algorithm. This algorithm allows for depth-wise enrichment of the neural network approximation, in a manner that can achieve good practical performance by first training a shallow network and then adapting. We prove universal approximation of the associated neural network class for $L^2_v$ functions on compact sets. Our overall framework allows for constructive means to detect appropriate breadth and depth, and related compact parametrizations of neural networks, significantly reducing the need for architectural hyperparameter tuning. Numerical experiments for parametric PDE problems and a 3D CFD wing design optimization parametric map demonstrate that the proposed methodology can achieve remarkably high accuracy for limited training data, and outperformed other neural network strategies we compared against.

42 ENGINEERING↗

Modeling Phenological Controls on Carbon Dynamics in Dryland Sagebrush Ecosystems

Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbonfluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle,we must also be able to understand and model its effects on dryland vegetation. Dynamic Global VegetationModels (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduceseasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to bothtotal productivity and competitive interactions with other plants, errors in how process-based models representphenology hinder our ability to predict climate change impacts. This may be particularly problematic in drylandecosystems where many species have developed a complex phenology in response to seasonal variability in bothmoisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure ofexisting phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) andgross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized modelparameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, wemodified the model to include a semi-deciduous phenology type and added flexibility to the representation ofgrass phenology. While optimizing parameters reduced model bias, the largest gains in model performance wereassociated with the development of our new representation of phenology. This modified model was able to bettercapture seasonal patterns of both leaf area index (R2=0.75) and gross primary productivity (R2=0.84), thoughits ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new modelalso resulted in a more realistic outcome of modeled competition between grass and shrubs. These findingsdemonstrate the importance of improving how DGVMs represent phenology in order to accurately forecastclimate change impacts in dryland ecosystems.

Ecosystem model Phenology Parameter optimization E↗