Onward but not always upward: individualistic elevational shifts of tree species in subtropical montane forests

Ongoing global climate change is driving widespread shifts in species distributions. Trends show frequent upwards shifts of treelines, but information on changes in montane forest below the treeline and in the tropics and subtropics is limited, despite the importance of these areas for biodiversity and ecosystem function. Patterns of species shifts in tropical and subtropical regions are likely to be more complex and individualistic than global averages suggest due to high species diversity and strong influence of competition, alongside direct climatic limitations on distributions. To address the question of how subtropical montane tree species are likely to move as climate changes, we used an extensive national forest inventory to estimate distribution shifts of 75 tree species in Taiwan by comparing the optimum elevation and range edges of adults and juveniles within species. Overall there was a significant difference in optimum elevation of adults and juveniles. Life stage mismatches suggested upward shifts in 35% of species but downward shifts of over half (56%), while 8% appeared stable. Upward elevation shifts were disproportionately common in high elevation species, whilst mid to low elevation species suggested greater variation in shift direction. Whilst previous research on mountain forest range shifts has been dominated by work addressing changes in treeline position, we show that although high elevation species shift up, below the treeline species may shift individualistically, heralding widespread changes in forest communities over coming decades. The wide variation of responses indicated is likely driven by individual species responses to interacting environmental factors such as competition, topography and anthropogenic influences across the broad range of forest types investigated. As global environmental changes continue, more detailed understanding of tree range shifts across a wide spectrum of forests will allow us to prepare for the implications of such changes for biodiversity, ecosystem function and dependent human populations.     

Anthropogenic‐driven rapid shifts in tree distribution lead to increased dominance of broadleaf species

Over the past century, major shifts in the geographic distribution of tree species have occurred in response to changes in land use and climate. We analyse species distribution and abundance from about 33 000 forest inventory plots in Spain sampled twice over a period of 10–12 years. We show a dominance of range contraction (extinction), and demographic decline over range expansion (colonization), with seven of 11 species exhibiting extinction downhill of their distribution. Contrary to expectations, these dynamics are not always consistent with climate warming over the study period, but result from legacies in forest structure due to past land use change and fire occurrence. We find that these changes have led to the expansion of broadleaf species (i.e. family Fagaceae) over areas formerly dominated by conifer species (i.e. family Pinaceae), due to the greater capacity of the former to respond to most disturbances and their higher competitive ability. This recent and rapid transition from conifers to broadleaves has important implications in forest dynamics and ecosystem services they provide. The finding raises the question as to whether the increasing dominance of relatively drought‐sensitive broadleaf species will diminish resilience of Mediterranean forests to very likely drier conditions in the future.     

Linking species thermal tolerance to elevational range shifts in upland dung beetles

Climate warming has been proposed as the main cause of the recent range shifts seen in many species. Although species' thermal tolerances are thought to play a key role in determining responses to climate change, especially in ectotherms, empirical evidence is still limited. We investigate the connection between species' thermal tolerances, elevational range and shifts in the lower elevational limit of dung beetle species (Coleoptera, Aphodiidea) in an upland region in the northwest of England. We measured thermal tolerances in the laboratory, and used current and historical distribution data to test specific hypotheses about the area's three dominant species, particularly the species most likely to suffer from warming: Agollinus lapponum. We found marked differences between species in their minimum and maximum thermal tolerance and in their elevational range and patterns of abundance. Overall, differences in thermal limits among species matched the abundance patterns along the elevation gradient expected if distributions were constrained by climate. Agollinus lapponum abundance increased with elevation and this species showed lower maximum and minimum thermal limits than Acrossus depressus, for which abundance declined with elevation. Consistent with lower tolerance to high temperature, we recorded an uphill retreat of the low elevation limit of A. lapponum (177 m over 57 yr) in line with the increase in summer temperature observed in the region over the same period. Moreover, this species has been replaced at low and mid‐elevations by the other two warm‐tolerant species (A. depressus and Agrilinus ater). Our results provide empirical evidence that species' thermal tolerance constrains elevational ranges and contributes to explain the observed responses to climate warming. A mechanistic understanding of how climate change directly affects species, such as the one presented here, will provide a robust base to inform predictions of how individual species and whole assemblages may change in the future.     

Potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum

