Natural montane peatlands generally are not a significant source of nitrous oxide (N₂O) due to environment limitation including phosphorus (P) scarcity and temperature lowness. Phosphorus enrichment and warming caused by global change are altering this limitation, likely to increase the source function of N₂O. However, the combined effects of P addition and warming on N₂O fluxes and biotic/abiotic factor in peatlands are still uncertain. To address this, we investigated the long-term (12 yrs) effects of P addition (5 and 10 kg ha^-1 yr^-1) and its interaction with warming on N₂O fluxes in a peatland. The results showed that although long-term P addition did not significantly affect the source/sink function of N₂O in the peatland, it stimulated enzyme activities, promoted peat decomposition. However, warming amplified the effect of P addition to increase N₂O emissions by stimulating enzyme activities and changing soil stoichiometry, even turned the peatland into a significant source of N₂O with an emission of approximate 100 g m^-2 during the growing season. Our study suggests that P enrichment under the current background of global warming will enhance the possibility of strong N₂O emissions in montane peatlands, which may increase the risk to further aggravate global warming.

The spatiotemporal variations of rainfall erosivity in the Hengduan Mountains, charac-terized by rugged terrain and high potential soil erosion risks, over the past 30 years was exam-ined. The changing trends of rainfall erosivity for 2025-2040 was also be investigated under the comprehensive scenario of moderate socio-economic development (SSP2-4.5) combined with me-dium-low radiative forcing, using four global climate models (GCMs) based on CMIP6. The results indicated that: (1) The annual distribution of rainfall erosivity in the Hengduan Mountains exhib-ited significant seasonal variations, with the order of erosivity being summer > autumn > spring > winter on a seasonal scale. (2) Over the past 30 years, there has been a slight decrease in annual precipitation and a slight increase in rainfall erosivity, with periodic extreme values occurring every 6-8 years. (3) Rainfall erosivity showed a decreasing gradient from southeast to northwest in terms of spatial distribution. There was a significant positive correlation between rainfall ero-sivity and precipitation, while a significant negative correlation existed with elevation in the ver-tical direction. Moreover, there was an increasing trend of rainfall erosivity in the northeastern part of the Hengduan Mountains and a decreasing trend in the southern region. (4) Under the joint driving forces of increased precipitation and erosive rainfall events, rainfall erosivity in the future is expected to significantly increase, posing a more severe risk of soil erosion in the Heng-duan Mountains.

A new subspecies of Sideroxylon mascatense (A.DC.) T.D.Penn. in the Western Hajar Mountains, differing from the type subspecies in the colour of the mature fruit, is described. The two colour forms are genetically isolated, showing that these two are distinct taxa.

This study was conducted to investigate the status of regional linkage in tourism development in in the Midlands and Northern Mountains of Vietnam. The data was collected from a survey of 755 people, including officials from State management bodies in charge of tourism, officials and staffs at tourism resorts, tourism firms, tourism scientists and tourists. In addition, we conducted 10 group discussions, interviewed 30 State tourism agency officials and tourism firms in the Midland and Mountainous provinces of Vietnam. The results show that tourism development in Vietnam in general and the Northwest region in particular is extremely fragmented, not yet forming a regional linkage; Regional and national tourism development programs are just formalistic. The main cause of the situation is the limited regional integration policy in Vietnam, the lack of appropriate regional governance mechanisms and inactive participation of the private sector in regional integration. Based on the findings, we propose a tourism sector linkage model; besides, policy implications are given for fulfilling the linkage policy in the Midlands and Northern Mountains area.

Improved understanding of climate-growth relationships of multi-species is fundamental to understand and predict response of forest growth to future climate change. Forests are mainly composed of conifers in Northwestern Yunnan Plateau, but variations of growth response to climates among the species are not well understood. To detect growth response of multiple species to climate change, we developed residual chronologies of four major conifers, i.e. Abies georgei, Picea likiangensis, Pinus densata and Larix potaninii at upper distributional limits in Shika Snow Mountain. By using dendroclimatology method, we analyzed correlations between the residual chronologies and climate variables. The results showed that conifer radial growth was influenced by both temperature and precipitation in Shika Snow Mountain. Previous November temperature, previous July mean maximum temperature (T_max) and current June precipitation were the common climatic factors, which had consistent influences on radial growth of four species. Temperature in previous post growing season (September–October) and current growing season (June-August), and precipitation in previous August were the common climatic factors, which had divergent impacts on four species radial growth. Current May T_max and early growing season (April-May) precipitation showed positive and negative influences on growth of P. likiangensis, respectively. Temperature in current post growing season positively affected growth of A. georgei. According to the prediction of climate models and our understanding in growth response of four species to climate variables, we may understand growth response to climate change at species level. It is difficult to predict future forest growth in the study area, since future climate change might cause both increases or decreases for four species and indirect effects of climate change on forest should be considered.

