Land surface conductance linked to precipitation: Co-evolution of vegetation and climate in Earth system models.

Vegetation and precipitation are known to fundamentally influence each other. However, this interdependence is not fully represented in climate models because the characteristics of land surface (canopy) conductance to water vapor and CO2 are determined independently of precipitation. Working within a coupled atmosphere and land modelling framework (CAM6/CLM5; coupled Community Atmosphere Model v6/Community Land Model v5), we have developed a new theoretical approach to characterizing land surface conductance by explicitly linking its dynamic properties to local precipitation, a robust proxy for moisture available to vegetation. This will enable regional surface conductance characteristics to shift fluidly with climate change in simulations, consistent with general principles of co-evolution of vegetation and climate. Testing within the CAM6/CLM5 framework shows that climate simulations incorporating the new theory outperform current default configurations across several error metrics for core output variables when measured against observational data. In climate simulations for the end of this century the new, adaptive stomatal conductance scheme provides a revised prognosis for average and extreme temperatures over several large regions, with increased primary productivity through central and east Asia, and higher rainfall through North Africa and the Middle East. The new projections also reveal more frequent heatwaves than originally estimated for the south-eastern US and sub-Saharan Africa but less frequent heatwaves across east Europe and northeast Asia. These developments have implications for evaluating food security and risks from extreme temperatures in areas that are vulnerable to climate change.

Large model structural uncertainty in global projections of urban heat waves.

Urban heat waves (UHWs) are strongly associated with socioeconomic impacts. Here, we use an urban climate emulator combined with large ensemble global climate simulations to show that, at the urban scale a large proportion of the variability results from the model structural uncertainty in projecting UHWs in the coming decades under climate change. Omission of this uncertainty would considerably underestimate the risk of UHW. Results show that, for cities in four high-stake regions - the Great Lakes of North America, Southern Europe, Central India, and North China - a virtually unlikely (0.01% probability) UHW projected by single-model ensembles is estimated by our model with probabilities of 23.73%, 4.24%, 1.56%, and 14.76% respectively in 2061-2070 under a high-emission scenario. Our findings suggest that for urban-scale extremes, policymakers and stakeholders will have to plan for larger uncertainties than what a single model predicts if decisions are informed based on urban climate simulations.

Beyond Static Benchmarking: Using Experimental Manipulations to Evaluate Land Model Assumptions.

Land models are often used to simulate terrestrial responses to future environmental changes, but these models are not commonly evaluated with data from experimental manipulations. Results from experimental manipulations can identify and evaluate model assumptions that are consistent with appropriate ecosystem responses to future environmental change. We conducted simulations using three coupled carbon-nitrogen versions of the Community Land Model (CLM, versions 4, 4.5, and-the newly developed-5), and compared the simulated response to nitrogen (N) and atmospheric carbon dioxide (CO2) enrichment with meta-analyses of observations from similar experimental manipulations. In control simulations, successive versions of CLM showed a poleward increase in gross primary productivity and an overall bias reduction, compared to FLUXNET-MTE observations. Simulations with N and CO2 enrichment demonstrate that CLM transitioned from a model that exhibited strong nitrogen limitation of the terrestrial carbon cycle (CLM4) to a model that showed greater responsiveness to elevated concentrations of CO2 in the atmosphere (CLM5). Overall, CLM5 simulations showed better agreement with observed ecosystem responses to experimental N and CO2 enrichment than previous versions of the model. These simulations also exposed shortcomings in structural assumptions and parameterizations. Specifically, no version of CLM captures changes in plant physiology, allocation, and nutrient uptake that are likely important aspects of terrestrial ecosystems' responses to environmental change. These highlight priority areas that should be addressed in future model developments. Moving forward, incorporating results from experimental manipulations into model benchmarking tools that are used to evaluate model performance will help increase confidence in terrestrial carbon cycle projections.

Comparing optimal and empirical stomatal conductance models for application in Earth system models.

