Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling vegetation change”

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 181 records · Page 10

A Greening Future Elevates Flash Drought Risk in Northern Mid‐to‐High Latitudes

Flash droughts have become a growing concern, as they can emerge rapidly and increase the risk of crop failure. Although past studies have investigated the meteorological drivers and future changes of flash drought, why flash drought is more frequent over humid and vegetated regions remains underexplored. This study delves further into the mechanism by which vegetation regulates flash drought and its future change using observations from multiple data sets and large ensemble simulations from three Earth system models. On an interannual timescale, both observations and simulations show robust increases in flash drought frequency and a higher flash-to-sub-seasonal drought ratio during spring or antecedent conditions with dense vegetation, supporting the important role of vegetation in flash drought occurrence, especially in the northern mid-to-high latitudes. In the latter regions, the large ensemble simulations show robust increases in flash drought (e.g., 67% and 46% increases in Eastern U.S. and North Asia in 2050–2100 relative to 1950–2000 under the high emission scenario), where the growing season is lengthening. Although greening might suggest reduced drought stress, it drives precipitation-soil moisture-evapotranspiration decoupling by increasing evapotranspiration partitioning to transpiration. As transpiration can access deep soil water through the plant root system, its increased portion can weaken the constraints of concurrent precipitation on evapotranspiration, thus accelerating soil moisture depletion under high evaporative demand, driving a slow-to-rapid drought transition. How vegetation regulates flash drought by regulating surface moisture budget is supported by observations and simulations. Although warming supports early planting, agriculture may increasingly be threatened by surging flash drought risk.

Drought↗

Global Carbon Budget 1800-2018 Loss of Additional Sink Capacity and Present Transient Difference

