Ecology

RESEARCH BRIEFINGS
07 September 2022

A diverse range of marine predators — including tunas, billfishes and sharks — in the North Pacific Ocean cluster together in clockwise-rotating eddies, seemingly to hunt deep-ocean prey, which are unusually abundant there. This suggests that there is a relationship between the foraging opportunities of predators and the energetics of this marine biome. More

We used two agent-based simulation models to investigate the concepts of “cooperate for the spread” and “pay for the escape,” both were net logo models created by Dr. Susan Hanisch.Afterward, we modified the first model to represent the concept of sharing the dispersal costs. We used the second model without modifications. Instead, we assigned definite values of some parameters that highlight the pay for the escape strategy.First modelThe original model was entitled “Evolution and patchy resource”18. She first developed it for educational purposes. It illustrates the concepts of cooperator-cheater competition, natural selection, spatial structure mechanisms, multilevel selection, and founder effects.Changeable variables

Distance-resource-areas: the distance between the centers of the resource areas.

Size-resource areas: the size of resource areas as a radius in the number of patches.

Living costs: the costs that each agent has to deduct from energy per iteration for basic survival.

Mutation rate: The probability that offspring agents have different traits than their parents.

Evolution: the ability of agents to produce offspring.

Constant variables

The number of patches is 112 × 112 patches.

Carrying capacity per patch: Resource = 10, Agents = 1

The growth rate of the resource = 0.2

The resources on a patch regrow by a logistic growth function up to the carrying capacity: New resource level = current resource level + (Growth-Rate × current resource level) × (1 – (Current resource level/carrying capacity)).

The cost for producing offspring is ten subtracted units of energy.

The initial level of energy of agents is set at living costs.

Role of randomness

Agents are distributed randomly in resource areas at the beginning of a simulation.

Sustainable behavior is distributed randomly with a probability of percent sustainables among the initial agent population.

The order in which agents move and harvest within one iteration is random.

Agents move to a randomly selected patch if several patches fulfill the objectives.

The order in which agents produce offspring within one iteration is random.

Agents reproduce offspring with a probability of (0.0005 × Energy).

Agents place offspring on a randomly selected unoccupied neighboring patch.

Offspring mutate with a potential mutation rate.

Model processesIn each iteration, each agent moves around in random order. There are three likelihoods:

If there are no unoccupied patches in a two-patch radius, they stay on the current patch.

If there are unoccupied patches with resources amounting to more than living costs, the agents move to them.

If the resource amount is less than the living costs, the agents move randomly to other unoccupied patches.

The agents harvest the resources from separated patches to gain energy for metabolism and proliferation. If the energy level of any agent falls to zero, it dies. The cooperator type harvests half of the resource, while the greedy type consumes 99%.The living costs are deducted from the energy amount of the agent constantly everywhere all the time. This process occurs whether an agent moves within the patch, between the patches, or even not. Therefore, the model does not consider dispersal cost explicitly.If there is an unoccupied neighbor patch, the agent can reproduce with a probability of 0.0005 of his energy, place the offspring on the unoccupied neighbor patch, and then transfer ten units of the energy to his offspring.Resources regrow only on resource patches. When the resource amount is more than or equal to 0.1, then it regrows. When the resource is less than 0.1, its value is set to 0.1.Output diagrams and monitors

The average energy of agents: average energy levels of sustainable and greedy agents, resulting from resource harvest minus living costs and reproduction.

Trait frequencies: the relative frequencies of sustainable and greedy agents in the total population, resulting from mutations, different reproduction rates, and death.

Agent population: the absolute number of the total population size resulting from reproduction and death.

ModificationsIn the first modification, we added a different type of cost that agents only incur when they disperse from one patch to another (in-between the patches). It is the slider entitled “dispersal costs”.In the second modification, we added another sharing dispersal costs tool to reduce them by dividing their value by the number of included agents (flock-mates) in the identified range from the same type. It is the slider entitled “group-dispersal-range.” which is the flock mate’s areas as a radius in the number of patches. Therefore, changing the value of the group dispersal range will change the area around every agent. Accordingly, the number of its flock mates who share the dispersal costs also adjusts.The group dispersal range is not confined to greedy agents but applies to all agents. Therefore, it represents the case of the wild-type cooperators who can also cooperate for the spread. The group dispersal range also does not only target the agents in between patches. However, it counts the agents inside and outside the patches. For example, once an agent starts its dispersion with a determined range containing ten agents, four from another type, three non-dispersal agents from the same type that existed inside a patch, and three dispersal agents from the same type outside the patches. The dispersal costs for this agent will be divided by 6.Our assumption that non-dispersal agents at the pre-departure stage share dispersion costs with dispersal agents; seems justified because they reap mutual benefits by reducing kin competition inside patches if they promote the migrators. However, can agents remotely pay the dispersion costs? Yes. For instance, some bacterial species can trigger the migration of other species if located in their vicinity, even if the two bacterial colonies are separated by a barrier19,20 or if they are non-motile21. On the other hand, dispersion is an extended process with many factors, including escape from predators, suppression of host defense mechanisms, and production of biosurfactants to reduce surface tension to facilitate motility. Therefore, the agent’s contribution (inside/outside the patches) to support such factors is considered a shared dispersal cost.Finally, cheaters can arise within cooperator patches by mutation or immigration. Therefore, to investigate the efficacy of migration, the mutation rate value should be 0 to cancel its effect in the meta-population dynamics.Second modelThe model is entitled “Evolution, resources, monitoring, and punishment.”22 is a simulation of a population with four types of agents competing for the same resource. It demonstrates many concepts, such as kin selection, cooperation, selfishness, public good, monitoring, punishment, sharing the costs, positive/negative frequency-dependent selection, and multilevel selection. The four agent colors and types: (1) Red: greedy, non-punishing. (2) Orange: greedy, punishing. (3) Turquoise: sustainable, non-punishing. (4) Green: sustainable, punishing.Punishing agents can perceive other agents in their environment to some degree (perception accuracy) and react to their behavior. There are three kinds of punishment: Punishers can kill agents with greedy harvesting behavior, stop them from harvesting in the next iteration, or have them pay a penalty fee to their neighbors.Agents have a cost (energy) to pay for, both detection and punishment, so this behavior is altruistic. Punisher agents of one type share punishment costs equally.Changeable variables

Death rate: the probability that agents die independent of their energy level.

Carrying capacity: the maximum amount of resource units on a patch from 1 to 100.

Growth rate: the rate at which resources on patches regrow. The maximum sustainable yield is calculated based on the carrying capacity and growth rate.

Harvest-sustainable: the number of resource units harvested by sustainable agents.

Harvest-greedy: the number of resource units harvested by sustainable agents.

Perception accuracy: the probability with which punishing agents notice greedy agents.

Costs-perception: the costs in units of energy, punishing agents have to pay for perceiving other agents.

Costs-punishment: the costs as units of energy that punishing agents have to pay in each iteration to punish other agents. All punishing agents of an agent divide the costs of punishment.

Punishment: the kinds of punishing behavior that punishing agents perform.

Fine: if the kind of punishment is “pay fine”, the fine in energy units that punished agents have to pay (shared between all their neighbors).

Living costs and mutation rate: see the first model.

Constant variables

The number of patches: There are 60 × 60 patches in the world.

The initial energy level of agents is set at living costs + 1.

The initial number of resource units on a patch is set to the carrying capacity.

The resources on a patch regrow: see the first model.

Role of randomness* In addition to items in the first model.

Agents take on their traits (harvest preference and ability to notice and punish) randomly based on the probability of percent-sustainable and percent-punishers.

The order in which punishing agents notice greedy agents within one iteration is random.

Greedy agents are noticed by punishing agents with a probability of perception accuracy.

The order in which detected greedy agents are punished within one iteration is random.

Agents produce offspring with a probability of (0.001 × Energy).

Agents die with a probability of (death-rate).

Model processesIn each iteration, each agent attempts to harvest resources from the patches it is on and the eight neighboring patches until the harvest preference level is reached, except for the punished agent with the sanction (suspend harvest once), its harvest amount = 0 in the current iteration. If the amount of resources available is lower than the amount that the unpunished agent attempts to harvest. Then, the agent moves to a neighboring unoccupied patch with the most resources after losing one energy unit as a move cost.Punishers pay the costs of perceiving the greedy agents. The greedy neighbors have been noticed with the probability of perception accuracy. The agent lost an amount of energy as living costs. The agent dies with the likelihood of death rate or if the energy level falls to zero.If there is an unoccupied neighbor patch, the agent can reproduce with a probability of 0.001 of its energy, place the offspring on the unoccupied neighbor patch, and then transfer half of its energy to its offspring that mutate according to the probability of the mutation rate.Resources regrow on all patches. When the resource amount is more than or equal to 0.1, then it regrows. When the resource is less than 0.1, its value is set to 0.1.Output diagrams and monitors

Populations (% of carrying capacity): the state of the resource and the agent population in the world as a percentage of total carrying capacity resulting from resource harvesting behavior and resource regrowth, agent reproduction, and death.

Average harvest per iteration: the average harvested amounts of agents per iteration by trait, resulting from harvested resource units, minus costs for monitoring and punishing (for punishing agents), minus fines (for punished agents in case of punishment “Pay fine”)

The average energy of agents and trait frequencies: see the first model.

How does the model represent a conditional defector strategy?The model aims to highlight the role of kin selection and punishment mechanisms in supporting cooperation evolution against cheats. We did not need to modify the model but just thought about what the conditional defector should do to upside down the game. The answer was to pay for the escape.For instance, if the standard Harvest-greedy of a cheater (greedy, non-punishing) was 13 and the Perception-accuracy of its actual punishers was 75%. Now suppose this cheater faces troubles, and it cannot dominate. However, if it gives up some of its profit to become 12, to escape punishment, and to reduce the perception accuracy to 60%, it could dominate and take over the population.The conditional cheater can pay something and reduce its profit to escape punishment by reducing perception accuracy if there is a positive correlation between these two variables. Therefore, this model is appropriate if it can support/deny such a correlation. More

Taco-Taype, N. & Zúñiga-Dávila, D. Efecto de la inoculación de plantas de Tarwi con cepas de Bradyrhizobium spp. aisladas de un lupino silvestre, en condiciones de invernadero. Revista peruana de biología. 27, 35–42 (2022).Article 

Google Scholar 
Atchison, G. W. et al. Lost crops of the Incas: Origins of domestication of the Andean pulse crop tarwi Lupinus mutabilis. Am. J. Bot. 103, 1592–1606 (2016).CAS 
Article 

Google Scholar 
Peru Origins. Tarwi (Lupinus Mutabilis). https://peruorigins.com/tarwi/ (2022).Guilengue, N., Alves, S., Talhinhas, P. & Neves-Martins, J. Genetic and genomic diversity in a tarwi (Lupinus mutabilis Sweet) germplasm collection and adaptability to Mediterranean climate conditions. Agronomy 10, 21 (2020).Article 

Google Scholar 
Repo-Carrasco-Valencia, R., Basilio-Atencio, J., Luna-Mercado, G. I., Pilco-Quesada, S. & VidaurreRuiz, J. Andean ancient grains: Nutritional value and novel uses. Biol. Life Sci. Forum. https://doi.org/10.3390/blsf2021008015 (2022).Article 

Google Scholar 
Gulisano, A., Alves, S., Martins, J. N. & Trindade, L. M. Genetics and breeding of Lupinus mutabilis: An emerging protein crop. Front. Plant Sci. https://doi.org/10.3389/fpls.2019.01385 (2019).Article 

Google Scholar 
Chen, Y., She, Y., Zhang, R., Wang, J. & Zhang, X. Use of starch-based fat replacers in foods as a strategy to reduce dietary intake of fat and risk of metabolic diseases. Food Sci. Nutr. 8, 16–22 (2020).CAS 
Article 

Google Scholar 
Frick, K. M., Kamphuis, L. G., Siddique, K. H. M., Singh, K. B. & Foley, R. C. Quinolizidine alkaloid biosynthesis in lupins and prospects for grain quality improvement. Front. Plant Sci. https://doi.org/10.3389/fpls.2017.00087 (2017).Article 

Google Scholar 
Chirinos-Arias, M. C. Andean Lupin (Lupinus mutabilis Sweet) a plant with nutraceutical and medicinal potential. Revista Bio. Ciencias. 3, 163–172 (2015).
Google Scholar 
Wink, M. Chemical defense of leguminosae. Are quinolizidine alkaloids part of the antimicrobial defense system of lupins?. Zeitschrift für Naturforschung C. 39, 548–552 (1984).Article 

Google Scholar 
Hidalgo, M. et al. Evaluation of in vitro suceptibility to spartein in four strain of Mycobacterium tuberculosis. Rev. Peru Med Exp Salud Publica. 39, 77–82 (2022).Article 

Google Scholar 
Muñoz, E. B., Luna-Vital, D. A., Fornasini, M., Baldeón, M. E. & Gonzalez de Mejia, E. Gamma-conglutin peptides from Andean lupin legume (Lupinus mutabilis Sweet) enhanced glucose uptake and reduced gluconeogenesis in vitro. J. Funct. 45, 339–347 (2018).Article 

Google Scholar 
Bryant, L., Rangan, A. & Grafenauer, S. Lupins and health outcomes: A systematic literature review. Nutrients 14, 327 (2022).CAS 
Article 

Google Scholar 
Jacobsen, S. & Mujica, A. Geographical distribution of the Andean lupin (Lupinus mutabilis Sweet). Plant Genet. Resour. Newslett. 155, 1–8 (2008).
Google Scholar 
Antunez de Mayolo, S. Nutricion en el antiguo Peru. Banco Central de la Republica. Lima, Peru. 127 (1981).FAO. Perfiles nutricionales por paises: Peru. (ed. FAO) 36 p. (2000).UNICEF. Estado Mundial de la Infancia 2019 incluye a Perú entre las experiencias exitosas de lucha contra la desnutrición crónica infantile. https://www.unicef.org/peru/nota-de-prensa/estado-mundial-infancia-nutricion-alimentos-derechos-peru-experiencias-exitosas-desnutricion-cronica-infantil-reporte (2022).MINSA (Ministry of health – Peru). Situacion actual de la anemia. https://anemia.ins.gob.pe/situacion-actual-de-la-anemia-c1 (2022).WHO. Anemia. https://www.who.int/es/health-topics/anaemia#tab=tab_1 (2022).Galani, Y. J. H., Orfila, C. & Gong, Y. Y. A review of micronutrient deficiencies and analysis of maize contribution to nutrient requirements of women and children in Eastern and Southern Africa. Crit. Rev. Food Sci. Nutr. 62, 1568–1591 (2022).CAS 
Article 

Google Scholar 
White, P. J. & Martin, R. B. Biofortifying crops with essential mineral elements. Trends Plant Sci. 10, 586–593 (2005).Article 

Google Scholar 
White, P. J. & Martin, R. B. Biofortification of crops with seven mineral elements often lacking in human diets-iron, zinc, copper, calcium, magnesium, selenium and iodine. New Phytol. 182, 49–84 (2009).CAS 
Article 