Earth system models (ESMs) rely on the calculation of canopy conductance in land surface models (LSMs) to quantify the partitioning of land surface energy, water, and CO2 fluxes. This is achieved by scaling stomatal conductance, gw , determined from physiological models developed for leaves. Traditionally, models for gw have been semi-empirical, combining physiological functions with empirically determined calibration constants. More recently, optimization theory has been applied to model gw in LSMs under the premise that it has a stronger grounding in physiological theory and might ultimately lead to improved predictive accuracy. However, this premise has not been thoroughly tested. Using original field data from contrasting forest systems, we compare a widely used empirical type and a more recently developed optimization-type gw model, termed BB and MED, respectively. Overall, we find no difference between the two models when used to simulate gw from photosynthesis data, or leaf gas exchange from a coupled photosynthesis-conductance model, or gross primary productivity and evapotranspiration for a FLUXNET tower site with the CLM5 community LSM. Field measurements reveal that the key fitted parameters for BB and MED, g1B and g1M, exhibit strong species specificity in magnitude and sensitivity to CO2 , and CLM5 simulations reveal that failure to include this sensitivity can result in significant overestimates of evapotranspiration for high-CO2 scenarios. Further, we show that g1B and g1M can be determined from mean ci /ca (ratio of leaf intercellular to ambient CO2 concentration). Applying this relationship with ci /ca values derived from a leaf δ13 C database, we obtain a global distribution of g1B and g1M , and these values correlate significantly with mean annual precipitation. This provides a new methodology for global parameterization of the BB and MED models in LSMs, tied directly to leaf physiology but unconstrained by spatial boundaries separating designated biomes or plant functional types.

Classifying heatwaves: Developing health-based models to predict high-mortality versus moderate United States heatwaves.

Heatwaves are divided between moderate, more common heatwaves and rare "high-mortality" heatwaves that have extremely large health effects per day, which we define as heatwaves with a 20% or higher increase in mortality risk. Better projections of the expected frequency of and exposure to these separate types of heatwaves could help communities optimize heat mitigation and response plans and gauge the potential benefits of limiting climate change. Whether a heatwave is high-mortality or moderate could depend on multiple heatwave characteristics, including intensity, length, and timing. We created heatwave classification models using a heatwave training dataset created using recent (1987-2005) health and weather data from 82 large US urban communities. We built twenty potential classification models and used Monte Carlo cross-validations to evaluate these models. We ultimately identified several models that can adequately classify high-mortality heatwaves. These models can be used to project future trends in high-mortality heatwaves under different scenarios of a changing future (e.g., climate change, population change). Further, these models are novel in the way they allow exploration of different scenarios of adaptation to heat, as they include, as predictive variables, heatwave characteristics that are measured relative to a community's temperature distribution, allowing different adaptation scenarios to be explored by selecting alternative community temperature distributions. The three selected models have been placed on GitHub for use by other researchers, and we use them in a companion paper to project trends in high-mortality heatwaves under different climate, population, and adaptation scenarios.

Projected trends in high-mortality heatwaves under different scenarios of climate, population, and adaptation in 82 US communities.

Some rare heatwaves have extreme daily mortality impacts; moderate heatwaves have lower daily impacts but occur much more frequently at present and so account for large aggregated impacts. We applied health-based models to project trends in high-mortality heatwaves, including proportion of all heatwaves expected to be high-mortality, using the definition that a high-mortality heatwave increases mortality risk by ≥20 %. We projected these trends in 82 US communities in 2061-2080 under two scenarios of climate change (RCP4.5, RCP8.5), two scenarios of population change (SSP3, SSP5), and three scenarios of community adaptation to heat (none, lagged, on-pace) for large- and medium-ensemble versions of the National Center for Atmospheric Research's Community Earth System Model. More high-mortality heatwaves were expected compared to present under all scenarios except on-pace adaptation, and population exposure was expected to increase under all scenarios. At least seven more high-mortality heatwaves were expected in a twenty-year period in the 82 study communities under RCP8.5 than RCP4.5 when assuming no adaptation. However, high-mortality heatwaves were expected to remain <1 % of all heatwaves and heatwave exposure under all scenarios. Projections were most strongly influenced by the adaptation scenario- going from a scenario of on-pace to lagged adaptation or from lagged to no adaptation more than doubled the projected number of and exposure to high-mortality heatwaves. Based on our results, fewer high-mortality heatwaves are expected when following RCP4.5 versus RCP8.5 and under higher levels of adaptation, but high-mortality heatwaves are expected to remain a very small proportion of total heatwave exposure.

Strong contributions of local background climate to urban heat islands.

The urban heat island (UHI), a common phenomenon in which surface temperatures are higher in urban areas than in surrounding rural areas, represents one of the most significant human-induced changes to Earth's surface climate. Even though they are localized hotspots in the landscape, UHIs have a profound impact on the lives of urban residents, who comprise more than half of the world's population. A barrier to UHI mitigation is the lack of quantitative attribution of the various contributions to UHI intensity (expressed as the temperature difference between urban and rural areas, ΔT). A common perception is that reduction in evaporative cooling in urban land is the dominant driver of ΔT (ref. 5). Here we use a climate model to show that, for cities across North America, geographic variations in daytime ΔT are largely explained by variations in the efficiency with which urban and rural areas convect heat to the lower atmosphere. If urban areas are aerodynamically smoother than surrounding rural areas, urban heat dissipation is relatively less efficient and urban warming occurs (and vice versa). This convection effect depends on the local background climate, increasing daytime ΔT by 3.0 ± 0.3 kelvin (mean and standard error) in humid climates but decreasing ΔT by 1.5 ± 0.2 kelvin in dry climates. In the humid eastern United States, there is evidence of higher ΔT in drier years. These relationships imply that UHIs will exacerbate heatwave stress on human health in wet climates where high temperature effects are already compounded by high air humidity and in drier years when positive temperature anomalies may be reinforced by a precipitation-temperature feedback. Our results support albedo management as a viable means of reducing ΔT on large scales.

