Magnitude of urban heat islands largely explained by climate and population.

Urban heat islands (UHIs) exacerbate the risk of heat-related mortality associated with global climate change. The intensity of UHIs varies with population size and mean annual precipitation, but a unifying explanation for this variation is lacking, and there are no geographically targeted guidelines for heat mitigation. Here we analyse summertime differences between urban and rural surface temperatures (ΔTs) worldwide and find a nonlinear increase in ΔTs with precipitation that is controlled by water or energy limitations on evapotranspiration and that modulates the scaling of ΔTs with city size. We introduce a coarse-grained model that links population, background climate, and UHI intensity, and show that urban-rural differences in evapotranspiration and convection efficiency are the main determinants of warming. The direct implication of these nonlinearities is that mitigation strategies aimed at increasing green cover and albedo are more efficient in dry regions, whereas the challenge of cooling tropical cities will require innovative solutions.

Intensity and frequency of extreme novel epidemics.

Observational knowledge of the epidemic intensity, defined as the number of deaths divided by global population and epidemic duration, and of the rate of emergence of infectious disease outbreaks is necessary to test theory and models and to inform public health risk assessment by quantifying the probability of extreme pandemics such as COVID-19. Despite its significance, assembling and analyzing a comprehensive global historical record spanning a variety of diseases remains an unexplored task. A global dataset of historical epidemics from 1600 to present is here compiled and examined using novel statistical methods to estimate the yearly probability of occurrence of extreme epidemics. Historical observations covering four orders of magnitude of epidemic intensity follow a common probability distribution with a slowly decaying power-law tail (generalized Pareto distribution, asymptotic exponent = -0.71). The yearly number of epidemics varies ninefold and shows systematic trends. Yearly occurrence probabilities of extreme epidemics, Py, vary widely: Py of an event with the intensity of the "Spanish influenza" (1918 to 1920) varies between 0.27 and 1.9% from 1600 to present, while its mean recurrence time today is 400 y (95% CI: 332 to 489 y). The slow decay of probability with epidemic intensity implies that extreme epidemics are relatively likely, a property previously undetected due to short observational records and stationary analysis methods. Using recent estimates of the rate of increase in disease emergence from zoonotic reservoirs associated with environmental change, we estimate that the yearly probability of occurrence of extreme epidemics can increase up to threefold in the coming decades.

A kernel-modulated SIR model for Covid-19 contagious spread from county to continent.

The tempo-spatial patterns of Covid-19 infections are a result of nested personal, societal, and political decisions that involve complicated epidemiological dynamics across overlapping spatial scales. High infection "hotspots" interspersed within regions where infections remained sporadic were ubiquitous early in the outbreak, but the spatial signature of the infection evolved to affect most regions equally, albeit with distinct temporal patterns. The sparseness of Covid-19 infections in the United States was analyzed at scales spanning from 10 to 2,600 km (county to continental scale). Spatial evolution of Covid-19 cases in the United States followed multifractal scaling. A rapid increase in the spatial correlation was identified early in the outbreak (March to April). Then, the increase continued at a slower rate and approached the spatial correlation of human population. Instead of adopting agent-based models that require tracking of individuals, a kernel-modulated approach is developed to characterize the dynamic spreading of disease in a multifractal distributed susceptible population. Multiphase Covid-19 epidemics were reasonably reproduced by the proposed kernel-modulated susceptible-infectious-recovered (SIR) model. The work explained the fact that while the reproduction number was reduced due to nonpharmaceutical interventions (e.g., masks, social distancing, etc.), subsequent multiple epidemic waves still occurred; this was due to an increase in susceptible population flow following a relaxation of travel restrictions and corollary stay-at-home orders. This study provides an original interpretation of Covid-19 spread together with a pragmatic approach that can be imminently used to capture the spatial intermittency at all epidemiologically relevant scales while preserving the "disordered" spatial pattern of infectious cases.

Seasonal hysteresis of surface urban heat islands.