A large variety of remote sensing-based land use/land cover (LULC) products are currently available on national and global scales. This literature review and in-situ validation study evaluates the suitability of these products for local scale applications in the complex terrain of the Ethiopian mountains. For the review, 146 research papers were analyzed. Most studies (73%) have been published since 2013 and are based on individually computed maps. Not a single study relied on readily available LULC products. Nine readily available LULC products with 20, 30, 300, 500 and 1,000 m spatial resolution have been identified at national and global scales. To complement and extend this body of research, the recent (since 2013) LULC products were validated using 185 ground truth points collected in the Bale Mountains National Park between 1,500 and 4,385 m a.s.l. The results indicate a rather poor overall accuracy (<50%).

This paper assesses the use the Global Sustainable Tourism Criteria (GSTC) indicators for monitoring tourism along a mountainous route. The study was carried out along a mountain route in the Drakensberg Mountains of South Africa. A survey was carried out among tourism businesses, local communities and local government officials located along the Maluti Route, which passes though the eastern part of the Free State Province of South Africa. This study constitutes the first assessment in which the GSTC indicators are applied in the assessment of the sustainability of route tourism in South Africa. The study first evaluated the importance of each indicator by considering the level of application in tourism business organisations before examining the perceptions of local community members and local government officials regarding the applicability of the indicators in the assessment of tourism sustainability. The results suggest that there is a general appreciation of tourism sustainability among tourism business operators although shortcomings were evident on environmental and socio-economic indicators. In conclusion, the paper suggests ways through which tourism sustainability could be enhanced among different players along the Maluti Route.

Wild Chinese prickly ash with elevated antioxidants is a valuable genetic resource for Zanthoxylum bungeanum Maxim improvement. There are rich wild germplasm resources in the Qinling Mountains. In a study with wild germplasm resources from different altitudes and six cultivated varieties, the phenolic and flavonoid compounds were analyzed by high performance liquid chromatography (HPLC). The chromatograms of them were basically the same, although their chemical composition content was greatly different. The thirty samples were divided into three categories through the hierarchical clustering analysis. And catechin, hyperoside and quercitrin were considered to be key compound for the quality evaluation, by contrast, the wild samples with an altitude of 2300±50 m (Ⅳ group) had the highest content of key compounds, and showed stronger antioxidant activity and antibacterial ability, indicating that these wild samples could be used as an excellent breeding resource. This is the first time to evaluate the quality of wild Chinese prickly ash in different altitude areas of Qinling Mountains. These excellent wild germplasm resources provided substantial potential accessions for use directly in Chinese prickly ash breeding programs.

Background: The Hengduan Mountains are one of the most species–rich mountainous areas in the world. The origin and evolution of such a remarkable biodiversity are likely to be associated with geological or climatic dynamics, as well as taxon-specific biotic processes (e.g., hybridization, polyploidization, etc.). Here, we investigate the mechanisms fostering the diversification of the endemic Himalopsyche martynovi complex, a poorly known group of aquatic insects. Methods: We used multiple allelic datasets generated from 691 AHE loci to reconstruct species and RaxML phylogenetic trees. We selected the most reliable phylogenetic tree to perform network and gene flow analyses. Results: Phylogenetic reconstructions and network analysis identified three clades, including H. epikur, H. martynovi sensu stricto and H. cf. martynovi. Himalopsyche martynovi sensu stricto and H. cf. martynovi present an intermediate morphology between H. epikur and H. viteceki, the closest known relative to the H. martynovi–complex. The gene flow analysis revealed extensive gene flow among these lineages. Conclusion: Our results suggest that H. viteceki and H. epikur are likely to have contributed to the evolution of H. martynovi sensu stricto and H. cf. martynovi via gene flow, and thus, our study provides insights in the diversification process of a lesser–known ecological group, and hints at the potential role of gene flow in the emergence of biological novelty in the Hengduan Mountains.

