Search

Plant activity in semi-arid ecosystems is largely controlled by pulses of precipitation, making them particularly vulnerable to increased aridity that is expected with climate change. Simple bucket-model hydrology schemes in land surface models (LSMs) have had limited ability in accurately capturing semi-arid water stores and fluxes. Recent, more complex, LSM hydrology models have not been widely evaluated against semi-arid ecosystem in situ data. We hypothesize that the failure of older LSM versions to represent evapotranspiration, ET, in arid lands is because simple bucket models do not capture realistic fluctuations in upper-layer soil moisture. We therefore predict that including a discretized soil hydrology scheme based on a mechanistic description of moisture diffusion will result in an improvement in model ET when compared to data because the temporal variability of upper-layer soil moisture content better corresponds to that of precipitation inputs. To test this prediction, we compared ORCHIDEE LSM simulations from (1) a simple conceptual 2-layer bucket scheme with fixed hydraulic parameters and (2) an 11-layer discretized mechanistic scheme of moisture diffusion in unsaturated soil based on Richards equations, against daily and monthly soil moisture and ET observations, together with data-derived estimates of transpiration / evapotranspiration, T∕ET, ratios, from six semi-arid grass, shrub, and forest sites in the south-western USA. The 11-layer scheme also has modified calculations of surface runoff, water limitation, and resistance to bare soil evaporation, E, to be compatible with the more complex hydrology configuration. To diagnose remaining discrepancies in the 11-layer model, we tested two further configurations: (i) the addition of a term that captures bare soil evaporation resistance to dry soil; and (ii) reduced bare soil fractional vegetation cover. We found that the more mechanistic 11-layer model results in a better representation of the daily and monthly ET observations. We show that, as predicted, this is because of improved simulation of soil moisture in the upper layers of soil (top ∼ 10 cm). Some discrepancies between observed and modelled soil moisture and ET may allow us to prioritize future model development and the collection of additional data. Biases in winter and spring soil moisture at the forest sites could be explained by inaccurate soil moisture data during periods of soil freezing and/or underestimated snow forcing data. Although ET is generally well captured by the 11-layer model, modelled T∕ET ratios were generally lower than estimated values across all sites, particularly during the monsoon season. Adding a soil resistance term generally decreased simulated bare soil evaporation, E, and increased soil moisture content, thus increasing transpiration, T, and reducing the negative bias between modelled and estimated monsoon T∕ET ratios. This negative bias could also be accounted for at the low-elevation sites by decreasing the model bare soil fraction, thus increasing the amount of transpiring leaf area. However, adding the bare soil resistance term and decreasing the bare soil fraction both degraded the model fit to ET observations. Furthermore, remaining discrepancies in the timing of the transition from minimum T∕ET ratios during the hot, dry May–June period to high values at the start of the monsoon in July–August may also point towards incorrect modelling of leaf phenology and vegetation growth in response to monsoon rains. We conclude that a discretized soil hydrology scheme and associated developments improve estimates of ET by allowing the modelled upper-layer soil moisture to more closely match the pulse precipitation dynamics of these semi-arid ecosystems; however, the partitioning of T from E is not solved by this modification alone.

Many ecosystem processes that influence Earth system feedbacks – vegetation growth, water and nutrient cycling, disturbance regimes – are strongly influenced by multidecadal- to millennial-scale climate variations that cannot be directly observed. Paleoclimate records provide information about these variations, forming the basis of our understanding and modeling of them. Fossil pollen records are abundant in the NE US, but cannot simultaneously provide information about paleoclimate and past vegetation in a modeling context because this leads to circular logic. If pollen data are used to constrain past vegetation changes, then the remaining paleoclimate archives in the northeastern US (NE US) are quite limited. Nonetheless, a growing number of diverse reconstructions have been developed but have not yet been examined together. Here we conduct a systematic review, assessment, and comparison of paleotemperature and paleohydrological proxies from the NE US for the last 3000 years. Regional temperature reconstructions (primarily summer) show a long-term cooling trend (1000 BCE–1700 CE) consistent with hemispheric-scale reconstructions, while hydroclimate data show gradually wetter conditions through the present day. Multiple proxies suggest that a prolonged, widespread drought occurred between 550 and 750 CE. Dry conditions are also evident during the Medieval Climate Anomaly, which was warmer and drier than the Little Ice Age and drier than today. There is some evidence for an acceleration of the longer-term wetting trend in the NE US during the past century; coupled with an abrupt shift from decreasing to increasing temperatures in the past century, these changes could have wide-ranging implications for species distributions, ecosystem dynamics, and extreme weather events. More work is needed to gather paleoclimate data in the NE US to make inter-proxy comparisons and to improve estimates of uncertainty in reconstructions.

How carbon (C) is allocated to different plant tissues (leaves, stem, and roots) determines how long C remains in plant biomass and thus remains a central challenge for understanding the global C cycle. We used a diverse set of observations (AmeriFlux eddy covariance tower observations, biomass estimates from tree-ring data, and leaf area index (LAI) measurements) to compare C fluxes, pools, and LAI data with those predicted by a land surface model (LSM), the Community Land Model (CLM4.5). We ran CLM4.5 for nine temperate (including evergreen and deciduous) forests in North America between 1980 and 2013 using four different C allocation schemes: i. dynamic C allocation scheme (named "D-CLM4.5") with one dynamic allometric parameter, which allocates C to the stem and leaves to vary in time as a function of annual net primary production (NPP); ii. an alternative dynamic C allocation scheme (named "D-Litton"), where, similar to (i), C allocation is a dynamic function of annual NPP, but unlike (i) includes two dynamic allometric parameters involving allocation to leaves, stem, and coarse roots; iii.–iv. a fixed C allocation scheme with two variants, one representative of observations in evergreen (named "F-Evergreen") and the other of observations in deciduous forests (named "F-Deciduous"). D-CLM4.5 generally overestimated gross primary production (GPP) and ecosystem respiration, and underestimated net ecosystem exchange (NEE). In D-CLM4.5, initial aboveground biomass in 1980 was largely overestimated (between 10 527 and 12 897 g C m−2) for deciduous forests, whereas aboveground biomass accumulation through time (between 1980 and 2011) was highly underestimated (between 1222 and 7557 g C m−2) for both evergreen and deciduous sites due to a lower stem turnover rate in the sites than the one used in the model. D-CLM4.5 overestimated LAI in both evergreen and deciduous sites because the leaf C–LAI relationship in the model did not match the observed leaf C–LAI relationship at our sites. Although the four C allocation schemes gave similar results for aggregated C fluxes, they translated to important differences in long-term aboveground biomass accumulation and aboveground NPP. For deciduous forests, D-Litton gave more realistic Cstem ∕ Cleaf ratios and strongly reduced the overestimation of initial aboveground biomass and aboveground NPP for deciduous forests by D-CLM4.5. We identified key structural and parameterization deficits that need refinement to improve the accuracy of LSMs in the near future. These include changing how C is allocated in fixed and dynamic schemes based on data from current forest syntheses and different parameterization of allocation schemes for different forest types. Our results highlight the utility of using measurements of aboveground biomass to evaluate and constrain the C allocation scheme in LSMs, and suggest that stem turnover is overestimated by CLM4.5 for these AmeriFlux sites. Understanding the controls of turnover will be critical to improving long-term C processes in LSMs.