STUDY PROTOCOL

# Local indicators of climate change impacts described by indigenous peoples and local communities: Study protocol

Victoria Reyes-García<sup>1,2*</sup>, Santiago Álvarez-Fernández<sup>2</sup>, Petra Benyei<sup>2</sup>, David García-del-Amo<sup>2</sup>, André B. Junqueira<sup>2</sup>, Vanesse Labeyrie<sup>3,4</sup>, Xiaoyue Li<sup>2</sup>, Vincent Porcher<sup>1,2</sup>, Anna Porcuna-Ferrer<sup>1,2</sup>, Anna Schlingmann<sup>1,2</sup>, Ramin Soleymani<sup>2</sup>

1 Institució Catalana de Recerca i Estudis Avançats, Barcelona, Spain, 2 Institut de Ciència i Tecnologia Ambientals, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Barcelona, Spain, 3 Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Gestion des Ressources Renouvelables et Environnement, Montpellier, France, 4 Gestion des Ressources Renouvelables et Environnement, Université de Montpellier, Montpellier, France

* Victoria.reyes@uab.cat

<figure><img src="image_1.png" /></figure>

## Abstract

<figure><img src="image_2.png" /></figure>

### OPEN ACCESS

**Citation:** Reyes-García V, Álvarez-Fernández S, Benyei P, García-del-Amo D, Junqueira AB, Labeyrie V, et al. (2023) Local indicators of climate change impacts described by indigenous peoples and local communities: Study protocol. PLoS ONE 18(1): e0279847. https://doi.org/10.1371/journal.pone.0279847

**Editor:** Shalini Dhyani, CSIR National Environmental Engineering Research Institute (CSIR-NEERI), INDIA

**Received:** September 28, 2022

**Accepted:** December 15, 2022

**Published:** January 5, 2023

**Copyright:** © 2023 Reyes-García et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** No datasets were generated or analysed during the current study. All relevant data from this study will be made available upon study completion.

**Funding:** This work is funded by the European Research Council (ERC) under grant agreement No 771056-LICCI-ERC-2017-COG. This work contributes to the ICTA-UAB "María de Maeztu" Programme for Units of Excellence funded by the

### Introduction

In the quest to improve the understanding of climate change impacts on elements of the atmospheric, physical, and life systems, scientists are challenged by the scarcity and uneven distribution of grounded data. Through their long history of interaction with the environment, Indigenous Peoples and local communities have developed complex knowledge systems that allow them to detect impacts of climate change in the local environment. The study protocol presented here is designed 1) to inventory climate change impacts on the atmospheric, physical, and life systems based on local knowledge and 2) to test hypotheses on the global spatial, socioeconomic, and demographic distribution of reported impacts. The protocol has been developed within the framework of a project aiming to bring insights from Indigenous and local knowledge systems to climate research (https://licci.eu).

### Methods

Data collection uses a mixed-method approach and relies on the collaboration of a team of 50 trained partners working in sites where people's livelihood directly depend on nature. The data collection protocol consists of two steps. Step 1 includes the collection of secondary data (e.g., spatial and meteorological data) and site contextual information (e.g., village infrastructure, services). Step 1 also includes the use of 1) semi-structured interviews (n = 20-30/site) to document observations of environmental change and their drivers and 2) focus group discussions to identify consensus in the information gathered. Step 2 consist in the application of a household (n from 75 to 125) and individual survey (n from 125 to 175) using a standardized but locally adapted instrument. The survey includes information on 1) individual and household socio-demographic characteristics, 2) direct dependence on nature, 3) household's vulnerability, and 4) individual perceptions of climate change impacts. Survey data are entered in a specifically designed database.
Spanish Ministry of Science, Innovation and
Universities (CEX2019-000940-M).

**Competing interests:** The authors have declared
that no competing interests exist.

## Expected results

This protocol allows the systematic documentation and analysis of the patterned distribution of local indicators of climate change impacts across climate types and livelihood activities. Data collected with this protocol helps fill important gaps on local climate change impacts research and can provide tangible outcomes for local people who will be able to better reflect on how climate change impacts them.

## Introduction

There is overwhelming evidence that anthropogenic climate change has a global influence on elements of the atmospheric, physical, and life systems [1, 2], with direct effects on local livelihoods and cultures [3]. However, scientists have a more meagre understanding on how climate change differently impacts local physical and life systems across the globe [4]. Natural scientists acknowledge that downscaling global models to resolutions that are relevant for regional or local level policymaking (e.g., adaptation planning) is challenging due to the uneven data availability across regions and to the uncertainties introduced by downscaling techniques [5]. Moreover, uncertainty in modeling exercises is even larger due to the complex interactions between biophysical systems and human society, which may result in unexpected changes and synergetic impacts [1]. Scientists also argue that a better grasp of climate change impacts requires to couple our understanding of ecological and social dynamics, for which climate change research needs to fully take into consideration impacts on local socioeconomic systems [1, 3, 6, 7].

Different societies are differently affected by climate change, not only because climate change impacts vary along place-specific geographical and ecological characteristics, but also because climate change affects people and societies through specific pathways largely mediated by socio-cultural factors. For example, while sea-level rise is a climate-related phenomenon potentially affecting millions of people living close to sea level, its actual impacts depend not only on biophysical conditions (e.g., magnitude of tidal influences, overall island size and relief) but also on site-specific socioeconomic conditions (e.g., resources to cope with sea-level rise, dependence on fisheries). In sum, the scientific community largely agrees that there is a need to better understand localized climate change impacts and how such impacts are experienced. Yet, in this quest, scientists are challenged by the scarcity of grounded data, which has resulted in climate scientists calling for the exploration of new sources of information [2, 8].

Indigenous and local knowledge systems can contribute to our understanding of local climate change impacts [9–11]. Throughout the world, Indigenous Peoples and local communities (IPLC) with a long history of interaction with their environment have developed intricate and complex knowledge systems (including knowledge, technologies, and forms of organization) that allow them to detect not only changes in local weather and climatic variability [12, 13], but also the direct effects of such changes on the environment on which they depend [14]. Moreover, Indigenous and local knowledge systems are particularly suited to explain how climate change affects local social-ecological systems, livelihoods, and cultures, because of the long-term detailed observation of these systems [2, 15].