Google Scholar 
Waters, B. M. & Sankaran, R. P. Moving micronutrients from the soil to the seeds: genes and physiological processes from a biofortification perspective. Plant Sciences. 180, 562–574 (2011).CAS 
Article 

Google Scholar 
Brooker, R. W. et al. Improving intercropping: A synthesis of research in agronomy, plant physiology and ecology. New Phytol. 206, 107–117 (2015).Article 

Google Scholar 
Ducsay, L. et al. Possibility of selenium biofortification of winter wheat grain. Plant Soil Environ. 62, 379–383 (2016).CAS 
Article 

Google Scholar 
Kumar, S. & Pandey, G. Biofortification of pulses and legumes to enhance nutrition. Heliyon. https://doi.org/10.1016/j.heliyon.2020.e03682 (2020).Article 

Google Scholar 
Diehn, T. A. et al. Boron demanding tissues of Brassica napus express specific sets of functional Nodulin26-like Intrinsic Proteins and BOR 1 transporters. Plant J. 100, 68–82 (2019).CAS 
Article 

Google Scholar 
Jayalakshmi, V. A., Reddy, T. & Nagamadhuri, K. V. Genetic diversity and variability for protein and micro nutrients in advance breeding lines and chickpea varieties grown in Andhra Pradesh.”. Legume Res. Int. J. 42, 768–772 (2019).
Google Scholar 
Bouis, H. & Saltzman, A. Improving nutrition through biofortification: A review of evidence from HarvestPlus, 2003 through 2016. Glob Food Sec. 12, 49–58 (2017).Article 

Google Scholar 
Sanca, D. Composición nutricional de diez genotipos de lupino (L. mutabilis y L. albus) desamargados por proceso acuoso. Thesis. Universidad Nacional Agraria La Molina. (2015).Rodríguez, A. Evaluación “in vitro” de la actividad antibacteriana de los alcaloides del agua de desamargado del chocho (Lupinus mutrabilis Sweet). Thesis. Escuela Superior Politécnica de Chimborazo, Ecuador (2009).Villacres, E. et al. Germination, an effective process to in-crease the nutritional value and reduce non-nutritive factors of lupine grain (Lupinus mutabilis Sweet). Int. J. Food Sci. Nutr. Eng. 5, 163–168 (2015).
Google Scholar 
Villacres, E., Rubio, A., Egas, L., Segovia, G. Usos alternativos del chocho: Chocho (Lupinus mutabilis Sweet) alimento andino redescubierto. IOP publishing: repositorio. https://repositorio.iniap.gob.ec/handle/41000/298 (2006).Ortega-David, E. A., Rodríguez, A. D. & Burbano, A. Z. Caracterización de semillas de lupino (Lupinus mutabilis) sembrado en los Andes de Colombia. Acta Agronómica. 59, 111–118 (2010).
Google Scholar 
White, P. J. & Broadley, M. R. Physiological limits to zinc biofortification of edible crops. Front Plant Sci. 80, 1–11 (2011).
Google Scholar 
Zhao, F., Su, Y. H., Dunham, S. J. & Rakszegiet, M. Variation in mineral micronutrient concentrations in grain of wheat lines of diverse origin. J. Cereal Sci. 49, 290–295 (2009).CAS 
Article 

Google Scholar 
Uauy, C., Distelfeld, A., Fahima, T., Blechl, A. & Dubcovsky, J. A NAC Gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science 24, 1298–1301 (2006).ADS 
Article 

Google Scholar 
Shorrocks, V. M. The occurrence and correction of boron deficiency. Plant Soil 193, 121–148 (1997).CAS 
Article 

Google Scholar 
D’Imperio, M. et al. Boron biofortification of Portulaca oleracea L. through soilless cultivation for a new tailored crop. Agronomy. 10, 999–1013 (2020).Article 

Google Scholar 
Boyacioglu, O., Orenay-Boyacioglu, S., Yildirim, H. & Korkmaz, M. Boron intake, osteocalcin polymorphism and serum level in postmenopausal osteoporosis. J. Trace Elem. Med. Biol. 48, 52–56 (2018).CAS 
Article 

Google Scholar 
Oliveira Araújo, E., Ferreira Dos Santos, E. & Camacho Oliveira, M. A. Boron-zinc interaction in the absorption of micronutrients by cotton. Agronomía Colombiana. 36, 51–57 (2018).Article 

Google Scholar 
Squitti, R., Siotto, M. & Polimanti, R. Low-copper diet as a preventive strategy for Alzheimer’s disease. Neurobiol. Aging 2, 40–50 (2014).Article 

Google Scholar 
Schilsky, M.L. Management of Wilson Disease (A Pocket Guide), 1st ed.; Publisher: Humana Press, Farmington, CT, USA. 154–196 (2018).Martins, A. C. et al. Manganese in the diet: Bioaccessibility, adequate intake, and neurotoxicological effects. J. Agric. Food Chem. 46, 12893–12903 (2020).Article 

Google Scholar 
Falah, S. A. & Saja, N. M. Essential trace elements and their vital roles in human body. Indian J. Adv. Chem. Sci. 3, 127–136 (2017).
Google Scholar 
National institutes of health. Manganese. Fact Sheet for Health Professionals. IOP Publishing ods.od.nih.gov. https://ods.od.nih.gov/factsheets/Manganese-HealthProfessional/. (2021).Savadi, S. Molecular regulation of seed development and strategies for engineering seed size in crop plants. Plant Growth Regul. 84, 401–422 (2018).CAS 
Article 

Google Scholar 
Ge, L. et al. (2016) Increasing seed size and quality by manipulating BIG SEEDS1 in legume species. Proc Natl Acad Sci. 113, 12414–12419 (2016).CAS 
Article 

Google Scholar 
Zou, L. Effects of gradual and sudden heat stress on seed quality of Andean lupin, Lupinus mutabilis. Thesis. University of Helsinki. https://helda.helsinki.fi/handle/10138/16501 (2009).Buircell, B.J., Cowling, A.W. Genetic Resources in Lupins (eds. Gladstones, J.S., Atkins, C.A., Hamblin, J.) (United Kingdom: CAB International, 1998).Aguilar-Angulo, L. A. Evaluación del rendimiento de grano y capacidad simbiótica de once accesiones de tarwi (Lupinus mutabilis Sweet), bajo condiciones de Otuzco-La Libertad (Universidad Nacional Agraria La Molina, 2015).
Google Scholar 
De La Cruz, N. Caracterización fenotípica y de rendimiento preliminar de ecotipos de tarwi (Lupinus mutabilis sweet), bajo condiciones del Callejón de Huaylas – Ancash (Universidad Nacional Agraria la Molina, 2018).
Google Scholar 
Huisa, J. Evaluación del comportamiento agronómico de catorce accesiones del ensayo nacional de tarwi (Lupinus mutabilis sweet.) en el CIP Camacani Puno – Perú”. Thesis. Universidad Nacional Agraria la Moina (2018).Cayo, B. Evaluación del comportamiento agronómico de ocho genotipos selectos de tarwi (Lupinus mutabilis sweet) bajo condiciones del CIP. CAMACANI – UNA – PUNO. Thesis. Universidad Nacional del Altiplano (2020).Buircell, B.J., Cowling, A.W. Lupin. Lupinus spp. Promoting the conservation and use of underutilized and ne-glected crops (eds. Gladstones, J.S., Atkins, C.A., Hamblin, J.) (United Kingdom: CAB International, 1998).Plata, J. Comportamiento Agronómico de dos Variedades de tarwi (Lupinus mutabilis Sweet), bajo tres densidades de siembra en la comunidad Marka Hilata Carabuco (Universidad San Andres, 2016).
Google Scholar 
Mendoza, C. Rendimiento de ecotipos regionales y variedades de tarwi (Lupínus mutabilis Sweet.) en el valle del Mantaro, Jauja, Junín. Thesis. Universidad Nacional Agraria la Moina (2020).Aguilar, S. Sistemas de producción de Lupinus mutabilis Sweet ‘chocho’ en terrazas y laderas con fertilización fosfatada en Cajamarca. Dissertation. La Molina National Agrarian University (2011).Aquino, S. Sustentabilidad del cultivo de tarwi (Lupinus mutabilis sweet) en la zona altoandina del Valle del Mantaro (Universidad Nacional Agraria la Molina, 2018).
Google Scholar 
Barda, M. S., Chatzigeorgiou, T., Papadopoulos, G. K. & Bebeli, P. J. Agro-morphological evaluation of Lupinus mutabilis in two locations in greece and association with insect pollinators. Agriculture https://doi.org/10.3390/agriculture11030236 (2021).Article 

Google Scholar 
Herniter, I. A., Jia, Z. & Kusi, F. Market preferences for cowpea (Vigna unguiculata [L.] Walp) dry grain in Ghana. African J Ag Res. 14, 928–934 (2019).Article 

Google Scholar 
Dordas, C. Foliar boron application affects lint and seed yield and improves seed quality of cotton grown on calcareous soils. Nutr. Cycl. Agroecosyst. 76, 19–28 (2006).CAS 
Article 

Google Scholar 
Kristek, S. et al. Effect of various rates of boron on yield and quality of high-grade sugar beet varieties. Listy Cukrovarnické a Řepařské. 4, 146–150 (2018).
Google Scholar 
Thomas, C. L. et al. Root morphology and seed and leaf ionomic traits in a Brassica napus L. diversity panel show wide phenotypic variation and are characteristic of crop habit. BMC Plant Biol. 16, 214–232 (2016).CAS 
Article 

Google Scholar 
Dursun, A. et al. Effects of boron fertilizer on tomato, pepper and cucumber yields and chemical composition. Commun Soil Sci Plant Anal. 1, 1576–1593 (2010).Article 

Google Scholar 
Sotiropoulos, T. E., Therios, T. N., Dimassi, K. N., Bosabalidis, A. & Kofidis, G. Nutritional status, growth, CO2 assimilation, and leaf anatomical responses in two kiwifruit species under boron toxicity. J Plant Nutr. 25, 1249–1261 (2002).CAS 
Article 

Google Scholar 
Muccifora, S. & Bellani, L. Effects of copper on germination and reserve mobilization in Vicia sativa L. seeds. Environ. Pollut. 179, 68–74 (2013).CAS 
Article 

Google Scholar 
Kobraee, S. Effect of foliar fertilization with zinc and manganese sulfate on yield, dry matter accumulation, and zinc and manganese contents in leaf and seed of chickpea (Cicer arietinum). J. Appl. Biol. Biotechnol. 7, 20–28 (2019).CAS 

Google Scholar 
IBPGR (1981) Lupin descriptors. https://www.bioversityinternational.org/fileadmin/bioversity/publications/Web_version/103/ (1981).Zasoski, R. J. & Burau, R. G. A rapid nitric-perchloric acid digestion method for multi-element tissue analysis. Commun. Soil Sci. Plant Anal. 8, 425–436 (1997).Article 

Google Scholar 
Pereira, T., Coelho, C. M. M., Bogo, A., Guidolin, A. F. & Miquelluti, D. J. Diversity in common bean landraces from south Brazil. Acta Bot. Croat. 1, 79–92 (2009).
Google Scholar 
Pujar, M., Govindaraj, M., Gangaprasad, S., Kanatti, A. & Shivade, H. Genetic variation and diversity for grain iron, zinc, protein and agronomic traits in advanced breeding lines of pearl millet [Pennisetum glaucum (L.) R Br] for biofortification breeding. Genet. Resour. Crop Evol. 67, 2009–2022 (2020).CAS 
Article 

Google Scholar 
Lira, J. P. E. et al. Safflower genetic diversity based on agronomic characteristics in Mato Grosso state, Brazil, for a crop improvement program. Genet. Mol. Res. 1, 1–12 (2021).
Google Scholar 
de Sá, S. F. et al. Genetic diversity via REML-BLUP of ex situ conserved macauba [Acrocomia aculeata (Jacq.) Lodd. ex Mart.] ecotypes. Genet. Resour. Crop Evol. 68, 3193–3204 (2021).Article 

Google Scholar 
Kuru, R., Yilmaz, S., Tasli, P. N., Yarat, A. & Sahin, F. Boron content of some foods consumed in Istanbul, Turkey. Biol. Trace Elem. Res. 187, 1–8 (2019).CAS 
Article 

Google Scholar 
Shokunbi, O., Adepoju, O., Mojapelo, P., Ramaite, I. & Akinyele, I. Copper, manganese, iron and zinc contents of Nigerian foods and estimates of adult dietary intakes. J. Food Compos. Anal. 82, 103–245 (2019).Article 

Google Scholar 
Norwegian scientific committee for food and environment. Assessment of dietary intake of manganese in rela-tion to tolerable upper intake. IOP Publishing wkm. www.vkm.no. (2018).Gil, V., Guzmán, L. & Quintero, E. Caracterización de la variabilidad morfológica de un “genotipo local” de maíz y dos de sus selecciones. Centro Agrícola. 4, 79–83 (2004).
Google Scholar  More

Chaigneau, A., Gizolme, A. & Grados, C. Mesoscale eddies off Peru in altimeter records: identification algorithms and eddy spatio-temporal patterns. Prog. Oceanogr. 79, 106–119 (2008).ADS 
Article 

Google Scholar 
McGillicuddy, D. J. Jr et al. Influence of mesoscale eddies on new production in the Sargasso Sea. Nature 394, 263–266 (1998).ADS 
CAS 
Article 

Google Scholar 
Dufois, F. et al. Anticyclonic eddies are more productive than cyclonic eddies in subtropical gyres because of winter mixing. Sci. Adv. 2, 1–7 (2016).Article 

Google Scholar 
Godø, O. R. et al. Mesoscale eddies are oases for higher trophic marine life. PLoS ONE 7, e30161 (2012).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Chelton, D. B., Gaube, P., Schlax, M. G., Early, J. J. & Samelson, R. M. The influence of nonlinear mesoscale eddies on near-surface oceanic chlorophyll. Science 334, 328–333 (2011).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Sarmiento, J. L. et al. Response of ocean ecosystems to climate warming. Global Biogeochem. Cycles 18, GB3003 (2004).ADS 
Article 
CAS 

Google Scholar 
Bell, J. D. et al. Diversifying the use of tuna to improve food security and public health in Pacific Island countries and territories. Mar. Policy 51, 584–591 (2015).Article 

Google Scholar 
Della Penna, A. & Gaube, P. Mesoscale eddies structure mesopelagic communities. Front. Mar. Sci. 7, 454 (2020).ADS 
Article 

Google Scholar 
Braun, C. D. et al. The functional and ecological significance of deep diving by large marine predators. Ann. Rev. Mar. Sci. 14, 129–159 (2022).PubMed 
Article 

Google Scholar 
McGillicuddy, D. J. Jr Mechanisms of physical-biological-biogeochemical interaction at the oceanic mesoscale. Ann. Rev. Mar. Sci. 8, 125–159 (2016).PubMed 
Article 

Google Scholar 
Fennell, S. & Rose, G. Oceanographic influences on deep scattering layers across the North Atlantic. Deep-Sea Res. Part I Oceanogr. Res. Pap. 105, 132–141 (2015).ADS 
Article 

