A unified explanation for the morphology of raised peatlands.

Raised peatlands, or bogs, are gently mounded landforms that are composed entirely of organic matter1-4 and store the most carbon per area of any terrestrial ecosystem5. The shapes of bogs are critically important because their domed morphology4,6,7 accounts for much of the carbon that bogs store and determines how they will respond to interventions8,9 to stop greenhouse gas emissions and fires after anthropogenic drainage10-13. However, a general theory to infer the morphology of bogs is still lacking4,6,7. Here we show that an equation based on the processes universal to bogs explains their morphology across biomes, from Alaska, through the tropics, to New Zealand. In contrast to earlier models of bog morphology that attempted to describe only long-term equilibrium shapes4,6,7 and were, therefore, inapplicable to most bogs14-16, our approach makes no such assumption and makes it possible to infer full shapes of bogs from a sample of elevations, such as a single elevation transect. Our findings provide a foundation for quantitative inference about the morphology, hydrology and carbon storage of bogs through Earth's history, as well as a basis for planning natural climate solutions by rewetting damaged bogs around the world.

Climate velocities and species tracking in global mountain regions.

Mountain ranges contain high concentrations of endemic species and are indispensable refugia for lowland species that are facing anthropogenic climate change1,2. Forecasting biodiversity redistribution hinges on assessing whether species can track shifting isotherms as the climate warms3,4. However, a global analysis of the velocities of isotherm shifts along elevation gradients is hindered by the scarcity of weather stations in mountainous regions5. Here we address this issue by mapping the lapse rate of temperature (LRT) across mountain regions globally, both by using satellite data (SLRT) and by using the laws of thermodynamics to account for water vapour6 (that is, the moist adiabatic lapse rate (MALRT)). By dividing the rate of surface warming from 1971 to 2020 by either the SLRT or the MALRT, we provide maps of vertical isotherm shift velocities. We identify 17 mountain regions with exceptionally high vertical isotherm shift velocities (greater than 11.67 m per year for the SLRT; greater than 8.25 m per year for the MALRT), predominantly in dry areas but also in wet regions with shallow lapse rates; for example, northern Sumatra, the Brazilian highlands and southern Africa. By linking these velocities to the velocities of species range shifts, we report instances of close tracking in mountains with lower climate velocities. However, many species lag behind, suggesting that range shift dynamics would persist even if we managed to curb climate-change trajectories. Our findings are key for devising global conservation strategies, particularly in the 17 high-velocity mountain regions that we have identified.

Disappearing cities on US coasts.

The sea level along the US coastlines is projected to rise by 0.25-0.3 m by 2050, increasing the probability of more destructive flooding and inundation in major cities1-3. However, these impacts may be exacerbated by coastal subsidence-the sinking of coastal land areas4-a factor that is often underrepresented in coastal-management policies and long-term urban planning2,5. In this study, we combine high-resolution vertical land motion (that is, raising or lowering of land) and elevation datasets with projections of sea-level rise to quantify the potential inundated areas in 32 major US coastal cities. Here we show that, even when considering the current coastal-defence structures, further land area of between 1,006 and 1,389 km2 is threatened by relative sea-level rise by 2050, posing a threat to a population of 55,000-273,000 people and 31,000-171,000 properties. Our analysis shows that not accounting for spatially variable land subsidence within the cities may lead to inaccurate projections of expected exposure. These potential consequences show the scale of the adaptation challenge, which is not appreciated in most US coastal cities.

Ancient Plasmodium genomes shed light on the history of human malaria.

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

Precipitation regime changes in High Mountain Asia driven by cleaner air.