A dataset of two global carbon fluxes: the loss of additional sink capacity (LASC) and present versus transient difference (PTD). The loss of additional sink capacity is the loss of indirect anthropogenic sink capacity in land ecosystems caused by the loss of ecosystems that can increase carbon in response to changing environments, primarily forests. The present transient difference refers to the direct anthropogenic flux in land ecosystems (from land cover and land use change) and the difference in that flux when calculated assuming "present day" (actually more like 1980) ecosystem carbon stocks or when calculated assuming transient carbon stocks (i.e. caused by environmental change). In this dataset, the LASC and PTD fluxes were calculated using the TRENDY ensemble of Dynamic Global Vegetation Models (DGVMs; Sitch et al., 2024) following the methods described therein and in Obermeier et al. (2021). For further details see Obermeier et al. (2021), Friedlingstein et al. (2019), Sitch et al. (2024), and Walker et al. (2025). The dataset includes ensemble means and standard deviations for annual fluxes, fluxes cumulated annually over the whole time period, and cumulated since 1959. These data were originally collected to compare different approaches (DGVMs versus bookkeeping models) for estimating land use and land cover change (LULCC) emissions (Obermeier et al. 2021). More recently, and similarly, these data were used to harmonize the estimates of the direct and indirect anthropogenic fluxes in the global carbon cycle and to recalculate the global carbon budget (Walker et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

Expansion of high-latitude deciduous forests driven by interactions between climate warming and fire: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication. More information in the related references. The dataset is a model (ecosys) output used in a study that examined the potential future shifts in high-latitude vegetation types (PFTs) driven by changes in 21st century climate and fire. The study highlights the importance of mechanistically representing differences in plant functional traits on current and possible 21st century trajectories of Alaskan PFTs. This dataset package includes 13 files in pdf, csv, and nc formats.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Wildfire Risk Support via Satellite-Derived Vegetation Health and Land Surface Model Soil Moisture

Land Surface Model (LSM) and evaporative demand products provide advanced lead time to wildfire conditions that complement traditional fire indices and represent short-term changes that add context to overall, long-term drought conditions. Established fire indices typically use weather indicators (i.e., precipitation, temperature) and estimated dead fuel moisture to indirectly obtain land surface and sub-surface characterization. The convergence of LSM shallow and deep-layer soil moisture output and satellite-derived vegetation health combine to provide a tool for stakeholders to examine trends in the state of land surface conditions that can help assess wildfire threat. Satellite remote sensing data can constrain near-real time vegetation characteristics within the LSM and/or provide a derived stress index in order to help characterize the wildfire risk. In addition, a percentile product of soil moisture is derived from a comparison of current LSM conditions to the historical record of the LSM in order to put the current conditions in perspective relative to the season and geographic region. The 2015 season as well as the 2018 Camp Fire Complex in California were examined in terms of the changes in LSM soil moisture and vegetation states. Satellite vegetation health consistently showed decreases a month prior to wildfire initiation. Additionally, maximum changes in total column soil moisture corresponded with the greatest concentration of fire locations. While soil moisture deficits occurred in the shallow layers across northern California in 2018, significant deficits at all sub-surface levels were seen ahead of the Camp Fire event. This presentation will demonstrate the complementary value of LSM output and satellite measured vegetation health to diagnose short-term deficits in sub-surface soil moisture and the rapid decline in vegetation health which precedes large wildfire events.

Wildfire↗

Sources of atmospheric methane from wetlands

Through the use of a computerized Geographic Information Service (GIS), field data on methane emissions and water inundation levels have been combined with remotely sensed vegetation cover data. This permits a more precise extrapolation of point flux measurements to regional scale flux estimates. Since the GIS allows changes in environmental variables to be included in the model, the sensitivity of emissions on large scales to changes in these parameters may be studied. Here, methane flux measurements from the Everglades marsh complex are reported, giving the flux by vegetation class.

Bartlett, K. B.↗

AMSR2 Soil Moisture Downscaling Using Temperature and Vegetation Data

Soil moisture (SM) applications in terrestrial hydrology require higher spatial resolution soil moisture products than those provided by passive microwave remote sensing instruments (grid resolution of 9 km or larger). In this investigation, an innovative algorithm that uses visible/infrared remote sensing observations to downscale Advanced Microwave Scanning Radiometer 2 (AMSR2) coarse spatial resolution SM products was developed and implemented for use with data provided by the Advanced Microwave Scanning Radiometer 2 (AMSR2). The method is based on using the Normalized Difference Vegetation Index (NDVI) modulated relationships between day/night SM and temperature change at corresponding times. Land surface model output variables from the North America Land Data Assimilation System (NLDAS), remote sensing data from the Moderate-Resolution Imaging Spectroradiometer (MODIS), and Advanced Very High Resolution Radiometer (AVHRR) were used in this methodology. The functional relationships developed using NLDAS data at a grid resolution of 12.5 km were applied to downscale AMSR2 JAXA (Japan Aerospace Exploration Agency) SM product (25 km) using MODIS land surface temperature (LST) and NDVI observations (1 km) to produce the 1 km SM estimates. The downscaled SM estimates were validated by comparing them with ISMN (International Soil Moisture Network) in situ SM in the Black Bear-Red Rock watershed, central Oklahoma between 2015-2017. The overall statistical variables of the downscaled AMSR2 SM validation R2, slope, RMSE and bias, demonstrate good accuracy. The downscaled SM better characterized the spatial and temporal variability of SM at watershed scales than the original SM product.

AMSR2; passive microwave soil moisture; soil moist↗

The Economics of Farming Expansion in the Brazilian Cerrado Under Possible Effects of Climate Change

This analysis assesses the financial viability of legally investing in native Cerrado vegetation deforestation for crop production, considering climate change. The study uses data from twelve different crop models based on three different climate models to predict potential future crop yields in cleared land for growing soy and maize. The outcomes show that in many micro-regions, investments in clearing land for crop production would destroy economic value, that is, generate a negative net present value because of low/negative and volatile cashflows driven primarily by future yields as affected by climate. Our analysis was carried out based on present agricultural practices and technology. As climate changes, farmers may adapt their practices, which can lead to more resilient and productive crops, or grow different crops, which could provide better returns on investment in clearing land than the ones resulting from our analysis. Despite various uncertainties, farmers, policy makers and financial institutions should be aware of the climatic and financial risks associated with land clearing in Brazil, mainly in micro-regions in which all scenarios resulted in negative outcomes in the investment analysis.

agricultural expansion↗

An Observation‐Driven Approach to Improve Vegetation Phenology in a Global Land Surface Model

An empirical model calibration approach is presented that aims to approximate missing biosphere processes in a global land surface model without the need for substantial model structural changes. The strategy is implemented here using the NASA Catchment‐CN land surface model and Moderate Resolution Imaging Spectroradiometer (MODIS) observations of the fraction of absorbed photosynthetically active radiation (FPAR). Existing plant functional types (PFTs) of the Catchment‐CN model are divided into three subtypes, based on the bias between the model‐simulated and MODIS‐observed FPAR. Separate sets of vegetation parameters for each subtype are then calibrated at a small number of grid cells with homogeneous, single‐PFT land cover, using MODIS FPAR reference observations from 2003 to 2009. The effectiveness of the empirical approach at improving the realism of modeled vegetation dynamics is investigated with two global model simulations for the period 2010–2016, one using the newly calibrated parameter values and the other using the original values. Globally, the calibrated parameters reduce the root mean square error (RMSE) of the modeled FPAR with respect to MODIS by 0.029 (∼10%) on average. In some regions, substantially larger RMSE reductions are achieved. RMSE reductions are primarily driven by model bias reductions, with neutral effects on the temporal correlation skill. While the empirical approach is suitable for achieving consistent model improvements, it is shown to be sensitive to the characteristics of the model error, specifically a dominance of the bias component in the case of Catchment‐CN. Ultimately, more fundamental model structural changes may be required to achieve better improvements.

J. Kolassa↗

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities

Abstract As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

54 ENVIRONMENTAL SCIENCES↗

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗

An Improved Analysis of Forest Carbon Dynamics using Data Assimilation

There are two broad approaches to quantifying landscape C dynamics - by measuring changes in C stocks over time, or by measuring fluxes of C directly. However, these data may be patchy, and have gaps or biases. An alternative approach to generating C budgets has been to use process-based models, constructed to simulate the key processes involved in C exchange. However, the process of model building is arguably subjective, and parameters may be poorly defined. This paper demonstrates why data assimilation (DA) techniques - which combine stock and flux observations with a dynamic model - improve estimates of, and provide insights into, ecosystem carbon (C) exchanges. We use an ensemble Kalman filter (EnKF) to link a series of measurements with a simple box model of C transformations. Measurements were collected at a young ponderosa pine stand in central Oregon over a 3-year period, and include eddy flux and soil C02 efflux data, litterfall collections, stem surveys, root and soil cores, and leaf area index data. The simple C model is a mass balance model with nine unknown parameters, tracking changes in C storage among five pools; foliar, wood and fine root pools in vegetation, and also fresh litter and soil organic matter (SOM) plus coarse woody debris pools. We nested the EnKF within an optimization routine to generate estimates from the data of the unknown parameters and the five initial conditions for the pools. The efficacy of the DA process can be judged by comparing the probability distributions of estimates produced with the EnKF analysis vs. those produced with reduced data or model alone. Using the model alone, estimated net ecosystem exchange of C (NEE)= -251 f 197g Cm-2 over the 3 years, compared with an estimate of -419 f 29gCm-2 when all observations were assimilated into the model. The uncertainty on daily measurements of NEE via eddy fluxes was estimated at 0.5gCm-2 day-1, but the uncertainty on assimilated estimates averaged 0.47 g Cm-2 day-1, and only exceeded 0.5gC m-2 day-1 on days where neither eddy flux nor soil efflux data were available. In generating C budgets, the assimilation process reduced the uncertainties associated with using data or model alone and the forecasts of NEE were statistically unbiased estimates. The results of the analysis emphasize the importance of time series as constraints. Occasional, rare measurements of stocks have limited use in constraining the estimates of other components of the C cycle. Long time series are particularly crucial for improving the analysis of pools with long time constants, such as SOM, woody biomass, and woody debris. Long-running forest stem surveys, and tree ring data, offer a rich resource that could be assimilated to provide an important constraint on C cycling of slow pools. For extending estimates of NEE across regions, DA can play a further important role, by assimilating remote-sensing data into the analysis of C cycles. We show, via sensitivity analysis, how assimilating an estimate of photosynthesis - which might be provided indirectly by remotely sensed data - improves the analysis of NEE.

Williams, Mathew↗

Quantifying Carbon Cycle Extremes and Attributing Their Causes Under Climate and Land Use and Land Cover Change From 1850 to 2300

The increasing atmospheric carbon dioxide (CO 2 ) mole fraction affects global climate through radiative (trapping longwave radiation) and physiological effects (reduction of plant transpiration). We use the simulations of the Community Earth System Model (CESM1-BGC) forced with Representative Concentration Pathway 8.5 to investigate climate-vegetation feedbacks from 1850 to the year 2300. Human-induced land use and land cover change (LULCC), through biogeochemical and biogeophysical processes, alter the climate and modify photosynthetic activity. The changing characteristics of extreme anomalies in photosynthesis, referred to as carbon cycle extremes, increase the uncertainty of terrestrial ecosystems to act as a net carbon sink. However, the role of LULCC in altering carbon cycle extremes under business-as-usual (continuously rising) CO 2 emissions is unknown. Here we show that LULCC magnifies the intensity, frequency, and extent of carbon cycle extremes, resulting in a net reduction in expected photosynthetic activity in the future. We found that large temporally contiguous negative carbon cycle extremes are due to a persistent decrease in soil moisture, which is triggered by declines in precipitation. With LULCC and global warming, vegetation exhibits increased vulnerability to hot and dry environmental conditions, increasing the frequency of fire events and resulting in considerable losses in photosynthetic activity. While most regions show strengthening of negative carbon cycle extremes, a few locations show a weakening effect driven by declining vegetation cover or benign climate conditions for photosynthesis. Increasing hot, dry, and fire-driven carbon cycle extremes are essential for improving carbon cycle modeling and estimation of ecosystem responses to LULCC and rising CO 2 mole fractions. Moreover, large aberrations in vegetation productivity represent potential and growing threats to human lives, wildlife, and food security.

54 ENVIRONMENTAL SCIENCES↗

Sampling Soil for Characterization and Site Description

The sampling scheme for soil characterization within the GLOBE program is uniquely different from the sampling methods of the other protocols. The strategy is based on an understanding of the 5 soil forming factors (parent material, climate, biota, topography, and time) at each study site, and how each of these interact to produce a soil profile with unique characteristics and unique input and control into the atmospheric, biological, and hydrological systems. Soil profile characteristics, as opposed to soil moisture and temperature, vegetative growth, and atmospheric and hydrologic conditions, change very slowly, depending on the parameter being measured, ranging from seasonally to many thousands of years. Thus, soil information, including profile description and lab analysis, is collected only one time for each profile at a site. These data serve two purposes: 1) to supplement existing spatial information about soil profile characteristics across the landscape at local, regional, and global scales, and 2) to provide specific information within a given area about the basic substrate to which elements within the other protocols are linked. Because of the intimate link between soil properties and these other environmental elements, the static soil properties at a given site are needed to accurately interpret and understand the continually changing dynamics of soil moisture and temperature, vegetation growth and phenology, atmospheric conditions, and chemistry and turbidity in surface waters. Both the spatial and specific soil information can be used for modeling purposes to assess and make predictions about global change.

