An evaluation of freshwater monitoring programs in ILTER nodes and mountain national parks: identifying key variables to monitor global change effects

Identifying and quantifying global change impacts on biotic and abiotic components of ecosystems is critical to promote an effective adaptation that increases the success of conservation strategies. To achieve this goal, global and regional assessment efforts require certain degree of harmonization on local monitoring programs to establish relevant comparisons at different spatio-temporal scales. Otherwise, the lack of harmonization might hinder the detection and assessment on the effects of human impacts. In this work we have compiled information on freshwater monitoring programs located in areas of intensive research and conservation interest: International Long Term Ecological Research (ILTER) nodes and mountain National Parks. We aimed at evaluating the quality and robustness of these programs to assess the impact of global change, addressing from the worldwide to the European and Spanish national scale. Results highlighted that freshwater monitoring programs lack a common strategy to monitor these ecosystems. Even at the continental and national scales, contrasting strategies and level of detail have been historically applied. Water quality, habitat and biodiversity are more commonly monitored than community structure and ecosystem functioning. Monitoring efforts on the Spanish Mountain National parks indicated differences on the targeted aquatic ecosystems. Rivers and lakes received a higher attention, while mires were rarely considered. Our results provide evidence that greater efforts should be directed towards constructing a coordinated strategy to monitor freshwater ecosystems at national, continental, and global scales. This strategy should involve a shared backbone of biophysical and biogeochemical variables for each habitat type on agreed protocols that are implemented across regions and administrative borders. Achieving this will support a substantial advance on the ecological research to further delineate proper conservation strategies to face the challenges imposed by global change.


