Data for wetlandscapes and their changes around the world

. Geography and associated hydrological, hydroclimate and land use conditions and their changes determine


Introduction
Wetlands contribute more than 20 % of the total value of global ecosystem services (Costanza et al., 2014), while covering only a small percentage (4 %-9 %) of global land surface (Morganti et al., 2019;Zedler and Kercher, 2005;Mitsch and Gosselink, 2000). Wetlands are associated with a diverse range of functions, such as water quality remediation (e.g., Chalov et al., 2017;Quin et al., 2015), regulation of soil moisture and groundwater replenishment (e.g., Ameli and Creed, 2019;Golden et al., 2017), flood control (e.g., Quin and Destouni, 2018;Acreman and Holden, 2013), and biodiversity conservation (e.g., Cohen et al., 2016;Mitchell et al., 2008). Through these functions, wetlands can support regional sustainability  but are also one of the most vulnerable ecosystems globally . For instance, human land and/or water use developments (Destouni et al., 2013;Jaramillo and Destouni, 2015;Maneas et al., 2019) in combination with climate variability and change (Orth and Destouni, 2018;Seneviratne et al., 2006) affect large-scale water fluxes, with impacts on wetland functions and ecosystem services. These impacts extend over coupled systems of multiple wetlands and the associated total hydrological catchment that integrates these, referred to as a wetlandscape , with even well-intended actions towards various sustainable development goals potentially affecting wetland functions and services in different directions . As a consequence of various change impacts, wet-land areas are now suffering rapid and continued decline in different regions worldwide (Davidson et al., 2018;Davidson, 2014).
The scale mismatch between the existing large-scale studies of various landscape changes and the still mostly local wetland impact studies  creates an urgent need for comprehensive, science-based assessment of the interactions between large-scale drivers of change and large-scale wetland systems (Ameli and Creed, 2019;Creed et al., 2017). Adopting a wetlandscape perspective involves moving away from the individual wetland scale to consider the large-scale functioning of the hydrologically coupled system of multiple wetlands and their surrounding landscape. Assessments at these larger scales are needed to enable the formulation of scientific evidence-based guidance and strategies to protect wetlands under global change (Thorslund et al., 2018;Ameli and Creed, 2019). The conceptual framework on wetlandscapes was developed over 30 years ago, by Preston and Bedford (1988), but the dynamics and impacts of many large-scale drivers or functions at wetlandscape scales remain still largely uninvestigated and unknown, with the interactions between large-scale hydroclimatic variability and change and wetland dynamics still being largely underexplored at the wetlandscape scale . The combination of high wetland vulnerability and rapid largescale changes subject to major knowledge and data gaps highlights the need to synthesize and create datasets available for evaluating change effects and feedbacks at the scales of whole wetlandscapes.
To address this need and support large-scale studies of whole wetlandscapes in and across different parts of the world, we have created a novel database, named the Wetlandscape Change Information Database (WetCID), for 27 wetlandscapes around the world and their associated geographical, wetland, hydrology, hydroclimate and land-use conditions. WetCID consists of a survey-based collection of local information and data, combined with compilation and synthesis of gridded large-scale datasets for a range of relevant hydroclimatic and land-use variables.
The remainder of this paper is structured as follows: in Sect. 2, we describe the methodology used in collecting, processing and summarizing different datasets. In Sect. 3, we present WetCID summaries and sample figures and maps from different components of the underlying datasets in order to exemplify and highlight the potential of new insights that can be gained from using this database as well as its limitations. In Sect. 4, we discuss data availability and the format and structure of different files in WetCID. Based on the findings, we present some conclusions in Sect. 5.