Google Scholar 
Duffy, L. M. et al. Global trophic ecology of yellowfin, bigeye, and albacore tunas: understanding predation on micronekton communities at ocean-basin scales. Deep-Sea Res. Part II Topical Stud. Oceanogr. 140, 55–73 (2017).ADS 
Article 

Google Scholar 
Gaube, P. et al. Mesoscale eddies influence the movements of mature female white sharks in the Gulf Stream and Sargasso Sea. Sci. Rep. 8, 7363 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Braun, C. D., Gaube, P., Sinclair-Taylor, T. H., Skomal, G. B. & Thorrold, S. R. Mesoscale eddies release pelagic sharks from thermal constraints to foraging in the ocean twilight zone. Proc. Natl Acad. Sci. USA 116, 17187–17192 (2019).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Doyle, T. K. et al. Leatherback turtles satellite-tagged in European waters. Endanger. Species Res. 4, 23–31 (2008).Article 

Google Scholar 
Pauly, D. & Christensen, V. Primary production required to sustain global fisheries. Nature 374, 255–257 (1995).ADS 
CAS 
Article 

Google Scholar 
Lynham, J., Nikolaev, A., Raynor, J., Vilela, T. & Villaseñor-Derbez, J. C. Impact of two of the world’s largest protected areas on longline fishery catch rates. Nat. Commun. 11, 979 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Polovina, J. J., Abecassis, M., Howell, E. A. & Woodworth, P. Increases in the relative abundance of mid-trophic level fishes concurrent with declines in apex predators in the subtropical North Pacific, 1996-2006. Fish. Bull. 107, 523–531 (2009).
Google Scholar 
Royer, T. C. Ocean eddies generated by seamounts in the North Pacific. Science 199, 1063–1064 (1978).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Liu, Y. et al. Eddy analysis in the subtropical zonal band of the North Pacific Ocean. Deep-Sea Res. Part I Oceanogr. Res. Pap. 68, 54–67 (2012).ADS 
Article 

Google Scholar 
Bernstein, R. L. & White, W. B. Time and length scales of baroclinic eddies in the central North Pacific Ocean. J. Phys. Oceanogr. 4, 613–624 (1974).ADS 
Article 

Google Scholar 
Maunder, M. N. & Punt, A. E. Standardizing catch and effort data: a review of recent approaches. Fish. Res. 70, 141–159 (2004).Article 

Google Scholar 
Woodworth, P. A. et al. Eddies as offshore foraging grounds for melon-headed whales (Peponocephala electra). Mar. Mammal Sci. 28, 638–647 (2012).Article 

Google Scholar 
Gaube, P. et al. The use of mesoscale eddies by juvenile loggerhead sea turtles (Caretta caretta) in the southwestern Atlantic. PLoS ONE 12, e0172839 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Chambault, P. et al. Swirling in the ocean: immature loggerhead turtles seasonally target old anticyclonic eddies at the fringe of the North Atlantic Gyre. Prog. Oceanogr. 175, 345–358 (2019).ADS 
Article 

Google Scholar 
Gaube, P., McGillicuddy Jr, D., Chelton, D., Behrenfeld, M. & Strutton, P. Regional variations in the influence of mesoscale eddies on near-surface chlorophyll. J. Geophys. Res. Oceans 119, 8195–8220 (2014).Waga, H., Kirawake, T. & Ueno, H. Impacts of mesoscale eddies on phytoplankton size structure. Geophys. Res. Lett. 46, 13191–13198 (2019).ADS 
Article 

Google Scholar 
Irigoien, X. et al. Large mesopelagic fishes biomass and trophic efficiency in the open ocean. Nat. Commun. 5, 3271 (2014).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Chen, Y.-lL. et al. Biologically active warm-core anticyclonic eddies in the marginal seas of the western Pacific Ocean. Deep Sea Res. Part I 106, 68–84 (2015).CAS 
Article 

Google Scholar 
Harke, M. J. et al. Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre. Environ. Microbiol. 23, 4807–4822 (2021).CAS 
PubMed 
Article 

Google Scholar 
Hawco, N. J. et al. Iron depletion in the deep chlorophyll maximum: mesoscale eddies as natural iron fertilization experiments. Global Biogeochem. Cycles 35, e2021GB007112 (2021).ADS 
CAS 
Article 

Google Scholar 
Klevjer, T. A. et al. Large scale patterns in vertical distribution and behaviour of mesopelagic scattering layers. Sci. Rep. 6, 19873 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Behrenfeld, M. J. et al. Global satellite-observed daily vertical migrations of ocean animals. Nature 576, 257–261 (2019).CAS 
PubMed 
Article 

Google Scholar 
Madigan, D. J. et al. Water column structure defines vertical habitat of twelve pelagic predators in the South Atlantic. ICES J. Mar. Sci. 78, 867–883 (2021).Article 

Google Scholar 
Arostegui, M., Gaube, P. & Braun, C. Movement ecology and stenothermy of satellite-tagged shortbill spearfish (Tetrapturus angustirostris). Fish. Res. 215, 21–26 (2019).Article 

Google Scholar 
Lehodey, P., Senina, I. & Murtugudde, R. A spatial ecosystem and populations dynamics model (SEAPODYM)—modeling of tuna and tuna-like populations. Prog. Oceanogr. 78, 304–318 (2008).ADS 
Article 

Google Scholar 
Varghese, S. P., Somvanshi, V. S. & Dalvi, R. S. Diet composition, feeding niche partitioning and trophic organisation of large pelagic predatory fishes in the eastern Arabian Sea. Hydrobiologia 736, 99–114 (2014).CAS 
Article 

Google Scholar 
Ward, P. & Myers, R. A. Inferring the depth distribution of catchability for pelagic fishes and correcting for variations in the depth of longline fishing gear. Can. J. Fish. Aquat.Sci. 62, 1130–1142 (2005).Article 

Google Scholar 
Kai, E. T. et al. Top marine predators track Lagrangian coherent structures. Proc. Natl Acad. Sci. USA 106, 8245–8250 (2009).ADS 
CAS 
Article 

Google Scholar 
Lima, I. D., Olson, D. B. & Doney, S. C. Biological response to frontal dynamics and mesoscale variability in oligotrophic environments: biological production and community structure. J. Geophys. Res. Oceans 107, 25-1–25-21 (2002).Article 

Google Scholar 
Spall, S. A. & Richards, K. J. A numerical model of mesoscale frontal instabilities and plankton dynamics—I. model formulation and initial experiments. Deep-Sea Res. Part I Oceanogr. Res. Pap. 47, 1261–1301 (2000).ADS 
Article 

Google Scholar 
Siegelman, L., O’Toole, M., Flexas, M., Rivière, P. & Klein, P. Submesoscale ocean fronts act as biological hotspot for southern elephant seal. Sci. Rep. 9, 5588 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Lévy, M., Ferrari, R., Franks, P. J., Martin, A. P. & Rivière, P. Bringing physics to life at the submesoscale. Geophys. Res. Lett. https://doi.org/10.1029/2012GL052756 (2012).Article 

Google Scholar 
Guidi, L. et al. Does eddy-eddy interaction control surface phytoplankton distribution and carbon export in the North Pacific Subtropical Gyre? J. Geophys. Res. Biogeosciences https://doi.org/10.1029/2012JG001984 (2012).Article 

Google Scholar 
Chow, C. H., Cheah, W., Tai, J. H. & Liu, S. F. Anomalous wind triggered the largest phytoplankton bloom in the oligotrophic North Pacific Subtropical Gyre. Sci. Rep. 9, 15550 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Guo, M., Xiu, P., Chai, F. & Xue, H. Mesoscale and submesoscale contributions to high sea surface chlorophyll in subtropical gyres. Geophys. Res. Lett. 46, 13217–13226 (2019).ADS 
Article 

Google Scholar 
Klein, P. et al. Ocean-scale interactions from space. Earth Space Sci. 6, 795–817 (2019).ADS 
Article 

Google Scholar 
Martin, A. et al. The oceans’ twilight zone must be studied now, before it is too late. Nature 580, 26–28 (2020).ADS 
CAS 
PubMed 
Article 

Google Scholar 
St. John, M. A. et al. A dark hole in our understanding of marine ecosystems and their services: perspectives from the mesopelagic community. Front. Marine Sci. 3, 31 (2016).
Google Scholar 
Bigelow, K., Musyl, M. K., Poisson, F. & Kleiber, P. Pelagic longline gear depth and shoaling. Fish. Res. 77, 173–183 (2006).Article 

Google Scholar 
Brodziak, J. & Walsh, W. A. Model selection and multimodel inference for standardizing catch rates of bycatch species: a case study of oceanic whitetip shark in the Hawaii-based longline fishery. Can. J. Fish. Aquat.Sci. 70, 1723–1740 (2013).Article 

Google Scholar 
Woodworth-Jefcoats, P. A., Polovina, J. & Drazen, J. Synergy among oceanographic variability, fishery expansion, and longline catch composition in the central North Pacific Ocean. Fish. Bull. 116, 228–239 (2018).Article 

Google Scholar 
Boggs, C. H. Depth, capture time, and hooked longevity of longline-caught pelagic fish: timing bites of fish with chips. Fish. Bull. 90, 642–658 (1992).
Google Scholar 
Walsh, W. A. & Brodziak, J. Applications of Hawaii longline fishery observer and logbook data for stock assessment and fishery research. NOAA Tech. Memo. 57, 62 (2016).
Google Scholar 
Walsh, W. A. & Brodziak, J. Billfish CPUE standardization in the Hawaii longline fishery: model selection and multimodel inference. Fish. Res. 166, 151–162 (2015).Article 

Google Scholar 
Gilman, E., Chaloupka, M., Fitchett, M., Cantrell, D. L. & Merrifield, M. Ecological responses to blue water MPAs. PLoS ONE 15, e0235129 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Portner, E. J., Polovina, J. J. & Choy, C. A. Patterns in micronekton diversity across the North Pacific Subtropical Gyre observed from the diet of longnose lancetfish (Alepisaurus ferox). Deep-Sea Research Part I 125, 40–51 (2017).ADS 
Article 

Google Scholar 
Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378–400 (2017).Article 

Google Scholar 
Hartig, F. DHARMa: Residual diagnostics for hierarchical (multi-level/mixed) regression models. R package version 0.3.3.0 http://florianhartig.github.io/DHARMa/ (2020).Jackson, C. H. Multi-state models for panel data: the msm package for R. J. Stat. Softw. https://doi.org/10.18637/jss.v038.i08 (2011).Article 

Google Scholar 
Bates, D. et al. lme4: Linear mixed-effects models using ’Eigen’ and S4. R package version 1.1-25 https://github.com/lme4/lme4/ (2020).Lenth, R. et al. emmeans: Estimated marginal means, aka least-squares mean. R package version 1.7.2 https://github.com/rvlenth/emmeans (2022).R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2020); http://www.r-project.org/ More

In line with the results of present study, we suggest that the exploitation of forage resources by agro-pastoralists occurs in a non-random manner. The use of forage resources is guided by a series of functional characters related to palatability and nutritional value, which determine preferential use due to the better quality of resource. At the same time, we understand that forage uses are complex and multifactorial in nature, and regulated in a substantial way by seasonality and ecological factors (Fig. 5), such as the availability of plant resources and local diversity.Figure 5Diagrammatic representation for the effects of vegetation cover and seasonality on forage resource selection in Dry Forests. Image created with Microsoft Office 2019 PowerPoint (www.office.com).Full size imageThe differences of plant species cited between areas reveal the positive effect of vegetation cover on the use and knowledge of plants by agro-pastoralists. Our findings reveal that the greater number of plant species mentioned by agro-pastoralists in Area II is directly associated with greater availability of resources in this area, as long as we consider vegetation cover as availability of resources, which allows different species to be used throughout the year. On the other hand, in regions with low vegetation cover (Area I), the low availability of resources limits the use and knowledge of plants by residents, which can lead to greater pressure on a small set of available species. Such findings reinforce the importance of vegetation cover for ecosystem provision of goods and services to human populations that depend directly or indirectly on these services.The most represented families found in the present study have also been reported in several other ethnobotanical studies6,16,17,29, with emphasis on Fabaceae and Poaceae, which are recognized for their high forage potential, which derives, above all, from high palatability and nutritional value30. Simultaneously, citations mostly for native species reflect the importance and potential of Caatinga resources as important components of the ruminant diet11, both for the woody and herbaceous strata, corroborating the estimate in the literature that 70% of vegetation has potential use as forage31.The characteristic seasonality of vegetation, on the other hand, represents a limiting factor for forage productivity, culminating in high fluctuations in quality and availability, as well as changes in the dominance of different strata and composition of forage species throughout the seasons11,32. The seasonal distribution of species explains the similarity of seasons between areas, with a higher similarity percentage for the dry seasons, since there is less availability of resources to be exploited compared to the rainy season. In this context, the potentially used species are commonly accessible woody species in both areas. However, during the rainy season, the high availability of herbaceous plants regulates different uses (Fig. 4), but even so, they also exhibit relatively similar patterns, mainly due to the woody component that denotes the common demand by ruminants at the beginning of this season.The effect of climatic variables on vegetation use patterns was documented by16,17, both of which showed greater richness in the use of herbaceous forage during the rainy season, a finding that reflects the seasonal distribution—restriction to that season—and decrease in the qualitative character of annual species33. At the same time, it also reflects the greater number of unique species for the rainy season. However, when compared to woody strata, significant differences in terms of richness are not found because although the diversity of herbaceous species in the Caatinga is greater24, it is much less known than that of the tree-shrub stratum11.Agro-pastoralists even characterize animal preferences for herbaceous stratum, but as its diversity is immense and ephemeral, they claim to have limited ability to identify the species. The high abundance of resources in the rainy season also reduces the concern with forage use, which implies less attention to the species that are consumed. In contrast, woody species, due to multiple uses and greater availability over time, tend to be better known10,34, with a different effect in the dry season making the optimal foraging pattern in this period inherent to the knowledge of agro-pastoralists35.In addition, according to the ecological appearance hypothesis, there is a general tendency for less apparent species to be neglected by populations36. Some studies have corroborated the hypothesis within the context of forage use, with woody species being cited more and having more uses6,15. In addition, people tend to focus on resources whose supply is given continuously10, which may explain why woody species are well represented in both seasons.Security in the provisioning of ecosystem services is an essential component for local populations, and thus woody species are highly valued because they reflect predictability of use15,35. This can be a particularly influential criterion because perennial or late leaf deciduous species, such as Cynophalla flexuosa and Myracrodruon urundeva, had significant amounts of citations and perceptions employing high valuation, as represented by some statements by some interviewees: “É um refrigero na seca” (it is savage in the dry season), “É uma ração boa na seca” (it is a good food in the dry season).In turn, differences in richness of the species cited by the two areas corroborate our first hypothesis that populations inserted in environments with greater vegetation cover tend to cite more species. In line with these findings, considerable floristic dissimilarity was also found between the two areas, given the exclusivity of species. Such dissimilarity may suggest particularities in the vegetation attributes of each area, such as greater floristic diversity7,37,38.Since anthropic processes are irregularly distributed in space, variation in the provisioning of ecosystem services by vegetation also occurs, and influences different collection profiles39. On the other hand, areas with greater species richness have been shown to have greater use patterns6,7. The larger number of species cited as woody and native for Area II is, therefore, associated with greater general richness, as well as herbaceous species present in the rainy season. In contrast, common species are reflected in trends of similar foraging patterns, as well as the presence of common species between areas38. In addition to different levels of disturbance, differences in floristic composition between areas may also be due to edaphic variation40.Our second hypothesis was refuted because the difference in the richness of exotic species between the areas. Plausible explanations for this finding are that, in general, exotic herbaceous species are commonly used for forage in the semi-arid region of Brazil41. Herbaceous species comprise the primary component of the ruminant diet. However, in the midst of their occurrence restricted to the short rainy period, exotic species, mainly of Fabaceae and Poaceae, have been introduced to increase the forage availability, which currently represents an important attribute of forage resources in the Caatinga41,42,43. At the same time, and to also increase the availability of forage resources, the cultivation of species by agro-pastoralists may be common in their properties44, mainly exotics, such as Prosopis juliflora, that have high adaptive potential and governmental incentives45.Regarding use patterns, according to the data presented here it is possible to state that agro-pastoralists ’ experiences with herding activities provide an accumulation of a vast knowledge about forage resources15. This knowledge allows forage resources to be characterized by their potential according to a variety of criteria associated with seasonal variation and qualitative attributes, as commonly found by other studies14,15,16,17,37. Such criteria are often revealed by qualitative approaches that define the valuation perception of resources. Thus, nutritional value and palatability can be implicitly associated with the definitions of “É uma ração boa” (it is a good food), “o bicho gosta muito” (the animals like it very much) and “Rico em proteínas” (rich in protein).It should be added that the establishment of intrinsic relationships with resources allows a particular understanding at a high level of detail15,35, such as changes in palatability throughout development with descriptions including chemical17 and structural changes. Studies confirm that some Caatinga species vary in their chemical composition during leaf maturation, which influences nutritional quality17,46.In addition to revealing the domain of information, this body of knowledge allows maximizing forage use based on nutritional properties weighted by availability14,37. Nunes37 confirmed that the forage species selected by informants and the criteria they adopted coincided with nutritional values measured by the literature, and that, as also found in the present study, younger plants were recognized as highly appreciated by animals. This appreciation is due to the greater palatability of plant organs at this stage47. This is a matter of concern for the sustainability of the Caatinga, since direct or indirect grazing has compromised the regeneration process12 since younger individuals are clearly more sensitive to damage48.Also, considering the potential of Caatinga, we suggest that investment through government actions encourage the cultivation of native species to ensure the production of forage and, consequently, guarantee the sustainability of livestock activity and the ecosystem in question. More