Temporal dynamics of urban warming have been extensively studied at the diurnal scale, but the impact of background climate on the observed seasonality of surface urban heat islands (SUHIs) remains largely unexplored. On seasonal time scales, the intensity of urban-rural surface temperature differences ([Formula: see text]) exhibits distinctive hysteretic cycles whose shape and looping direction vary across climatic zones. These observations highlight possible delays underlying the dynamics of the coupled urban-biosphere system. However, a general argument explaining the observed hysteretic patterns remains elusive. A coarse-grained model of SUHI coupled with a stochastic soil water balance is developed to demonstrate that the time lags between radiation forcing, air temperature, and rainfall generate a rate-dependent hysteresis, explaining the observed seasonal variations of [Formula: see text] If solar radiation is in phase with water availability, summer conditions cause strong SUHI intensities due to high rural evaporative cooling. Conversely, cities in seasonally dry regions where evapotranspiration is out of phase with radiation show a summertime oasis effect controlled by background climate and vegetation properties. These seasonal patterns of warming and cooling have significant implications for heat mitigation strategies as urban green spaces can reduce [Formula: see text] during summertime, while potentially negative effects of albedo management during winter are mitigated by the seasonality of solar radiation.

Peak grain forecasts for the US High Plains amid withering waters.

Irrigated agriculture contributes 40% of total global food production. In the US High Plains, which produces more than 50 million tons per year of grain, as much as 90% of irrigation originates from groundwater resources, including the Ogallala aquifer. In parts of the High Plains, groundwater resources are being depleted so rapidly that they are considered nonrenewable, compromising food security. When groundwater becomes scarce, groundwater withdrawals peak, causing a subsequent peak in crop production. Previous descriptions of finite natural resource depletion have utilized the Hubbert curve. By coupling the dynamics of groundwater pumping, recharge, and crop production, Hubbert-like curves emerge, responding to the linked variations in groundwater pumping and grain production. On a state level, this approach predicted when groundwater withdrawal and grain production peaked and the lag between them. The lags increased with the adoption of efficient irrigation practices and higher recharge rates. Results indicate that, in Texas, withdrawals peaked in 1966, followed by a peak in grain production 9 y later. After better irrigation technologies were adopted, the lag increased to 15 y from 1997 to 2012. In Kansas, where these technologies were employed concurrently with the rise of irrigated grain production, this lag was predicted to be 24 y starting in 1994. In Nebraska, grain production is projected to continue rising through 2050 because of high recharge rates. While Texas and Nebraska had equal irrigated output in 1975, by 2050, it is projected that Nebraska will have almost 10 times the groundwater-based production of Texas.

Does growing atmospheric CO2 explain increasing carbon sink in a boreal coniferous forest?

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001-2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [ ± 42] gC m-2a-1 ), and NEP increased at rate +6.4-7.0 gC m-2a-1 (or ca. +2.5% a-1 ) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d a-1 , and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct CO2 fertilization effect diminishes when moving from leaf to ecosystem, and only 30-40% of the observed ecosystem GPP increase could be attributed to CO2 . The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing CO2 and LAI, explaining the apparent proportionality between observed GPP and CO2 trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone.

Reduced ecosystem resilience quantifies fine-scale heterogeneity in tropical forest mortality responses to drought.

Sensitivity of forest mortality to drought in carbon-dense tropical forests remains fraught with uncertainty, while extreme droughts are predicted to be more frequent and intense. Here, the potential of temporal autocorrelation of high-frequency variability in Landsat Enhanced Vegetation Index (EVI), an indicator of ecosystem resilience, to predict spatial and temporal variations of forest biomass mortality is evaluated against in situ census observations for 64 site-year combinations in Costa Rican tropical dry forests during the 2015 ENSO drought. Temporal autocorrelation, within the optimal moving window of 24 months, demonstrated robust predictive power for in situ mortality (leave-one-out cross-validation R2  = 0.54), which allows for estimates of annual biomass mortality patterns at 30 m resolution. Subsequent spatial analysis showed substantial fine-scale heterogeneity of forest mortality patterns, largely driven by drought intensity and ecosystem properties related to plant water use such as forest deciduousness and topography. Highly deciduous forest patches demonstrated much lower mortality sensitivity to drought stress than less deciduous forest patches after elevation was controlled. Our results highlight the potential of high-resolution remote sensing to "fingerprint" forest mortality and the significant role of ecosystem heterogeneity in forest biomass resistance to drought.