Many works have documented Indigenous Peoples and local communities’ first-hand observations of changes in social-ecological systems attributed to climate change [see 9, 10, 14 for reviews]. These works provide rich multi-site qualitative place-based information, but do not consider a common strategy to gather place-specific, yet comparable, knowledge from
different locations. Consequently, insights from local knowledge systems continue to lack transferability, integration, and scalability to climate research and policy [10]. To increase transferability, there is a need for a research tool able to draw correspondences between Indigenous and local knowledge's qualitative and interpretative nature into categories that allow generalization. To allow knowledge integration, there is a need for a tool that allows combining inputs from multi-site place-based research. Finally, to address scalability, there is a need for a community of practice that considers the need to ensure that placed-based information is effectively upscale to climate research [10].

The project Local Indicators of Climate Change Impacts (LICCI, https://licci.eu) was designed to address these challenges through 1) the development of a standardized protocol for the collection of cross-culturally comparable data on local indicators of climate change impacts and 2) the creation of a network of researchers who use this standard protocol. In this paper, we describe the research protocol used for collecting cross-culturally comparable information on direct, observable, localized climate change impacts on elements of the local atmospheric, physical, and life systems and covariates that might affect how this knowledge is distributed within and across societies. We define 'local indicators of climate change impacts' (hereafter LICCI) as first-hand observations of environmental change reported by people living in close interaction with a particular environment and attributed-totally or partially- to changes in elements of the atmospheric system (i.e., temperature, precipitation, winds, or seasonality). The overall framework of the project can be found in [10] and the specific tools used for data collection can be found in [16].

## Materials and methods

### Objectives, hypotheses, and research strategy

The LICCI study protocol aims to inventory and identify patterns in reports of local indicators of climate change impacts. Specifically, data collected with this protocol will be used to test hypotheses related to the globally patterned *i*) spatial, *ii*) socioeconomic, and *iii*) demographic distribution of LICCI. Table 1 summarizes the premises, provides examples, and presents the formal hypotheses that guide this work.

<table><caption>Table 1. Hypotheses that guide the protocol.</caption><thead><tr><th></th><th>Premise</th><th>Example</th><th>Hypothesis</th></tr></thead><tbody><tr><td rowspan="2">Spatial patterns</td><td rowspan="2">Climate change impacts are not uniform.</td><td>Climate change impacts vary spatially in occurrence and intensity [17].</td><td>H1a. Groups living in strongly climate-affected areas will provide higher number and more diverse reports of LICCI than groups living in less climate-affected areas.</td></tr><tr><td>Climate change impacts vary across biomes and climatic regions.</td><td>H1b. Groups living within the same climatic regions will report similar LICCI.</td></tr><tr><td rowspan="2">Socioeconomic patterns</td><td rowspan="2">Climate change impacts are not equally distributed across the social spectrum.</td><td>Climate change impacts affect more intensively people with direct dependence on natural resources [2].</td><td>H2a. The higher a household's direct dependence on natural resources, i) the more and ii) the higher severity of LICCI reported.</td></tr><tr><td>Climate change impacts more intensively people who already face other stressors, such as poverty or lack of access to natural resources [2, 18].</td><td>H2b. The higher a household's vulnerability, i) the more LICCI reported and ii) the higher their severity.</td></tr><tr><td rowspan="2">Demographic patterns</td><td rowspan="2">As other components of local knowledge systems [19], the distribution of knowledge on climate change impacts might be patterned by individual demographic characteristics.</td><td>Elders' reports of changes might have a wider temporal depth than youngsters' reports [20].</td><td>H3a. Elders will i) report more LICCI and ii) a higher severity of their impacts than youngsters.</td></tr><tr><td>Because their different roles and responsibilities [21], women and men might observe different climate change impacts.</td><td>H3b. In each site, the similarity on LICCI reports will be higher among women (or among men), than between women and men.</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0279847.t001
<table><thead><tr><th rowspan="3"></th><th colspan="6">Preparation</th><th colspan="6">Data collection</th><th colspan="6">Analysis</th></tr><tr><th colspan="3">Year 1</th><th colspan="3">Year 2</th><th colspan="3">Year 3</th><th colspan="3">Year 4</th><th colspan="3">Year 5</th></tr><tr><th>2018</th><th>2019</th><th></th><th>2018</th><th>2019</th><th></th><th>2020</th><th>2021</th><th></th><th>2022</th><th>2023</th><th></th><th>2022</th><th>2023</th><th></th></tr><tr><th></th><th>Q1</th><th>Q2</th><th>Q3</th><th>Q4</th><th>Q1</th><th>Q2</th><th>Q3</th><th>Q4</th><th>Q1</th><th>Q2</th><th>Q3</th><th>Q4</th><th>Q1</th><th>Q2</th><th>Q3</th><th>Q4</th></tr></thead><tbody><tr><td><strong>Preparation</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Recruit core team</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Develop instruments</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pilot test</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Recruit & train partners</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Data collection</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Collect field data (Steps 1 &2)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Analysis</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Data curation</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Data analysis</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Note:** Project duration 1st June 2018 to 31st May 2023.

**Fig 1. Gantt chart.** Development and application of the LICCI protocol. **Note:** Project duration 1st June 2018 to 31st May 2023.

https://doi.org/10.1371/journal.pone.0279847.g001

The LICCI research project started in June 2018 and has a 5-years duration, although the protocol could be used by independent researchers afterwards. Within the project frame, the design, implementation, and use of data collected through this protocol is organized in three phases (Fig 1): The preparation phase (June 2018-December 2019), the data collection phase (January 2020-December 2022), and the data analysis phase (June 2022- June 2023). The estimated time to complete data collection is two person/months, but the overall period of data collection is longer to accommodate different partners' needs. Due to calendar adjustments derived from the COVID-19 pandemic, data collection and data curation overlapped. The analysis of global data will start once all datasets are curated (December 2022).