Aim To identify potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum (LGM). Location Europe. Methods We developed a new approach combining past‐climate simulations with the concept of analogous climate space. Our index gives a continuous measure of the potential of a location to have acted as a source or a sink for species that have shifted their ranges since the LGM. High glacial source potential is indicated by LGM climatic conditions that are widespread now; high post‐glacial sink potential is indicated by current climatic conditions that were widespread at the LGM. The degree of isolation of source and sink areas was calculated as the median distance to areas with analogous climate conditions. Results We identified areas of high glacial source potential in the previously recognized refugial areas in the southern European peninsulas, but also in large areas in central‐western Europe. The most climatically isolated source areas were located in northern Spain, in north‐western Europe and in eastern Turkey. From here species would have had to cover substantial distances to find current climate conditions analogous to LGM conditions of these areas. Areas with high post‐glacial sink potential were mainly located in Fennoscandia and in central and south‐eastern Europe. Some of the most isolated sink areas were located in the Spanish highlands and around the Baltic Sea. Main conclusions Our species‐independent approach successfully identified previously recognized glacial refugial areas with high source potential for species range shifts in southern Europe and in addition highlighted other potential source areas in central Europe. This study offers new insights into how the distribution of past and current climatic conditions may have influenced past species range shifts and current large‐scale biodiversity patterns.     

Evolutionary consequences of climate‐induced range shifts in insects

Range shifts can rapidly create new areas of geographic overlap between formerly allopatric taxa and evidence is accumulating that this can affect species persistence. We review the emerging literature on the short‐ and long‐term consequences of these geographic range shifts. Specifically, we focus on the evolutionary consequences of novel species interactions in newly created sympatric areas by describing the potential (i) short‐term processes acting on reproductive barriers between species and (ii) long‐term consequences of range shifts on the stability of hybrid zones, introgression and ultimately speciation and extinction rates. Subsequently, we (iii) review the empirical literature on insects to evaluate which processes have been studied, and (iv) outline some areas that deserve increased attention in the future, namely the genomics of hybridisation and introgression, our ability to forecast range shifts and the impending threat from insect vectors and pests on biodiversity, human health and crop production. Our review shows that species interactions in de novo sympatric areas can be manifold, sometimes increasing and sometimes decreasing species diversity. A key issue that emerges is that climate‐induced hybridisations in insects are much more widespread than anticipated and that rising temperatures and increased anthropogenic disturbances are accelerating the process of species mixing. The existing evidence only shows the tip of the iceberg and we are likely to see many more cases of species mixing following range shifts in the near future.     

Upslope migration of Andean trees

Aim Climate change causes shifts in species distributions, or ‘migrations’. Despite the centrality of species distributions to biodiversity conservation, the demonstrated large migration of tropical plant species in response to climate change in the past, and the expected sensitivity of species distributions to modern climate change, no study has tested for modern species migrations in tropical plants. Here we conduct a first test of the hypothesis that increasing temperatures are causing tropical trees to migrate to cooler areas. Location Tropical Andes biodiversity hotspot, south‐eastern Peru, South America. Methods We use data from repeated (2003/04–2007/08) censuses of 14 1‐ha forest inventory plots spanning an elevational gradient from 950 to 3400 m in Manu National Park in south‐eastern Peru, to characterize changes in the elevational distributions of 38 Andean tree genera. We also analyse changes in the genus‐level composition of the inventory plots through time. Results We show that most tropical Andean tree genera shifted their mean distributions upslope over the study period and that the mean rate of migration is approximately 2.5–3.5 vertical metres upslope per year. Consistent with upward migrations we also find increasing abundances of tree genera previously distributed at lower elevations in the majority of study plots. Main conclusions These findings are in accord with the a priori hypothesis of upward shifts in species ranges due to elevated temperatures, and are potentially the first documented evidence of present‐day climate‐driven migrations in a tropical plant community. The observed mean rate of change is less than predicted from the temperature increases for the region, possibly due to the influence of changes in moisture or non‐climatic factors such as substrate, species interactions, lags in tree community response and/or dispersal limitations. Whatever the cause(s), continued slower‐than‐expected migration of tropical Andean trees would indicate a limited ability to respond to increased temperatures, which may lead to increased extinction risks with further climate change.     