Data acquisition
In compilation of WetCID for the 27 wetlandscapes, we employed three sources of primary data. These were (1) local site survey data, depicting general characteristics of each wetlandscape (catchment) and its geographical characteristics (including shapefiles for its spatial extent) and associated hydrological, climate and land-use conditions and their observed or perceived changes; (2) gridded historical data time series of monthly precipitation and temperature from Climate Research Unit Time-Series version 4.02 (CRU_TS4.02; Harris et al., 2014); and (3) historical data of annual land cover and its changes from the NOAA-HYDE dataset provided by the NOAA's National Climate Data Center (Jain et al., 2013;Meiyappan and Jain, 2012).
The survey for local site data was given to researchers within the Global Wetland Ecohydrological Network (GWEN) (http://www.gwennetwork.se/, last access: 9 May 2020). The GWEN researchers responding to the survey specified the relevant wetlandscape extent (total hydrological catchment with wetlands) and provided boundaries in GIS format for the 27 wetlandscapes, located as shown in Fig. 1. Information and data of all three types (local survey-based, hydroclimate, land use) were collected and synthesized for each of these wetlandscapes from all three sources (1-3). In addition to the local survey information, data on hydroclimate and land-use variables were thus also compiled from the global datasets in both gridded and aggregated form for each wetlandscape, as described further in the following.

Site information surveys
A questionnaire for collecting local site knowledge and information on the availability and accessibility of local data was developed during a GWEN workshop held in Santa Marta, Colombia, on 24-28 April 2018. The questionnaire was sent out by email after the workshop to all participating GWEN researchers. The researchers responding to it related their answers to a specific wetlandscape in which they had active research.
The questionnaire comprised two main parts. Part 1 contained general questions about the geography, climate, hydrology, and wetland-relevant human activities and changes in the wetlandscapes. Part 2 focused on the availability and accessibility of local site data, structured into "Hydroclimate", "Land use" and "Other" data (see templates in the database files for a full outline of the questionnaire). The collective knowledge obtained on conditions and changes in the 27 wetlandscapes and on data availability-accessibility is summarized in Sect. 3.1.
To complement this local knowledge and information basis, we further extracted and synthesized data for the 27 wetlandscapes from relevant global hydroclimate and land-use datasets as described below.

Hydroclimate data
The temperature and precipitation data taken from the CRU_TS4.02 global datasets (Harris et al., 2014) covered a 30-year period , to be consistent with the time span of existing global land-use change data. CRU_TS4.02 provides hydroclimate data with spatial resolution of 0.5 • × 0.5 • and at a monthly temporal scale. In preparing temperature and precipitation datasets for each wetlandscape, the gridded data within the area of the wetlandscape were extracted from the global datasets and also spatially aggregated over that area, based on area-weighted averaging over the grid cells covered by the shapefile of each wetlandscape (catchment). This provided wetlandscape-specific data time series for each variable at each grid cell and aggregated over the whole wetlandscape. To facilitate analyses at different spatial resolutions, both the gridded and the aggregated time series were included in WetCID for each of the 27 wetlandscapes.
In addition to the gridded and aggregated data time series, period-specific temperature and precipitation changes were also calculated for each wetlandscape by dividing the total 30-year time span of the collected data into the two 15year periods, 1981-1995(Per1) and 1996-2010. Such period-specific change quantification can facilitate relatively simple and straightforward analysis of how these hydroclimatic changes correlate with and may have driven other wetlandscape changes (e.g., in runoff, evapotranspiration, wetland area) between the same time periods (Destouni et al., 2013;Jaramillo andDestouni, 2014, 2015). Absolute and  Peel et al., 2007), with the number of wetlandscapes extended from those included in similar GWEN-site mapping by Thorslund et al. (2017). The site numbering is in order of latitude from north to south, covering a latitude range from 70 • N to 25 • S. relative (%) changes between these periods (AbsChng and RelChng, respectively) were calculated from the mean annual values of temperature and precipitation during Per1 and Per2 as where Var Per1 and Var Per2 are average temperature (in C • ) or precipitation (in mm yr −1 ) over Per1 (1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995) and Per2 (1996Per2 ( -2010, respectively. Equation (1) was applied to both temperature and precipitation data to calculate their absolute changes in each wetlandscape, while Eq.
(2) was only applied to precipitation data to calculate the corresponding percentage change in precipitation.

Land-use data
The NOAA-HYDE dataset was used to estimate land uses and their changes in each wetlandscape. NOAA-HYDE estimates annual changes in land-cover area over the global land mass, starting from a base map for the year 1765. The estimations follow a predefined pathway, determined by relevant land-use and management datasets (cropland, pastureland, urbanization, timber harvesting), to obtain forest area distributions close to satellite-based estimates of forests in recent years (Meiyappan and Jain, 2012). NOAA-HYDE data cover the period 1770-2010 with yearly temporal resolution and spatial resolution of 0.5 • × 0.5 • , from which data for the period 1981-2010 were used for the development of this database, consistent with the hydroclimate data period described above.
The NOAA-HYDE land-cover maps show the percentage of grid cell area containing 28 different land-cover types (see Table A1 in Appendix A). In this study, we reclassified these 28 land-cover types into 10 distinct land-cover types, namely urban, shrubland, grassland, pastureland, cropland, forest, water, desert, tundra and savannah, by combining similar land-cover classes (see Table A1). As done for the hydroclimate data, the gridded land-use data were also spatially aggregated over each wetlandscape based on the area-weighted averaging method (with weights of specific land-cover area in each grid cell relative to total wetlandscape area). This provided a wetlandscape-specific data time series of annual land use and land cover for each of the reclassified 10 landcover types. The final WetCID files comprised gridded time series data on absolute grid cell area (in km 2 ) covered by each land-cover type, time series data in percentage of grid cell area covered by each land-cover type, and aggregated absolute and percentage time series data for each wetlandscape area.
Analogous to the hydroclimatic changes, period-specific change quantification can facilitate relatively simple and straightforward analysis of how different types of land-use changes between time periods correlate with and may have driven associated wetlandscape changes (Destouni et al., 2013;Jaramillo and Destouni, 2015). Equation (1) was therefore also used to calculate absolute change in the area of each land-cover type (km 2 ) within each wetlandscape between Per1 (1981Per1 ( -1995 and Per2 (1996Per2 ( -2010. In the land-use case, Var Per1 and Var Per2 represent annual average area cov-ered by a land-cover type within each wetlandscape during Per1 and Per2, respectively. Furthermore, the corresponding change in relative land-cover area (ChngRel in % of total wetlandscape area) was calculated as where Area C is the total wetlandscape (catchment) area (in km 2 ) and Var Per1 and Var Per2 are the annual average areas covered by each land-cover type in the wetlandscape during Per1 and Per2, respectively. Table 1 summarizes some general geographical, climate and wetland type information provided by GWEN researchers in the survey information forms. Each site represents either an individual wetland or a wetlandscape (e.g., a catchment) including multiple wetlands. The country, main climate zone and wetland area relative to total wetlandscape (catchment) area are also given for each site in Table 1. Moreover, a summary of the availability-accessibility of local data on the wetlands, hydrology, climate and land uses as well as the wetlandscape (catchment) area in each of the 27 wetlandscapes is also shown in Fig. 2. The variables of evapotranspiration and soil moisture were revealed as having large data gaps (red color in Fig. 2), indicating an overall need to use other data sources (e.g., gridded global data products) for quantifying these variables and associated processes. The synthesized survey dataset also contains information about different types of wetland, hydroclimatic and/or landuse changes observed or perceived to have occurred in the 27 investigated wetlandscapes (Fig. 3). Substantial changes are reported for most of these wetlandscapes, but a few sites have no known changes (e.g., in the arid Moroccan site 17) or have important knowledge gaps regarding changes (e.g., in the cold Swedish sites 2 and 5, even though availability to at least some data is relatively good there). The information on local data availability-accessibility (Fig. 2) and observed or perceived change occurrence (Fig. 3) summarized and structured in WetCID can guide further study directions and support identification of key needs for complementary new local data and/or use of additional large-scale (regionalglobal) gridded data. Furthermore, the wetlandscapes of Wet-CID are located in different regions of the world, with seven sites in northern Europe (sites 1-7), seven in the Amazon and Caribbean region (sites 20 and 23-27), four in North America (sites 10, 15, 18 and 19), three in the Middle East (sites 12, 13 and 16), two in the Mediterranean region (sites 11 and 14), two in Siberia (sites 8 and 9), and two more in other parts of the world (northern Africa and eastern Asia). As such, regional patterns and characteristics can be identified and regional strategies developed, e.g., to enhance availability of data and information and determine further research needed to bridge region-specific knowledge gaps and decide on relevant management plans for each region's wetland ecosystems. Such regional characterizations and assessments can be initialized with the current version of WetCID and further updated as more data for already-included and possible additional regional wetlandscapes become available in future database versions.

Hydroclimatic data
Data for long-term average temperature and precipitation conditions, and changes in these between Per1 (1981Per1 ( -1995 and Per2 (1996Per2 ( -2010 at the 27 wetlandscape sites, are presented in Fig. 4. The horizontal axis in the diagrams shows the wetlandscape site numbers in order of their latitude from north to south, covering the latitude range from 70 • N to 25 • S. The increase in average temperature and precipitation with decreasing latitude (Fig. 4a-b) illustrates that the wetlandscapes also cover a wide range of hydroclimate conditions, from low to high temperature and precipitation values (see also Fig. 1). Temperature has increased over almost all wetlandscapes and considerably more so in the more northern and colder areas than in the warmer areas around and south of the Equator (Fig. 4a-b). In contrast, precipitation changes are relatively small, varying around zero, in the more northern, colder and drier areas, while precipitation has mostly increased in the warmer and also wetter areas around and south of the Equator (Fig. 4c-e). Overall, the changes in mean annual temperature range from zero to +1 • C, while the changes in precipitation range from −70 to +170 mm yr −1 , with the Iranian site (12) (Lake Urmia catchment) exhibiting the greatest increase in temperature (+1 • C) and the greatest relative decrease in precipitation (−17 %).
Figures 5 and 6 exemplify gridded variability and change data for temperature and precipitation over the Volga (no. 9) and the León-Atrato (no. 23) wetlandscapes. The data time series of wetlandscape-aggregated annual average temperature and precipitation in these wetlandscapes (Fig. 5) exemplify such data prepared and included in WetCID for all 27 wetlandscapes. These two wetlandscapes were chosen for data exemplification because they represent different hydroclimatic conditions, with Volga being cold and dry while León-Atrato is warm and wet (Fig. 5), as well as having different sizes, with Volga being the largest (1 360 000 km 2 ) and León-Atrato (2344 km 2 ) one of the smallest studied wetlandscapes. The data for these examples (Fig. 5) are consistent with corresponding data implications across the dif-    ferent wetlandscapes over the world (Fig. 4) in indicating an overall positive (warmer-wetter) spatial correlation between long-term average temperature and precipitation. Temporally, however, the recent changes in these variables imply a negative correlation (towards warmer and mostly drier conditions) for the Volga wetlandscape ( Fig. 6a and c) as for several other northern wetlandscapes in WetCID (Fig. 4). In contrast, a positive correlation (towards mostly warmer and wetter conditions) is implied by the recent temporal changes in the León-Atrato wetlandscape (Fig. 6b, d), as one of the most southern wetlandscapes in WetCID (Fig. 4). Such spatiotemporal sign shifts and dipole emergence in temperatureprecipitation correlations have been noted in other recent studies of long-term variations and short-term changes of hydroclimate over Europe (Charpentier Ljungqvist et al., 2019). WetCID can facilitate further studies of these correlation conditions for and across the different wetlandscapes around the world. The data for the Volga and León-Atrato examples also emphasize that wetlandscapes can have very different area extents (spatial scales), with potentially important implications for the spatial resolution (Fig. 6) and related usefulness of data provided in WetCID. For example, the Volga wetlandscape includes 982 grid cells with complete or par-tial coverage in the hydroclimate datasets, while the León-Atrato wetlandscape only includes four of such grid cells. Most of the available global datasets from climate and earth system models have coarser spatial resolution than the size of most individual wetlands. Thus, model data for individual wetlands are subject to high uncertainty, whereas data aggregated over whole wetlandscapes have greater potential for accuracy (Bring et al., 2015), highlighting the need for considering the whole-wetlandscape scales in assessments of how wetland systems interact with hydroclimate and landuse changes.

Land-use data
The aggregated and gridded land-use data in WetCID can also be used for different types of whole-wetlandscape analyses. Figure 7 summarizes the data for the long-term average relative area of each land-cover type (% of total wetlandscape area), and associated absolute area changes (km 2 ) and changes in relative area coverage (% of total wetlandscape area), for different land-cover types across the 27 wetlandscapes. The data reveal, for example, the high percentage of forest area in wetlandscapes at high latitudes and in the tropics, while relative cropland area increases towards the  temperate regions (Fig. 7a). Figure 7 also summarizes the different types of land-cover transformations, for example from "forest" into "cropland and pastureland" in the tropical Mekong wetlandscape (21), "pastureland" into "grassland" in the temperate Irish wetlandscape (7) and into "cropland" in the borderline cold-dry Iranian wetlandscape of the dramatically shrinking Lake Urmia (12) (Khazaei et al., 2019), "shrubland" into "cropland" in the borderline temperate Iranian wetlandscape (13), and "cropland" into "shrubland" in the warm temperate Greek wetlandscape (14).
The data time series of different land-cover types and their changes between Per1 (1981-1995) and Per2 (1996-2010 show, for example, forest and (decreasing) cropland, followed by pastureland and grassland, to be dominant in the large Volga wetlandscape, while forest, pastureland and (decreasing) cropland areas dominate the small León-Atrato wetlandscape (Fig. 8). Gridded maps of land-cover area changes in these wetlandscape examples (Figs. 9-10) again demonstrate large spatial-resolution differences, with potentially important implications for the usefulness of land-use datasets for wetlandscapes of smaller scales. For example, in the most northern Swedish-Arctic wetlandscape (1), grassland is obtained as the second most dominant land-cover type after tundra (Fig. 7a), which is not normally seen in this northern Arctic region.
These three summary documents (all in Excel) were created from responses obtained in the main survey of  GWEN researchers (see survey template and structure in WetCID files). Summary document A contains summarized site-specific information on the wetlands, hydrology, climate and land uses in each of for the 27 wetlandscapes. Summary documents B and C contain local knowledge relating to the availability-accessibility (or lack) of land-use and hydroclimatic data, respectively, for each of the 27 wetlandscapes.
-Folder 2: gridded land-use and hydroclimatic datasets (NetCDF database files). In WetCID, there is a separate NetCDF file for each wetlandscape that contains a complete set of gridded hydroclimate and land-use data time series for the closest rectangular window around the catchment polygon of the wetlandscape. The gridded hydroclimate datasets were created by subsetting the CRU_TS4.02 original global datasets over the area of each wetlandscape (catchment). The gridded landuse dataset for each wetlandscape (catchment) was created by first reclassifying the land-cover types and then subsetting the global gridded data. All these gridded data time series are saved in separate NetCDF files for each wetlandscape, which is an appropriate file type for storing gridded data. Each NetCDF file contains 18 variables, including hydroclimate, land-cover and some auxiliary variables. Appendix B presents the table of general attributes (Table B1) and information and explanations of all 18 variables included in the NetCDF database files (Table B2). Sample MATLAB and R codes for reading and extracting data from the NetCDF files are also provided in Appendix C.
-Folder 3: aggregated land-use and hydroclimate data (Excel databases). The time series of land-use and hydroclimatic data aggregated over each wetlandscape (catchment) were created from the gridded datasets (NetCDF files) and stored as Excel files for each wetlandscape. The Excel file for each wetlandscape contains three sheets: (1) annual time series of covered area by each land-cover type (in km 2 ), (2) time series of annual relative area (%) occupied by each land-cover type, and (3) time series of monthly temperature ( • C) and precipitation (mm per month) data.
To perform any spatial analysis of the wetlandscapes, one needs to have access to the shapefile and polygons of the wetlandscape (catchment) and wetlands within it. These shapefiles were provided by the GWEN researchers and can be downloaded from WetCID files.
The data for land-cover changes are provided for all considered land-use variables.

Conclusions
The presented new database (WetCID) combines surveybased local information and knowledge with gridded largescale hydroclimate and land-use datasets for 27 wetlandscapes around the world. The gridded datasets contain 30year time series of mean monthly precipitation and temperature along with annual average land uses and their changes over this time period for each wetlandscape. Wet-CID can support site assessments; cross-regional compar-isons; and scenario analyses of the roles and impacts of various land-use, hydroclimatic and wetland conditions and their changes on whole-wetlandscape functions and associated ecosystem services. The information on local data availability-accessibility and observed or perceived change occurrence summarized and structured in WetCID can guide further study directions and support identification of key needs for complementary new local data and/or use of additional regional-global gridded datasets. The gridded large-scale hydroclimatic and land-use data included in WetCID have been derived using open-access data sources and processed with open-source tools, while WetCID has been designed so that more data can readily be added to it. The site-specific usefulness of different included data varies for wetlandscapes of different scales, but WetCID can be updated with small time investment as new datasets become available or current datasets are expanded or refined.