## Preparation

The preparation phase includes the development and testing of data collection instruments, the site selection and partners' recruitment, and partners' training.

**Development and testing of data collection instruments.** The original development of the data collection instruments was done by the authors. The overall protocol strategy was conceived to rely on a mixed-method approach and the protocol is divided in two steps. Step 1 focuses on the collection of secondary and primary background information, including a list of observations of environmental changes and its drivers, and uses focus group discussions to identify consensus in the information gathered. Step 2 includes a survey to households and individuals to gather data for hypotheses testing.

To improve the overall clarity of the protocol, minimize participant's fatigue, and avoid sensitive issues that might cause discomfort, we pre-tested all methods in the field. Instruments to be used in Step 1 (i.e., semi-structured interviews and FGDs) were tested with Akha smallholders communities in China, with Betsimisaraka agroforestry farmers in Madagascar, and with Sereer agriculturarists in Senegal. Such testing resulted in 1) a better focus on environmental (vs. global) change, 2) the refinement of our system to code observations of environmental change into indicators of climate change; and 3) a larger consideration of drivers of environmental change other that climate change. Instruments to be used in Step 2 (i.e., individual and household surveys) were tested with Tsimane' Indigenous Peoples in the Bolivian Amazon and Daasanach agro-pastoralists in Kenya. The second testing also resulted in adjustments of
<figure><img src="image_3.png" /><figcaption>**Fig 2. Selected study sites.** Sites are displayed across the main different climate types defined by Koeppen-Geiger [22]. Figure elaborated by the authors. Original map from the Natural Earth data set in the public domain (available from https://www.naturalearthdata.com), elaborated with the R libraries maptools [23] and sp [24].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0279847.g002

the instrument, including removing sensitive or repetitive questions, or adding questions to capture factors not considering, and reordering or rewording questions to improve the flow.

**Site selection and partners' recruitment.** In 2019, the author team worked in the selection of sites where the protocol would be applied. Fig 2 shows the geographical location of the study sites where the protocol has been applied.

The following criteria guided site selection:

1. **Data availability:** World coverage of weather stations leaves large areas of the world with poor, incomplete, or unreliable weather records, and the literature on local knowledge of climate change impacts is also spatially biased [9, 10]. Hence, we prioritized the selection of sites where instrumental data are deficient and no/few studies on local knowledge of climate change impacts had been conducted.

2. **Climate type:** To test hypotheses *H*1a and *H*1b, on LICCI spatial distribution, site selection was done to cover the five main different climate types defined by Koeppen-Geiger [22]: tropical, arid, temperate, continental, and polar/cold.

3. **Predominant livelihood activity:** To test hypotheses *H*2a and *H*2b, about the socioeconomic distribution of LICCI, site selection was done to include sites mainly depending on four livelihood activities (i.e., agriculture, fishing, pastoralism, and foraging).

4. **Partnership feasibility:** Site selection also included logistical considerations derived from the need to establish an extended network of partners responsible for data collection. To recruit partners, a call was widely circulated. Selection criteria included research experience and previous relations with the proposed study site. The call encouraged South-based researchers to apply.

**Partners' training.** To minimize interviewer and coder biases during the application of the protocol, all partners participated in a training workshop. Originally, we conducted three one-week face-to-face workshops during 2019 in Barcelona with 48 partners (S1 File). In March 2021, we conducted one additional two-day online workshop to train a set of 10 new
partners to substitute partners who had to abandon the project because of the COVID-19
pandemic.

During the workshops, partners were exposed to the rationale of the project and received
detailed explanations concerning the application of all the instruments of data collection. Part-
ners also had the opportunity to discuss their schedule for data collection and other practicali-
ties. The workshop included discussions on the code of ethics for the inclusion of Indigenous
knowledge in research and the elaboration of a "Local Knowledge Research Agreement" to be
discussed and negotiated with communities. The training workshop was also used to discuss
issues related to data ownership and sharing (S1 File). Training also included safety consider-
ations related to 1) personal health (e.g., medical contacts, health insurance coverage, vaccina-
tion) and safety (e.g., avoiding places locally considered dangerous, avoiding openly exhibiting
valuables, driving at night), 2) preparing for emergencies (e.g., allying with a local institution
or a trusted local leader, registering with local authorities), and 3) avoiding conflict (e.g., select-
ing non-conflictive villages). After the training workshop, partners signed a contract with the
host institution (Universitat Autònoma de Barcelona) detailing the terms of the agreement
(S1 File).

**Ethical considerations.** The research protocol was approved by the Ethics Committee of the Universitat Autonoma de Barcelona (CEEAH 4781) and the LICCI project follows the European Research Council ethical guidelines. Ethical considerations include obtaining Free, Prior and Informed Consent (FPIC) and following Safety and Security guidelines and guidelines for ethical conduct discussed during partners training. All procedures and documents were reviewed by an external and independent Ethics Advisor (Prof. Dr. Michael Schönhuth), who provided continuous reports to the project funders and feedback to the team.

The protocol specifies that before the onset of data collection, all partners need to obtain
Free, Prior and Informed Consent (FPIC) from the organizations representing the communi-
ties where they aim to work, the villages and the people we work with (S1 File). In the first visit
to a village, partners held a meeting to present the research and ask for consent to stay in the
village. In the community meetings, partners presented detailed information on the objectives
and scopes of the study, the participation of subjects, and the costs and benefits associated with
participation, and obtain written consent to participate. In this meeting, partners were also
instructed to hold an open discussion leading to a "community engagement protocol" in
which the participants had the opportunity to ask in which way they want the information to
be returned, or other requirements. Finally, partners also asked for FPIC from each individual.

**Data collection**

In this section, we describe the sampling strategy to select villages within a specific site and the data collection methods for Steps 1 and 2. The logic of the protocol and the relation between instruments is depicted in Fig 3.