Maximizing leaf carbon gain in varying saline conditions: An optimization model with dynamic mesophyll conductance.

While the adverse effects of elevated salinity levels on leaf gas exchange in many crops are not in dispute, representing such effects on leaf photosynthetic rates (A) continues to draw research attention. Here, an optimization model for stomatal conductance (gc ) that maximizes A while accounting for mesophyll conductance (gm ) was used to interpret new leaf gas exchange measurements collected for five irrigation water salinity levels. A function between chloroplastic CO2 concentration (cc ) and intercellular CO2 concentration (ci ) modified by salinity stress to estimate gm was proposed. Results showed that with increased salinity, the estimated gm and maximum photosynthetic capacity were both reduced, whereas the marginal water use efficiency λ increased linearly. Adjustments of gm , λ and photosynthetic capacity were shown to be consistent with a large corpus of drought-stress experiments. The inferred model parameters were then used to evaluate the combined effects of elevated salinity and atmospheric CO2 concentration (ca ) on leaf gas exchange. For a given salinity level, increasing ca increased A linearly, but these increases were accompanied by mild reductions in gc and transpiration. The ca level needed to ameliorate A reductions due to increased salinity is also discussed using the aforementioned model calculations.

Detecting forest response to droughts with global observations of vegetation water content.

Droughts in a warming climate have become more common and more extreme, making understanding forest responses to water stress increasingly pressing. Analysis of water stress in trees has long focused on water potential in xylem and leaves, which influences stomatal closure and water flow through the soil-plant-atmosphere continuum. At the same time, changes of vegetation water content (VWC) are linked to a range of tree responses, including fluxes of water and carbon, mortality, flammability, and more. Unlike water potential, which requires demanding in situ measurements, VWC can be retrieved from remote sensing measurements, particularly at microwave frequencies using radar and radiometry. Here, we highlight key frontiers through which VWC has the potential to significantly increase our understanding of forest responses to water stress. To validate remote sensing observations of VWC at landscape scale and to better relate them to data assimilation model parameters, we introduce an ecosystem-scale analog of the pressure-volume curve, the non-linear relationship between average leaf or branch water potential and water content commonly used in plant hydraulics. The sources of variability in these ecosystem-scale pressure-volume curves and their relationship to forest response to water stress are discussed. We further show to what extent diel, seasonal, and decadal dynamics of VWC reflect variations in different processes relating the tree response to water stress. VWC can also be used for inferring belowground conditions-which are difficult to impossible to observe directly. Lastly, we discuss how a dedicated geostationary spaceborne observational system for VWC, when combined with existing datasets, can capture diel and seasonal water dynamics to advance the science and applications of global forest vulnerability to future droughts.

Universal Return to Isotropy of Inhomogeneous Atmospheric Boundary Layer Turbulence.

A recalcitrant problem in the physics of turbulence is the representation of the tendency of large-scale anisotropic eddies to redistribute their energy content with decreasing scales, a phenomenon referred to as return to isotropy. An unprecedented dataset of atmospheric turbulence measurements covering flat to mountainous terrain, stratification spanning convective to very stable conditions, surface roughness ranging over several orders of magnitude, and Reynolds numbers that far exceed the limits of direct numerical simulations and laboratory experiments was assembled for the first time and used to explore the scalewise return to isotropy. The multiple routes to energy equipartitioning among velocity components are shown to be universal once the initial anisotropy at large scales, linked to turbulence generation, is accounted for.

A Multiscale Approach to Timescale Analysis: Isolating Diel Signals from Solute Concentration Time Series.