This study proposed a protocol to evaluate the behavior of geese reared outdoors in agroforestry systems. A data collection form (i.e., BOF) was developed and validated both in relation to its reliability and its validity. In this context, moreover, ABMs useful for a welfare assessment protocol could be defined, and changes in the behavior of geese due to daily time and environmental context could be identified.Behavioral observations, based on the capture of the major changes in an animal’s body language17, are used daily in the assessment of animal health and welfare. Body language is a type of dynamic expression of the interactions among conspecifics or between animals and their environment. Behavioral changes can happen quickly or as subtle shifts not easily detectable18. Indeed, especially in the case of direct observation in the field, it becomes difficult to identify each behavioral variation. Furthermore, the on-farm use of the BOF proposed in the present study involved focal subgroup sampling, as ten geese were simultaneously observed, which may increase the difficulties. Indirect observation by videos, which allow the review of a certain action several times and the focal-animal approach, is a useful tool to partially overcome these issues and thus improve the accuracy of observation. The validation process of the BOF adopted in this study, therefore, included the definition of both its interobserver reliability and correlation with indirect observations.In this study, the direct observations in the field were performed by both an expert (i.e., main observer) and an inexperienced trained observer. As expected, the main observer was able to detect a higher frequency of behaviors, especially the rarer ones. For example, the inexperienced observer did not report any examples of allo-grooming, squawking, wagging tail, stretching, or panting behavior. However, the two observers showed excellent interobserver reliability (ICC  > 0.75). Major agreements were found for walking, roosting, and foraging. Accordingly, several studies have shown that observers with little experience can also provide a valuable contribution in observational research19,20. Overall, these results support the reliability of the BOF even if the observer’s experience helps him or her to better grasp rarer behaviors, as these behaviors could play an important role as welfare indicators.In the last two decades, important technological developments have occurred in the livestock sector. The use of sensors, cameras, and other devices can generate objective information about individual behavior, thereby allowing its evaluation in large observation areas and for large groups of animals and resulting in the better detection of natural animal behavior. Thus, in our study, the data collected by a video recording system (Noldus XT) were used as a gold standard measure to define the criterion validity of the BOF. Our results indicated excellent agreement between direct and indirect observations, supporting the BOF criterion validity. A poor correlation was only found for 2 variables (i.e., squawking and wagging tail), which were more difficult to collect by direct observation. The use of the BOF involved the simultaneous observation of 10 animals, but the geese had a synchronized behavior and moved in groups within the grazing area. This greatly facilitated focal subgroup sampling, allowed all animals to always be under observation, and could explain the high correlation between the two observation methods. However, the comparison between the observations collected in the field by the main observer and those recorded using the computerized system confirmed the greater accuracy of the latter. The analysis of the video in continuous with the use of some tools, such as the zoom or slow-motion functions, and the focal-animal sampling provided an easier identification of some behaviors and, in general, greater accuracy. Due to its nonintrusive approach, video recording has become a common practice for behavior assessment21, but it can be expensive and time-consuming. On the other hand, direct observations made by the BOF were valid and less expensive, suggesting that it could be a feasible tool with which to evaluate the welfare principle of Appropriate behavior. As recommended for welfare assessment protocols22, the BOF ethogram included indicators of both positive and negative states; however, it would be necessary to integrate it with behavioral tests and other ABMs evaluating the human-animal relationship.As mentioned above, there is no standardized geese behavior ethogram. Thus, to verify the content validity of the BOF, its behavior variables were analyzed through a PCA. The 4 extracted PCs could represent the broad behavioral dimensions of geese. In particular, the geese’s activity reported in PC1 was characterized by locomotor, foraging, and exploratory behaviors, with opposite signs with respect to roosting. The positive correlation between explorative and grazing activities and their negative correlation with static behaviors has been widely demonstrated in chickens. Chicken genotypes characterized by low exploratory aptitude exhibited low kinetic behaviors but a high frequency of roost and rest behaviors23. Göransson et al.24 showed that 50% of the observed birds exhibited sitting behavior, whereas less than 10% performed foraging activity.PC2 included all the variables that characterized the geese’s social aspects, including both positive and negative interactions. Usually, greylag geese live in a large flock because the offspring remain with their parents for an entire year. Such groups are characterized by complex relationships based on social interactions25. The formation of a group is characterized by agonistic behaviors such as fighting, pecking, and threatening, as well as submissive behaviors such as avoiding contact, crouching, and escaping26 to establish a hierarchical order. After this phase, a tolerance status develops, and birds maintain their social interactions through the use of body postures and vocalizations. Accordingly, the variables reported in PC2 were related not only to aggressive behaviors but also to geese’s vocalization and posture, which probably helped to maintain flock stability. Therefore, a higher PC2 score could indicate the need to establish and maintain a hierarchical order within the group, resulting in high social interactions.PC3 reported comfort and body care behaviors. The opportunity to spend a lot of time on body care, which should also include access to water for bathing, is of paramount importance with regard to fulfilling the biological requirements of geese27. Thus, a higher loading of this PC means that animals showed a good degree of both welfare and adaptability. In our study, a high frequency of self-cleaning and wing flapping behaviors was recorded, and the geese often took advantage of the water tub. In contrast, a very low frequency of aggression behaviors was observed, suggesting that the groups of geese were quite stable and that the animals felt safe in the environment in which they were rear. These findings confirm that agroforestry has a favorable impact on bird welfare by allowing the display of the full range of behavior, improving the animals’ comfort28.PC4 was mainly represented by the neck forward behavior. This position only occasionally represents an attack behavior and is not utilized during the establishment of hierarchical order but when it is necessary to maintain and reinforce the order inside the group. Furthermore, a goose that assumes this posture often does so while continuing another activity29. The neck forward behavior was positively associated with the stretching behavior. Stretching is usually categorized as a comfort behavior for broilers30, but it could also be used when the animal needs to relax stress-related tension in their muscles31,32 or as an adaptive strategy for dealing with unknown contexts33. Neck forward and stretching were eventually considered social avoidance behaviors, although they could be ambivalent and thus require further study, case-by-case assessment, and perhaps a better description in the ethogram.Finally, some interesting results emerged regarding the comparison of geese’s behavior during the morning and afternoon and between the two different agroforestry systems. In particular, geese showed a higher frequency of active behaviors such as walking, foraging, drinking, neck forward, and feeding during the morning compared to the afternoon. All of these behaviors suggest that geese concentrate their grazing and exploration activities during the morning. When and where to move is crucial for the food search and to avoid both predators and adverse climate conditions34. Cartoni Mancinelli et al.35 included exploratory attitude, walking, and eating grass activities in a multifactorial score as important parameters to consider to evaluate the adaptability of different organically reared chicken genotypes. Thus, exploratory and kinetic behaviors are fundamental, especially in animals reared outdoors. Moreover, the positive correlation between walking and grazing behaviors is widely known36,37. In contrast, during the afternoon, geese showed higher frequencies of static behaviors such as resting, roosting, and self-grooming, suggesting that geese are more dedicated to comfort and body care activities during this time. These trials were performed in the hottest season; thus, the geese’s behavioral differences during the day could also depend on the fact that animals preferred to carry out active behaviors during the cooler hours (morning), while in the hottest hours (afternoon), they engaged in static activities. Active behaviors cause an increase in metabolism and body temperature38, whereas static behavior, such as roosting, is considered adaptative behavior to promote heat dissipation31,39.This could also explain why higher frequencies of walking and foraging and lower frequencies of static behaviors were found in the orchard system than in the vineyard system. Studies carried out on chickens have reported that, among different pasture enrichments, the presence of trees promotes walking animal activity compared with crop inclusion40,41. The cover provided by trees made the animals feel protected from predators and provided shade during the hottest part of the day40, thereby stimulating the animals to explore all the available space in the pen. Accordingly, geese reared in the apple orchard ingested more grass than those reared in a vineyard36. However, there were no differences between the two systems for social behaviors. Moreover, the highest frequency of roosting and self-cleanliness behaviors was recorded in the vineyard, suggesting that this space offered a comfortable environment and that both systems seem respectful of the biological needs and welfare of the geese.The behavioral assessment protocol proposed in this study involving the BOF ethogram was feasible, low-cost, fast, and responsive both over time and between housing systems. It could thus be used for the assessment of Appropriate behavior in a welfare assessment protocol for geese reared in outdoor or free-range systems, although it lacks indicators of the human-animal relationship, such as avoidance distance or handling tests; such a scoring system should be developed. Regarding the specific behaviors in the two agroforestry systems, it should also be noted that they are difficult to generalize, as the characteristics of the plants, the environment, and management could have influenced these traits. Specifically, the behaviors could have been affected by the temperatures; therefore, further trials at different altitudes, seasons (i.e., autumn and winter), and climate are necessary for external validation. More

Editorial Committee of Vegetation Map China. Vegetation Map of China (1:1000 000) (Geology Press, 2007).Fu, B. et al. Current condition and protection strategies of Qinghai–Tibet Plateau ecological security barrier. Bull. Chin. Acad. Sci. 36, 1298–1306 (2021).
Google Scholar 
Zhang, Y. et al. Spatial and temporal variability in the net primary production (NPP) of alpine grassland on Tibetan Plateau from 1982 to 2009. Acta Geogr. Sin. 68, 1197–1211 (2013).
Google Scholar 
Bao, C. & Liu, R. Spatiotemporal evolution of the urban system in the Tibetan Plateau. J. Geoinf. Sci. 21, 1330–1340 (2019).
Google Scholar 
Miehe, G. et al. The Kobresia pygmaea ecosystem of the Tibetan highlands — origin, functioning and degradation of the world’s largest pastoral alpine ecosystem: Kobresia pastures of Tibet. Sci. Total Environ. 648, 754–771 (2019). This work describes features of K. pygmaea grassland and reveals that overstocking has caused pasture degradation and soil deterioration.Article 

Google Scholar 
Liu, Y. et al. Grassland dynamics in responses to climate variation and human activities in China from 2000 to 2013. Sci. Total Environ. 690, 27–39 (2019).Article 

Google Scholar 
Cao, J. et al. Grassland degradation on the Qinghai–Tibetan Plateau: reevaluation of causative factors. Rangel. Ecol. Manag. 72, 988–995 (2019).Article 

Google Scholar 
Zhao, X. Restoration and Sustainable Management of Degradaded Grassland in the Three Rivers Headwater (Science Press, 2011).Gao, Q. Exploration and Study on Eoclogical Revelization Fuatures in Qiangtang Plateau (China Agriculture Press, 2015).Gu, X. et al. Soil extractable organic C and N contents, methanotrophic activity under warming and degradation in a Tibetan alpine meadow. Agric. Ecosyst. Environ. 278, 6–14 (2019).Article 

Google Scholar 
Li, Y. et al. Changes of soil microbial community under different degraded gradients of alpine meadow. Agric. Ecosyst. Environ. 222, 213–222 (2016).Article 

Google Scholar 
Wang, W., Wang, Q. & Wang, H. The effect of land management on plant community composition, species diversity, and productivity of alpine Kobersia steppe meadow. Ecol. Res. 21, 181–187 (2005).Article 

Google Scholar 
Xu, H., Wang, X. & Zhang, X. Alpine grasslands response to climatic factors and anthropogenic activities on the Tibetan Plateau from 2000 to 2012. Ecol. Eng. 92, 251–259 (2016).Article 

Google Scholar 
Yu, L., Tang, L., Wei, D., Mei, M. & Zhou, H. Characteristics and causes of changes of alpine grassland productivity in the source region of Yellow River. Int. Conf. Geoinformatics https://doi.org/10.1109/GEOINFORMATICS.2010.5567879 (2010).Article 

Google Scholar 
Yang, Y. et al. Responses of the functional structure of soil microbial community to livestock grazing in the Tibetan alpine grassland. Glob. Chang. Biol. 19, 637–648 (2013). This work shows that soil microbial community functional structure is very sensitive to livestock grazing.Article 

Google Scholar 
Gao, Y. Z. et al. Belowground net primary productivity and biomass allocation of a grassland in Inner Mongolia is affected by grazing intensity. Plant Soil 307, 41–50 (2008).Article 

Google Scholar 
Dlamini, A. P. et al. Controlling factors of sheet erosion under degraded grasslands in the sloping lands of KwaZulu-Natal, South Africa. Agric. Water Manag. 98, 1711–1718 (2011).Article 

Google Scholar 
Niemandt, C. & Greve, M. Fragmentation metric proxies provide insights into historical biodiversity loss in critically endangered grassland. Agric. Ecosyst. Environ. 235, 172–181 (2016).Article 