High Mountain Asia (HMA) has experienced a spatial imbalance in water resources in recent decades, partly because of a dipolar pattern of precipitation changes known as South Drying-North Wetting1. These changes can be influenced by both human activities and internal climate variability2,3. Although climate projections indicate a future widespread wetting trend over HMA1,4, the timing and mechanism of the transition from a dipolar to a monopolar pattern remain unknown. Here we demonstrate that the observed dipolar precipitation change in HMA during summer is primarily driven by westerly- and monsoon-associated precipitation patterns. The weakening of the Asian westerly jet, caused by the uneven emission of anthropogenic aerosols, favoured a dipolar precipitation trend from 1951 to 2020. Moreover, the phase transition of the Interdecadal Pacific Oscillation induces an out-of-phase precipitation change between the core region of the South Asian monsoon and southeastern HMA. Under medium- or high-emission scenarios, corresponding to a global warming of 0.6-1.1 °C compared with the present, the dipolar pattern is projected to shift to a monopolar wetting trend in the 2040s. This shift in precipitation patterns is mainly attributed to the intensified jet stream resulting from reduced emissions of anthropogenic aerosols. These findings underscore the importance of considering the impact of aerosol emission reduction in future social planning by policymakers.

Oceanic climate changes threaten the sustainability of Asia's water tower.

Water resources sustainability in High Mountain Asia (HMA) surrounding the Tibetan Plateau (TP)-known as Asia's water tower-has triggered widespread concerns because HMA protects millions of people against water stress1,2. However, the mechanisms behind the heterogeneous trends observed in terrestrial water storage (TWS) over the TP remain poorly understood. Here we use a Lagrangian particle dispersion model and satellite observations to attribute about 1 Gt of monthly TWS decline in the southern TP during 2003-2016 to westerlies-carried deficit in precipitation minus evaporation (PME) from the southeast North Atlantic. We further show that HMA blocks the propagation of PME deficit into the central TP, causing a monthly TWS increase by about 0.5 Gt. Furthermore, warming-induced snow and glacial melt as well as drying-induced TWS depletion in HMA weaken the blocking of HMA's mountains, causing persistent northward expansion of the TP's TWS deficit since 2009. Future projections under two emissions scenarios verified by satellite observations during 2020-2021 indicate that, by the end of the twenty-first century, up to 84% (for scenario SSP245) and 97% (for scenario SSP585) of the TP could be afflicted by TWS deficits. Our findings indicate a trajectory towards unsustainable water systems in HMA that could exacerbate downstream water stress.

Less extreme and earlier outbursts of ice-dammed lakes since 1900.

Episodic failures of ice-dammed lakes have produced some of the largest floods in history, with disastrous consequences for communities in high mountains1-7. Yet, estimating changes in the activity of ice-dam failures through time remains controversial because of inconsistent regional flood databases. Here, by collating 1,569 ice-dam failures in six major mountain regions, we systematically assess trends in peak discharge, volume, annual timing and source elevation between 1900 and 2021. We show that extreme peak flows and volumes (10 per cent highest) have declined by about an order of magnitude over this period in five of the six regions, whereas median flood discharges have fallen less or have remained unchanged. Ice-dam floods worldwide today originate at higher elevations and happen about six weeks earlier in the year than in 1900. Individual ice-dammed lakes with repeated outbursts show similar negative trends in magnitude and earlier occurrence, although with only moderate correlation to glacier thinning8. We anticipate that ice dams will continue to fail in the near future, even as glaciers thin and recede. Yet widespread deglaciation, projected for nearly all regions by the end of the twenty-first century9, may bring most outburst activity to a halt.

Trait-based sensitivity of large mammals to a catastrophic tropical cyclone.

Extreme weather events perturb ecosystems and increasingly threaten biodiversity1. Ecologists emphasize the need to forecast and mitigate the impacts of these events, which requires knowledge of how risk is distributed among species and environments. However, the scale and unpredictability of extreme events complicate risk assessment1-4-especially for large animals (megafauna), which are ecologically important and disproportionately threatened but are wide-ranging and difficult to monitor5. Traits such as body size, dispersal ability and habitat affiliation are hypothesized to determine the vulnerability of animals to natural hazards1,6,7. Yet it has rarely been possible to test these hypotheses or, more generally, to link the short-term and long-term ecological effects of weather-related disturbance8,9. Here we show how large herbivores and carnivores in Mozambique responded to Intense Tropical Cyclone Idai, the deadliest storm on record in Africa, across scales ranging from individual decisions in the hours after landfall to changes in community composition nearly 2 years later. Animals responded behaviourally to rising floodwaters by moving upslope and shifting their diets. Body size and habitat association independently predicted population-level impacts: five of the smallest and most lowland-affiliated herbivore species declined by an average of 28% in the 20 months after landfall, while four of the largest and most upland-affiliated species increased by an average of 26%. We attribute the sensitivity of small-bodied species to their limited mobility and physiological constraints, which restricted their ability to avoid the flood and endure subsequent reductions in the quantity and quality of food. Our results identify general traits that govern animal responses to severe weather, which may help to inform wildlife conservation in a volatile climate.