**Sampling villages.** In each site, one partner is responsible for the collection of data across 3–5 villages, defined as the lowest administrative unit in an area and normally under the jurisdiction of a village leader/council. Selected villages should be representative and relatively homogeneous in terms of the environmental and socio-cultural conditions of the site. Villages with particular conditions (e.g., villages with high donor intervention) are to be avoided. To reduce complexity in logistics, targeted villages should have between 20 and 500 households. Villages with more than 500 households should be broken into smaller units and sample one or a few of them.

**Step 1.** In Step 1, partners collect secondary and primary background information from different sources and using different sampling strategies.
<figure><img src="image_4.png" /><figcaption>**Fig 3.** Diagram of data collection protocol displaying the information to be collected, the data source or sample, and the method of data collection. The protocol includes two steps. Step 1 includes the collection of secondary background information, site contextual information, observations of environmental change, and consensual indicators of climate change impacts. Step 2 consists the collection of household and individual characteristics through a survey instrument.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0279847.g003

**Background information** to be obtained includes general information (e.g., historical events, protected areas limits) and meteorological data. Background information could be obtained before conducting fieldwork from secondary sources (e.g., literature review, public records, etc.), or while in the field through participant observations and conversation with local officials/knowledgeable individuals. Partners should also obtain geo-referenced data on the location of 1) the study site, villages, neighboring protected areas, and officially recognized Indigenous territories; 2) the closest weather station, airport, and/or harbor/port; and 3) the closest market town and administrative center. Partners will obtain 20-years (minimum) series of data on temperature and precipitation at monthly resolution (minimum) from weather stations within 200 km (maximum) from the study site.

**Site contextual information** includes the collection of information on local livelihoods and dependence on the natural environment and a timeline of events that are important to the community. Through field observations, secondary sources, and interviews to local experts, partners will collect information about 1) local livelihood activities, including the timing (i.e., seasonal calendar), location, and distribution of these activities, 2) the timing of important events in the recent history of the community (i.e., dated back 80–100 years), and 3) a list of household assets with market value (i.e., fishing nets, machetes etc.), which capture variation in ownership of market assets across households.

**Observations of environmental change and their drivers:** Partners will use semi-structured interviews to obtain a list of observations of environmental changes. To recruit participants for these interviews, partners will use “quota sampling”. To capture diversity in perceptions of climate change impacts, each partner should select 20–30 participants aiming for a minimum of three interviews per quota across gender (men and women), age (young, a middle aged, and elder), and livelihoods (site dependent). To select people within each quota, partners will rely on key-informants.
Interviews start with the question “Compared to when you were young, what changes in the environment have you noticed?” Additional questions can be then directed to stimulate interviewees to report changes in (a) the atmospheric (e.g., weather/seasons, temperature, rainfall and snowfall, wind, storms); (b) physical (e.g., soil, river, streams), and (c) life system (e.g., wild animals, wild plants, crops, pastures). For each observation of change, partners will also ask about the direction of change (e.g., increase/decrease, earlier/later) and the driver of the change (i.e., “why do you think that this happens? Or “What do you think is the cause of that?”). If the answer to this question is another environmental change, the partner should continue asking about drivers until the respondent could not identify any further driver of the reported observations of change. Note that, to avoid inducing answers, at this stage, the protocol focuses on ‘environmental changes’ in a broad sense and not on ‘climate change’. The attribution to climate change (and/or other drivers) is done *a posteriori*, combining the information on observations of change with that on drivers of change.

**Consensual indicators of climate change impacts.** Each partner should organize 3–5 focus group discussions (FGD) in the study site to discuss information from semi-structured interviews. Between 4 and 12 participants should participate in each discussion. To capture the site diversity in terms of livelihood activities (e.g., farmers, fishers, herders), age, and gender, and to avoid potential conflicts or power imbalances, participants will be selected through convenience sampling.

To prepare for FGD, partners should first classify observations of environmental changes. Following the classification proposed by [14, see also S1 File], partners should start by grouping verbatim observations depicting the same environmental change into ‘indicators’ of change. The classification used differentiates between indicators refereeing to changes in elements of i) the atmospheric (e.g., temperature, rain), ii) the physical (e.g., water and soil temperature, sea-level rise) and iii) the life systems (e.g., the morphology, abundance, or distribution of wild and managed/cultivated plant and animal species). For each indicator, partners should also compile information on all the drivers leading to that change. In FGD, partners should discuss the resulting list of indicators of change and their drivers. In particular, partners should ask FGD participants whether they agree on the reported changes and their drivers, focusing on potentially unclear or contradictory information.

Using reports from semi-structured interviews and FGD, partners should also code the level of group agreement in three categories: ‘agreement’ i) when there is consistency in the report of the indicator and its drivers in semi-structured interviews or ii) when FGD participants agreed with the indicator and its driver during the discussion; ‘inconclusive’ when observations referring to that indicator are mentioned by less than 10% of people participating in semi-structured interviews and the indicator is not discussed in FGD; and ‘disagreement’ when observations referring to the same indicator report different directions or drivers and/or there is no agreement on the indicator and its drivers during FGD.

The consensual list of LICCI for a site should include all indicators for which there is agreement and that are reportedly driven (directly or through cascading effects) by changes in elements of the atmospheric system. The participatory, collective, and iterative nature of the process ensures that the final list of LICCI in a site reflects the group *social memory*.

**Step 2.** Step 2 consists of the application of a structurally identical survey, but with site-specific questions derived from information collected in Step 1. The survey includes questions directed to the household and questions directed to the individuals (see [16] for the exact survey questions). To select households (between 75 and 125), partners should use *simple random sampling*, drawing households from a local census. Household questions can be answered by anyone who is considered a household head, defined as a person who–by formal or informal rules– can take decisions on household’s labor and income. Within each household, partners should use *convenience quota sampling* among household heads to independently answer
individual questions. Individuals (between 125 and 175) should be approximately evenly distributed across gender and age categories.

The survey will include four sections.