Google Scholar 
Kang, S. C. et al. Review of climate and cryospheric change in the Tibetan Plateau. Environ. Res. Lett. 5, 15101–15101 (2010).Article 

Google Scholar 
Shen, H. et al. Effects of simulated N deposition on photosynthesis and productivity of key plants from different functional groups of alpine meadow on Qinghai–Tibetan Plateau. Environ. Pollut. 251, 731–737 (2019).Article 

Google Scholar 
Yu, G. R. et al. Stabilization of atmospheric nitrogen deposition in China over the past decade. Nat. Geosci. 12, 424 (2019).Article 

Google Scholar 
Lu, C. & Tian, H. Spatial and temporal patterns of nitrogen deposition in China: synthesis of observational data. J. Geophys. Res. 112, D22S05 (2007).
Google Scholar 
National Bureau of Statistics of China. China Statistics Yearbook (China Statistics Press, 2020).Mo, X. Sustainable livestock carring capacity and overgrazing rate of grassland over Qinghai–Tibet plateau since 1980. Natl Tibetan Plateau Data Center https://doi.org/10.11888/Socioeco.tpdc.270347 (2020).Article 

Google Scholar 
Tian, Y. Y., Jiang, G. H., Zhou, D. Y. & Li, G. Y. Systematically addressing the heterogeneity in the response of ecosystem services to agricultural modernization, industrialization and urbanization in the Qinghai–Tibetan Plateau from 2000 to 2018. J. Clean. Prod. https://doi.org/10.1016/j.jclepro.2020.125323 (2021).Article 

Google Scholar 
Yao, Y. et al. Spatiotemporal pattern of gross primary productivity and its covariation with climate in China over the last thirty years. Glob. Chang. Biol. 24, 184–196 (2018).Article 

Google Scholar 
Li, L. et al. Current challenges in distinguishing climatic and anthropogenic contributions to alpine grassland variation on the Tibetan Plateau. Ecol. Evol. 8, 5949–5963 (2018). This work finds that large inconsistencies exist in distinguishing the respective contribution of climatic and anthropogenic forces to grassland dynamics.Article 

Google Scholar 
Zhong, L., Ma, Y. M., Xue, Y. K. & Piao, S. L. Climate change trends and impacts on vegetation greening over the Tibetan Plateau. J. Geophys. Res. Atmos. 124, 7540–7552 (2019). This work demonstrates that the general increasing trends in vegetation density and greening of the QTP are mainly caused by climate factors, using satellite-derived climate and vegetation data from 1999 to 2014.Article 

Google Scholar 
Yang, K. & He, J. China meteorological forcing dataset (1979–2018). Natl Tibetan Plateau Data Center https://doi.org/10.11888/AtmosphericPhysics.tpe.249369.file (2019).Article 

Google Scholar 
Xiong, Q. et al. Monitoring the impact of climate change and human activities on grassland vegetation dynamics in the northeastern Qinghai–Tibet Plateau of China during 2000–2015. J. Arid. Land 11, 637–651 (2019).Article 

Google Scholar 
Pan, T., Zou, X. T., Liu, Y. J., Wu, S. H. & He, G. M. Contributions of climatic and non-climatic drivers to grassland variations on the Tibetan Plateau. Ecol. Eng. 108, 307–317 (2017).Article 

Google Scholar 
Hou, X. 1:1 Million vegetation map of China (National Tibetan Plateau Data Center, 2019).Peng, S. S. et al. Recent change of vegetation growth trend in China. Environ. Res. Lett. 6, 044027 (2011).Article 

Google Scholar 
Yuan, W. et al. Increased atmospheric vapor pressure deficit reduces global vegetation growth. Sci. Adv. 5, eaax1396 (2019).Article 

Google Scholar 
Zhu, Z. C. et al. Greening of the earth and its drivers. Nat. Clim. Chang. 6, 791 (2016).Article 

Google Scholar 
Vermote, E. et al. NOAA climate data record (CDR) of normalized difference vegetation index (NDVI), version 4. NOAA Natl Cent. Environ. Inf. https://doi.org/10.7289/V5PZ56R6 (2014).Article 

Google Scholar 
Chen, H. et al. Attribution analyses of changes in alpine grasslands on the Qinghai–Tibetan Plateau. Chin. Sci. Bull. 65, 2406–2418 (2020). This work demonstrates that human activities play an increasingly important role in the restoration of degraded grasslands.Article 

Google Scholar 
Shen, M. et al. Evaporative cooling over the Tibetan Plateau induced by vegetation growth. Proc. Natl Acad. Sci. USA 112, 9299–9304 (2015).Article 

Google Scholar 
Cai, D. et al. Vegetation dynamics on the Tibetan Plateau (1982–2006): an attribution by ecohydrological diagnostics. J. Clim. 28, 4576–4584 (2015).Article 

Google Scholar 
Zhou, W. et al. Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 83, 303–313 (2017).Article 

Google Scholar 
Liu, Z. et al. Patterns of plant species diversity along an altitudinal gradient and its effect on above-ground biomass in alpine meadows in Qinghai–Tibet Plateau. Biodivers. Sci. 23, 451–462 (2015).Article 

Google Scholar 
Harris, R. B. Rangeland degradation on the Qinghai–Tibetan Plateau: a review of the evidence of its magnitude and causes. J. Arid. Environ. 74, 1–12 (2010).Article 

Google Scholar 
Lu, S. et al. Basic characteristics of Stipa sareptana var. krylovii communities in China. Chin. J. Plant. Ecol. 44, 1087–1094 (2020).Article 

Google Scholar 
Qiao, X. et al. Distribution, community characteristics and classification of Stipa tianschanica var. klemenzii steppe in China. Chin. J. Plant. Ecol. 41, 231–237 (2017).Article 

Google Scholar 
Qiao, X., Guo, K., Zhao, L., Yang, Y. & Zhao, H. Distribution, community characteristics and classification of Stipa basiplumosa steppe on Tibetan Plateau. Geogr. Res. 36, 2432–2440 (2017).
Google Scholar 
Qiao, X., Guo, K., Zhao, L., Wang, Z. & Liu, C. Community characteristics of Stipa bungeana alliance in China. Chin. J. Plant Ecol. 44, 986–994 (2020).Article 

Google Scholar 
Zhu, Y., Qiao, X., Guo, K., Xu, R. & Zhao, L. Distribution, community characteristics and classification of Stipa tianschanica var. gobica steppe in China. Chin. J. Plant Ecol. 42, 785–792 (2018).Article 

Google Scholar 
Li, X. R., Jia, X. H. & Dong, G. R. Influence of desertification on vegetation pattern variations in the cold semi-arid grasslands of Qinghai–Tibet plateau, north-west China. J. Arid. Environ. 64, 505–522 (2006).Article 

Google Scholar 
Tang, L. et al. Changes in vegetation composition and plant diversity with rangeland degradation in the alpine region of Qinghai–Tibet Plateau. Rangel. J. 37, 107–115 (2015).Article 

Google Scholar 
Zhou, X. Chinese Kobresia Meadow (Science Press, 2001).Wang, B. Z. et al. Potential distribution patterns of Stipa bungeana in China and the major factors influencing distribution. Acta Prataculturae Sinica 28, 3–13 (2019).
Google Scholar 
Sun, H., Li, W., Zhang, M. & Han, Y. A comprehensive scientific expedition to the Qinghai–Tibet Plateau. Resour. Sci. 8, 22–30 (1986).
Google Scholar 
Zhu, F. X. et al. Spatiotemporal variations of annual shallow soil temperature on the Tibetan Plateau during 1983–2013. Clim. Dyn. 51, 2209–2227 (2018).Article 

Google Scholar 
Chen, L. T. et al. Changes of carbon stocks in alpine grassland soils from 2002 to 2011 on the Tibetan Plateau and their climatic causes. Geoderma 288, 166–174 (2017).Article 

Google Scholar 
Ding, J. et al. Decadal soil carbon accumulation across Tibetan permafrost regions. Nat. Geosci. 10, 420–424 (2017).Article 

Google Scholar 
Tian, L. M. et al. Variations in soil nutrient availability across Tibetan grassland from the 1980s to 2010s. Geoderma 338, 197–205 (2019).Article 

Google Scholar 
Pepin, N. et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 5, 424–430 (2015).Article 

Google Scholar 
Chen, H. et al. The impacts of climate change and human activities on biogeochemical cycles on the Qinghai–Tibetan Plateau. Glob. Chang. Biol. 19, 2940–2955 (2013). This work suggests that warming enhanced NPP and soil respiration but many uncertainties remain.Article 

Google Scholar 
Shen, M. G. et al. Plant phenological responses to climate change on the Tibetan Plateau: research status and challenges. Natl Sci. Rev. 2, 454–467 (2015).Article 

Google Scholar 
Liu, X. D., Yin, Z. Y., Shao, X. M. & Qin, N. S. Temporal trends and variability of daily maximum and minimum, extreme temperature events, and growing season length over the eastern and central Tibetan Plateau during 1961–2003. J. Geophys. Res. Atmos. https://doi.org/10.1029/2005jd006915 (2006).Article 

Google Scholar 
Yang, K. et al. Recent climate changes over the Tibetan Plateau and their impacts on energy and water cycle: a review. Glob. Planet. Change 112, 79–91 (2014). This work reviews the main spatio-temporal characteristics of climate change on the QTP.Article 

Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Experimental warming, not grazing, decreases rangeland quality on the Tibetan Plateau. Ecol. Appl. 17, 541–557 (2007).Article 

Google Scholar 
Li, C. et al. Productivity and quality of alpine grassland vary with soil water availability under experimental warming. Front. Plant. Sci. 9, 1790 (2018).Article 

Google Scholar 
Peng, A. H. et al. Plant community responses to warming modified by soil moisture in the Tibetan Plateau. Arct. Antarct. Alp. Res. 52, 60–69 (2020).Article 

Google Scholar 
Li, F. et al. Leaf area rather than photosynthetic rate determines the response of ecosystem productivity to experimental warming in an alpine steppe. J. Geophys. Res. Biogeosci. 124, 2277–2287 (2019).Article 

Google Scholar 
Zong, N. et al. Responses of ecosystem CO2 fluxes to short-term experimental warming and nitrogen enrichment in an alpine meadow, northern Tibet Plateau. Sci. World J. 2013, 415318 (2013).Article 

Google Scholar 
Chen, Q., Niu, B., Hu, Y., Luo, T. & Zhang, G. Warming and increased precipitation indirectly affect the composition and turnover of labile-fraction soil organic matter by directly affecting vegetation and microorganisms. Sci. Total Environ. 714, 136787 (2020).Article 

Google Scholar 
Wang, X. X. et al. Effects of short-term and long-term warming on soil nutrients, microbial biomass and enzyme activities in an alpine meadow on the Qinghai–Tibet Plateau of China. Soil. Biol. Biochem. 76, 140–142 (2014).Article 

Google Scholar 
Jiang, L. L. et al. Plant organic N uptake maintains species dominance under long-term warming. Plant Soil 433, 243–255 (2018).Article 

Google Scholar 
Li, N. et al. Short-term effects of temperature enhancement on community structure and biomass of alpine meadow in the Qinghai–Tibet Plateau. Acta Ecol. Sin. 31, 895–905 (2011).
Google Scholar 
Jiang, Y., Fan, M. & Zhang, Y. Effect of short-term warming on plant community features of alpine meadow in northern Tibet. Chin. J. Ecol. 36, 616–622 (2017).
Google Scholar 
Wang, S. et al. Effects of warming and grazing on soil N availability, species composition, and ANPP in an alpine meadow. Ecology 93, 2365–2376 (2012). This work shows the effects of asymmetric warming and moderate grazing on plant composition, diversity, productivity and their relationships.Article 

Google Scholar 
Liu, P. et al. Ambient climate determines the directional trend of community stability under warming and grazing. Glob. Change Biol. 27, 5198–5210 (2021). This work finds that the negative effect of warming on plant diversity disappears with experimental duration, and ambient climate modulates the effects of warming and grazing on productivity stability.Article 

Google Scholar 
Zhang, B. et al. Responses of soil microbial communities to experimental warming in alpine grasslands on the Qinghai–Tibet Plateau. PLoS ONE 9, e103859 (2014).Article 

Google Scholar 
Chen, X. et al. Effects of warming and nitrogen fertilization on GHG flux in the permafrost region of an alpine meadow. Atmos. Environ. 157, 111–124 (2017).Article 

Google Scholar 
Zhang, Y. et al. Effects of grazing and climate warming on plant diversity, productivity and living state in the alpine rangelands and cultivated grasslands of the Qinghai–Tibetan Plateau. Rangel. J. 37, 57–65 (2015).Article 

Google Scholar 
Quan, Q. et al. High-level rather than low-level warming destabilizes plant community biomass production. J. Ecol. 109, 1607–1617 (2021).Article 

Google Scholar 
Wang, X. et al. Response of greenhouse gases emission fluxes to long-term warming in alpine meadow of northern Tibet. Chin. J. Agrometeorol. 39, 152–161 (2018).
Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Experimental warming causes large and rapid species loss, dampened by simulated grazing, on the Tibetan Plateau. Ecol. Lett. 7, 1170–1179 (2004).Article 

Google Scholar 
Zhang, C. H. et al. Recovery of plant species diversity during long-term experimental warming of a species-rich alpine meadow community on the Qinghai–Tibet Plateau. Biol. Conserv. 213, 218–224 (2017).Article 

Google Scholar 
Li, X. et al. Responses of biotic interactions of dominant and subordinate species to decadal warming and simulated rotational grazing in Tibetan alpine meadow. Sci. China Life Sci. 61, 849–859 (2018).Article 

Google Scholar 
Klein, J. A., Harte, J. & Zhao, X. Q. Dynamic and complex microclimate responses to warming and grazing manipulations. Glob. Chang. Biol. 11, 1440–1451 (2005).Article 

Google Scholar 
Chen, J. et al. Warming effects on ecosystem carbon fluxes are modulated by plant functional types. Ecosystems 20, 515–526 (2017).Article 

Google Scholar 
Zhang, Y. Q. & Welker, J. M. Tibetan alpine tundra responses to simulated changes in climate: aboveground biomass and community responses. Arct. Alp. Res. 28, 203–209 (1996).Article 

Google Scholar 
Liu, H. et al. Shifting plant species composition in response to climate change stabilizes grassland primary production. Proc. Natl Acad. Sci. USA 115, 4051–4056 (2018). This work demonstrates that shifting plant species composition in response to climate change may have stabilized primary production in this high-elevation ecosystem, but also causes a shift from above-ground to below-ground productivity.Article 

Google Scholar 
Ganjurjav, H. et al. Differential response of alpine steppe and alpine meadow to climate warming in the central Qinghai–Tibetan Plateau. Agric. For. Meteorol. 223, 233–240 (2016).Article 

Google Scholar 
Jiang, F., Wei, X., Kang, B. & Shao, X. Effects of warming on alpine meadow diversity and primary productivity. Acta Agrestia Sin. 27, 298–305 (2019).
Google Scholar 
Zong, N. et al. Responses of plant community structure and species composition to warming and N addition in an alpine meadow, northern Tibetan Plateau, China. Chin. J. Appl. Ecol. 27, 3739–3748 (2016).
Google Scholar 
Peng, F. et al. Warming-induced shift towards forbs and grasses and its relation to the carbon sequestration in an alpine meadow. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/aa6508 (2017).Article 