Hypoxia Inducible Factor pathway proteins in high-altitude mammals.

Humans and other mammals inhabit hypoxic high-altitude locales. In many of these species, genes under positive selection include ones in the Hypoxia Inducible Factor (HIF) pathway. One is PHD2 (EGLN1), which encodes for a key oxygen sensor. Another is HIF2A (EPAS1), which encodes for a PHD2-regulated transcription factor. Recent studies have provided insights into mechanisms for these high-altitude alleles. These studies have (i) shown that selection can occur on nonconserved, unstructured regions of proteins, (ii) revealed that high altitude-associated amino acid substitutions can have differential effects on protein-protein interactions, (iii) provided evidence for convergent evolution by different molecular mechanisms, and (iv) suggested that mutations in different genes can complement one another to produce a set of adaptive phenotypes.

Topography of mountain belts controlled by rheology and surface processes.

It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process5,6 and mantle-scale tectonic model7. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state8, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.

An oxygen-sensing mechanism for angiosperm adaptation to altitude.

Flowering plants (angiosperms) can grow at extreme altitudes, and have been observed growing as high as 6,400 metres above sea level1,2; however, the molecular mechanisms that enable plant adaptation specifically to altitude are unknown. One distinguishing feature of increasing altitude is a reduction in the partial pressure of oxygen (pO2). Here we investigated the relationship between altitude and oxygen sensing in relation to chlorophyll biosynthesis-which requires molecular oxygen3-and hypoxia-related gene expression. We show that in etiolated seedlings of angiosperm species, steady-state levels of the phototoxic chlorophyll precursor protochlorophyllide are influenced by sensing of atmospheric oxygen concentration. In Arabidopsis thaliana, this is mediated by the PLANT CYSTEINE OXIDASE (PCO) N-degron pathway substrates GROUP VII ETHYLENE RESPONSE FACTOR transcription factors (ERFVIIs). ERFVIIs positively regulate expression of FLUORESCENT IN BLUE LIGHT (FLU), which represses the first committed step of chlorophyll biosynthesis, forming an inactivation complex with tetrapyrrole synthesis enzymes that are negatively regulated by ERFVIIs, thereby suppressing protochlorophyllide. In natural populations representing diverse angiosperm clades, we find oxygen-dependent altitudinal clines for steady-state levels of protochlorophyllide, expression of inactivation complex components and hypoxia-related genes. Finally, A. thaliana accessions from contrasting altitudes display altitude-dependent ERFVII activity and accumulation. We thus identify a mechanism for genetic adaptation to absolute altitude through alteration of the sensitivity of the oxygen-sensing system.

Global fine-scale changes in ambient NO2 during COVID-19 lockdowns.

Nitrogen dioxide (NO2) is an important contributor to air pollution and can adversely affect human health1-9. A decrease in NO2 concentrations has been reported as a result of lockdown measures to reduce the spread of COVID-1910-20. Questions remain, however, regarding the relationship of satellite-derived atmospheric column NO2 data with health-relevant ambient ground-level concentrations, and the representativeness of limited ground-based monitoring data for global assessment. Here we derive spatially resolved, global ground-level NO2 concentrations from NO2 column densities observed by the TROPOMI satellite instrument at sufficiently fine resolution (approximately one kilometre) to allow assessment of individual cities during COVID-19 lockdowns in 2020 compared to 2019. We apply these estimates to quantify NO2 changes in more than 200 cities, including 65 cities without available ground monitoring, largely in lower-income regions. Mean country-level population-weighted NO2 concentrations are 29% ± 3% lower in countries with strict lockdown conditions than in those without. Relative to long-term trends, NO2 decreases during COVID-19 lockdowns exceed recent Ozone Monitoring Instrument (OMI)-derived year-to-year decreases from emission controls, comparable to 15 ± 4 years of reductions globally. Our case studies indicate that the sensitivity of NO2 to lockdowns varies by country and emissions sector, demonstrating the critical need for spatially resolved observational information provided by these satellite-derived surface concentration estimates.