Distributional migrations,expansions, and contractions of tropical plant species as revealed in dated herbarium records

Species are predicted to respond to global warming through ‘cold‐ward’ shifts in their geographic distributions due to encroachment into newly suitable habitats and/or dieback in areas that become climatically unsuitable. I conduct one of the first‐ever tests of this hypothesis for tropical plant species. I test for changes in the thermal distributions of 239 South American tropical plant species using dated herbarium records for specimens collected between 1970 and 2009. Supporting a priori predictions, many species (59%) exhibit some evidence of significant cold‐ward range shifts even after correcting for collection biases. Over 1/3 of species (35%) show significant cold‐ward movement in their hot thermal limits (mean rate of change = 0.022 °C yr^?1). Most of these species (85%; 30% of all study species) show no corresponding shift in their cold thermal limits. These unbalanced changes in the species’ thermal range limits may indicate species that are experiencing dieback due to their intolerance of rising temperatures coupled with an inability to expand into newly climatically suitable habitats. On the other hand, 25% of species show significant cold‐ward shifts in their cold thermal range limits (mean rate of change = 0.003 °C yr^?1), but 80% of these species (20% of all study species) show no corresponding shift in their hot thermal range limits. In these cases, the unbalanced shifts may indicate species that are able to ‘benefit’ under global warming, at least temporally, by both tolerating rising temperatures and expanding into new suitable habitat. An important ancillary result of this study is that the number of species exhibiting significant range shifts was greatly influenced by shifting collector biases. This highlights the need to account for biases when analyzing natural history records or other long‐term records.     

Future climate‐driven shifts in distribution of Calanus finmarchicus

Calanus finmarchicus is a key‐structural species of the North Atlantic polar biome. The species plays an important trophic role in subpolar and polar ecosystems as a grazer of phytoplankton and as a prey for higher trophic levels such as the larval stages of many fish species. Here, we used a recently developed ecological niche model to assess the ecological niche (sensu Hutchinson) of C. finmarchicus and characterize its spatial distribution. This model explained about 65% of the total variance of the observed spatial distribution inferred from an independent dataset (data of the continuous plankton recorder survey). Comparisons with other types of models (structured population and ecophysiological models) revealed a clear similarity between modeled spatial distributions at the scale of the North Atlantic. Contemporary models coupled with future projections indicated a progressive reduction of the spatial habitat of the species at the southern edge and a more pronounced one in the Georges Bank, the Scotian Shelf and the North Sea and a potential increase in abundance at the northern edge of its spatial distribution, especially in the Barents Sea. These major changes will probably lead to a major alteration of the trophodynamics of North Atlantic ecosystems affecting the trophodynamics and the biological carbon pump.     

Patterns of climate‐induced density shifts of species: poleward shifts faster in northern boreal birds than in southern birds

Climate change has been shown to cause poleward range shifts of species. These shifts are typically demonstrated using presence–absence data, which can mask the potential changes in the abundance of species. Moreover, changes in the mean centre of weighted density of species are seldom examined, and comparisons between these two methods are even rarer. Here, we studied the change in the mean weighted latitude of density (MWLD) of 94 bird species in Finland, northern Europe, using data covering a north–south gradient of over 1000 km from the 1970s to the 2010s. The MWLD shifted northward on average 1.26 km yr^?1, and this shift was significantly stronger in northern species compared to southern species. These shifts can be related to climate warming during the study period, because the annual temperature had increased more in northern Finland (by 1.7 °C) than in southern Finland (by 1.4 °C), although direct causal links cannot be shown. Density shifts of species distributed over the whole country did not differ from shifts in species situated on the edge of the species range in southern and northern species. This means that density shifts occur both in the core and on the edge of species distribution. The species‐specific comparison of MWLD values with corresponding changes in the mean weighted latitude using presence–absence atlas data (MWL) revealed that the MWLD moved more slowly than the MWL in the atlas data in the southern species examined, but more rapidly in the northern species. Our findings highlight that population densities are also moving rapidly towards the poles and the use of presence–absence data can mask the shift of population densities. We encourage use of abundance data in studies considering the effects of climate change on biodiversity.     

Temperature‐driven shifts in a host‐parasite interaction drive nonlinear changes in disease risk