Google Scholar 
Dorji, T. et al. Grazing and spring snow counteract the effects of warming on an alpine plant community in Tibet through effects on the dominant species. Agric. For. Meteorol. 263, 188–197 (2018).Article 

Google Scholar 
Xue, X., Peng, F., You, Q., Xu, M. & Dong, S. Belowground carbon responses to experimental warming regulated by soil moisture change in an alpine ecosystem of the Qinghai–Tibet Plateau. Ecol. Evol. 5, 4063–4078 (2015).Article 

Google Scholar 
Jing, X. et al. No temperature acclimation of soil extracellular enzymes to experimental warming in an alpine grassland ecosystem on the Tibetan Plateau. Biogeochemistry 117, 39–54 (2014).Article 

Google Scholar 
Yu, C. Q., Shen, Z. X., Zhang, X. Z., Sun, W. & Fu, G. Response of soil C and N, dissolved organic C and N, and inorganic N to short-term experimental warming in an alpine meadow on the Tibetan Plateau. Sci. World J. 2014, 152576 (2014).
Google Scholar 
Zhang, Y. et al. Simulated warming enhances the responses of microbial N transformations to reactive N input in a Tibetan alpine meadow. Environ. Int. 141, 105795 (2020).Article 

Google Scholar 
Jia, J. et al. Climate warming alters subsoil but not topsoil carbon dynamics in alpine grassland. Glob. Chang. Biol. 25, 4383–4393 (2019).Article 

Google Scholar 
Ding, X. L. et al. Warming increases microbial residue contribution to soil organic carbon in an alpine meadow. Soil. Biol. Biochem. 135, 13–19 (2019).Article 

Google Scholar 
Guan, S. et al. Climate warming impacts on soil organic carbon fractions and aggregate stability in a Tibetan alpine meadow. Soil. Biol. Biochem. 116, 224–236 (2018).Article 

Google Scholar 
Rui, Y. C. et al. Warming and grazing affect soil labile carbon and nitrogen pools differently in an alpine meadow of the Qinghai–Tibet Plateau in China. J. Soils Sediment. 11, 903–914 (2011).Article 

Google Scholar 
Heng, T., Wu, J., Xie, S. & Wu, M. The responses of soil C and N, microbial biomass C or N under alpine meadow of Qinghai–Tibet Plateau to the change of temperature and precipitation. Chin. Agric. Sci. Bull. 27, 425–430 (2011).
Google Scholar 
Li, N., Wang, G., Yang, Y., Gao, Y. & Liu, G. Plant production, and carbon and nitrogen source pools, are strongly intensified by experimental warming in alpine ecosystems in the Qinghai–Tibet Plateau. Soil. Biol. Biochem. 43, 942–953 (2011).Article 

Google Scholar 
Zhao, J. X. et al. Increased precipitation offsets the negative effect of warming on plant biomass and ecosystem respiration in a Tibetan alpine steppe. Agric. For. Meteorol. https://doi.org/10.1016/j.agrformet.2019.107761 (2019). This work shows that increased precipitation offsets the negative effect of warming on plant biomass and ecosystem respiration in a Tibetan alpine steppe.Article 

Google Scholar 
Wu, H. et al. Effects of increased precipitation combined with nitrogen addition and increased temperature on methane fluxes in alpine meadows of the Tibetan Plateau. Sci. Total Environ. 705, 135818 (2020).Article 

Google Scholar 
Shi, F. S., Chen, H., Chen, H. F., Wu, Y. & Wu, N. The combined effects of warming and drying suppress CO2 and N2O emission rates in an alpine meadow of the eastern Tibetan Plateau. Ecol. Res. 27, 725–733 (2012).Article 

Google Scholar 
Fu, G., Zhang, H. R. & Sun, W. Response of plant production to growing/non-growing season asymmetric warming in an alpine meadow of the northern Tibetan Plateau. Sci. Total Environ. 650, 2666–2673 (2019).Article 

Google Scholar 
Xiong, Q. L. et al. Warming and nitrogen deposition are interactive in shaping surface soil microbial communities near the alpine timberline zone on the eastern Qinghai–Tibet Plateau, southwestern China. Appl. Soil. Ecol. 101, 72–83 (2016).Article 

Google Scholar 
Wang, C. et al. Responses of plant leaf traits to simulated rainfall changes in alpine region. Acta Ecol. Sin. 41, 1–13 (2021).Article 

Google Scholar 
Zhang, K. et al. Effects of short-term warming and altered precipitation on soil microbial communities in alpine grassland of the Tibetan Plateau. Front. Microbiol. 7, 1032 (2016).
Google Scholar 
Evans, R. D. & Ehleringer, J. R. Water and nitrogen dynamics in an arid woodland. Oecologia 99, 233–242 (1994).Article 

Google Scholar 
Swap, R. J., Aranibar, J. N., Dowty, P. R., Gilhooly, W. P. III & Macko, S. A. Natural abundance of 13C and 15N in C3 and C4 vegetation of southern Africa: patterns and implications. Glob. Chang. Biol. 10, 350–358 (2004).Article 

Google Scholar 
Jia, Y. et al. Spatial and decadal variations in inorganic nitrogen wet deposition in China induced by human activity. Sci. Rep. 4, 3763–3763 (2014).Article 

Google Scholar 
Liu, Y. W., Xu, R., Wang, Y. S., Pan, Y. P. & Piao, S. L. Wet deposition of atmospheric inorganic nitrogen at five remote sites in the Tibetan Plateau. Atmos. Chem. Phys. 15, 11683–11700 (2015).Article 

Google Scholar 
Wang, W. et al. Atmospheric nitrogen deposition to a southeast Tibetan forest ecosystem. Atmosphere https://doi.org/10.3390/atmos11121331 (2020).Article 

Google Scholar 
Zou, X. et al. Ice-core based assessment of nitrogen deposition in the central Tibetan Plateau over the last millennium. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2021.152692 (2022).Article 

Google Scholar 
Aerts, R., Wallen, B. & Malmer, N. Growth-limiting nutrients in sphagnum-dominated bogs subject to low and high amospheric nitrogen supply. J. Ecol. 80, 131–140 (1992).Article 

Google Scholar 
Bai, Y. F. et al. Tradeoffs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning: evidence from Inner Mongolia grasslands. Glob. Chang. Biol. 16, 358–372 (2010).Article 

Google Scholar 
Du, Y. Population statistics of Qinghai–Tibet Plateau (1952–2016) (National Tibetan Plateau Data Center, 2019).Zhang, Y. J., Zhang, X. Q., Wang, X. Y., Liu, N. & Kan, H. M. Establishing the carrying capacity of the grasslands of China: a review. Rangel. J. 36, 1–9 (2014).Article 

Google Scholar 
Bardgett, R. D. et al. Combatting global grassland degradation. Nat. Rev. Earth Environ. 2, 720–735 (2021). This work shows that socio-ecological solutions are needed to combat degradation and promote restoration.Article 

Google Scholar 
Liu, M. et al. Effects of rotational and continuous overgrazing on newly assimilated C allocation. Biol. Fertil. Soils 57, 193–202 (2021).Article 

Google Scholar 
Yang, X. X. et al. Different responses of soil element contents and their stoichiometry (C:N:P) to yak grazing and Tibetan sheep grazing in an alpine grassland on the eastern Qinghai–Tibetan Plateau. Agric. Ecosyst. Environ. https://doi.org/10.1016/j.agee.2019.106628 (2019).Article 

Google Scholar 
Lin, B., Zhao, X. R., Zheng, Y., Qi, S. & Liu, X. Z. Effect of grazing intensity on protozoan community, microbial biomass, and enzyme activity in an alpine meadow on the Tibetan Plateau. J. Soils Sediment. 17, 2752–2762 (2017).Article 

Google Scholar 
Ma, W. M., Ding, K. Y. & Li, Z. W. Comparison of soil carbon and nitrogen stocks at grazing-excluded and yak grazed alpine meadow sites in Qinghai–Tibetan Plateau, China. Ecol. Eng. 87, 203–211 (2016).Article 

Google Scholar 
Li, W. et al. Effects of grazing regime on vegetation structure, productivity, soil quality, carbon and nitrogen storage of alpine meadow on the Qinghai–Tibetan Plateau. Ecol. Eng. 98, 123–133 (2017).Article 

Google Scholar 
Sun, J. et al. Effects of grazing regimes on plant traits and soil nutrients in an alpine steppe, northern Tibetan Plateau. PLoS ONE 9, e108821 (2014).Article 

Google Scholar 
Niu, K. C., He, J. S. & Lechowicz, M. J. Grazing-induced shifts in community functional composition and soil nutrient availability in Tibetan alpine meadows. J. Appl. Ecol. 53, 1554–1564 (2016).Article 

Google Scholar 
Luan, J. W. et al. Different grazing removal exclosures effects on soil C stocks among alpine ecosystems in east Qinghai–Tibet Plateau. Ecol. Eng. 64, 262–268 (2014).Article 

Google Scholar 
Wei, D. et al. Responses of CO2, CH4 and N2O fluxes to livestock exclosure in an alpine steppe on the Tibetan Plateau, China. Plant Soil 359, 45–55 (2012).Article 

Google Scholar 
Shen, H. et al. Grazing enhances plant photosynthetic capacity by altering soil nitrogen in alpine grasslands on the Qinghai–Tibetan Plateau. Agric. Ecosyst. Environ. 280, 161–168 (2019).Article 

Google Scholar 
Jiang, L. et al. Grazing modifies inorganic and organic nitrogen uptake by coexisting plant species in alpine grassland. Biol. Fertil. Soils 52, 211–221 (2016).Article 

Google Scholar 
Sun, Y., Schleuss, P. M., Pausch, J., Xu, X. L. & Kuzyakov, Y. Nitrogen pools and cycles in Tibetan Kobresia pastures depending on grazing. Biol. Fertil. Soils 54, 569–581 (2018).Article 

Google Scholar 
Chen, B. et al. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai–Tibet Plateau. Agric. For. Meteorol. 189-190, 11–18 (2014).Article 

Google Scholar 
Wang, Z. Q. et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai–Tibet Plateau, in China. Ecol. Inform. 33, 32–44 (2016).Article 

Google Scholar 
Huang, K. et al. The influences of climate change and human activities on vegetation dynamics in the Qinghai–Tibet Plateau. Remote Sens. 8, 876 (2016).Article 

Google Scholar 
Li, L. et al. Increasing sensitivity of alpine grasslands to climate variability along an elevational gradient on the Qinghai–Tibet Plateau. Sci. Total Environ. 678, 21–29 (2019).Article 

Google Scholar 
Wang, Z. et al. Vegetation expansion on the Tibetan Plateau and its relationship with climate change. Remote. Sens. https://doi.org/10.3390/rs12244150 (2020).Article 

Google Scholar 
Wu, J. et al. Disentangling climatic and anthropogenic contributions to nonlinear dynamics of alpine grassland productivity on the Qinghai–Tibetan Plateau. J. Environ. Manag. 281, 111875 (2021).Article 

Google Scholar 
Fu, G., Shen, Z. X. & Zhang, X. Z. Increased precipitation has stronger effects on plant production of an alpine meadow than does experimental warming in the northern Tibetan Plateau. Agric. For. Meteorol. 249, 11–21 (2018).Article 

Google Scholar 
Hu, Y. et al. Effect of increasing precipitation and warming on microbial community in Tibetan alpine steppe. Environ. Res. 189, 109917 (2020).Article 

Google Scholar 
Ma, Z. et al. Climate warming reduces the temporal stability of plant community biomass production. Nat. Commun. 8, 15378 (2017).Article 

Google Scholar 
Bai, W., Xi, J. & Wang, G. Effects of short-term warming and nitrogen addition on CO2 emission during growing season in an alpine swamp meadow ecosystem of Qinghai–Tibetan Plateau. Chin. J. Ecol. 38, 927–936 (2019).
Google Scholar 
Bai, W., Wang, G. X., Xi, J. Y., Liu, Y. W. & Yin, P. S. Short-term responses of ecosystem respiration to warming and nitrogen addition in an alpine swamp meadow. Eur. J. Soil Biol. 92, 16–23 (2019).Article 

Google Scholar 
Ge, Y. et al. Effects of warming and nitrogen addition on plant community structure and species diversity of alpine meadow in northern Tibet. Ecol. Environ. Sci. 28, 2185–2191 (2019).
Google Scholar 
Zong, N. et al. Effects of warming and nitrogen addition on community production and biomass allocation in an alpine meadow. Chin. J. Appl. Ecol. 29, 59–67 (2018).
Google Scholar 
Zhu, X. X. et al. Effects of warming, grazing/cutting and nitrogen fertilization on greenhouse gas fluxes during growing seasons in an alpine meadow on the Tibetan Plateau. Agric. For. Meteorol. 214, 506–514 (2015).Article 

Google Scholar 
Fu, G. et al. Clipping alters the response of biomass production to experimental warming: a case study in an alpine meadow on the Tibetan Plateau, China. J. Mt. Sci. 12, 935–942 (2015).Article 

Google Scholar 
Chen, S. P. et al. Plant diversity enhances productivity and soil carbon storage. Proc. Natl Acad. Sci. USA. 115, 4027–4032 (2018).Article 

Google Scholar 
Wu, J. S. et al. Effects of livestock exclusion and climate change on aboveground biomass accumulation in alpine pastures across the northern Tibetan Plateau. Chin. Sci. Bull. 59, 4332–4340 (2014).Article 

Google Scholar 
Sun, J., Cheng, G. W., Li, W. P., Sha, Y. K. & Yang, Y. C. On the variation of NDVI with the principal climatic elements in the Tibetan Plateau. Remote Sens. 5, 1894–1911 (2013).Article 

Google Scholar 
Sun, J. et al. Reconsidering the efficiency of grazing exclusion using fences on the Tibetan Plateau. Sci. Bull. 65, 1405–1414 (2020). This work finds that fencing enclosures lead to some negative impacts, such as hindering wildlife movement.Article 

Google Scholar 
Yu, C. et al. Grazing exclusion to recover degraded alpine pastures needs scientific assessments across the northern Tibetan Plateau. Sustainability https://doi.org/10.3390/su8111162 (2016).Article 

Google Scholar 
Wu, J. & Wang, X. Effect of enclosure ages on community characters and biomas of the degraded alpine steppe at the northern Tibet. Acta Agrestia Sin. 25, 261–266 (2017).
Google Scholar 
Zhao, J. X., Luo, T. X., Li, R. C., Li, X. & Tian, L. H. Grazing effect on growing season ecosystem respiration and its temperature sensitivity in alpine grasslands along a large altitudinal gradient on the central Tibetan Plateau. Agric. For. Meteorol. 218, 114–121 (2016).Article 

Google Scholar 
Deng, L., Sweeney, S. & Shangguan, Z. P. Grassland responses to grazing disturbance: plant diversity changes with grazing intensity in a desert steppe. Grass Forage Sci. 69, 524–533 (2014).Article 