Recent decline in the global land evapotranspiration trend due to limited moisture supply.

More than half of the solar energy absorbed by land surfaces is currently used to evaporate water. Climate change is expected to intensify the hydrological cycle and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land−a key diagnostic criterion of the effects of climate change and variability−remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network, meteorological and remote-sensing observations, and a machine-learning algorithm. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate.

Terrestrial gross primary production (GPP) is the largest global CO(2) flux driving several ecosystem functions. We provide an observation-based estimate of this flux at 123 +/- 8 petagrams of carbon per year (Pg C year(-1)) using eddy covariance flux data and various diagnostic models. Tropical forests and savannahs account for 60%. GPP over 40% of the vegetated land is associated with precipitation. State-of-the-art process-oriented biosphere models used for climate predictions exhibit a large between-model variation of GPP's latitudinal patterns and show higher spatial correlations between GPP and precipitation, suggesting the existence of missing processes or feedback mechanisms which attenuate the vegetation response to climate. Our estimates of spatially distributed GPP and its covariation with climate can help improve coupled climate-carbon cycle process models.

Surfactant-mediated modulation of hydrophobic interactions in associative polymer solutions containing cyclodextrin.

The ability of nonionic surfactants to modulate the rheological characteristics of comblike hydrophobically modified associative polymer solutions containing cyclodextrin (CD) is examined. Addition of either alpha- or beta-CD to these polymers results in a marked decrease in solution viscosity and viscoelastic properties because of the encapsulation of the polymer hydrophobes by CD. Nonionic surfactants, introduced to such a system, alter the hydrophobic interactions by competing with the polymer hydrophobes for complexation with the CDs. In this regard, nonylphenol ethoxylates (NPe) with different ethylene oxide chain lengths, which determine the hydrophilic-lipophilic balance (HLB) of the surfactant, are used. Our results reveal that the extent and rate of recovery of zero shear viscosity as well as dynamic moduli are strongly influenced by the type of CD (alpha versus beta) as well as the HLB of the surfactant. For polymer solutions containing alpha-CD, recovery is observed solely in the presence of a low-HLB surfactant (NP6 and NP8). Additionally, in the case of NP6, the viscosity increases monotonically above the original hydrophobically modified alkali-soluble emulsion viscosity with surfactant addition, whereas for that of a higher HLB surfactant (NP15), the viscosity shows no appreciable change. In the case of beta-CD, on the other hand, a complete recovery as well as further enhancement of rheological properties is achieved using the above surfactants. However, the trends of viscosity modulation are dissimilar for different surfactants as we witness a monotonous increase in the case of NP6 while a maximum in viscosity is observed in the presence of NP8 and NP15. The contrasting observations can be explained in terms of varying affinities of alpha- and beta-CDs to bind with NP surfactants and the existence of different micellar structures in solution as governed by surfactant HLB. These results are further confirmed by UV/vis spectroscopy and cloud point measurements.

The importance of land-cover change in simulating future climates.

Adding the effects of changes in land cover to the A2 and B1 transient climate simulations described in the Special Report on Emissions Scenarios (SRES) by the Intergovernmental Panel on Climate Change leads to significantly different regional climates in 2100 as compared with climates resulting from atmospheric SRES forcings alone. Agricultural expansion in the A2 scenario results in significant additional warming over the Amazon and cooling of the upper air column and nearby oceans. These and other influences on the Hadley and monsoon circulations affect extratropical climates. Agricultural expansion in the mid-latitudes produces cooling and decreases in the mean daily temperature range over many areas. The A2 scenario results in more significant change, often of opposite sign, than does the B1 scenario.

Regions of strong coupling between soil moisture and precipitation.

Previous estimates of land-atmosphere interaction (the impact of soil moisture on precipitation) have been limited by a lack of observational data and by the model dependence of computational estimates. To counter the second limitation, a dozen climate-modeling groups have recently performed the same highly controlled numerical experiment as part of a coordinated comparison project. This allows a multimodel estimation of the regions on Earth where precipitation is affected by soil moisture anomalies during Northern Hemisphere summer. Potential benefits of this estimation may include improved seasonal rainfall forecasts.