1. *Individual and household socio-demographic characteristics* including household composition and information from household members (i.e., level of formal education, occupation, migratory status), history of family belonging to the region, access to meteorological information, and individual's subjective wellbeing and climate change awareness.

2. *Individual and household dependence on natural resources* will be proxied using the 'pebble distribution method' [25, 26], in which respondents will be asked to distribute a given number of points ('pebbles') across different livelihood activities practiced locally, placing pebbles based on the time invested (individual level) and income obtained (household level) from each activity. Livelihood activities include natural resource-related activities as well as activities unrelated to natural resources.

3. *Household's vulnerability:* We draw on the sustainable rural livelihoods approach [27] and proxy vulnerability as a factor of (lack of) access to the five capitals (i.e., financial, physical, human, social, and natural capital) that allow households to pursue their livelihood strategies. To assess *financial capital*, we ask questions about income from different sources, savings, and credit. To assess *physical capital*, we ask about ownership of 10 market assets, housing material, and access to communal assets (e.g., access to roads or local infrastructure). To assess *human capital*, we ask about household composition, formal education, number of languages spoken, and local knowledge. To assess *social capital*, we ask about membership to associations, social relations inside and outside the village, and perceptions of trust and cooperativeness. Finally, to assess access to *natural capital*, we ask about access to water and local food, land and livestock ownership, and access to communal resources (e.g., common forest, pastures, fisheries) for both sustenance and income generation.

4. *Individual perceptions of climate change impacts:* Partners will use the site-specific list of consensual LICCI to randomly select 15 LICCI to be included in the survey. For each LICCI, we will ask survey respondents to report a) whether they have observed the change (i.e., *indicator perceived*), 2) whether that change has any noticeable (positive or negative) impact on their livelihood (i.e., *impactful indicator perceived*), and 3) the *severity of the impact*. Adapting a method developed in medical research [28], to measure the severity of a LICCI, we will ask whether the impact affected the person 1) a lot, 2) a little, 3) not at all.

As some survey questions are site-specific a draft version of the survey should be pre-tested with 10 informants. The draft version of the survey should include 30 LICCI from the consensuated list and 15 commercial items displaying variation. Results from the test will be used to select the assets and LICCI to be used in the final survey to be applied in a site.

## Data analysis

*To test hypotheses on spatial patterns*, we will rely on site-specific lists of LICCI and geographical information. We will elaborate a map displaying all LICCI in the inventory and use open-source platforms to visualize the geo-referenced database of reported impacts and use spatial matching techniques to search for patterns across different *i*) climate-affected areas [29] (*H1a*) and *ii*) Köppen climatic areas* [17] (*H1b*).

*To test hypotheses on socioeconomic patterns*, we will use survey data. We will aggregate responses to survey questions to create three dependent (i.e., *number of indicators perceived, number of impactful indicators perceived, severity of impacts*) and two explanatory variables
(i.e., *direct dependence on natural resources* and *vulnerability*) to be used in multivariate regressions. The *number of indicators* and *impactful indicators perceived* correspond to the sum of LICCI from those included in our survey instrument reported by the interviewed person. The score of *perceived severity* will be a composite measure created from data on the perceived severity of the various LICCI included in the survey. The number of pebbles allocated to different livelihood activities will be used as a measure of the importance of that activity to the respondent and the number of pebbles allocated to natural resource-related activities will be used as a measure of household *dependence on natural resources*. We will create an index of household *vulnerability* by normalizing the measures of the different indicators of household's relative access to the five asset types defined and aggregating them into a single measure. To estimate how the *number of indicators*, *impactful indicators* and *severity of impacts perceived* vary across households depending on their direct dependence on natural resources (*H2a*) and vulnerability (*H2b*), we will use regression analysis. Specifically, we will model each of the three dependent variables as a function of *direct dependence on natural resources* and *vulnerability* while controlling for additional variables that might affect the studied relation (e.g., site). As data to test H2 are naturally hierarchical in structure, with households nested within villages, and villages nested within climatic areas, we will consider the multilevel structure of model variance by fitting mixed-effect models. For the statistical analysis, we will use the package *nlme* [30] in R statistical software version 4.2.1.

*To test hypotheses on demographic patterns*, we will use the exact same dependent variables and a similar analytical approach. In this analysis, we will use as main explanatory variables the *age* or *position in the life cycle* (i.e., young, adult, or elder), as aggregated measures of age might improve the reliability of estimates (*H3a*), and the *sex* of the informant (*H3b*).

## Data management

Partners submit data through a web application specifically developed for this study protocol and accessible for registered partners (S1 File). The application stores data in the browser of the device used, which can then be submitted to a particular website following a set of guidelines (S1 File). To strengthen intercoder reliability and minimize data coding and submission errors, during the data collection period (2020–2022), authors closely worked with partners to clarify doubts in data collection, coding, and submission.

The LICCI project participates in the Horizon 2020 Open Research Data pilot and follows the ‘FAIR data principles’, for which it will make the project’s research data findable, accessible, interoperable, and re-usable. Procedures for the management of data, from collection to use after the project ends are outlined in the project’s Data Management Plan (DMP), which includes data and metadata standards for use of data, policies for access and sharing, policies for re-use and distribution, and plans for data archiving and preservation.

Data collected with this protocol will only be made open if study participants agree on it. Results from Step 1 of the project will be made available only in aggregated form. Anonymized survey data will be made available with the publication of research results, also in an aggregated format. All data will be stored on and available from Harvard's Dataverse. To describe the underlying data, metadata with sufficient detail will be created to be intelligible for other users. The LICCI metadata will follow a generalized scheme in the Dataverse including, among others, title, creator(s) and contact person(s); date; version; location; contributor, e.g., funding body, including the grant agreement number; data format(s); keywords; identifiers (DOI and url); access rights (license(s)), and other labels, including labels on the particularities for access and use of cultural heritage data that is digitally circulating outside community contexts to respect local data sovereignty [31]. The LICCI team applies an embargo period of 24
months from data submission to data release to ensure that the partner collecting the data has a first right to publication. In June 2024, the data will be made public on Dataverse after a thor- ough anonymization process to ensure no sensitive information is released, considering data sensitivity from an individual and a cultural perspective.