Google Scholar 
Yuan, Z., Epstein, H. & Li, G. Grazing exclusion did not affect soil properties in alpine meadows in the Tibetan permafrost region. Ecol. Eng. https://doi.org/10.1016/j.ecoleng.2019.105657 (2020).Article 

Google Scholar 
Zhang, W. et al. Effects of banning grazing and delaying grazing on species diversity and biomass of alpine meadow in northern Tibet. J. Agric. Sci. Technol. 15, 143–149 (2013).
Google Scholar 
Miao, F., Guo, Y., Miao, P., Guo, Z. & Shen, Y. The northeast edge of Qinghai–Tibet Plateau area of alpine meadow community characteristics respond to nurture. Acta Prataculture Sin. 21, 11–16 (2012).
Google Scholar 
Lu, X. et al. Short-term grazing exclusion has no impact on soil properties and nutrients of degraded alpine grassland in Tibet, China. Solid Earth 6, 1195–1205 (2015).Article 

Google Scholar 
Gao, Y. H., Zeng, X. Y., Schumann, M. & Chen, H. Effectiveness of exclosures on restoration of degraded alpine meadow in the eastern Tibetan Plateau. Arid. Land. Res. Manag. 25, 164–175 (2011).Article 

Google Scholar 
Yao, X. X. et al. Effects of long term fencing on biomass, coverage, density, biodiversity and nutritional values of vegetation community in an alpine meadow of the Qinghai–Tibet Plateau. Ecol. Eng. 130, 80–93 (2019).Article 

Google Scholar 
Chen, W., Chang, H. & Liu, R. Fractal features of soil particle size distributions and their implications for indicating enclosure management in a semiarid grassland ecosystem. Pol. J. Ecol. 68, 132–144 (2020).
Google Scholar 
Smith, D., King, R. & Allen, B. L. Impacts of exclusion fencing on target and non-target fauna: a global review. Biol. Rev. 95, 1590–1606 (2020).Article 

Google Scholar 
Zhang, Y. et al. Assessment of effectiveness of nature reserves on the Tibetan Plateau based on net primary production and the large sample comparison method. J. Geogr. Sci. 26, 27–44 (2016).Article 

Google Scholar 
Hu, J. Research on the status quo and problems of natural reserve construction in Qinghai–Tibet Plateau. Environ. Dev. 32, 204–206 (2020).
Google Scholar 
Shao, Q., Fan, J., Liu, J., Cao, W. & Liu, L. Target-based assessment on effects of first-stage ecological conservation and restoration project in three-river source region, China and policy recommendations. Bull. Chin. Acad. Sci. 32, 35–44 (2017).
Google Scholar 
Liu, F. & Zeng, Y. N. Spatial–temporal change in vegetation net primary productivity and its response to climate and human activities in Qinghai Plateau in the past 16 years. Acta Ecol. Sin. 39, 1528–1540 (2019).
Google Scholar 
Zhang, Y., Wu, X., Qi, W., Li, S. & Bai, W. Characteristics and protection effectiveness of nature reserves on the Tibetan Plateau, China. Resources. Science 37, 1455–1464 (2015).
Google Scholar 
Buckley, M. C. & Crone, E. E. Negative off-site impacts of ecological restoration: understanding and addressing the conflict. Conserv. Biol. 22, 1118–1124 (2008).Article 

Google Scholar 
Cao, S. X. & Zhang, J. Political risks arising from the impacts of large-scale afforestation on water resources of the Tibetan Plateau. Gondwana Res. 28, 898–903 (2015).Article 

Google Scholar 
Li, Y. & Li, W. Why “Balance of Forage and Livestock” system failed to reach sustainable grassland utilization. China Agric. Univ. J. Soc. Sci. Ed. 29, 124–131 (2012).
Google Scholar 
Du, S. A Study on the Satisfaction Degree of Herdsmen’s Income and Grassland Ecological Compensation Policy. Master thesis, Lanzhou Univ. (2019).Yu, H., Wang, G., Yang, Y. & Lü, Y. Concept of grassland green carrying capacity and its application framework in national park. Acta Ecol. Sin. 40, 7248–7254 (2020).
Google Scholar 
Deng, Y. & Li, C. The investigation and research about the Farmland Retirement and Environment Project in the Yangtze River headwaters area. Ecol. Econ. 2, 77–80 (2006).
Google Scholar 
Zhou, Q. et al. Analysis on the relationship between grassland area and forage-livestock balance in Qinghai–Tibet Plateau. Chin. J. Grassl. 41, 110–117 (2019).
Google Scholar 
Li, Y. et al. Awareness and reaction of herdsmen to the policy of returning grazing land to grasslands in the Changtang Plateau,Tibet. Pratacultural Sci. 30, 788–794 (2013).
Google Scholar 
Fan, J. et al. Third pole national park group construction is scientific choice for implementing strategy of major function zoning and green development in Tibet, China. Bull. Chin. Acad. Sci. 32, 932–944 (2017).
Google Scholar 
Xu, Z., Cheng, S. & Gao, L. Impacts of herders sedentarization on regional spatial heterogeneity and grassland ecosystem change in pastoral area. J. Arid. Land. Resour. Environ. 31, 8–13 (2017).
Google Scholar 
Ptackova, J. Sedentarisation of Tibetan nomads in China: implementation of the Nomadic settlement project in the Tibetan Amdo area; Qinghai and Sichuan Provinces. Pastoralism https://doi.org/10.1186/2041-7136-1-4 (2011).Article 

Google Scholar 
Weber, K. T. & Horst, S. Desertification and livestock grazing: the roles of sedentarization, mobility and rest. Pastoralism https://doi.org/10.1186/2041-7136-1-19 (2011).Article 

Google Scholar 
Zhang, J. et al. Ecological consequence of nomad settlement policy in the pasture area of Qinghai–Tibetan Plateau: from plant and soil perspectives. J. Environ. Manage. https://doi.org/10.1016/j.jenvman.2020.110114 (2020).Article 

Google Scholar 
Li, C. X., de Jong, R., Schmid, B., Wulf, H. & Schaepman, M. E. Spatial variation of human influences on grassland biomass on the Qinghai–Tibetan Plateau. Sci. Total Environ. 665, 678–689 (2019).Article 

Google Scholar 
Kuang, W. Dataset of Urban Distribution, Urban Population and Built-up Area in Tibetan Plateau (2000–2015) (National Tibetan Plateau Data Center, 2021).Tian, L. & Chen, J. Urban expansion inferenced by ecosystem production on the Qinghai–Tibet plateau. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/ac3178 (2022).Article 

Google Scholar 
Liu, Y. & Lu, C. Quantifying grass coverage trends to identify the hot plots of grassland degradation in the Tibetan Plateau during 2000–2019. Int. J. Environ. Res. Public Health https://doi.org/10.3390/ijerph18020416 (2021).Article 

Google Scholar 
Tang, L. et al. Warming counteracts grazing effects on the functional structure of the soil microbial community in a Tibetan grassland. Soil. Biol. Biochem. 134, 113–121 (2019).Article 

Google Scholar 
Zhong, L. Tourism Planning Case in Tibetan Plateau (China Tourism Press, 2018).La, M. Discussion of the coordinated development of tourism development and ecological Environment Protection in Tibetan. Soc. Sci. Res. 6, 118–120, (2013).
Google Scholar 
Zhuang, M. et al. Opportunities for household energy on the Qinghai–Tibet Plateau in line with United Nations’ Sustainable Development Goals. Renew. Sustain. Energy Rev. https://doi.org/10.1016/j.rser.2021.110982 (2021).Article 

Google Scholar 
Ruess, R. W. & Mcnaughton, S. J. Grazing and the dynamics of nutrient and energy regulated microbial processes in the serengeti grasslands. Oikos 49, 101–110 (1987).Article 

Google Scholar 
Li, M. et al. Changes in plant species richness distribution in Tibetan alpine grasslands under different precipitation scenarios. Glob. Ecol. Conserv. 21, 13 (2020).
Google Scholar 
Wang, Z. et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai–Tibet Plateau, in China. Ecol. Inform. 33, 32–44 (2016).Article 

Google Scholar 
Muñoz Sabater, J. ERA5-Land monthly averaged data from 1981 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS), https://doi.org/10.24381/cds.68d2bb30 (2019).Abatzoglou, J. T., Dobrowski, S. Z., Parks, S. A. & Hegewisch, K. C. Terraclimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 5, 170191 (2018).Article 

Google Scholar  More

Metcalf, C. J. & Pavard, S. Why evolutionary biologists should be demographers. Trends Ecol. Evol. 22, 205–212 (2007).PubMed 
Article 

Google Scholar 
French, J. C., Riris, P., Fernandez-Lopez de Pablo, J., Lozano, S. & Silva, F. A manifesto for palaeodemography in the twenty-first century. Philos. Trans. R. Soc. Lond. B Biol. Sci. 376, 20190707 (2021).PubMed 
Article 

Google Scholar 
French, J. C. Demography and the Palaeolithic archaeological record. J. Archaeol. Method Th. 23, 150–199 (2016).Article 

Google Scholar 
Henrich, J. Demography and cultural evolution: how adaptive cultural processes can produce maladaptive losses – The Tasmanian case. Am. Antiquity 69, 197–214 (2004).Article 

Google Scholar 
Powell, A., Shennan, S. & Thomas, M. G. Late Pleistocene demography and the appearance of modern human behavior. Science 324, 1298–1301 (2009).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Shennan, S. Demography and cultural innovation: a model and its implications for the emergence of modern human culture. Camb. Archaeol. J. 11, 5–16 (2001).Article 

Google Scholar 
Jorgensen, E. K. The palaeodemographic and environmental dynamics of prehistoric Arctic Norway: an overview of human-climate covariation. Quat. Int. 549, 36–51 (2020).Article 

Google Scholar 
Jorgensen, E. K. & Riede, F. Convergent catastrophes and the termination of the Arctic Norwegian Stone Age: a multi-proxy assessment of the demographic and adaptive responses of mid-Holocene collectors to biophysical forcing. Holocene 29, 1782–1800 (2019).ADS 
Article 

Google Scholar 
Riede, F. Lateglacial and Postglacial Pioneers in Northern Europe (Archaeopress, 2014).Tallavaara, M. & Seppä, H. Did the mid-Holocene environmental changes cause the boom and bust of hunter-gatherer population size in eastern Fennoscandia? Holocene 22, 215–225 (2011).ADS 
Article 

Google Scholar 
Kavanagh, P. H. et al. Hindcasting global population densities reveals forces enabling the origin of agriculture. Nat. Hum. Behav. 2, 478–484 (2018).PubMed 
Article 

Google Scholar 
Tallavaara, M., Luoto, M., Korhonen, N., Jarvinen, H. & Seppa, H. Human population dynamics in Europe over the Last Glacial Maximum. Proc. Natl Acad. Sci. USA 112, 8232–8237 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bliege Bird, R. & Codding, B. F. Promise and peril of ecological and evolutionary modelling using cross-cultural datasets. Nat. Ecol. Evol. 6, 1–3 (2021).Hamilton, M. J. & Tallavaara, M. Statistical inference, scale and noise in comparative anthropology. Nat. Ecol. Evol. 6, 122 (2022).PubMed 
Article 

Google Scholar 
Gurven, M. D. & Davison, R. J. Periodic catastrophes over human evolutionary history are necessary to explain the forager population paradox. Proc. Natl Acad. Sci. USA 116, 12758–12766 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tallavaara, M. & Jorgensen, E. K. Why are population growth rate estimates of past and present hunter-gatherers so different? Philos. T R Soc. B 376, 20190708 (2021).Blackman, F. F. Optima and limiting factors. With two diagrams in the text. Ann. Bot. Lond. 19, 281–296 (1905).Article 

Google Scholar 
Maier, A. et al. Cultural evolution and environmental change in Central Europe between 40 and 15 ka. Quat. Int. 581-582, 225–240 (2021).Article 

Google Scholar 
Zhu, D., Galbraith, E. D., Reyes-Garcia, V. & Ciais, P. Global hunter-gatherer population densities constrained by influence of seasonality on diet composition. Nat. Ecol. Evol. 5, 1536 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
Binford, L. R. Archaeology as anthropology. Am. Antiquity 28, 217–225 (1962).Article 

Google Scholar 
Lowe, J. J. et al. Synchronisation of palaeoenvironmental events in the North Atlantic region during the Last Termination: a revised protocol recommended by the INTIMATE group. Quat. Sci. Rev. 27, 6–17 (2008).ADS 
Article 

Google Scholar 
Bocquet-Appel, J. P., Demars, P. Y., Noiret, L. & Dobrowsky, D. Estimates of upper Palaeolithic meta-population size in Europe from archaeological data. J. Archaeol. Sci. 32, 1656–1668 (2005).Article 

Google Scholar 
Fort, J., Pujol, T. & Cavalli-Sforza, L. L. Palaeolithic populations and waves of advance (Human range expansions). Camb. Archaeol. J. 14, 53–61 (2004).Article 

Google Scholar 
Schmidt, I. et al. Approaching prehistoric demography: proxies, scales and scope of the Cologne Protocol in European contexts. Philos. Trans. R. Soc. Lond. B Biol. Sci. 376, 20190714 (2021).PubMed 
Article 

Google Scholar 
de Pablo, J. F. L. et al. Palaeodemographic modelling supports a population bottleneck during the Pleistocene-Holocene transition in Iberia. Nat. Commun. 10, 1872 (2019).Binford, L. R. Constructing Frames of Reference: An Analytical Method for Archaeological Theory Building Using Ethnographic and Environmental Data Sets. (Univ. California Press, 2019).Johnson, A. L. Exploring adaptive variation among hunter-gatherers with Binford’s frames of reference. J. Archaeol. Res. 22, 1–42 (2014).Article 

Google Scholar 
Graham, M. H. Confronting multicollinearity in ecological multiple regression. Ecology 84, 2809–2815 (2003).Article 

Google Scholar 
Tallavaara, M., Eronen, J. T. & Luoto, M. Productivity, biodiversity, and pathogens influence the global hunter-gatherer population density. Proc. Natl Acad. Sci. USA 115, 1232–1237 (2018).CAS 
PubMed 
Article 

Google Scholar 
Cade, B. S. & Noon, B. R. A gentle introduction to quantile regression for ecologists. Front Ecol. Environ. 1, 412–420 (2003).Article 

Google Scholar 
Cade, B. S., Terrell, J. W. & Schroeder, R. L. Estimating effects of limiting factors with regression quantiles. Ecology 80, 311–323 (1999).Article 

Google Scholar 
Burman, P., Chow, E. & Nolan, D. A cross-validatory method for dependent data. Biometrika 81, 351–358 (1994).MathSciNet 
MATH 
Article 

Google Scholar 
Burke, K. D. et al. Differing climatic mechanisms control transient and accumulated vegetation novelty in Europe and eastern North America. Philos. Trans. R. Soc. Lond. B Biol. Sci. 374, 20190218 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Currie, D. J. Energy and large-scale patterns of animal-species and plant-species richness. Am. Nat. 137, 27–49 (1991).Article 