Conservation scientists recognize that additional protected areas are needed to maintain biological diversity and ecological processes. As regional conservation planners embark on recommending additional areas for protection in formal conservation reserves, it is important to evaluate candidate lands for their role in building a resilient protected areas system of the future. Here, we evaluate North Carolina’s Mountain Treasures with respect to their (i) ecological integrity, (ii) role in connecting existing core protected areas, (iii) potential to diversify the ecosystem representation of reserves, and (iv) role in maintaining hotspots of biologically-rich areas not well protected. Mountain Treasures represent a citizen inventory of roadless areas and serve as candidates for elevated levels of conservation protection on U.S. federal lands. We compared Mountain Treasures to other candidate lands throughout the country to evaluate their potential national significance. While the Mountain Treasures tended to be more impacted by human modifications than other roadless areas, they are as important as other roadless areas with respect to their role in connecting existing protected areas and diversifying representation of ecosystems in conservation reserves. However, Mountain Treasures tended to have a much higher biodiversity priority index than other roadless areas leading to an overall higher composite score compared to other roadless areas. Our analysis serves as an example of how using broad-scale datasets can help conservation planners assess the national significance of local areas.

South rear-edge populations of widely distributed temperate and boreal tree species such as birches (Betula pubescens and B. pendula) are considered particularly vulnerable to climate warming, and at the same time they constitute genetic reservoirs of drought-adapted ecotypes. Here, we compared radial growth patterns and responses to climate, river or dam flows and a drought index of rear-edge (southernmost) populations (Toledo Mountains, central-southern Spain) with populations located in northern Spain of B. pubescens and B. pendula. Then, we per-formed a comparative analysis across Europe of B. pendula populations. The main climatic con-straint of birch growth was high summer water deficit, although the effect of local hydrological conditions was particularly important in rear-edge populations. We found declining growth trends in rear-edge stands since the early 21st century, related to decreasing water availability and increasing aridity. Our results also suggested distinct growth patterns and climate-growth associations of B. pendula across Europe that show how populations further south and in warmer locations were more sensitive to the drought. In conclusion, drought-induced growth decline can be exacerbated by local human land uses leading to reduced river inflow, thus endangering birch populations at their southern distribution limit. Protection of endangered rear-edge birch populations requires adequate management of local water resources.

Snow is a vital environmental parameter and dynamically responsive to climate change, particularly in mountainous regions. Snow cover can be monitored at variable spatial scales using Earth Observation (EO) data. Long-lasting remote sensing missions enable the generation of multi-decadal time series and thus the detection of long-term trends. However, there have been few attempts to use these to model future snow cover dynamics. In this study, we therefore explore the potential of such time series to forecast the Snow Line Elevation (SLE) in the European Alps. We generate monthly SLE time series from the entire Landsat archive (1985-2021) in 43 Alpine catchments. Positive long-term SLE change rates are detected, with the highest rates (5-8 m/y) in the Western and Central Alps. We utilize this SLE dataset to implement and evaluate seven uni-variate time series modeling and forecasting approaches. The best results were achieved by Random Forests, with a Nash-Sutcliffe efficiency (NSE) of 0.79 and a Mean Absolut Error (MAE) of 258 m, Telescope (0.76, 268 m), and seasonal ARIMA (0.75, 270 m). Since the model performance varies strongly with the input data, we developed a Combined forecast based on the best performing methods in each catchment. This approach was then used to forecast the SLE for the years 2022-2029. In the majority of the catchments the shift of the forecast median SLE level retained the sign of the long-term trend. In cases where a deviating SLE dynamic is forecast a discussion based on the unique properties of the catchment and past SLE dynamics is required. In the future, we expect major improvements in our SLE forecasting efforts by including external predictor variables in a multi-variate modeling approach.