Levine, Elissa↗

Evaluating models of climate and forest vegetation

Understanding how the biosphere may respond to increasing trace gas concentrations in the atmosphere requires models that contain vegetation responses to regional climate. Most of the processes ecologists study in forests, including trophic interactions, nutrient cycling, and disturbance regimes, and vital components of the world economy, such as forest products and agriculture, will be influenced in potentially unexpected ways by changing climate. These vegetation changes affect climate in the following ways: changing C, N, and S pools; trace gases; albedo; and water balance. The complexity of the indirect interactions among variables that depend on climate, together with the range of different space/time scales that best describe these processes, make the problems of modeling and prediction enormously difficult. These problems of predicting vegetation response to climate warming and potential ways of testing model predictions are the subjects of this chapter.

Clark, James S.↗

Application of remote sensing in regional scale estimates of vegetation carbon budgets: The Belfix project

A concept for coupling the remote sensing derived fraction of the absorbed photosynthetic active radiation (FAPAR) with a functional ecosystem model was developed. The study was named the Belfix procedure. The quantification of changes in carbon dynamics at the ecosystem level is a key issue in studies of global climatic change effects at the vegetation atmosphere interface. An operational procedure, for the determination of carbon fluxes at the regional scale (Belgian territory), is presented. The approach allows for the determination of the sink function of vegetation for carbon (dioxide). The phyto- and litter mass, photosynthetic assimilation, autotroph and heterotroph carbon fluxes and net ecosystem exchange (NEE) of carbon, were evaluated. The results suggest that a single solution can be obtained for ecosystem rates and states, applying an iterative procedure, based on minimizing the change in maximal seasonal green phytomass in function of yearly FAPAR temporal profiles. Total phytomass values obtained are in close range with those obtained by ground sampling.