Mechanical forcing of the North American monsoon by orography.

A band of intense rainfall extends more than 1,000 km along Mexico's west coast during Northern Hemisphere summer, constituting the core of the North American monsoon1,2. As in other tropical monsoons, this rainfall maximum is commonly thought to be thermally forced by emission of heat from land and elevated terrain into the overlying atmosphere3-5, but a clear understanding of the fundamental mechanism governing this monsoon is lacking. Here we show that the core North American monsoon is generated when Mexico's Sierra Madre mountains deflect the extratropical jet stream towards the Equator, mechanically forcing eastward, upslope flow that lifts warm and moist air to produce convective rainfall. These findings are based on analyses of dynamic and thermodynamic structures in observations, global climate model integrations and adiabatic stationary wave solutions. Land surface heat fluxes do precondition the atmosphere for convection, particularly in summer afternoons, but these heat fluxes alone are insufficient for producing the observed rainfall maximum. Our results indicate that the core North American monsoon should be understood as convectively enhanced orographic rainfall in a mechanically forced stationary wave, not as a classic, thermally forced tropical monsoon. This has implications for the response of the North American monsoon to past and future global climate change, making trends in jet stream interactions with orography of central importance.

Widespread six degrees Celsius cooling on land during the Last Glacial Maximum.

The magnitude of global cooling during the Last Glacial Maximum (LGM, the coldest multimillennial interval of the last glacial period) is an important constraint for evaluating estimates of Earth's climate sensitivity1,2. Reliable LGM temperatures come from high-latitude ice cores3,4, but substantial disagreement exists between proxy records in the low latitudes1,5-8, where quantitative low-elevation records on land are scarce. Filling this data gap, noble gases in ancient groundwater record past land surface temperatures through a direct physical relationship that is rooted in their temperature-dependent solubility in water9,10. Dissolved noble gases are suitable tracers of LGM temperature because of their complete insensitivity to biological and chemical processes and the ubiquity of LGM-aged groundwater around the globe11,12. However, although several individual noble gas studies have found substantial tropical LGM cooling13-16, they have used different methodologies and provide limited spatial coverage. Here we use noble gases in groundwater to show that the low-altitude, low-to-mid-latitude land surface (45 degrees south to 35 degrees north) cooled by 5.8 ± 0.6 degrees Celsius (mean ± 95% confidence interval) during the LGM. Our analysis includes four decades of groundwater noble gas data from six continents, along with new records from the tropics, all of which were interpreted using the same physical framework. Our land-based result broadly supports a recent reconstruction based on marine proxy data assimilation1 that suggested greater climate sensitivity than previous estimates5-7.

Extreme elevational migration spurred cryptic speciation in giant hummingbirds.

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Massive seasonal high-altitude migrations of nocturnal insects above the agricultural plains of East China.

Long-distance migrations of insects contribute to ecosystem functioning but also have important economic impacts when the migrants are pests or provide ecosystem services. We combined radar monitoring, aerial sampling, and searchlight trapping, to quantify the annual pattern of nocturnal insect migration above the densely populated agricultural lands of East China. A total of ~9.3 trillion nocturnal insect migrants (15,000 t of biomass), predominantly Lepidoptera, Hemiptera, and Diptera, including many crop pests and disease vectors, fly at heights up to 1 km above this 600 km-wide region every year. Larger migrants (>10 mg) exhibited seasonal reversal of movement directions, comprising northward expansion during spring and summer, followed by southward movements during fall. This north-south transfer was not balanced, however, with southward movement in fall 0.66× that of northward movement in spring and summer. Spring and summer migrations were strongest when the wind had a northward component, while in fall, stronger movements occurred on winds that allowed movement with a southward component; heading directions of larger insects were generally close to the track direction. These findings indicate adaptations leading to movement in seasonally favorable directions. We compare our results from China with similar studies in Europe and North America and conclude that ecological patterns and behavioral adaptations are similar across the Northern Hemisphere. The predominance of pests among these nocturnal migrants has severe implications for food security and grower prosperity throughout this heavily populated region, and knowledge of their migrations is potentially valuable for forecasting pest impacts and planning timely management actions.