Temperature-index modeling is used to determine the magnitude of temperature depression on the Blanca Massif, Colorado, required to maintain steady-state mass balances of nine reconstructed glaciers at their extent during the Last Glacial Maximum (LGM). The mean temperature depression thus determined is ~8.6 +0.7/–0.9 °C where the uncertainties account for those inherent in the glacier reconstructions, in model parameters (e.g. melt factors), and possible modest changes in LGM precipitation. Associated equilibrium-line altitudes (ELAs) exhibit a statistically significant directional dependency being lower toward the north and east. Under the assumption that regional temperature change was uniform, required changes in precipitation vary systematically – also exhibiting a directional dependency coinciding with that in ELAs – and indicate increases (over modern) occurred on the eastern side of the massif while decreases occurred on the western side. This disparity represents a strengthening of a precipitation asymmetry, particularly winter precipitation, that exists today. The modern precipitation asymmetry may be a consequence of snow being blown over to the eastern side of the massif (advective transport) by southwesterly flow. Intensification of this flow during the LGM would have enhanced advection, and augmented snow accumulation on glaciers, thus explaining the lower ELAs and increased precipitation on that side of the massif.

After World War II, the economic recovery of Western Europe implied a swift economic transition for all regions, including the area of the Alps, although affecting various parts at different pace and stages. The resulting out-migration led to population decline in some mountain valleys and regions already since the 1950s. A similar out-migration movement began in China after its rural reform started in the 1970s. The effect was in some cases even more significant than in the Alps, with the first village being deserted in the 1980s. Current estimations report of about 380,000 abandoned rural villages in China between 2000 and 2016, particularly in its mountain regions. While lower population densities might alleviate the pressures on ecology and contribute to environmental benefits, these movements aggravate a spiraling-down process of local economies and culture. In the Alps many regions facing challenges of out-migration and economic changes elaborated agritourism schemes that provided both economic incentives and stability to involved mountain farmers, and continuation of local land management systems. In contrast, in China hardly any comparable trends of rural tourism developed. However, in recent years China's interest for tourism-oriented farm diversification increased and a range of rural tourism and agricultural tourism initiatives emerged. This paper focuses on the analysis of successful initiatives, problems and development prospects in the Alps and China's rural areas, redefining agritourism as a systematic integrated activity. Agritourism might therefore be assessed as a core element of the future sustainable development of the Alps and the Chinese countryside.

In the future, the pursuit of high-quality economic development and a focus on ecological environmental protection in China will inevitably result in significant conflicts between land use and ecological land use. The challenge lies in achieving sustainable high-quality development while simultaneously protecting the ecological environment, optimizing the land use structure, and promoting a harmonious relationship between humans and the land. These challenges are faced by all regions. Land use conflicts primarily occur in peri-urban areas characterized by prominent economic development and urban agglomeration. Previous studies have mainly focused on analyzing the effects of land use on habitat quality during intense urbanization. However, it is important to recognize that land pressure encompasses economic, ecological, and social aspects. To gain a comprehensive understanding of the spatial conflict of land use and the impact on habitat quality in Ankang, a city that has been advocating ecological protection for the past two decades, this study aims to objectively analyze the spatial trends in land use changes in such cities. Additionally, it aims to provide insights for the harmonious development of land use in eco-region-oriented cities. Using the SSP-RCP scenarios provided by CMIP6, this paper employs a system analysis method, PLUS model, InVEST model, and land use conflict measurement model to dynamically simulate the future habitat quality and spatial conflict patterns of land use in Ankang City. The study explores the spatial coupling effect of both factors under different scenarios. The results indicate the following:(1) Under different future shared socio-economic path scenarios, land use intensity and land conflict levels follow the order of SSP585 (high forcing scenario), SSP370 (medium to high forcing scenario), SSP245 (medium forcing scenario), and SSP126 (low forcing scenario), with intensity and conflict decreasing accordingly. (2) The overall development trend in Ankang City reveals an intensification of land use conflicts and a decrease in habitat quality. The expansion rate of construction land is increasing and exhibiting aggregation, while agricultural land area is expanding and forest land area is continuously decreasing. (3) Land use intensity exhibits a significant positive correlation with land conflict levels, while land conflict levels demonstrate a significant negative correlation with habitat quality. These findings suggest that land use has had some impact on the ecological environment, with indications of habitat degradation. Even in Ankang, where ecological development is highly valued, the city will gradually face conflicts between ecological protection and economic development in future scenarios. The study highlights that Ankang's future development space will be constrained within the context of environmental protection, leading to greater land use conflicts in urban and surrounding areas. Consequently, the quality of habitats will inevitably decline. Therefore, it is crucial to allocate sufficient space for economic development while simultaneously prioritizing ecological protection. This approach will ensure a healthy economic development trajectory and foster a harmonious relationship between humans and the land.