Veroustraete, Frank↗

The impact of plant‐derived fire management prescriptions on fire‐responsive bird species

Abstract In fire‐prone regions, the occurrence of some faunal species is contingent on the presence of resources that arise through post‐fire plant succession. Through planned burning, managers can alter resource availability and aim to provide the conditions required to promote biodiversity. Understanding how species occurrence changes at different spatial and temporal scales after fire is essential to achieve this goal. However, many fire prescriptions are guided primarily by the responses of fire‐sensitive plants when setting tolerable fire intervals. This approach assumes that maintaining floristic diversity will satisfy the requirements of fauna. We surveyed bird species in two semi‐arid vegetation types across an environmental gradient in south‐eastern Australia. We conducted four surveys at each of 253 sites across a 75‐year chronosequence of time since fire and used generalized additive mixed models to examine changes in the occurrence of birds in response to time since fire. Model predictions were compared to plant‐derived fire prescriptions currently guiding fire management in the region. Time since fire was a significant predictor for 18 of 28 species modeled, in at least one vegetation type, over a gradient of 1.3° of latitude. We detected considerable variation in the responses of some species, both between vegetation types and geographically within a vegetation type. Our evaluation of plant‐derived fire prescriptions suggests that the intervals considered acceptable for maintaining floristic diversity may not be sustainable for populations of birds requiring longer unburnt vegetation, with 6 of the 12 species assessed attaining a mean occurrence probability of 20.3% by the minimum tolerable fire interval, and 57.3% by the maximum tolerable fire interval, in their respective vegetation types. Our findings highlight the potential vulnerability of fire‐responsive bird species if fire prescriptions are applied in a manner that fails to account for the slow development of habitat resources needed by some species, and the variation detected within the responses of species. This highlights the need for species‐specific data collected at an appropriate spatial scale to inform management plans.