Genetic polymorphisms and expression of Rhesus blood group RHCE are associated with 2,3-bisphosphoglycerate in humans at high altitude.

Red blood cell (RBC) metabolic reprogramming upon exposure to high altitude contributes to physiological human adaptations to hypoxia, a multifaceted process critical to health and disease. To delve into the molecular underpinnings of this phenomenon, first, we performed a multi-omics analysis of RBCs from six lowlanders after exposure to high-altitude hypoxia, with longitudinal sampling at baseline, upon ascent to 5,100 m and descent to sea level. Results highlighted an association between erythrocyte levels of 2,3-bisphosphoglycerate (BPG), an allosteric regulator of hemoglobin that favors oxygen off-loading in the face of hypoxia, and expression levels of the Rhesus blood group RHCE protein. We then expanded on these findings by measuring BPG in RBCs from 13,091 blood donors from the Recipient Epidemiology and Donor Evaluation Study. These data informed a genome-wide association study using BPG levels as a quantitative trait, which identified genetic polymorphisms in the region coding for the Rhesus blood group RHCE as critical determinants of BPG levels in erythrocytes from healthy human volunteers. Mechanistically, we suggest that the Rh group complex, which participates in the exchange of ammonium with the extracellular compartment, may contribute to intracellular alkalinization, thus favoring BPG mutase activity.

Importance and vulnerability of the world's water towers.

Mountains are the water towers of the world, supplying a substantial part of both natural and anthropogenic water demands1,2. They are highly sensitive and prone to climate change3,4, yet their importance and vulnerability have not been quantified at the global scale. Here we present a global water tower index (WTI), which ranks all water towers in terms of their water-supplying role and the downstream dependence of ecosystems and society. For each water tower, we assess its vulnerability related to water stress, governance, hydropolitical tension and future climatic and socio-economic changes. We conclude that the most important (highest WTI) water towers are also among the most vulnerable, and that climatic and socio-economic changes will affect them profoundly. This could negatively impact 1.9 billion people living in (0.3 billion) or directly downstream of (1.6 billion) mountainous areas. Immediate action is required to safeguard the future of the world's most important and vulnerable water towers.

Evidence of human occupation in Mexico around the Last Glacial Maximum.

The initial colonization of the Americas remains a highly debated topic1, and the exact timing of the first arrivals is unknown. The earliest archaeological record of Mexico-which holds a key geographical position in the Americas-is poorly known and understudied. Historically, the region has remained on the periphery of research focused on the first American populations2. However, recent investigations provide reliable evidence of a human presence in the northwest region of Mexico3,4, the Chiapas Highlands5, Central Mexico6 and the Caribbean coast7-9 during the Late Pleistocene and Early Holocene epochs. Here we present results of recent excavations at Chiquihuite Cave-a high-altitude site in central-northern Mexico-that corroborate previous findings in the Americas10-17of cultural evidence that dates to the Last Glacial Maximum (26,500-19,000 years ago)18, and which push back dates for human dispersal to the region possibly as early as 33,000-31,000 years ago. The site yielded about 1,900 stone artefacts within a 3-m-deep stratified sequence, revealing a previously unknown lithic industry that underwent only minor changes over millennia. More than 50 radiocarbon and luminescence dates provide chronological control, and genetic, palaeoenvironmental and chemical data document the changing environments in which the occupants lived. Our results provide new evidence for the antiquity of humans in the Americas, illustrate the cultural diversity of the earliest dispersal groups (which predate those of the Clovis culture) and open new directions of research.