Climate change may shift the timing and consequences of interspecific interactions, including those important to disease spread. Because hosts and pathogens may respond differentially to climate shifts, however, predicting the net effects on disease patterns remains challenging. Here, we used field data to guide a series of laboratory experiments that systematically evaluated the effects of temperature on the full infection process, including survival, penetration, establishment, persistence, and virulence of a highly pathogenic trematode (Ribeiroia ondatrae), and the development and survival of its amphibian host. Our results revealed nonlinearities in pathology as a function of temperature, which likely resulted from changes in both host and parasite processes. Both hosts and parasites responded strongly to temperature; hosts accelerated development while parasites showed enhanced host penetration but reduced establishment (encystment) and survival outside the host. While there were no differences in host survival among treatments, we observed a mid‐temperature peak in parasite‐induced deformities (63% at 20 °C), with the lowest frequency of deformities (12%) occurring at the highest temperature (26 °C). This nonlinear effect could result from temperature‐driven changes in parasite burden owing to shifts in host penetration and/or clearance, reductions in host vulnerability owing to faster development, or both. Furthermore, despite strong temperature‐driven changes in parasite penetration, survival, and establishment, the opposing nature of these effects lead to no difference in tadpole parasite burdens shortly after infection. These findings suggest that temperature‐driven changes to the disease process may not be easily observable from comparison of parasite burdens alone, but multi‐tiered experiments quantifying the responses of hosts, parasites and their interactions can enhance our ability to predict temperature‐driven changes to disease risk. Climate‐driven changes to disease patterns will therefore depend on underlying shifts in host and parasite development rates and the timing of their interactions.     

Accelerating regrowth of temperate‐maritime forests due to environmental change

To understand how environmental changes have influenced forest productivity, stemwood biomass (B) dynamics were analyzed at 1267 permanent inventory plots, covering a combined 209 ha area of unmanaged temperate‐maritime forest in southwest British Columbia, Canada. Net stemwood production (ΔB) was derived from periodic remeasurements of B collected over a 40‐year measurement period (1959–1998) in stands ranging from 20 to 150 years old. Comparison between the integrated age response of net stemwood production, ΔB(A), and the age response of stemwood biomass, B(A), suggested a 58 ± 11% increase in ΔB between the first 40 years of the chronosequence period (1859–1898) and the measurement period. To estimate extrinsic forcing on ΔB, several different candidate models were developed to remove variation explained by intrinsic factors. All models exhibited temporal bias, with positive trends in (observed minus predicted) residual ΔB ranging between of 0.40 and 0.64% yr^?1. Applying the same methods to stemwood growth (G) indicated residual increases ranging from 0.43 and 0.67% yr^?1. Higher trend estimates corresponded with models that included site index (SI) as a predictor, which may reflect exaggeration of the age‐decline in SI tables. Choosing a model that excluded SI, suggested that ΔB increased by 0.40 ± 0.18% yr^?1, while G increased by 0.43 ± 0.12% yr^?1 over the measurement period. Residual G was significantly correlated with atmospheric carbon dioxide (CO₂), temperature (T), and climate moisture index (CMI). However, models driven with climate and CO₂, alone, could not simultaneously explain long‐term and measurement‐period trends without additional representation of indirect effects, perhaps reflecting compound interest on direct physiological responses to environmental change. Evidence of accelerating forest regrowth highlights the value of permanent inventories to detect and understand systematic changes in forest productivity caused by environmental change.     

Response of Sierra Nevada forests to projected climate–wildfire interactions

Climate influences forests directly and indirectly through disturbance. The interaction of climate change and increasing area burned has the potential to alter forest composition and community assembly. However, the overall forest response is likely to be influenced by species‐specific responses to environmental change and the scale of change in overstory species cover. In this study, we sought to quantify how projected changes in climate and large wildfire size would alter forest communities and carbon (C) dynamics, irrespective of competition from nontree species and potential changes in other fire regimes, across the Sierra Nevada, USA. We used a species‐specific, spatially explicit forest landscape model (LANDIS‐II) to evaluate forest response to climate–wildfire interactions under historical (baseline) climate and climate projections from three climate models (GFDL, CCSM3, and CNRM) forced by a medium–high emission scenario (A2) in combination with corresponding climate‐specific large wildfire projections. By late century, we found modest changes in the spatial distribution of dominant species by biomass relative to baseline, but extensive changes in recruitment distribution. Although forest recruitment declined across much of the Sierra, we found that projected climate and wildfire favored the recruitment of more drought‐tolerant species over less drought‐tolerant species relative to baseline, and this change was greatest at mid‐elevations. We also found that projected climate and wildfire decreased tree species richness across a large proportion of the study area and transitioned more area to a C source, which reduced landscape‐level C sequestration potential. Our study, although a conservative estimate, suggests that by late century, forest community distributions may not change as intact units as predicted by biome‐based modeling, but are likely to trend toward simplified community composition as communities gradually disaggregate and the least tolerant species are no longer able to establish. The potential exists for substantial community composition change and forest simplification beyond this century.     