Solute concentration time series reflect hydrological and biological drivers through various frequencies, phases, and amplitudes of change. Untangling these signals facilitates the understanding of dynamic ecosystem conditions and transient water quality issues. A case in point is the inference of biogeochemical processes from diel solute concentration variations. This analysis requires approaches capable of isolating subtle diel signals from background variability at other scales. Conventional time series analyses typically assume stationary or deterministic background variability; however, most rivers do not respect such niceties. We developed a time-series filtering method that uses empirical mode decomposition to decompose a measured solute concentration time series into intrinsic mode frequencies. Based on externally supplied mechanistic knowledge, we then filter these modes by periodicity, phase, and coherence with neighboring days. This method is tested on three synthetic series that incorporate environmental variability and sensor noise and on a year of 15 min sampled concentration time series from three hydrologically and ecologically distinct rivers in the eastern United States. The proposed method successfully isolated signals in the measured data sets that corresponded with variability in gross primary productivity. The strength the diel signal isolated through this method was smaller compared to the true signal in the synthetic series; however, uncertainty analysis showed that the process-model-based estimates derived from these signals were similar to other inference methods. This signal decomposition method retains information that can be used for further process modeling while making different assumptions about the data than Fourier and wavelet analyses.

Population agglomeration is a harbinger of the spatial complexity of COVID-19.

The spatial template over which COVID-19 infections operate is a result of nested societal decisions involving complex political and epidemiological processes at a broad range of spatial scales. It is characterized by 'hotspots' of high infections interspersed within regions where infections are sporadic to absent. In this work, the sparseness of COVID-19 infections and their time variations were analyzed across the US at scales ranging from 10 km (county scale) to 2600 km (continental scale). It was found that COVID-19 cases are multi-scaling with a multifractality kernel that monotonically approached that of the underlying population. The spatial correlation of infections between counties increased rapidly in March 2020; that rise continued but at a slower pace subsequently, trending towards the spatial correlation of the population agglomeration. This shows that the disease had already spread across the USA in early March such that travel restriction thereafter (starting on March 15th 2020) had minor impact on the subsequent spatial propagation of COVID-19. The ramifications of targeted interventions on spatial patterns of new infections were explored using the epidemiological susceptible-infectious-recovered (SIR) model mapped onto the population agglomeration template. These revealed that re-opening rural areas would have a smaller impact on the spread and evolution of the disease than re-opening urban (dense) centers which would disturb the system for months. This study provided a novel way for interpreting the spatial spread of COVID-19, along with a practical approach (multifractals/SIR/spectral slope) that could be employed to capture the variability and intermittency at all scales while maintaining the spatial structure.

Increasing atmospheric humidity and CO2 concentration alleviate forest mortality risk.

Climate-induced forest mortality is being increasingly observed throughout the globe. Alarmingly, it is expected to exacerbate under climate change due to shifting precipitation patterns and rising air temperature. However, the impact of concomitant changes in atmospheric humidity and CO2 concentration through their influence on stomatal kinetics remains a subject of debate and inquiry. By using a dynamic soil-plant-atmosphere model, mortality risks associated with hydraulic failure and stomatal closure for 13 temperate and tropical forest biomes across the globe are analyzed. The mortality risk is evaluated in response to both individual and combined changes in precipitation amounts and their seasonal distribution, mean air temperature, specific humidity, and atmospheric CO2 concentration. Model results show that the risk is predicted to significantly increase due to changes in precipitation and air temperature regime for the period 2050-2069. However, this increase may largely get alleviated by concurrent increases in atmospheric specific humidity and CO2 concentration. The increase in mortality risk is expected to be higher for needleleaf forests than for broadleaf forests, as a result of disparity in hydraulic traits. These findings will facilitate decisions about intervention and management of different forest types under changing climate.

Xylem-phloem hydraulic coupling explains multiple osmoregulatory responses to salt stress.