Discussion

This section discusses 1) the main opportunities of using this protocol and the potential impact of data collected using it and 2) the main challenges of using this protocol and the limitations of data collected using it.

Opportunities and potential impact

Researchers are challenged to provide robust and implementable information on climate
change impacts, particularly at the local scale. The protocol presented here provides an oppor-
tunity to address this challenge for several reasons. First, the inventory of LICCI will improve
our understanding of the localized responses of the physical and life systems to climate change.
Importantly, such data refers to remote areas, difficult to extensively sample with traditional
field methods. The data collected through this protocol can complement research based in
other sources of data (e.g., weather stations, remote sensing, biological surveys) to help define
more precisely the linkages between climate change and the way it is reflected in local systems.
The potential to fill data gaps regarding how climate change impacts the biophysical world
directly answers the call of the Intergovernmental Panel on Climate Change (IPCC) to develop
the evidence base for the potential contribution of local knowledge to climate research [32].

Second, data collected through this protocol will improve our understanding of the impacts of climate change on local livelihoods, from changes in cropping patterns and damages related to extreme events [33], to disease outbreaks or to the increase in conflicts over dwindling resources [34]. While responding to climate change demands adjusting to weather and climate changes; identifying risks, making decisions on how to respond, and implementing such deci- sions are all mediated by socioeconomic and cultural factors [3]. The information collected through this protocol will expand our knowledge on climate change impacts on local liveli- hoods and thus improve our understanding of the human dimensions of climate change and guide policies attempting to mitigate its immediate impacts.

Third, the involvement of local communities in data collection holds the potential to increase local agency for adaptation [35]. Local people will actively participate in identifying climate change impacts and the resulting information will be reported back to them. This pro- cess increases the likelihood that local people take ownership over the process and results, potentially empowering them in crafting, determining, and adopting viable strategies for cop- ing with and adapting to climate change immediate impacts on their livelihood.

Challenges and limitations

The use of the protocol faces three main challenges and one important limitation. First, the use
of the protocol faces the challenge to mobilize Indigenous and local knowledge in a structured
and organized way, while -at the same time- defying the prevailing technical approach to docu-
ment information from other knowledge systems which seldom captures their intricacies [15].
To overcome the challenge, we selected partners with previous or planned long-term contact
with the study site and train them to collect information that reflects group social memory.
We recommend that further applications of this protocol continue to rely on researchers with
long-term involvement with local communities or on researchers-local communities’ partner-
ships to limit misinterpreting Indigenous and local knowledge.
The second challenge is to adequately capture the temporal dimension of change, as neglecting it might lead to conflate climate change and weather variability. For example, in a case study derived from this work [33], researchers found that Tuareg participants reported the impact of rain and temperature irregularities and severe drought events on their pastoral and semi-pastoral livelihoods. Paradoxically, they did not explicitly report decadal trends in the frequency of extreme events. The differential perception of climate change impacts across time scales can have important implications for undertaking climate change adaptation measures. To address the challenge, we designed the protocol to capture information encapsulated in the group's *social memory*, or the cumulative and shared knowledge encompassing long-term, intergenerational observations of immersive experience in a particular place [36]. Since social memory provides a more accurate temporal assessment than individual observations, we recommend that future work continues developing research based on social (rather than individual) information.

A third challenge of the use of this protocol refers to the need to be adapted to different context and social-ecological systems. To overcome the challenge, we put large emphasis in designing a flexible protocol. This strategy also increases the diversity of outcomes from the project, as reflected in the first case studies published by partners who have used data collected through Step 1 [33, 34, 37]. We also invested in the training of partners to understand the concepts behind the questions, so they could use their experience to locally adapt the protocol. For example, due to the COVID-19 pandemic, in many sites was not possible to conduct FGD to maintain social distancing. Several partners modified the protocol to validate the consensual nature of indicators collected in semi-structured interviews through additional interviews to other informants. Given the large carbon impact of our training strategy [38], we have created a number of written and online training materials, as well as a Massive Open Online Course (MOOC) about climate change and Indigenous peoples and local communities (S1 File). We recommend the use of these online materials to guide and train researchers aiming to use this protocol in the future. The use of the online training materials will continue to ensure comparability across datasets, while minimizing carbon impact.

The main limitation of this research relates to the ability to discern whether reported changes in elements of the atmospheric, physical, and life systems can be really attributed to climate change, as complexity in such systems and confounding drivers of change (e.g., geological processes, land-use change, or human demography shifts) make attribution difficult [8]. Indeed, researchers increasingly recognize the synergistic effects of many drivers of change [1]. While such problem is common to any research on climate change impacts, it seems to be aggravated when assessing localized impacts, the focus of this work. We cannot fully tackle this limitation, but we notice that for many Indigenous Peoples and local communities the effect of different drivers of change are actually intertwined.

## Supporting information

**S1 File. Online resources associated to the LICCI project.**
(DOCX)

**S2 File. Online resources associated to the LICCI project.**
(DOCX)

## Acknowledgments

We thank all people who have and who will participate in this study for their time and insights. We thank L. Calvet-Mir, C. Drouin, E. Poncela for their help in managing the project and all LICCI partners and field assistants for their insights while testing the protocol.
## Author Contributions

**Conceptualization:** Victoria Reyes-García.

**Data curation:** Santiago Álvarez-Fernández, Ramin Soleymani.

**Funding acquisition:** Victoria Reyes-García.

**Investigation:** André B. Junqueira, Xiaoyue Li, Vincent Porcher, Anna Porcuna-Ferrer.

**Methodology:** Victoria Reyes-García, Petra Benyei, David García-del-Amo, André B. Junqueira, Vanesse Labeyrie, Xiaoyue Li, Vincent Porcher, Anna Porcuna-Ferrer, Anna Schlingmann, Ramin Soleymani.