Introduction
At the end of the XX century a group of scientists set up the World Climate Research Programme (WCRP) with the aim of determining whether the world climate patterns were changing, and whether the causes were natural (e.g., glaciation-interglaciation periods) or aggravated as a consequence of human activities (see WCRP 1993). After decades of debate, the publication of the first part of the IPCC Sixth Assessment Report (IPCC 2021) has confirmed the unequivocal influence of human action on climate change, showing evidence of intensification and spreading across the planet. In addition, other programs, such as the International Human Dimensions Programme stated that human activity was also the cause of global changes on environmental conditions and on species distribution patterns, coining the term "global change" (Steffen et al. 2001). Global change effects pose a serious threat for the health and functioning of the biosphere (Hoffmann 2021;Malhi et al. 2020) and have important implications for human societies, including the availability of food, energy and water, the increase of environmental hazards, and the occurrence of new pathogens and diseases, among others (Steffen et al. 2004). Understanding how biotic and abiotic components of ecosystems respond to human related changes (e.g., new pollutants, land use and land cover change, climate change, dams, etc.) is paramount to find adaptation strategies for future scenarios (United Nations 2018) and design effective conservation strategies. To achieve this goal, monitoring programs need to be designed and implemented at different spatial and temporal scales to measure and detect abiotic and biotic changes (Parr et al. 2002;Wohner et al. 2019). Several proposals and programs have been developed worldwide for this purpose, with the most relevant being: • The International Long Term Ecological Research Network (ILTER): composed of hundreds of research sites (i.e., nodes) covering a wide array of ecosystems (marine, continental, freshwater and terrestrial). ILTER focuses on node-based long-term research and monitoring to identify and understand environmental changes across the planet. This network was established in 1993. ILTER nodes have increased drastically in the last decade, reaching 800 nodes in 44 countries in 2016 (Dirnböck et al. 2019). • The Geo Biodiversity Observation Network (GEO BON): a global social network and community created to coordinate efforts to design and implement national and regional biodiversity monitoring programs with the final aim of supporting decision-makers (Navarro et al. 2017). • ICP-Waters: The International Co-operative Programme on Assessment and Monitoring of the Effects of Air Pollution on Rivers and Lakes (Kvaeven et al. 2001): Established in 1985, ICP-Waters collects hydrochemichal and biological data from more than 500 sites in more than 20 countries across Europe and North America.
ILTER network could be considered as the largest and most powerful international initiative to understand environmental change across the globe attending to the broad variety of environmental variables that are regularly measured in most of the nodes (climatic, hydrological, biodiversity, socio-economic, etc.). Freshwater ecosystems are particularly vulnerable to global change (Woodward et al. 2010) and hence, they represent priority ecosystems that should be monitored. In this regard, this globally distributed network of long-term biomonitoring and research nodes monitors multiple ecological components that are not currently considered simultaneously in other freshwater monitoring programs (e.g., Water Framework Directive in Europe, the Clean Water Act in the U.S., the National Water Act in South Africa, the National Water Management Strategy of Australia and New Zealand or the Water Act in Canada). For instance, in the European context, the Water Framework Directive (WFD) has been criticized for the overrepresentation of monitoring programs focused on the effect of nutrient or organic pollution (Kelly et al. 2016), while other potential impacts related to global change, such as hydrological and morphological pressures, had much less attention (Carvalho et al. 2019). Moreover, ILTER network covers the monitoring schemes of freshwater ecosystems (lakes, lagoons, big rivers, small streams, wetlands, etc.) of a large variety of characteristics, while the ecosystem typologies covered by the WFD monitoring network is confined in such a way that many water systems are not included, e.g., some wetlands such as mires, river basins < 10 km 2 , or tributaries with less than 15 km length. Another strength of ILTER network is related to the metadata documentation system. In this regard, the Drupal Ecological Information Management System (DEIMS; see below) allows to easily search different information and data about nodes, staff and published results (Haase et al. 2016) although some limitations to data accessibility have been underlined (Ondei et al. 2018).
In the context of freshwater ecosystems, those located in mountain areas are recognized as very sensitive and vulnerable to changes in their environmental conditions (Beniston and Stoffel 2013). Mountain areas and their associated ecosystems have been recognized as "sentinels for global change", as their physical and biological features show patterns of changes that are often more readily identifiable than in any other geographical settings of the globe (Gobi et al. 2014;Zamora et al. 2016). On one hand, altitudinal gradients facilitate space for time climatic change understanding (Michalet et al. 2014) and, on the other, they often have a low human impact due to remoteness that facilitates recording trace impacts (Kollmair et al. 2005). Moreover, the altitudinal gradients characteristic of mountain areas influence other key variables (e.g., temperature or slope) creating unique environments which shape biological communities (Körner 2007). All these characteristics accentuate the fragility and susceptibility of mountain freshwater ecosystems to changes on temperature and precipitation regimes associated with global warming ) and changes in land use and land cover (Álvarez-Martínez et al. 2014;Paduel et al. 2016;Pérez-Silos et al. 2021).
On the socioeconomic dimension, mountain aquatic ecosystems are important for human societies as providers of essential natural resources and ecosystem services (energy, food, water, biodiversity, etc.;Barquín et al. 2015;Beniston and Stoffel 2013). As such, they have been historically considered conservation priority areas (Catalan et al. 2017) and about 17% of mountain areas worldwide have a protection figure (e.g., National Park; NP). This percentage supposes the 32% of the continental surface of the world´s protected areas (IUCN Mountain Specialist Group). Mountain systems are protected by environmental legislation worldwide, and different projects have been developed for monitoring the effects of global change on mountain areas at global (e.g., UNESCO GLOCHAMORE Project) and regional scales (e.g., The Eastern Rivers and Mountains Network in the USA; Save the Blue Heart of Europe in the Balkans region).
Despite the efforts carried out to generate common frameworks for the monitoring of freshwater ecosystems, there are still many gaps that need to be filled (e.g. Schmeller et al. 2017). In many cases, data compiled by monitoring programs do not comply with the minimum required design to properly evaluate the effects of global change (Wade 2006). Some of the shortcomings of these programs include: (1) the low periodicity of measurements; (2) the use of contrasting methodologies, and (3) the lack of a common set of essential variables/indicators to identify and quantify the effects of global change (e.g., climate or land use and land cover).
The main objective of the present study is to review and assess the potential of monitoring programs to evaluate the effects of global change on streams, lakes and mires at three spatial scales; (1) worldwide, (2) regional scale (Europe) and (3) local scale (Spain). At the worldwide scale we have analyzed the monitoring programs of 157 river and 150 lakes on the ILTER network; at the regional scale we have analyzed the monitoring programs provided by 28 European mountain NP; at the local scale we have analysed the monitoring programs in the 5 Spanish Iberian mountain NP. Beyond the general objective, this study aims to: (1) identify the set of biotic and abiotic variables measured on different freshwater ecosystems within different programs, and (2) to evaluate the quality of these programs to monitor global change effects according to a set of key variables.