Dispersal through stream networks: modelling climate‐driven range expansions of fishes

Aim To incorporate dispersal through stream networks into models predicting the future distribution of a native, freshwater fish given climate change scenarios. Location Sweden. Methods We used logistic regression to fit climate and habitat data to observed pike (Esox lucius Linnaeus) distributions in 13,476 lakes. We used GIS to map dispersal pathways through streams. Lakes either (1) contained pike or were downstream from pike lakes, (2) were upstream from pike lakes, but downstream from natural dispersal barriers, or (3) were isolated from streams or were upstream from natural dispersal barriers. We then used climate projections to model future distributions of pike and compared our results with and without including dispersal. Results Given climate and habitat, pike were predicted present in all of 99,249 Swedish lakes by 2100. After accounting for dispersal barriers, we only predicted pike presence in 31,538 lakes. Dispersal barriers most strongly limited pike invasion in mountainous regions, but low connectivity also characterized some relatively flat regions. Main conclusions The dendritic network structure of streams and interconnected lakes makes a two‐dimensional representation of the landscape unsuitable for predicting range shifts of many freshwater organisms. If dispersal through stream networks is not accounted for, predictions of future fish distributions in a warmer climate might grossly overestimate range expansions of warm and cool‐water fishes and underestimate range contractions of cold‐water fishes. Dispersal through stream networks can be modelled in any region for which a digital elevation model and species occurrence data are available.     

Climate‐driven shifts in the distribution of koala‐browse species from the Last Interglacial to the near future

The koala's Phascolarctos cinereus distribution is currently restricted to eastern and south‐eastern Australia. However, fossil records dating from 70 ± 4 ka (ka = 10³ yr) from south‐western Australia and the Nullarbor Plain are evidence of subpopulation extinctions in the southwest at least after the Last Interglacial (~128–116 ka). We hypothesize that koala sub‐population extinctions resulted from the eastward retraction of the koala's main browse species in response to unsuitable climatic conditions. We further posit a general reduction in the distribution of main koala‐browse trees in the near future in response climate change. We modelled 60 koala‐browse species and constructed a set of correlative species distribution models for five time periods: Last Interglacial (~128–116 ka), Last Glacial Maximum (~23–19 ka), Mid‐Holocene (~7–5 ka), present (interpolations of observed data, representative of 1960–1990), and 2070. We based our projections on five hindcasts and one forecast of climatic variables extracted from WorldClim based on two general circulation models (considering the most pessimistic scenario of high greenhouse‐gas emissions) and topsoil clay fraction. We used 17 dates of koala fossil specimens identified as reliable from 70 (± 4) to 535 (± 49) ka, with the last appearance of koalas at 70 ka in the southwest. The main simulated koala‐browse species were at their greatest modelled extent of suitability during the Last Glacial Maximum, with the greatest loss of koala habitat occurring between the Mid‐Holocene and the present. We predict a similar habitat loss between the present and 2070. The spatial patterns of habitat change support our hypothesis that koala extinctions in the southwest, Nullarbor Plain and central South Australia resulted from the eastward retraction of the dominant koala‐browse species in response to long‐term climate changes. Future climate patterns will likely increase the extinction risk of koalas in their remaining eastern ranges.     

Predicting when climate‐driven phenotypic change affects population dynamics