Little is known about how differences in water availability within the “super humid” tropics can influence the distribution of understory plant species and the composition of understory plant communities. We investigated variation in the physiological drought tolerances of understory plants and plant communities across a large elevation and precipitation gradient. We established 58 understory plots along a gradient of 400 – 3600 m asl elevation and 1000 – 6000 mm yr-1 rainfall in and around Manu National Park in southeastern Peru. Within the plots, we sampled all understory woody plants and measured three metrics of physiological leaf drought tolerance - turgor loss point (TLP), cuticular conductance (Gmin), and solute leakage (SL) - and assessed how the community-level means of these traits related to mean annual precipitation (MAP) and elevation (in the study gradient temperature decreases linearly and vapor pressure deficit increases monotonically with elevation). We did not find any correlations between the three metrics of leaf drought tolerance, suggesting that they represent independent strategies for coping with low water availability. Despite being widely used metrics of leaf drought tolerance, neither TLP nor Gmin showed any significant relationships with elevation or MAP. In contrast, SL, which has only recently been developed for use in ecological field studies, increased significantly at higher precipitations and at lower elevations (i.e., plants in colder and drier habitats have lower SL, indicating greater drought tolerances). Our results illustrate that differences in water availability may play a strong role in structuring plant communities even in in the wet tropics and highlight the potential for SL assays to be an efficient and effective tool for measuring drought tolerances in the field.

SOC stock (SOCS) changes caused by land use changes are still unclear, and understanding this response is essential for many environmental policies and land management practices. In this study, we investigated the temporal-spatial and vertical distribution characteristics of SOCS on the Western Sichuan Plateau (WSP) using the sparrow search algorithm-random forest regression (SSA-RFR) models with excellent model applicability and accuracy. The temporal-spatial varia-tions in the SOCS were modeled using 1080 soil samples and a set of nine environmental covariates. We analyzed the effect of land use changes on the SOCS on the WSP. The total SOCS increased by 18.03 Tg C from 1990 to 2020. The results of this study confirmed that the SOCS in the study area has increased significantly since 2010, with an increase of 27.88 Tg C compared to the total SOCS in 2010. We found that the spatial distribution of the SOCS increased from southeast to northwest, and the vertical distribution of the SOCS in the study area decreased with increasing soil depth. Forest and grassland are the main sources of SOCS, the total SOCS in the forest and grassland accounted for 37.53 and 59.39% of the total SOC pool in 2020, respectively. Expansion of the wetlands, forest, and grassland areas could increase SOCS in the study area.

Monitoring snowmelt dynamics in mountains is crucial to understand water releases downstream. Sentinel-1 (S-1) synthetic aperture radar (SAR) has become one of the most widely used techniques to achieve this aim due to its high frequency of acquisitions and all-weather capability. This work aims to understand the possibilities of S-1 SAR imagery to capture snowmelt dynamics and related changes in streamflow response in semiarid mountains. The results proved that S-1 SAR imagery was able not only to capture the final spring melting but also all melting cycles that commonly appear throughout the year in these types of environments. The general change detection approach to identify wet snow was adapted for these regions using as reference the average S-1 SAR image from the previous summer, and a threshold of -3.00 dB. In addition, four different type of melting runoff onsets depending on physical snow condition were identified. When translating that at the catchment scale, distributed melting runoff onset maps were defined to better understand the spatiotemporal evolution of melting dynamics. Finally, a linear connection between melting dynamics and streamflow was found for long-lasting melting cycles, with a determination coefficient (R2) ranging from 0.62 to 0.83 and an average delay between the melting onset and streamflow peak of about 21 days.