Worldwide review
We used the DEIMS and dataset registry (DEIMS-SDR) to obtain a global list of ILTER nodes that are included in their Rivers (ILTER-Rivers) or Lakes (ILTER-Lakes) monitoring programs. DEIMS-SDR is a worldwide catalogue of in situ observation or experimentation facilities covering all biomes, allowing data associated with each site to be documented and accessible to the public and particularly to scientists (for more details see: Wohner et al. 2019; https://deims.org/). It should be noted this data hub does not contain raw data, but contains metadata associated to the type of variables monitored in each ILTER node. Moreover, the database is rather asymmetrical in the level of detail when describing sampling methods or specific variables monitored in each ILTER node.
The database was accessed on the 17 of February 2020 to retrieve information from sites that had an "active" status (i.e., the site is currently being operated) and "freshwater rivers" or "freshwater lakes" were annotated as their GEO-BON biome (Fig. 1). Then, we queried which variables were being monitored in each of the selected list of LTER sites. We filtered out variables that correspond to other biomes (e.g., forest or soil variables) except for atmospheric variables. Then, we classified the obtained list of variables on 5 types (atmospheric, water quality, physical habitat, biological and landscape) and on 32 and 34 subtypes for rivers and lakes, respectively (see Supplementary material Table S1). We only took into account for subsequent analyses the subtype level dropping those variables that were too general and not different to the type level (e.g., "ecosystem parameter" or "atmospheric parameter").

Monitoring programs at the european level
To assess the monitoring programs carried out on aquatic ecosystems located within European Mountain National Parks (EMNP), we visited the official websites of 51 European countries (including Turkey, Azerbaijan, Armenia and Georgia; Russia was not included) to list all the EMNP. We selected those with altitudes higher than 800 m a.s.l. Following these criteria we identified 113 EMNP, belonging to 27 different countries. Then, we contacted the selected EMNP via e-mail. Moreover, we also contacted different national organisms/ agencies responsible for the management of these EMNP (e.g., Environment Ministries), EUROPARC (www.europarc.org) and the Dinarides Parks Net (www.parksdinarides.org). E-mails were sent in October 2019, requesting information on the existing monitoring programs to assess the effect of global change on aquatic ecosystems. Specifically, we asked about the variables being measured, the periodicity of measurements and the type of aquatic ecosystem monitored (lake, river, mire, etc.).

Fig. 1
Worldwide and European distribution of active LTER sites that monitor river (black circles; 157) and lake ecosystems (red circles; 150) included in the DEIMS-SDR database

Working partnership in the Spanish Mountain National Parks
At the most detailed scale, we assessed active or recent monitoring programs carried out in aquatic ecosystem, including rivers, lakes and mires, located within the Iberian Spanish Mountain National Parks (SMNP) network. This network consists in five SMNP ( We first compiled the information (research and technical reports) developed in these SMNP (2000-2018 period) from the repositories available in the organism that coordinates the management of all the Spanish National Parks (Autonomous Organism of National Parks; Spanish Environmental Ministry). We identified 20 projects and programs dealing with monitoring programs in the five SMNP, which are lead by 22 principal researchers and 12 different institutions. Based on this information we designed a query form to contact via email with researchers and institutions involved in these projects and programs, obtaining an affirmative answer and interest in participating in this consulting process from 18 of them. We further obtained 10 answered online questionnaires where we required information about the spatial and temporal design of the monitoring programs previously identified. Attending to the reviewed bibliography and the questionnaires, we organized the workshop "Synthesis of the monitoring programs of global change in aquatic ecosystems SMNP" on the 4 of December 2017, held in the headquarters of the Autonomous Organism of National Parks in Madrid, Spain. During the workshop, 12 short presentations took place to present the initiatives currently established in different SMNP to monitor global change in aquatic ecosystems, as well as the most relevant results obtained to date. The meeting ended with a round table of synthesis and analysis of the presented results.
The review of this information allowed us to identify a set of 40 environmental and ecological variables (Table 1; Supplementary material Tables S2 to S5) that would be potentially useful to assess and monitor the effect of global change on aquatic ecosystems on SMNP. These 40 variables were classified into 6 different types according to the environmental/ecological attribute they addressed: Climate and Hydrology (10); Water Quality (7); Physical Habitat (9); composition (6) and structure (6) of biological communities and Ecosystem functioning (4). In addition, monitoring programs were evaluated individually attending to the targeted ecosystems: rivers, lakes and mires. We also considered monitoring programs focused on hydrometeorological variables because they are intrinsically related to freshwater ecosystems and global change.

Evaluation of Monitoring Programs
We generated an assessment system based on scores (0-4) for the different relevant characteristics of a monitoring program ( Table 2). We provided individual assessments for the monitoring programs in each one of the considered domains (hydrometeorology, rivers, lakes and mires), and taking into account the 5 types of variables described in Table 1 (except for cartography, see below). Specifically, we considered 5 relevant elements for each monitoring program as quality indicators: (1) Cartography; (2) Spatial design of the monitoring program; (3) Length of the historical series; (4) Frequency of the measurement and (5) Percentage of considered variables in relation to the maximum potential variables for each type. These quality elements have been largely highlighted as relevant to design good monitoring programs (e.g., Bartram and Balance 1986;Downes 2002). The cartography was evaluated based on the existence of cartography for river, lakes and mires in each SMNP, while the assessment was not segregated by the type of variable. The spatial monitoring design and the percentage of considered variables was evaluated for each type of variable in each domain. The duration and frequency of the measurements were evaluated independently for each variable (Supplementary material Tables S2 to S5) and the final assessment value was calculated as the average of the values of independent variables by type of variable).

Worldwide review
The LTER-River dataset included 157 nodes and 214 variables to monitor river ecosystems, while LTER-Lake dataset had 150 nodes and 204 variables (see Supplementary material, Table S1). All these monitoring variables could be lumped into 31 major subtypes and 5 major types (Atmospheric, Landscape, Water chemistry, Physical habitat and Biological; Table S1). Only 21% of the LTER-River and 18% of the LTER-Lake nodes measured variables simultaneously in the 5 major variable types, while an extra 25% for river and 20% for lakes measured variables in at least 4 of the 5 major variable types. Water temperature and depth where the most represented in monitoring programs in both, LTER-River and Lake (above 55% for river and 45% for lakes nodes; Fig. 3). Then, the most common variables Table 1 Variables and type of variables identified as key indicators of global change in monitoring programs of the five SMNP. We indicated the domains/types of ecosystems in which they can be potentially used; (X) is used to indicate groundwater dependence TYPE VARIABLES Hydro-meteorology Rivers Lakes Mires Climate and Hydrology Air temperature, pressure and humidity in LTER-River were water conductivity (55% of sites), suspended solids-turbidity (50%), nutrients (carbon: 50% or nitrogen: 48%) and pH (47%; Fig. 3 A). The biological variables most used in monitoring LTER-River where those related to the diversity of the biological communities (43%), the structure and composition of these communities (35%; mainly benthic macroinvertebrates, 10% of sites) and the water chlorophyll concentration (35%). LTER-Lake presented a slightly lower proportion of nodes in which water chemistry was monitored through conductivity (45%), suspended solids-turbidity (45%), nutrients (carbon: 41% or nitrogen: 37%) and pH (33%; Fig. 3B). Although diversity, community and alkalinity variables ranked a bit higher on importance in lakes than in rivers, their frequency of use in monitoring programs was similar to rivers (45%, 37% and 37%, respectively; Fig. 3).
Remarkably, all physical habitat variables (except water depth), phosphorous and oxygen concentration among the water quality variables, and genetic and ecosystem functioning among the biotic variables, are the most underrepresented on LTER-River and LTER-Lake (less than 22% of nodes). Finally, only benthic invertebrate communities and water chlorophyll a in rivers and lakes have a reasonable number of sites to construct regional cover assessments, while many other biotic attributes are only collected from ILTER nodes within Europe (see Figure S1 in the supplementary material). Finally, it is important to remark that the DEIMS-SDR dataset did not always include enough detail about the type of variables, organisms/community monitored or the sampling protocols, what made difficult to reach a higher level of detail in our analyses.

Monitoring european mountain aquatic ecosystems
We contacted 113 EMNP (Fig. 4) of which 28 (25%) responded. From these 28, 17 (61%) confirmed that, currently, they are not carrying out any monitoring program to assess the effect of global change on their aquatic ecosystems. The positive responses (11; 39%) were from ( Fig. 4): • Spain. 5 MNP with monitoring programs on different aquatic ecosystems. A detailed description of these programs is included in Sect. 3.3. • France. 4 MNP with a monitoring program on lakes (Sentinel Lakes Network; www. lacs-sentinelles.fr). This program analyses: water chemistry (2 times/year), chlorophyll concentration and the structure and composition of phytoplankton and zooplankton communities. In this case, the environmental conditions and biological communities of these 4 MNP are analyzed following the same protocol, developing a common data base and a standardized long-term monitoring program. • Switzerland. The only MNP in Switzerland has a complete monitoring program on rivers, lakes and reservoirs. This program includes the analysis of water quantity and quality, the biological communities (invertebrates, fishes, diatoms, etc.), and the physical elements of the ecosystems (e.g., substrate composition and mobilization). For a more detailed description, see Schanz et al. (2012). • Poland. The Bieszczady MNP has 4 automatic stations monitoring hydrometeorological variables (rainfall, temperature, wind, humidity, etc.; www.imgw.pl). This MNP also has a program of biological surveillance (diatoms, fishes and invertebrates).

Evaluation of monitoring programs at the Spanish Mountain National Parks
The type and quality of cartographic information describing the different type of ecosystems was very similar across the five SMNP. Cartography of rivers consisted of GIS databases representing the hydrographic network with associated environmental information (AEI) in the five SMNP. Differences in the river cartographic information were related to the amount and quality of this AEI. The most basic information included topographic variables (e.g., catchment area, reach slope, elevation), while the most complete AEI database provided geological, land uses and/or climatic information. Cartography of lakes was essentially composed of GIS shapefiles indicating the geographical position and lake perimeter, while in the case of mires, most of the cartographic information consisted only in digitalized maps indicating the location of some of the largest mires. The assessment of the hydrometeorological information also evidenced a good degree of agreement in the extent and quality of the monitoring programs across the SMNP. All meteorological monitoring programs incorporated air temperature, atmospheric pressure and precipitation, while SN, SG and AESM considered the complete set of potential parameters (e.g., additionally, solar radiation, wind speed, soil temperature and humidity, and snowfall; Table 3), although these variables were not included in all the weather stations comprising each network. None of the SMNP considered a specific spatial design to locate the meteorological stations while all of them had series longer than ten years with daily frequency. In the case of hydrology, we also observed similar results, where all the SMNP had daily measurements of discharge and water table height, but without a specific spatial design. It must be pointed out the exception represented by PE where a specific control-impact monitoring design was conducted to locate several flow gauges matching with biological monitoring sites (Table 3; see Álvarez-Cabria et al. 2019 for details). Only SN provided measurements of permafrost, but time-series showed low frequency and duration.
The assessment of the monitoring programs conducted in the 3 types of aquatic ecosystems (rivers, lakes and mires) evidenced that rivers were monitored more intensively than lakes and, especially, than mires. The results highlighted that in all the SMNP the monitoring programs considered water quality, physical habitat and the composition of biological communities (Table 4). Nonetheless, the number of variables covered and the quality of data varied between SMNP. In general, results highlighted that most of the monitoring programs of these three types of variables were designed following the reference condition approach, with the exception of PE (all variables) and SG (Biological communities), where a controlimpact design was established. In addition, most of the series had at least between 5 and 10 years and presented a yearly or higher sampling frequency (Table 4). Moreover, all of these monitoring programs considered more than 50% of the identified variables included in the water quality and physical habitat types. By contrast, we observed that a smaller effort has been made to monitor the composition and structure of communities, and especially ecosystem functioning in rivers ( Table 4).
Monitoring programs conducted in lakes were especially intense in SN, SG and AESM where the assessment of water quality, physical habitat, community composition and community structure were considered (Table 5). Specifically, monitoring programs in these 3 SMNP measured 75-100% of all the potential variables for water quality and physical habitat and 50-75% for community composition and community structure. Moreover, programs showed a measurement frequency of at least 1 campaign per year, and the time length of the data series was above 10 years in most of the cases. The monitoring design rarely considered a design strategy beyond the reference condition and in many cases, there was no reference to the spatial design. The monitoring program of AESM represented an exception to this rule, and a before after control impact design was implemented in a set of selected lakes. Monitoring ecosystem functioning was very infrequent in most of the SMNP except for SN and AESM. For instance, SN monitored primary production while AESM monitored primary, secondary and fish production (Table 5).
Finally, only 2 SMNP monitored mire ecosystems. In general, these programs were less complete and had lower quality than those in rivers and lakes. AESM had the most complete mire monitoring program. It included 3 types of variables (water quality, physical habitat and biological community composition) but only the 25-50%, 50-75% and 25-50% of all the potential variables were used within each variable type, respectively (Table 6). Mires monitoring programs have been usually deployed without a clear spatial design and in most cases they are based on the reference condition approach (physical habitat and community composition in AESM). Most of the time series were shorter than 10 years and the frequency does not reach one measure per year, with the exception of water table regime which is measured once a year ( Table 6).

Discussion
This study identified the variables most commonly used in different monitoring programs of freshwater ecosystems at different geographical scales. Looking on river and lake nodes of the ILTER network and on the different MNP considered at the European and National (Spain) spatial scales, we can assert that monitoring of water quality variables is the most widespread approach worldwide (e.g., pH, nutrients, conductivity, etc.). Regarding biological variables, river invertebrates and lake phytoplankton communities were the most common biological assemblages monitored. On the other hand, physical variables of aquatic ecosystems (except water depth) and those related to the adjacent terrestrial landscape are scarcely represented in these monitoring programs, reducing their reliability to evaluate the effect of global change on freshwater ecosystems. The detailed assessment of the monitoring programs conduced in SMNP revealed that their design should be reconsidered, while important discrepancies have also been observed in relation to the quality of the programs and targeted types of variables. This contrasted with the monitoring of meteorological and hydrological variables. In this case, most of the monitoring programs across SMNP notably agreed. Nonetheless, owing to the importance of snowfall and snow accumulation in mountains areas and the potential impact of global change, a larger effort should be done to capture the changes in snow-related variables. Finally, maintaining long term monitoring networks (minimum of 10 years) is a crucial prerequisite to understand the interannual variability and trends of aquatic ecosystems and the potential influence of global change. Thus, long term monitoring designed under a BACI or regional-scale trend analysis will allow a better analysis of the temporal evolution of the affected and control sites. This is paramount to discriminate if the observed variation is generated by the effect of an impact/alteration, or whether the observed changes are more related to the natural ecosystem variability (e.g. dry or wet years).

Worldwide review
The global scale review highlighted that the most widespread variables considered in river and lake monitoring were those related to water quality (mainly: water temperature, conductivity and suspended solids -turbidity) and water depth, while other physical variables were only reported from less than 20% of the ILTER nodes (e.g., water velocity, lake volume, or other physical habitat characteristics). The higher frequency use of water quality variables is not surprising given four main reasons: (1) they respond to changes in the aquatic environ-  to assess the divergence between natural and impacted locations (e.g., using reference conditions, CI or BACI designs). Surprisingly, biological variables were represented only in 40% of the LTER nodes, especially biodiversity and community variables and to a lesser degree, population variables (less than 30% of sites). These results corroborate the findings of Pereira et al. (2013) and Schmeller et al. (2017), who highlighted the lack of harmonized monitoring programs for delivering adequate data on biodiversity change at the global scale. Our results also agreed with Turak et al. (2016), who observed a generalized use of information on freshwater species populations and biological communities for measuring changes in global freshwater biodiversity. By contrast, we have observed that ecosystem functioning variables are the less common biological variables used to monitor changes in freshwater ecosystems.
Quite surprisingly, less than 50% and 40% of the river and lake ILTER nodes, respectively, included more than 3 groups of variables simultaneously. Changes on water or aquatic biodiversity conditions can be associated to environmental changes induced by natural, human or a combination of both elements. Consequently, in many cases the monitoring of water characteristics or biodiversity does not provide enough information to identify the "source of the problem". In this regard, producing hydro-meteorological (e.g., air temperature, discharge, surface water or groundwater levels) or landscape level information is key to better understand cause-effect relationships (Ervinia et al. 2019).

European Mountain National Parks
Only 4 countries confirmed the implementation of any kind of monitoring program in at least one of their MNP (Spain, France, Poland and Switzerland), while EMNP from the other 9 countries confirmed the absence of monitoring programs. Existing monitoring programs from the 4 EU countries presented contrasting approaches. For example, France assessed several EMNP with the same methods and periodicity. However, the French program only monitored mountain lakes and did not include springs, streams or mires. By contrast, Spanish monitoring programs (SMNP) included rivers, lakes, and mires but they were not monitored consistently across all MNPs, showing different intensities and approaches across locations. These differences complicate comparisons of datasets and analyses of long term trends (see Sect. 4.3 for further details). Differences between both countries, Spain and France, could be related to the politic administration of NP. In France all the NP (9) are managed by the public body Parcs Nationaux de France, while in Spain each NP (15) is managed by the regional Government where the Park is located, and therefore, the monitoring programs might differ according to the regional administration. The case of Switzerland is different, with only one MNP (0,4% of country area) concentrating all the monitoring and conservation nationwide efforts. In this regard, Switzerland has implemented a complete monitoring program in the running and standing waters of this MNP. This program takes information from 2 catchments with different water sources (glacier-fed and precipitationfed) from 2001 to present, including measurements and analysis on water physico-chemical characteristics, biological community structure and composition, atmospheric and weather conditions and variables related to landscape (see Schanz et al. 2012). All this information, compiled in a long-term database, constitutes one of the best in Europe to identify consequences of global change in mountain areas.

The Spanish Mountain National Parks
Monitoring the effects of global change in mountain aquatic ecosystems requires accurate cartographic information (Barquín et al. 2015). The assessment of the available cartographic databases for the five SMNP indicated that there was only sufficient information for rivers, while these databases seemed incomplete for lakes and mires. Moreover, even for rivers, it must be pointed out that a digital hydrographic network derived from a digital elevation model is not sufficient to establish appropriate spatial designs for monitoring. In this case, specific catchment environmental information from different digital sources should be appended to the digital river network. In addition, González-Ferreras and Barquín (2017) evidenced that available river networks still lack a complete representation of the spatial distribution of river channels. In this regard, improving mapping tools and coupling this digital cartography with analytical information is well needed, not only to understand the interactions between large scale environmental modifications and local alterations, but also to avoid significant errors in the water budget calculations and resulting flows (Egüen et al. 2012). Nowadays, global data sets and maps at reasonable resolution are available worldwide (e.g., HydroSHEDS, HydroRIVERS and HydroATLAS (Linke et al. 2019) and Hydro-LAKES (Messager et al. 2016) for rivers and lakes, respectively. The use of these databases would allow initially filling some information gaps until more specific cartographic products could be developed (but see Benda et al. 2016).
This study highlights the lack of a common national strategy to monitor the biological and ecological impacts in mountain rivers, lakes and mires in Spain. This uncoordinated strategy is prone to generating weaker conclusions in comparison to the possibility of generating a multiple and simultaneous program with a diversity of well-structured datasets (Tydecks et al. 2019). In addition, we found a general weakness related to the lack of adequate spatio-temporal designs, which would be needed to statistically discern changes that can be attributed to global change, natural variability or others (Downes 2002). In most cases the best design does not go beyond the reference condition, while Control-Impact and Before-After/Control-Impact designs would help to discriminate changes caused by global change to those caused by natural variability .
Most of the installed meteorological stations in the 5 SMNP (14) measure basic climatic variables such as temperature, precipitation, pressure or air humidity. However, the small number of stations measuring snow-related variables contrasts with the shifts in mountain snowpacks and glacier retreats that have been observed in the last decades in mountain areas (Kneib et al. 2020). Changes on the timing of snow-melting phenology will also affect runoff and the hydrological regime (Beniston and Stoffel 2013), which might further influence mountain freshwater ecosystems. In this context, additional monitoring efforts are paramount to provide insight into the impacts of the snowpack dynamics on mountain freshwater ecosystems. A remarkable example is the Global Monitoring System for Snow in SN (www.uco.es/dfh/snowmed;Polo et al. 2019), with additional instrumentation beyond the standard weather stations to measure snowfall, snow cover area, and specific components of the energy balance.
Most river and lake monitoring programs in the five SMNP provided water quality and physical habitat information, and biological data rely mostly on invertebrates community composition (Supplementary material Table S3), ignoring, in other cases, community structure. In this regard, early warning signals on biological communities because of hydrological or thermal regimen changes might be more easily detected as changes in community structure (Carlisle and Hawkins 1998;Whiterod et al. 2015) than as sudden changes such as species replacement. Moreover, the SMNP monitoring programs revealed an underrepresentation of ecological functioning variables. Functioning indicators present some advantages over structural ecosystem components since change in ecosystem process rates can be immediately linked to certain levels of environmental change (Palmer and Ruhi, 2019). For instance, the proliferation of reliable and relatively inexpensive sensors for monitoring dissolved oxygen, and easy-to-use software for the calculation of gross primary production (GPP) and Ecosystem Respiration (ER; Engel et al. 2019) is underpinning a revolution in the use of this kind of indicator (Bernhardt et al. 2018;Smits et al. 2021). These two factors enable scientists to estimate ecosystem metabolism more frequently and in greater detail.
Our study also reveals large differences in the extent and quality of the monitoring programs between rivers, lakes and mires. Results highlighted that monitoring programs in mires were the poorest within the five SMNP. By contrast, our results evidenced a much larger effort to monitoring lakes and rivers than mires in the five SMNP. The rich monitoring networks on lakes existing in three of the SMNP should be remarked. Although mountain lakes are exceptionally susceptible to the effects of global change (Thompson et al. 2009;Alcocer et al. 2021;Zamora and Oliva, 2022s), they appeared underrepresented in monitoring networks many times, given the large amount of resources and effort needed to survey remote lakes (Schaeffer et al. 2013). Moreover, the resources needed for monitoring continuously remote mountain lakes might prevent maintaining the required temporal and spatial resolution of relevant phenomena to capture the effects of global change. On the other hand, there are now larger chances to overcome these drawbacks by Remote Sensing (RS) techniques, as the potentiality of detecting key water quality-, biota-, and functioningrelated indicators in lakes is increasing substantially (Dörnhöfer and Oppelt 2016) or by using in-situ sensing systems (Porter et al. 2009). New opportunities raised from the use of RS and novel optical sensors for lakes, might also include the improvement and extend of the monitoring for mires. However, several characteristics of mire ecosystems might challenge the use of these technologies for their surveillance (Gallant 2015).

Recommendations for regional and global assessments
Lessons learnt from the different analysis carried out in this study could help illustrating how the harmonization of local monitoring programs could be improved for comparisons at regional or global scales. This harmonization does not imply that freshwater monitoring programs must converge on a unique monitoring design. The complexity of biodiversity and the idiosyncrasy of each geographical and socioeconomic context make this task unfeasible. Nonetheless, the international scientific community, in coordination with public institutions, should consider a greater effort to agree on a core of ecosystem characteristics and essential variables to be effective in understanding the effects of global change in freshwater ecosystems. In this regard, the comparison of the ILTER nodes information has actually shown that most of the metadata associated to the monitoring programs allows only for the identification of high-level information. Thus, moving beyond the state of the art is rather difficult nowadays. In this regard, we advise that freshwater monitoring programs substantially improve their metadata (e.g. type of community surveyed, habitat sampled, methods and equipment, site selection criteria, survey frequency, date/period of surveys, etc.) so as to find coherences and inconsistencies when trying to build regional or global assessments.
On the other hand, the analysis of monitoring programs at the most local scale and detailed information source (i.e. the SMNP) has identified a key set of variables that could be used as a start to assess the effects of global change in freshwater ecosystems (see Table 1). Thus, one of the first steps towards harmonization is to find a consensus on which of these variables should be consistently measured across aquatic ecosystems at local scales and which might be the targeted biological communities (a bottom-up approach within regions). On a second round, the methodological approach should be clearly defined so as to find a common agreement on the core steps and on how sampling design and temporal frequency should be defined, establishing the minimum conditions for each variable to be comparable across localities. On a third step, the design of the monitoring program (BACI, trend analyses, reference condition approach, etc.) might have importance at the local level (depending on objectives and resources), but is less relevant for regional or global assessments as long as natural and impacted sites are both monitored. For example, a control-impact design with summer surveys on an annual basis using electrofishing and identifying fishes to the lowest possible taxonomic level could be a good start for monitoring trends in fish communities in rivers. This will help to reach a minimum consensus within a region so that the obtained data will be comparable across localities. On a final step, this information could actually be used to inform Essential Biodiversity Variables that summarise the state or trend of the targeted populations, communities or ecosystems for a given region (Pereria et al. 2013). Although, simple in theory these steps have proven to be rather difficult to follow in practice. In this regard, mechanisms for integration of data bases provided for each protected area should be improved so that they can be part of national and international data repositories (e.g., Freshwater Information Platform; Schmidt-Kloiber et al. 2019). Moreover, these efforts should come along with new funding schemes that improve and allow a sustained scientific collaboration at regional and global levels.

Conclusions and future improvements
This is one of the first attempts at accessing the quality of monitoring programs for measuring the effects of global change in freshwater ecosystems. There has been notable progress towards improving monitoring programs in the last decades, however there is still a lack of agreement on what, how and when to measure in order to effectively oversee the effect of global change in freshwater ecosystems. According to our initial expectations, we have found important contrasts among monitoring programs worldwide and, surprisingly, also within those enclosed under the same protection figure in a given country (i.e., SMNP).
Our analysis at the national Spanish scale has evidenced important shortcomings in the monitoring programs that should be solved to increase our understanding on the complex effects of global change and to provide with a stronger body of evidence on current trends and directions at a national level. This evidence is paramount to further design policies and concentrate efforts to design effective adaptation and mitigation actions (e.g., development of blue and green infrastructures networks). In this sense, we have found that specific monitoring designs that allow the application of well-established statistical approaches (e.g., control-impact), grounded in a complete, sound and environmentally-attributed cartography is essential to comply with these objectives. In addition, inclusion of new biological and functional indicators that capture the progressive rates of global change, rather than abrupt changes, are needed within monitoring programs.
All these improvements necessarily lie on a stronger collaboration and coordination between government agencies, research institutions, international scientific initiatives and other stakeholders to maximize efforts in supporting coordinated and joint assessments of freshwater ecosystems. Improving and harmonizing monitoring is paramount to understand the processes, asses the impacts and face the challenges imposed by global change in freshwater ecosystems. This requires the implementation of long-term data bases from different sites, like proposed by ILTER, the standardization of the data collection in the field, like proposed by GEO BON (e.g., Essential Biodiversity Variables; Pereira et al. 2013;Schmeller et al. 2017), as well as an appropriate periodicity of measurements (daily, monthly, annual, etc.; Haase et al. 2018) and a proper statistical design (Downes 2002). In addition, only through a harmonization process built with the consensus of the localities involved in monitoring aquatic ecosystems, can sufficient coherence and consistency be obtained to build regional or global aquatic biodiversity assessments. This effort must therefore grow from the local scale and gain in regional or global vision as new challenges are incorporated to reach higher scale assessments.
Finally, it is remarkable the advances of recent techniques and tools to further improve the monitoring and assessment of ecosystems. This is the case of RS and the programs and services developed for the use/application of this information, like the Copernicus Land Monitoring Service in Europe (https://www.copernicus.eu/en) joined to the open data cube tools needed to manage this information (Giuliani et al. 2019). Remote sensing has long been recognized as having the potential to complement conventional approaches to freshwater monitoring (Bukata 2013), increasingly being used as a complementary source of information to in situ monitoring networks. Regarding in situ monitoring, a continuous or quasi-continuous monitoring of different limnological parameters is possible using autonomous multiprobes and temperature dataloggers. These devices have large data storage capacity, high battery life and sensors which are capable of working unattended from weeks to months (Granados et al. 2020). On the other hand, there has been an increased interest in the use of new molecular techniques to monitor and manage ecosystems (Erickson et al. 2019), which complement traditional approaches to define the structure and composition of biological communities. Even though our results showed an underrepresentation of genetic variables in the ILTER network and they are rarely included as an objective in biodiversity monitoring programs (Turak et al. 2016), the rapid evolution of these techniques (e.g., environmental DNA), could produce a radical change in this tendency (Rees et al. 2014;Pawlosky et al. 2018). Funding This study was carried out with the financial support of the EU LIFE programme, as part of the project LIFE-DIVAQUA (Proyecto LIFE18 NAT/ES/000121). The information showed in this publication only reflects the point of view of the authors. The European Comission and CINEA cannot be held responsible for any use which be made of the information contained here. Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.

Data Availability
The datasets used to review the monitoring programs at the Global scale conducted in the ILTER nodes that are included in Rivers (ILTER-Rivers) or Lakes (ILTER-Lakes) are available at the DEIMS dataset registry (DEIMS-SDR) (for more details see: Wohner et al. 2019; https://deims.org/). The datasets used to review the monitoring programs at the European scale conducted in the EMNP are available by requesting them via email to the corresponding EMNP. The datasets used to review the monitoring programs at the National scale conducted in the five SMNP are mostly included in this article and it supplementary information file. In addition, meteorological data recorded by the Spanish Global Change Monitoring Network in available at the Meteorological data download web application (https://www.miteco.gob.es/es/redparques-nacionales/red-seguimiento/acceso-datos.aspx). Further and more detailed information is available from the corresponding author on reasonable request.

Declarations
Conflict of interest All authors declare that they have no conflicts of interest and no relevant financial or nonfinancial interests to declare that are relevant to the content of this article.

Consents
All authors consent to participate in this work. All authors agreed with the content and consent to publish the research presented in this work. Unfortunately, the author Carmen Zamora Muñoz died on March 2022.
Other statements Animal Research, Plant Reproducibility and Clinical Trials Registration are not applicable in the research presented in this work.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.