Species' responses to climate change are variable and diverse, yet our understanding of how different responses (e.g. physiological, behavioural, demographic) relate and how they affect the parameters most relevant for conservation (e.g. population persistence) is lacking. Despite this, studies that observe changes in one type of response typically assume that effects on population dynamics will occur, perhaps fallaciously. We use a hierarchical framework to explain and test when impacts of climate on traits (e.g. phenology) affect demographic rates (e.g. reproduction) and in turn population dynamics. Using this conceptual framework, we distinguish four mechanisms that can prevent lower‐level responses from impacting population dynamics. Testable hypotheses were identified from the literature that suggest life‐history and ecological characteristics which could predict when these mechanisms are likely to be important. A quantitative example on birds illustrates how, even with limited data and without fully‐parameterized population models, new insights can be gained; differences among species in the impacts of climate‐driven phenological changes on population growth were not explained by the number of broods or density dependence. Our approach helps to predict the types of species in which climate sensitivities of phenotypic traits have strong demographic and population consequences, which is crucial for conservation prioritization of data‐deficient species.     

Incorporating spatial autocorrelation into species distribution models alters forecasts of climate‐mediated range shifts

Species distribution models (SDMs) are widely used to forecast changes in the spatial distributions of species and communities in response to climate change. However, spatial autocorrelation (SA) is rarely accounted for in these models, despite its ubiquity in broad‐scale ecological data. While spatial autocorrelation in model residuals is known to result in biased parameter estimates and the inflation of type I errors, the influence of unmodeled SA on species' range forecasts is poorly understood. Here we quantify how accounting for SA in SDMs influences the magnitude of range shift forecasts produced by SDMs for multiple climate change scenarios. SDMs were fitted to simulated data with a known autocorrelation structure, and to field observations of three mangrove communities from northern Australia displaying strong spatial autocorrelation. Three modeling approaches were implemented: environment‐only models (most frequently applied in species' range forecasts), and two approaches that incorporate SA; autologistic models and residuals autocovariate (RAC) models. Differences in forecasts among modeling approaches and climate scenarios were quantified. While all model predictions at the current time closely matched that of the actual current distribution of the mangrove communities, under the climate change scenarios environment‐only models forecast substantially greater range shifts than models incorporating SA. Furthermore, the magnitude of these differences intensified with increasing increments of climate change across the scenarios. When models do not account for SA, forecasts of species' range shifts indicate more extreme impacts of climate change, compared to models that explicitly account for SA. Therefore, where biological or population processes induce substantial autocorrelation in the distribution of organisms, and this is not modeled, model predictions will be inaccurate. These results have global importance for conservation efforts as inaccurate forecasts lead to ineffective prioritization of conservation activities and potentially to avoidable species extinctions.     

Fire severity mediates climate‐driven shifts in understorey community composition of black spruce stands of interior Alaska

Question: How do pre‐fire conditions (community composition and environmental characteristics) and climate‐driven disturbance characteristics (fire severity) affect post‐fire community composition in black spruce stands? Location: Northern boreal forest, interior Alaska. Methods: We compared plant community composition and environmental stand characteristics in 14 black spruce stands before and after multiple, naturally occurring wildfires. We used a combination of vegetation table sorting, univariate (ANOVA, paired t‐tests), and multivariate (detrended correspondence analysis) statistics to determine the impact of fire severity and site moisture on community composition, dominant species and growth forms. Results: Severe wildfires caused a 50% reduction in number of plant species in our study sites. The largest species loss, and therefore the greatest change in species composition, occurred in severely burned sites. This was due mostly to loss of non‐vascular species (mosses and lichens) and evergreen shrubs. New species recruited most abundantly to severely burned sites, contributing to high species turnover on these sites. As well as the strong effect of fire severity, pre‐fire and post‐fire mineral soil pH had an effect on post‐fire vegetation patterns, suggesting a legacy effect of site acidity. In contrast, pre‐fire site moisture, which was a strong determinant of pre‐fire community composition, showed no relationship with post‐fire community composition. Site moisture was altered by fire, due to changes in permafrost, and therefore post‐fire site moisture overrode pre‐fire site moisture as a strong correlate. Conclusions: In the rapidly warming climate of interior Alaska, changes in fire severity had more effect on post‐fire community composition than did environmental factors (moisture and pH) that govern landscape patterns of unburned vegetation. This suggests that climate change effects on future community composition of black spruce forests may be mediated more strongly by fire severity than by current landscape patterns. Hence, models that represent the effects of climate change on boreal forests could improve their accuracy by including dynamic responses to fire disturbance.