Makdissi, Rhys↗

Development of A Composite Drought Indicator for Operational Drought Monitoring in the MENA Region

This paper presents the composite drought indicator (CDI) that Jordanian, Lebanese, Moroccan, and Tunisian government agencies now produce monthly to support operational drought management decision making, and it describes their iterative co-development processes. The CDI is primarily intended to monitor agricultural and ecological drought on a seasonal time scale. It uses remote sensing and modelled data inputs, and it reflects anomalies in precipitation, vegetation, soil moisture, and evapotranspiration. Following quantitative and qualitative validation assessments, engagements with policymakers, and consideration of agencies’ technical and institutional capabilities and constraints, we made changes to CDI input data, modelling procedures, and integration to tailor the system for each national context. We summarize validation results, drought modelling challenges and how we overcame them through CDI improvements, and we describe the monthly CDI production process and outputs. Finally, we synthesize procedural and technical aspects of CDI development and reflect on the constraints we faced as well as trade-offs made to optimize the CDI for operational monitoring to support policy decision-making—including aspects of salience, credibility, and legitimacy—within each national context.

Karim Bergaoui↗

Predicting spring green-up across diverse North American grasslands

Vegetation phenology influences many ecosystem and climate processes, such as carbon uptake and energy and water cycles. Thus, understanding drivers of vegetation phenology is crucial for predicting current and future impacts of climate change on ecological systems. Existing models can accurately predict the date of spring green-up in temperate forests but tend to perform poorly in grassland systems. We hypothesize this is because most do not incorporate water availability, a primary limiting factor for grassland plants. In this study, we used long-term datasets of digital imagery from the PhenoCam Network of 43 diverse North American grassland sites (195 site-years) to test existing spring phenology models, as well as develop several new models that incorporate precipitation or soil moisture (53 models). As a result, most of the new models performed substantially better, with the best model requiring sufficient accumulated precipitation followed by warm temperatures to trigger spring onset (root mean square error, RMSE, between predicted and observed dates = 16.0 days). Importantly, the best model performed well across all grassland types using a single set of parameters, from temperate to arid grasslands. Since plants are adapted to their local climates, model performance was further improved when parameters were independently optimized for four separate climate regions (RMSE = 10.4 days). Therefore, both sufficient precipitation and temperature are required for grassland green-up, but optimal thresholds vary by region. Running the top model with projected climate data (representative concentration pathway 8.5) suggests that, depending on the climate region, spring onset will occur up to 12 days earlier within 100 years in temperature-limited sites, but the trend is unclear for precipitation-limited sites (3.5 ± 8.0 days later). This new phenology model improves our ability to understand and predict grassland dynamics, with implications for both current and future ecosystem processes related to carbon and water cycling.

54 ENVIRONMENTAL SCIENCES↗