The altitude of the habitat is one of the important regulators of species survival. Kaçkar Mountains National Park is located in the Eastern Black Sea region of Turkey. This is the first study on the benthic diatom flora of the high mountain lakes in Kaçkar Mountains National Park, which is situated between 2782 and 3075 m a.s.l. A total of 84 diatom species were identified from benthic communities of 15 habitats in summer (19 July, 28 August) and autumn (10 September) months of 2020. The genus Pinnularia (13 species) formed the basis of the taxonomic list followed by Eunotia (5 species), Navicula (5 species), and Frustulia (4 species) genera, respectively. The waters in all the studied lakes were fresh, low-saline, with low-alkaline or circumneutral pH and organically uncontaminated as evidenced by prevailed bioindicator groups. Statistical methods and comparative floristic results confirm the role of the lake altitude for the diatom species distribution. Species richness was higher in lakes with lower altitude, but a statistical approach also revealed a tendency for it to increase with altitude if species richness increased. Further studies will be needed to continue exploring this pattern. In the first round of analysis, diatom species groups were categorized according to IUCN categories, revealing one potentially threatened and one rare species. To protect studied high mountain lakes, their ecological conditions must be constantly monitored in the Kaçkar Mountains National Park.

Snow-covered regions are the main source of reflection of incident shortwave radiation (ISR) on Earth’s surface. The deposition of light-absorbing particles (LAPs) on these regions increases the capacity of snow to absorb ISR and decreases surface snow albedo (SSA), which intensifies the radiative forcing leading to accelerated snowmelt and modifications of the hydrologic cycle. In this work we investigate changes in SSA and radiative forcing (RF) induced by LAPs in the Upper Aconcagua River Basin (Chilean central Andes) using remote sensing satellite data (MODIS), in-situ spectral SSA measurements, and the ISR (Chilean Solar Explorer platform) during the austral-winter months (May to August) for the 2004-2016 period. To estimate the changes in SSA and RF, we define two spectral ranges: i) an enclosed range (Ecr) between 841-876 nm, which isolates effects of Black-Carbon, an important LAP derived from anthropogenic activities, and ii) a broadband range (Bbr) between 300-2500 nm. Our results show that percent variations in SSA in the Ecr are higher than in the Bbr, regardless of the total amount of radiation received, which may be attributed to the presence of LAPs as these particles have a greater impact on SSA at wavelengths in the Ecr band than in the Bbr band.

In both managed and unmanaged forests, termites are functionally important members of the dead-wood-associated (saproxylic) insect community. However, little is known about regional-scale environmental drivers of geographic distributions of termite species, and how these environmental factors impact co-occurrence among congeneric species. Here we focus on the southern Appalachian Mountains—a well-known center of endemism for forest biota—and use Ecological Niche Modeling (ENM) to examine the distributions of three species of Reticulitermes termites (i.e., R. flavipes, R. virginicus, and R. malletei). To overcome deficiencies in public databases, ENMs were underpinned by field-collected high-resolution occurrence records coupled with molecular taxonomic species identification. Spatial overlap among areas of predicted occurrence of each species was mapped, and aspects of niche similarity were quantified. We also identified environmental factors that most strongly contribute to among-species differences in occupancy. Overall, we found that R. flavipes and R. virginicus showed significant niche divergence, which was primarily driven by dry-season precipitation. Also, all three species were most likely to co-occur in the mid-latitudes of the study area (i.e., northern Alabama and Georgia, eastern Tennessee and western North Carolina), which is an area of considerable topographic complexity. This work provides important baseline information for follow-up studies of local-scale drivers of these species’ distributions. It also identifies specific geographic areas where future assessments of the frequency of true syntopy vs. micro-allopatry, and associated interspecific competitive interactions, should be focused.

Iturbide is located in the Northeast of Mexico, it has a rich native and exotic flora, however, there are no ethnobotanical records, therefore, it requires attention in the documentation of traditional practices and uses of its botanical resources. In 2021, twelve field trips were carried out, applying 110 semi-structured interviews. Plant samples were collected, identified and deposited in an herbarium. We used the Chi-square test to compare ethnobotanical uses with respect to others reported in Mexico. To determine the cultural importance, three ethnobotanical indices were applied (UVI, ICF and FL). We recorded 250 species with ethnobotanical uses associated with 121 genera and 83 families, including 140 native and 110 exotic species. The most common plant families were Asteraceae, Lamiaceae, and Fabaceae. The main uses were ornamental, medicinal and food. The species with the highest UVI values ​​were Lepidium peruvianum, Ocimum basilicum and Rosmarinus officinale. The multifunctionality of the native and exotic flora demonstrates the extensive knowledge associated with botanical resources. Examples, the role of ornamental plants, with a direct impact on human well-being, the resilience of healers and traditional inhabitants by using different species for the treatment of various ailments. or indigenous edible plants in the daily diet.