Google Scholar 
Franklin, J. Mapping Species Distributions: Spatial Inference and Prediction (Cambridge Univ. Press, 2010).Harcourt, A. Human Biogeography (Univ. California Press, 2012).Marlowe, F. W. Hunter-gatherers and human evolution. Evol. Anthropol. 14, 54–67 (2005).Article 

Google Scholar 
Belovsky, G. E. An optimal foraging-based model of hunter-gatherer population-dynamics. J. Anthropol. Archaeol. 7, 329–372 (1988).Article 

Google Scholar 
Williams, J. W. & Jackson, S. T. Novel climates, no-analog communities, and ecological surprises. Front. Ecol. Environ. 5, 475–482 (2007).Article 

Google Scholar 
Ohlemuller, R. Climate. Running out of climate space. Science 334, 613–614 (2011).ADS 
PubMed 
Article 

Google Scholar 
Williams, J. W., Jackson, S. T. & Kutzbach, J. E. Projected distributions of novel and disappearing climates by 2100 AD. Proc. Natl Acad. Sci. USA 104, 5738–5742 (2007).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Warren, D. L., Glor, R. E. & Turelli, M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution 62, 2868–2883 (2008).PubMed 
Article 

Google Scholar 
Broennimann, O. et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Glob. Ecol. Biogeogr. 21, 481–497 (2012).Article 

Google Scholar 
Warren, D. L., Cardillo, M., Rosauer, D. F. & Bolnick, D. I. Mistaking geography for biology: inferring processes from species distributions. Trends Ecol. Evol. 29, 572–580 (2014).PubMed 
Article 

Google Scholar 
Wobst, H. M. The archaeo-ethnology of hunter-gatherers or the tyranny of the ethnographic record in archaeology. Am Antiquity 43, 303–309 (1978).Maier, A. et al. Demographic estimates of hunter-gatherers during the Last Glacial Maximum in Europe against the background of palaeoenvironmental data. Quat. Int. 425, 49–61 (2016).Article 

Google Scholar 
Riede, F. Oxford Handbook of the Archaeology and Anthropology of Hunter-Gatherers (Oxford Univ. Press, 2014).Jochim, M., Herhahn, C. & Starr, H. The Magdalenian colonization of southern Germany. Am. Anthropol. 101, 129–142 (1999).Article 

Google Scholar 
Arts, N. & Deeben, J. On the Northwestern Border of Late Magdalenian Territory: Ecology and Archaeology of Early Late Glacial Band Societies in Northwestern Europe. In Late Glacial in Central Europe. Culture and Environment. (eds Burdukiewicz, J. M. & Kobusiewicz, M.) (Polska Akademia Nauk, Warszawa 1987).Maier, A. Population and settlement dynamics from the Gravettian to the Magdalenian. Mitteilungen der Ges. f.ür. Urgesch. 26, 83–101 (2017).
Google Scholar 
Maier, A., Liebermann, C. & Pfeifer, S. J. Beyond the Alps and Tatra Mountains-the 20-14 ka repopulation of the northern mid-latitudes as inferred from palimpsests deciphered with keys from Western and Central Europe. J. Paleolit. Archaeol. 3, 398–452 (2020).Article 

Google Scholar 
Gamble, C., Davies, W., Pettitt, P. & Richards, M., Climate change. and evolving human diversity in Europe during the last glacial. Philos. Trans. R. Soc. Lond. B Biol. Sci. 359, 243–253 (2004).PubMed 
PubMed Central 
Article 

Google Scholar 
Housley, R. A., Gamble, C. S., Street, M. & Pettitt, P. Proceedings of the Prehistoric Society. (Cambridge Univ. Press).Bellwood, P. S. First Farmers: the Origins of Agricultural Societies. (Blackwell, Oxford 2005).d’Errico, F. et al. The origin and evolution of sewing technologies in Eurasia and North America. J. Hum. Evol. 125, 71–86 (2018).PubMed 
Article 

Google Scholar 
Moseler, F. Brandstrukturen im späten Magdalénien: Betrieb, Nutzung und Funktion (Verlag des Römisch-Germanischen Zentralmuseums, 2020).Simova, I. & Storch, D. The enigma of terrestrial primary productivity: measurements, models, scales and the diversity-productivity relationship. Ecography 40, 239–252 (2017).Rosenzweig, M. L. Net primary productivity of terrestrial communities – prediction from climatological data. Am. Nat. 102, 67 (1968).Article 

Google Scholar 
Jensen, H. J. & Møberg, T. Et røgeri fra ældre stenalder ved Bølling Sø? Midtjyske Fortaellinger 2007, 51–62 (2008).Holst, D. Hazelnut economy of early Holocene hunter-gatherers: a case study from Mesolithic Duvensee, northern Germany. J. Archaeol. Sci. 37, 2871–2880 (2010).Article 

Google Scholar 
Boethius, A. Something rotten in Scandinavia: the world’s earliest evidence of fermentation. J. Archaeol. Sci. 66, 169–180 (2016).Article 

Google Scholar 
Dyson‐Hudson, R. & Smith, E. A. Human territoriality: an ecological reassessment. Am. Anthropol. 80, 21–41 (1978).Article 

Google Scholar 
Finlayson, C. The water optimisation hypothesis and the human occupation of the mid-latitude belt in the Pleistocene. Quat. Int 300, 22–31 (2013).Article 

Google Scholar 
Laland, K. N. & Brown, G. R. Niche construction, human behavior, and the adaptive-lag hypothesis. Evol. Anthropol. 15, 95–104 (2006).Article 

Google Scholar 
Laland, K. N. & O’Brien, M. J. Niche construction theory and archaeology. J. Archaeol. Method Th. 17, 303–322 (2010).Article 

Google Scholar 
Riede, F. Handbook of Evolutionary Research in Archaeology (Springer, 2019).Jöris, O. & Terberger, T. Zur Rekonstruktion eines Zeltes mit Trapezförmigem Grundriss am Magdalénien-Fundplatz Gönnersdorf/Mittelrhein: Eine» Quadratur des Kreises «? Arch.äologisches Korrespondenzblatt 31, 163–172 (2001).
Google Scholar 
Salomon, H., Vignaud, C., Lahlil, S. & Menguy, N. Solutrean and Magdalenian ferruginous rocks heat-treatment: accidental and/or deliberate action? J. Archaeol. Sci. 55, 100–112 (2015).CAS 
Article 

Google Scholar 
Nakazawa, Y., Straus, L. G., Gonzalez-Morales, M. R., Solana, D. C. & Saiz, J. C. On stone-boiling technology in the Upper Paleolithic: behavioral implications from an Early Magdalenian hearth in El Miron Cave, Cantabria, Spain. J. Archaeol. Sci. 36, 684–693 (2009).Article 

Google Scholar 
Pedersen, J., Maier, A. & Riede, F. A punctuated model for the colonisation of the Late Glacial margins of northern Europe by Hamburgian hunter-gatherers. Quart.är. 65, 85–104 (2018).
Google Scholar 
Whallon, R. Social networks and information: non-“utilitarian” mobility among hunter-gatherers. J. Anthropol. Archaeol. 25, 259–270 (2006).Article 

Google Scholar 
Leal Filho, W. et al. Impacts of climate change to African indigenous communities and examples of adaptation responses. Nat. Commun. 12, 6224 (2021).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Heitz, C. F., Hinz, M., Laabs, J. & Hafner, A. Mobility as resilience capacity in northern Alpine Neolithic settlement communities. Archaeol. Rev. Camb. 36, 75–106 (2021).
Google Scholar 
Riede, F., Oetelaar, G. A. & VanderHoek, R. From crisis to collapse in hunter-gatherer societies. A comparative investigation of the cultural impacts of three large volcanic eruptions on past hunter-gatherers. Crisis to Collapse–The Archaeology of Social Breakdown. Louvain-la-Neuve: UCL Presses Universitaires De Louvian 23–39 (2017).Halstead, P., O’Shea, J. & O’Shea, J. M. Bad Year Economics: Cultural Responses to Risk and Uncertainty. (Cambridge Univ. Press, 2004).Brovkin, V. et al. Past abrupt changes, tipping points and cascading impacts in the Earth system. Nat. Geosci. 14, 550–558 (2021).ADS 
CAS 
Article 

Google Scholar 
Burke, A. et al. The archaeology of climate change: the case for cultural diversity. Proc Natl Acad Sci USA 118, e2108537118 (2021).Binford, L. R. Willow smoke and dogs tails – Hunter-gatherer settlement systems and archaeological site formation. Am. Antiquity 45, 4–20 (1980).Article 

Google Scholar 
Birdsell, J. B. Some environmental and cultural factors influencing the structuring of Australian aboriginal populations. Am. Nat. 87, 171–207 (1953).Article 

Google Scholar 
Kelly, R. L. The Lifeways of Hunter-Gatherers: The Foraging Spectrum (Cambridge Univ. Press, 2013).Penington, R. Hunter-gatherer demography. In Hunter-Gatherers: An Interdisciplinary Perspective. (eds. Panter-Brick, C., Layton, R. H. & Rowley-Conwy, P.) (Cambridge University Press, Cambridge, 2001).Wobst, H. M. Locational relationships in Paleolithic society. J. Hum. Evol. 5, 49–58 (1976).Article 

Google Scholar 
Richards, M. P., Pettitt, P. B., Stiner, M. C. & Trinkaus, E. Stable isotope evidence for increasing dietary breadth in the European mid-Upper Paleolithic. Proc. Natl Acad. Sci. USA 98, 6528–6532 (2001).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Drucker, D. & Bocherens, H. Carbon and nitrogen stable isotopes as tracers of change in diet breadth during Middle and Upper Palaeolithic in Europe. Int J. Osteoarchaeol. 14, 162–177 (2004).Article 

Google Scholar 
Kretschmer, I. Demographische Untersuchungen zu Bevölkerungsdichten, Mobilität und Landnutzungsmustern im späten Jungpaläolithikum (Verlag Marie Leidorf GmbH, 2015).Langley, M. C. & Street, M. Long range inland-coastal networks during the Late Magdalenian: evidence for individual acquisition of marine resources at Andernach-Martinsberg, German Central Rhineland. J. Hum. Evol. 64, 457–465 (2013).PubMed 
Article 

Google Scholar 
Lanczont, M. et al. Late Glacial environment and human settlement of the Central Western Carpathians: a case study of the Nowa Biala 1 open-air site (Podhale Region, southern Poland). Quat. Int 512, 113–132 (2019).Article 

Google Scholar 
Cziesla, E. Robbenjagd in Brandenburg? Gedanken zur Verwendung großer Widerhakenspitzen. Ethnographisch-archaologische Z. 48, 1–48 (2007).
Google Scholar 
Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1‐km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

Google Scholar 
Le Cook, B. & Manning, W. G. Thinking beyond the mean: a practical guide for using quantile regression methods for health services research. Shanghai Arch. Psychiatry 25, 55 (2013).PubMed 
PubMed Central 

Google Scholar 
Yee, T. W. & Mitchell, N. D. Generalized additive-models in plant ecology. J. Veg. Sci. 2, 587–602 (1991).Article 

Google Scholar 
Guisan, A., Edwards, T. C. & Hastie, T. Generalized linear and generalized additive models in studies of species distributions: setting the scene. Ecol. Model 157, 89–100 (2002).Article 

Google Scholar 
Fewster, R. M., Buckland, S. T., Siriwardena, G. M., Baillie, S. R. & Wilson, J. D. Analysis of population trends for farmland birds using generalized additive models. Ecology 81, 1970–1984 (2000).Article 

Google Scholar 
Drexler, M. & Ainsworth, C. H. Generalized additive models used to predict species abundance in the Gulf of Mexico: an ecosystem modeling tool. PLos ONE 8, e64458 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Moisen, G. G. & Frescino, T. S. Comparing five modelling techniques for predicting forest characteristics. Ecol. Model 157, 209–225 (2002).Article 

Google Scholar 
Wood, S. N. Generalized Additive Models: An Introduction with R (CRC Press, 2006).Zuur, A. F. A Beginner’s Guide to Generalized Additive Models with R (Highland Statistics Limited, 2012).Team, R. C. R: a language and environment for statistical computing. (2013).Wood, S. N. Generalized Additive Models: An Introduction with R (CRC Press, 2017).Fasiolo, M., Wood, S. N., Zaffran, M., Nedellec, R. & Goude, Y. Fast calibrated additive quantile regression. J. Am. Stat. Assoc. 116, 1402–1412 (2021).MathSciNet 
CAS 
Article 

Google Scholar 
AlejoOrdonez/PaleoPopDen: (Version NatCommV0) [Computer software]. Zenodo. https://doi.org/10.5281/ZENODO.6962693 (2022).Liu, Z. et al. Transient simulation of last deglaciation with a new mechanism for Bolling-Allerod warming. Science 325, 310–314 (2009).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lorenz, D. J., Nieto-Lugilde, D., Blois, J. L., Fitzpatrick, M. C. & Williams, J. W. Downscaled and debiased climate simulations for North America from 21,000 years ago to 2100AD. Sci. Data 3, 160048 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Peltier, W. R., Argus, D. & Drummond, R. Space geodesy constrains ice age terminal deglaciation: the global ICE‐6G_C (VM5a) model. J. Geophys. Res. Solid Earth 120, 450–487 (2015).ADS 
Article 

Google Scholar 
Vermeersch, P. M. European population changes during the Marine Isotope Stages 2 and 3. Quat. Int 137, 77–85 (2005).Article 

Google Scholar 
Gamble, C., Davies, W., Pettitt, P., Hazelwood, L. & Richards, M. The archaeological and genetic foundations of the European population during the late glacial: Implications for ‘agricultural thinking’. Camb. Archaeol. J. 15, 193–223 (2005).Article 

Google Scholar 
Steele, J. Radiocarbon dates as data: quantitative strategies for estimating colonization front speeds and event densities. J. Archaeol. Sci. 37, 2017–2030 (2010).Article 

Google Scholar 
Shennan, S. et al. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat. Commun. 4, 2486 (2013).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Surovell, T. A., Finley, J. B., Smith, G. M., Brantingham, P. J. & Kelly, R. Correcting temporal frequency distributions for taphonomic bias. J. Archaeol. Sci. 36, 1715–1724 (2009).Article 

Google Scholar 
Williams, A. N. The use of summed radiocarbon probability distributions in archaeology: a review of methods. J. Archaeol. Sci. 39, 578–589 (2012).Article 

Google Scholar 
Kelly, R. L., Surovell, T. A., Shuman, B. N. & Smith, G. M. A continuous climatic impact on Holocene human population in the Rocky Mountains. Proc. Natl Acad. Sci. USA 110, 443–447 (2013).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Hijmans, R. J. et al. Package ‘raster’. R package 734, (2015).Lewin-Koh, N. J. et al. Package ‘maptools’. Internet: http://cran.r-project.org/web/packages/maptools/maptools.pdf (2012).Thornthwaite, C. W. An approach toward a rational classification of climate. Geogr. Rev. 38, 55–94 (1948).Article 

Google Scholar 
Lieth, H.Primary Productivity of the Biosphere (Springer, 1975). More