Salinity is known to affect plant productivity by limiting leaf-level carbon exchange, root water uptake, and carbohydrates transport in the phloem. However, the mechanisms through which plants respond to salt exposure by adjusting leaf gas exchange and xylem-phloem flow are still mostly unexplored. A physically based model coupling xylem, leaf, and phloem flows is here developed to explain different osmoregulation patterns across species. Hydraulic coupling is controlled by leaf water potential, ψl , and determined under four different maximization hypotheses: water uptake (1), carbon assimilation (2), sucrose transport (3), or (4) profit function - i.e. carbon gain minus hydraulic risk. All four hypotheses assume that finite transpiration occurs, providing a further constraint on ψl . With increasing salinity, the model captures different transpiration patterns observed in halophytes (nonmonotonic) and glycophytes (monotonically decreasing) by reproducing the species-specific strength of xylem-leaf-phloem coupling. Salt tolerance thus emerges as plant's capability of differentiating between salt- and drought-induced hydraulic risk, and to regulate internal flows and osmolytes accordingly. Results are shown to be consistent across optimization schemes (1-3) for both halophytes and glycophytes. In halophytes, however, profit-maximization (4) predicts systematically higher ψl than (1-3), pointing to the need of an updated definition of hydraulic cost for halophytes under saline conditions.

The effect of plant water storage on water fluxes within the coupled soil-plant system.

In addition to buffering plants from water stress during severe droughts, plant water storage (PWS) alters many features of the spatio-temporal dynamics of water movement in the soil-plant system. How PWS impacts water dynamics and drought resilience is explored using a multi-layer porous media model. The model numerically resolves soil-plant hydrodynamics by coupling them to leaf-level gas exchange and soil-root interfacial layers. Novel features of the model are the considerations of a coordinated relationship between stomatal aperture variation and whole-system hydraulics and of the effects of PWS and nocturnal transpiration (Fe,night) on hydraulic redistribution (HR) in the soil. The model results suggest that daytime PWS usage and Fe,night generate a residual water potential gradient (Δψp,night) along the plant vascular system overnight. This Δψp,night represents a non-negligible competing sink strength that diminishes the significance of HR. Considering the co-occurrence of PWS usage and HR during a single extended dry-down, a wide range of plant attributes and environmental/soil conditions selected to enhance or suppress plant drought resilience is discussed. When compared with HR, model calculations suggest that increased root water influx into plant conducting-tissues overnight maintains a more favorable water status at the leaf, thereby delaying the onset of drought stress.

Do the energy fluxes and surface conductance of boreal coniferous forests in Europe scale with leaf area?

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (Gs ), water- and light-use efficiency and surface-atmosphere coupling of European boreal coniferous forests was explored using eddy-covariance (EC) energy and CO2 fluxes. The observed scaling relations were then explained using a biophysical multilayer soil-vegetation-atmosphere transfer model as well as by a bulk Gs representation. The LAI variations significantly alter radiation regime, within-canopy microclimate, sink/source distributions of CO2 , H2 O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem-scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry-canopy evapotranspiration (ET) was reasonably 'conservative' over the studied LAI range 0.5-7.0 m2 m-2 . Both ET and Gs experienced a minimum in the LAI range 1-2 m2 m-2 caused by opposing nonproportional response of stomatally controlled transpiration and 'free' forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m2 m-2 ) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m2 m-2 ). This finding emphasizes the significance of stand-replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light-saturated water-use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests.

Increasing water use efficiency along the C3 to C4 evolutionary pathway: a stomatal optimization perspective.

C4 photosynthesis evolved independently numerous times, probably in response to declining atmospheric CO2 concentrations, but also to high temperatures and aridity, which enhance water losses through transpiration. Here, the environmental factors controlling stomatal behaviour of leaf-level carbon and water exchange were examined across the evolutionary continuum from C3 to C4 photosynthesis at current (400 µmol mol(-1)) and low (280 µmol mol(-1)) atmospheric CO2 conditions. To this aim, a stomatal optimization model was further developed to describe the evolutionary continuum from C3 to C4 species within a unified framework. Data on C3, three categories of C3-C4 intermediates, and C4 Flaveria species were used to parameterize the stomatal model, including parameters for the marginal water use efficiency and the efficiency of the CO2-concentrating mechanism (or C4 pump); these two parameters are interpreted as traits reflecting the stomatal and photosynthetic adjustments during the C3 to C4 transformation. Neither the marginal water use efficiency nor the C4 pump strength changed significantly from C3 to early C3-C4 intermediate stages, but both traits significantly increased between early C3-C4 intermediates and the C4-like intermediates with an operational C4 cycle. At low CO2, net photosynthetic rates showed continuous increases from a C3 state, across the intermediates and towards C4 photosynthesis, but only C4-like intermediates and C4 species (with an operational C4 cycle) had higher water use efficiencies than C3 Flaveria. The results demonstrate that both the marginal water use efficiency and the C4 pump strength increase in C4 Flaveria to improve their photosynthesis and water use efficiency compared with C3 species. These findings emphasize that the advantage of the early intermediate stages is predominantly carbon based, not water related.