**Project administration:** Victoria Reyes-García, Xiaoyue Li.

**Software:** Ramin Soleymani.

**Supervision:** Victoria Reyes-García.

**Validation:** Petra Benyei, David García-del-Amo, André B. Junqueira, Xiaoyue Li.

**Visualization:** Vincent Porcher.

**Writing – original draft:** Victoria Reyes-García.

**Writing – review & editing:** Victoria Reyes-García, Santiago Álvarez-Fernández, Petra Benyei, David García-del-Amo, André B. Junqueira, Vanesse Labeyrie, Xiaoyue Li, Vincent Porcher, Anna Porcuna-Ferrer, Anna Schlingmann, Ramin Soleymani.

## References

1. Pörtner H-O, Scholes B, Agard J, Archer E, Arneth A, Bai X, et al. Scientific outcome of the IPBES-IPCC co-sponsored workshop on biodiversity and climate change. 2021 [cited 4 Apr 2022]. https://doi.org/10.5281/ZENODO.4659159

2. IPCC. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Pörtner H-O, Roberts DC, Tignor M, Poloczanska ES, Mintenbeck K, Alegría A, et al., editors. Cambridge: Cambridge University Press; 2022.

3. Adger WN, Barnett J, Brown K, Marshall N, O'Brien K. Cultural dimensions of climate change impacts and adaptation. Nat Clim Chang. 2013; 3: 112–117. https://doi.org/10.1038/nclimate1666

4. Stott PA, Gillett NP, Hegerl GC, Karoly DJ, Stone DA, Zhang X, et al. Detection and attribution of climate change: a regional perspective. Wiley Interdiscip Rev Clim Chang. 2010; 1: 192–211. https://doi.org/10.1002/wcc.34

5. Zou L, Zhou T, Li L, Zhang J, Zou L, Zhou T, et al. East China Summer Rainfall Variability of 1958–2000: Dynamical Downscaling with a Variable-Resolution AGCM. J Clim. 2010; 23: 6394–6408. https://doi.org/10.1175/2010JCLI3689.1

6. Crate SA. Climate and Culture: Anthropology in the Era of Contemporary Climate Change. Annu Rev Anthropol. 2011; 40: 175–194. https://doi.org/10.1146/annurev.anthro.012809.104925

7. Arneth A, Shin YJ, Leadley P, Rondinini C, Bukvareva E, Kolb M, et al. Post-2020 biodiversity targets need to embrace climate change. Proc Natl Acad Sci U S A. 2020; 117: 30882–30891. https://doi.org/10.1073/pnas.2009584117 PMID: 33288709

8. Rosenzweig C, Neofotis P. Detection and attribution of anthropogenic climate change impacts. Wiley Interdiscip Rev Clim Chang. 2013; 4: 121–150. https://doi.org/10.1002/wcc.209

9. Savo V, Lepofsky D, Benner JP, Kohfeld KE, Bailey J, Lertzman K. Observations of climate change among subsistence-oriented communities around the world. Nat Clim Chang. 2016; 6: 462–473. https://doi.org/10.1038/nclimate2958

10. Reyes-García V, García-del-Amo D, Benyei P, Fernández-Llamazares Á, Gravani K, Junqueira AB, et al. A collaborative approach to bring insights from local observations of climate change impacts into global climate change research. Curr Opin Environ Sustain. 2019;39. https://doi.org/10.1016/j.cosust.2019.04.007
11. Orlove B., Dawson N., Sherpa P., Adelekan I., Alangui W., Carmona R., et al. ICSM CHC White Paper I: Intangible Cultural Heritage, Diverse Knowledge Systems and Climate Change: Contribution of Knowledge Systems Group I to the International Co-Sponsored Meeting on Culture, Heritage and Climate Change. Charenton-le-Pont & Paris, France: ICOMOS & ICSM CHC, 2022. ISBN: 978-2-918086-71-0. URL: https://openarchive.icomos.org/id/eprint/2717/

12. Orlove BS, Chiang JCH, Cane MA. Forecasting Andean rainfall and crop yield from the influence of El Niño on Pleiades visibility. Nature. 2000; 403: 68–71. https://doi.org/10.1038/47456 PMID: 10638752

13. Fernández-Llamazares Á, Méndez-López ME, Díaz-Reviriego I, McBride MF, Pyhälä A, Rosell-Melé A, et al. Links between media communication and local perceptions of climate change in an indigenous society. Clim Change. 2015;131. https://doi.org/10.1007/s10584-015-1381-7 PMID: 26166919

14. Reyes-García V, Fernández-Llamazares Á, Guèze M, Garcés A, Mallo M, Vila-Gómez M, et al. Local indicators of climate change: the potential contribution of local knowledge to climate research. Wiley Interdiscip Rev Clim Chang. 2016; 7: 109–124. https://doi.org/10.1002/wcc.374 PMID: 27642368

15. Ford JD, Cameron L, Rubis J, Maillet M, Nakashima D, Willox AC, et al. Including indigenous knowledge and experience in IPCC assessment reports. Nat Clim Chang. 2016; 6: 349–353. https://doi.org/10.1038/nclimate2954

16. Reyes-García V, García del Amo D, Benyei P, Junqueira, André B, Labeyrie V, Li X, et al. Local Indicators of Climate Change Impacts. Data collection protocol. Figshare; 2020. Available: https://doi.org/10.6084/m9.figshare.11513511.v3

17. Seddon AWR, Macias-Fauria M, Long PR, Benz D, Willis KJ. Sensitivity of global terrestrial ecosystems to climate variability. Nat 2016 5317593. 2016; 531: 229–232. https://doi.org/10.1038/nature16986 PMID: 26886790

18. Adger WN. Vulnerability. Glob Environ Chang. 2006; 16: 268–281. https://doi.org/10.1016/J.GLOENVCHA.2006.02.006

19. Reyes-García V, Guèze M, Díaz-Reviriego I, Duda R, Fernández-Llamazares Á, Gallois S, et al. The adaptive nature of culture. A cross-cultural analysis of the returns of local environmental knowledge in three indigenous societies. Curr Anthropol. 2016; 57: 761–784. https://doi.org/10.1086/689307 PMID: 28104924

20. Fernández-Llamazares Á, Díaz-Reviriego I, Luz AC, Cabeza M, Pyhälä A, Reyes-García V. Rapid ecosystem change challenges the adaptive capacity of local environmental knowledge. Glob Environ Chang. 2015;31. https://doi.org/10.1016/j.gloenvcha.2015.02.001 PMID: 26097291

21. Díaz-Reviriego I, Fernández-Llamazares Á, Salpeteur M, Howard PL, Reyes-García V. Gendered medicinal plant knowledge contributions to adaptive capacity and health sovereignty in Amazonia. Ambio. 2016; 45: 263–275. https://doi.org/10.1007/s13280-016-0826-1 PMID: 27878530

22. Beck HE, Zimmermann NE, McVicar TR, et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Scientific Data. 2018 Oct; 5:180214. https://doi.org/10.1038/sdata.2018.214 PMID: 30375988; PMCID: PMC6207062.

23. Bivand R, Lewin-Koh N. _maptools: Tools for Handling Spatial Objects_. R package version 1.1–4, 2022. https://CRAN.R-project.org/package=maptools

24. Bivand Roger S., Pebesma Edzer, Virgilio Gomez-Rubio, 2013. Applied spatial data analysis with R, Second edition. Springer, NY. https://asdar-book.org/

25. Colfer CJP, Brocklesby MA, Diaw C, Etuge P, Günter M, Harwell E, et al. The Grab bag: supplementary methods for assessing human well-being. Grab bag Suppl methods Assess Hum well-being. Bogor, Indonesia: Center for International Forestry Research (CIFOR); 1999. Report No.: 6. https://doi.org/10.17528/CIFOR/000766

26. Lynam T, Jong W De, Sheil D, Kusumanto T, Evans K. A review of tools for incorporating community knowledge, preferences, and values into decision making in natural resources management. Ecol Soc. 2007; 12: 5. Available: http://www.ecologyandsociety.org/vol12/iss1/art5/

27. Scoones I. Sustainable Rural Livelihoods: A Framework for Analysis. Brighton; 1998. Report No.: 72.

28. Glanz K, Bishop DB. The Role of Behavioral Science Theory in Development and Implementation of Public Health Interventions. http://dx.doi.org/101146/annurev.pubhealth012809103604. 2010;31: 399–418. https://doi.org/10.1146/ANNUREV.PUBLHEALTH.012809.103604

29. Harris I, Osborn TJ, Jones P, Lister D. Version 4 of the CRU TS monthly high-resolution gridded multi-variate climate dataset. Sci Data 2020 71. 2020; 7: 1–18. https://doi.org/10.1038/s41597-020-0453-3 PMID: 32246091

30. Pinheiro J, Bates D, R Core Team. Linear and Nonlinear Mixed Effects Models [R package nlme version 3.1–159]. Comprehensive R Archive Network (CRAN); 9 Aug 2022 [cited 16 Sep 2022]. Available: https://cran.r-project.org/package=nlme
31. Reyes-García V, Tofighi-Niaki A, Austin BJ, Benyei P, Danielsen F, Fernández-Llamazares Á, et al. Data Sovereignty in Community-Based Environmental Monitoring: Toward Equitable Environmental Data Governance. Bioscience. 2022; 72: 714–717. https://doi.org/10.1093/biosci/biac048 PMID: 35923191

32. IPCC. Climate Change 2014: Impacts, Adaptation, and Vulnerability Part B: Regional Aspects. Barros VR, Field C.B., Dokken D.J., Mastrandrea M.D., Mach K.J., MC T.E. Bilir, Ebi K.L., Estrada Y.O., Genova R.C., Girma B., Kissel E.S., SM A.N. Levy, Mastrandrea P.R. and LLW, editors. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge U.K.; New York: Cambridge University Press; 2014. https://doi.org/10.1007/s13398-014-0173-7.2

33. Miara MD, Negadi M, Tabak S, Bendif H, Dahmani W, Hammou MA, et al. Climate Change Impacts Can Be Differentially Perceived Across Time Scales: A Study Among the Tuareg of the Algerian Sahara. Geohealth. 2022; 6: e2022GH000620. https://doi.org/10.1029/2022GH000620 PMID: 36330077

34. Junqueira AB, Fernández-Llamazares Á, Torrents-Ticó M, Haira PL, Nasak JG, Burgas D, et al. Interactions between Climate Change and Infrastructure Projects in Changing Water Resources: An Ethnobiological Perspective from the Daasanach, Kenya. J Ethnobiol. 2021;41. https://doi.org/10.2993/0278-0771-41.3.331 PMID: 35692568

35. Danielsen F, Marlynn M. Mendoza, Anson Tagtag, Alviola PA, Balete DS, Jensen AE, et al. Increasing Conservation Management Action by Involving Local People in Natural Resource Monitoring. Ambio. 2007; 36: 566–570. https://doi.org/10.1579/0044-7447(2007)36[566:icmabi]2.0.co;2 PMID: 18074893

36. Turner NJ, Berkes F. Coming to Understanding: Developing Conservation through Incremental Learning in the Pacific Northwest. Hum Ecol 2006 344. 2006; 34: 495–513. https://doi.org/10.1007/S10745-006-9042-0

37. Blanch-Ramirez J, Calvet-Mir L, Aceituno-Mata L, Benyei P. Climate change in the Catalan Pyrenees intersects with socioeconomic factors to shape crop diversity and management. Agron Sustain Dev. 2022; 42. https://doi.org/10.1007/s13593-022-00806-3 PMID: 36059570

38. Reyes-García V, Graf L, Junqueira AB, Madrid C. Decarbonizing the academic sector: Lessons from an international research project. J Clean Prod. 2022; 368: 133174. https://doi.org/10.1016/J.JCLEPRO.2022.133174