On the complementary relationship between marginal nitrogen and water-use efficiencies among Pinus taeda leaves grown under ambient and CO2-enriched environments.

BACKGROUND AND AIMS: Water and nitrogen (N) are two limiting resources for biomass production of terrestrial vegetation. Water losses in transpiration (E) can be decreased by reducing leaf stomatal conductance (g(s)) at the expense of lowering CO(2) uptake (A), resulting in increased water-use efficiency. However, with more N available, higher allocation of N to photosynthetic proteins improves A so that N-use efficiency is reduced when g(s) declines. Hence, a trade-off is expected between these two resource-use efficiencies. In this study it is hypothesized that when foliar concentration (N) varies on time scales much longer than g(s), an explicit complementary relationship between the marginal water- and N-use efficiency emerges. Furthermore, a shift in this relationship is anticipated with increasing atmospheric CO(2) concentration (c(a)). METHODS: Optimization theory is employed to quantify interactions between resource-use efficiencies under elevated c(a) and soil N amendments. The analyses are based on marginal water- and N-use efficiencies, λ = (∂A/∂g(s))/(∂E/∂g(s)) and η = ∂A/∂N, respectively. The relationship between the two efficiencies and related variation in intercellular CO(2) concentration (c(i)) were examined using A/c(i) curves and foliar N measured on Pinus taeda needles collected at various canopy locations at the Duke Forest Free Air CO(2) Enrichment experiment (North Carolina, USA). KEY RESULTS: Optimality theory allowed the definition of a novel, explicit relationship between two intrinsic leaf-scale properties where η is complementary to the square-root of λ. The data support the model predictions that elevated c(a) increased η and λ, and at given c(a) and needle age-class, the two quantities varied among needles in an approximately complementary manner. CONCLUSIONS: The derived analytical expressions can be employed in scaling-up carbon, water and N fluxes from leaf to ecosystem, but also to derive transpiration estimates from those of η, and assist in predicting how increasing c(a) influences ecosystem water use.

Spread of North American wind-dispersed trees in future environments.

Despite ample research, understanding plant spread and predicting their ability to track projected climate changes remain a formidable challenge to be confronted. We modelled the spread of North American wind-dispersed trees in current and future (c. 2060) conditions, accounting for variation in 10 key dispersal, demographic and environmental factors affecting population spread. Predicted spread rates vary substantially among 12 study species, primarily due to inter-specific variation in maturation age, fecundity and seed terminal velocity. Future spread is predicted to be faster if atmospheric CO(2) enrichment would increase fecundity and advance maturation, irrespective of the projected changes in mean surface windspeed. Yet, for only a few species, predicted wind-driven spread will match future climate changes, conditioned on seed abscission occurring only in strong winds and environmental conditions favouring high survival of the farthest-dispersed seeds. Because such conditions are unlikely, North American wind-dispersed trees are expected to lag behind the projected climate range shift.

Mean velocity profile in a sheared and thermally stratified atmospheric boundary layer.

A stability correction function φ(m)(ζ) that accounts for distortions to the logarithmic mean velocity profile (MVP) in the lower atmosphere caused by thermal stratification was proposed by Monin and Obukhov in the 1950s using dimensional analysis. Its universal character was established from many field experiments. However, theories that describe the canonical shape of φ(m)(ζ) are still lacking. A previous link between the spectrum of turbulence and the MVP is expanded here to include the effects of thermal stratification on the turbulent kinetic energy dissipation rate and eddy-size anisotropy. The resulting theory provides a novel explanation for the power-law exponents and coefficients already reported for φ(m)(ζ) from numerous field experiments.