Assessment of the potential responses of ecosystem services to anthropogenic threats in the Eten wetland, Peru

ABSTRACT Introduction: Coastal wetlands are currently threatened by human drivers, such as agriculture, infrastructure development, and urban sprawl. Pressures on these ecosystems disturb their morphology and biogeochemical cycles, resulting in the degradation of ecosystem services. However, little has been done to understand the wetland response to identify proper conservation strategies. Along the Peruvian coast, wetlands present a diversity of landscapes that face similar threats and pressures; however, the ecosystem response in each one may be different.Objective: This study aims to assess the environmental impacts on ecosystem services based on the understanding of geomorphic features and the status of the Eten coastal wetland (Peru).Methods: The methodology combines the application of open-source GIS tools and the collection of field data to characterize the geomorphic settings and to analyze the changes in environmental parameters. Then, the main threats and pressures on the Eten wetland are defined and related to impacts on ecosystem services using a cause-effect model.Results: Overall, the results indicated that the river plays a vital role in defining the wetland landscape and functions. The biological diversity of aquatic habitats is disturbed by hydraulic structures and agricultural activities. Current land use affects supporting and regulation services, such as water regulation.


Introduction
Worldwide, coastal wetlands are considered as one of the most valuable ecosystems due to their provision of important benefits to humans, such as carbon sequestration, storm protection, flood control, and waste processing (Bai et al. 2019;Camacho-Valdez et al. 2013). However, the quality of these wetlands is declining globally because of inadequate knowledge to manage changes in ecosystem structure and functionality of services and goods (Barbier 2013;Li et al. 2014). Over recent years, due to anthropogenic impacts (e.g., livestock, agriculture, water abstractions, urban expansion, and climate change), a widespread loss and ecological degradation of wet areas has occurred (Darrah et al. 2019;Rojas et al. 2019). Currently, changes in climate and anthropogenic activities alter the hydrology, biogeochemical cycles, and other processes that sustain wetlands and their ability of adaptation and resilience to external drivers (Burns et al. 2012;Elliott et al. 2010;Mitsch and Gossilink 2000;Morris et al. 2002;Wagner and Breil 2013). Consequently, the degradation of coastal wetlands results in the loss of habitats, changes in water flow and the increase in pollutants that reach the sea (Stedman and Dahl 2008). However, little is known about the past processes and landform changes in wetlands that control their development, and characteristics (Tooth and McCarthy 2007).
The Ramsar Convention designed a set of inventory sheets to prioritize areas for conservation around the world. Peru has only 13 RAMSAR sites, but 92 wetlands have been identified only along the Peruvian coastline (ProNaturaleza 2010). These ecosystems provide a variety of services and benefits, such as biodiversity, water regulation, carbon sequestration, storm protection, and cultural values. For instance, Peruvian wetlands are part of the Pacific Flyway of waterbirds (Perillo et al. 2009), playing an important role as stopover and breeding sites during birds' migration (Martini et al. 2009). Herein, the Humboldt current recirculates plankton and nutrient-rich waters along the coastline (Bakun and Weeks 2008), causing coastal areas to become the hotspots of shorebirds due to their high productivity (Butler, Davidson, and Morrison 2001). Even though the Peruvian Ministry of Environment (MINAM, in Spanish) created the National Wetland Strategy to foster the conservation and sustainable use of wetlands through the prevention, reduction, and mitigation of ecosystem degradation, currently many of the Peruvian coastal wetlands do not have any legal protection or tailored-conservation management plans. Previous research has been done to only focus on inventories and conservation status without considering an integrated landscape assessment (Angulo-Pratolongo, Schulenberg, and Puse-Fernández 2010; Aponte and Cano 2013;INRENA, ProNaturaleza, and WWF 1998;Pulido Capurro and Bermúdez Díaz 2018).
Nowadays, Peruvian coastal wetlands face many threats such water and soil pollution, agricultural development, overexploitation of resources, and urbanization sprawl. For instance, the Ramsar-site mangroves of Tumbes, located in the northwest of Peru and supporting over 120 birds and 93 fish species, is currently threatened by the construction of extensive aquaculture ponds (e.g., prawn and shrimp farming) and overharvesting of timber (Spalding, Kainuma, and Collins 2010). In the case of the Santa Rosa wetland on the Peruvian central coast north of Lima, extensive water (from surface and ground) withdrawal for irrigation purposes, livestock keeping, and wastewater release are leading to the degradation of the wetland and the loss of diversity (Aponte Ubillús and Ramírez Huaroto 2011;Ramirez, Aponte, and Cano 2010;Verones et al. 2012). The Ramsar-site Pantanos de Villa, located on the south-central coast of Lima, has also been affected by fires due to open burning of trash and campfires (Wyld, Audisio, and Poston 2015), reduced water availability, eutrophication, as well as urban sprawl and informal land management (Aponte, Ramirez, and Lertora 2018).
Located in northern Peru, the Eten wetland presents a diverse landscape composed by a variety of aquatic and terrestrial habitats which are defined by the fresh and brackish waters´ dynamics, that also supports the bird diversity in the region (Angulo-Pratolongo, Schulenberg, and Puse-Fernández 2010;Tabilo et al. 2016;Zavalaga et al. 2009). However, little has been done to protect this unique environment and the ecosystem services that it provides. This study develops an approach to assess the environmental impacts on ecosystem services based on the understanding of geomorphic features and current status of the Eten coastal wetland.
The goal of this research is to demonstrate that the characterization of both the main functions and the current conditions of the ecosystem can inform about the drivers and impacts. As a result, a tool for decisionmaking that provides a basic analysis of the impacts and the status of the wetland´s ecosystem is developed. This tool aims to be versatile and provide a rapid diagnosis and data for decision support for national and international organizations working on wetland restauration, management, or conservation. In the Peruvian region, prospective institutions that can use this tool include the MINAM, the Peruvian National Protected Areas Service (SERNANP), the Peruvian National Water Authority (ANA), and local governments.

Study area
This study focuses on Eten wetland, located on the Northern Peruvian coastline ( Figure 1A), in the lowland of the Chancay-Lambayeque basin. The Chancay-Lambayeque basin is a bi-regional watershed and is considered one of the main Peruvian Pacific catchments because of its large area (5,555 km 2 ) and extensive agricultural development (e.g., maize, beet, squash, and onion). Therefore, the basin supports many river-related infrastructures, such as reservoirs, intakes, tunnels, and irrigation channels (Autoridad Nacional del Agua 2010; 2011). As illustrated in Figure 1A, the Eten wetland is limited by the Reque river and the Pacific Ocean. Here, the geographic setting of the wetland is supported by the connectivity between the river and the ocean, creating an estuarine environment.
Seasonally connected to the ocean, the Reque river supports the agricultural productivity of the region. As in many coastal areas in Peru, agriculture represents the main economic activity in the study area, followed by livestock farming. Water for irrigation is supplied by the Reque river and groundwater wells. Moreover, the Eten wetland constitutes the main source of raw material for local people. As a heritage handed down from the Moche civilization (about 100 to 700 AD), local villagers elaborate handicrafts using local materials such as clay and vegetable fiber, particularly from the totora (Typha angustifolia), carrizo (Bacharis saliccifolia), junco (Scirpus limensis), and saltmarsh bulrush (Scirpus maritimus) (Angelakis et al. 2012;Aponte 2017). Elaborated products are bags, straw hats, purses, rugs, paintings, musical instruments, ceramics, jewelry, and textiles (Ministerio de Comercio Exterior y Turismo 2009). Furthermore, the Eten wetland is part of the birdwatching route "Northern Route for Bird Observation," an initiative created to foster ecotourism in Peru (Williams, Coppin, and Álvarez Alonso 2005). The current human demand for this wetland makes this an ideal study area for evaluating the impacts on ecosystem services on Peruvian coastal wetlands.

Methods
This section provides an outline of the research methods used to assess environmental impacts in Eten coastal wetland based on a hydrogeomorphic and environmental characterization ( Figure 2). First, a hydro-geomorphic classification was applied to categorize the Eten wetland under the Brinson (1993) methodology and, according to the class, the expected ecosystem services were identified following the meth-odology described by Malekmohammadi and Jahanishakib (2017). Second, an environmental characterization was performed using open-source tools and field measurements to determine the status of wetland components. Third, a DPSIR (Drivers, Pressures, States, Impacts, and Responses) was applied to assess possible impacts and responses in ecosystem services (Kristensen 2004;Malekmohammadi and Jahanishakib 2017;Müller and Burkhard 2012). The results of this assessment provide information to analyze cause-effect relationships between environmental changes and impacts on ecosystem services. (1) Hydro-geomorphic classification and identification of expected ecosystem services and human benefits The hydro-geomorphic approach assesses the functionality of a wetland ecosystem by analyzing the physical, chemical, and biological interactions of the structural components of the ecosystem and the surrounding landscape (Engineer Research and Development Center 2010). According to Brinson (1993), this characterization is based on: 1) the geomorphic setting or wetland geomorphology, 2) sources of water, and 3) hydrodynamics. According to these three factors, wetlands can be classified as: Depressional, Riverine, Flats (mineral and organic), Fringe (tidal/estuarine and lacustrine), and Slope. In this study, this approach is applied to the Peruvian coastal wetland Eten. Following the methodology described by Malekmohammadi and Jahanishakib (2017), the Eten wetland ecosystem functions and services were then identified based on previous research, literature review, and the guidelines described by the International Union for Conservation of Nature -IUCN (Dugan 1992) and the Millennium Ecosystem Assessment (MEA 2005) due to the lack of local data.

River morphology
The characterization of the Reque river morphology was performed to describe the lateral connectivity with the Eten wetland. The planform dynamics was described using past and recent aerial photographs, GIS-based tools, such as QGIS, and open-source satellite imagery of Landsat and Google Earth Pro. Past aerial photographs of 20-cm resolution from 1949 were acquired from the Peruvian National Aerial Photography Service, and then digitalized and georeferenced. To characterize the current landscape, photogrammetry was also applied to record recent aerial imagery from August 2019. These data were processed and georeferenced to obtain a high-resolution imagery of 20-cm resolution of sections along the Reque river ( Figure 1, items B, C, and D), and a Digital Elevation Model (DEM). The 2019 river morphology was described and characterized along the stream using variables such as sinuosity and bar patterns. Imagery from Google Earth Pro was also used to define agricultural, urban, and dune areas. Moreover, a hillshade of azimuth 45° was obtained from the DEM to characterize geomorphic characteristics such as paleochannels or "footprints" of past dynamics.

Land Cover and Land Use analysis
Landsat satellite imagery of 5-m resolution of , 2016, and aerial photographs for the 1949s environment were processed to assess the changes in land use within the Eten wetland. To compare changes in land cover with respect to agricultural and urban development, three classes were used: water, agricultural, and built-up areas. Built-up or urban area includes all type of buildings and roads; agricultural area consists of crops and agricultural lands; and water bodies comprise oxbow lakes, ponds, and the river. Therefore, croplands were manually digitalized per month and then grouped per year. Both croplands and urban areas were considered only within the Eten estuary region.

Sediment grain size and composition
Sediment material was sampled between July 2017 and August 2019 from different sites along the Reque River and within the Eten wetland (Figure 1.A) to characterize the sedimentary environment in the aquatic habitats, and to evaluate changes in sediment grain size distribution along the stream. The analysis also included the identification of the environmental conditions in which sediments were deposited. Sampling started from the Reque town down to the river mouth in which cores of 40-cm depth and surface sediment samples of 500 gr were obtained from (see Figure 1A): shoreline (8 samples), cropland (15 samples), marshes (15 samples) riverbanks (19 samples), and river bedforms (5 samples).
In the laboratory, samples were mixed and quartered to uniformly distribute the sediment material. Quartered samples were weighted, and then they were dried at a temperature of 110°C ± 5°C for about 12 hours. Then, dried samples were also weighted and washed using a #230 sieve to remove sediment less than 200 μm size. Samples were one more time dried, nested, and agitated to obtain the weight of passing and retained material to finally plot the particle size distribution and material composition.

Physical-chemical parameter of water quality
In July 2017, water samples were collected once from different locations along the Reque River and within the Eten wetland to assess the influence of freshwater inputs in the wetland and the influence of agriculture and urban sprawl on water quality. Three samples were collected along the Reque river from Cascajales town to the river delta, and three samples from lakes within the Eten wetland ( Figure 1A), following the Environmental Protection Agency (EPA) guidelines. Physical-chemical parameters included pH, Dissolved Oxygen (DO), Chemical Oxygen Demand (COD), Biochemical Oxygen Demand (BOD) conductivity, and temperature.

Aquatic biota: estuarine macroinvertebrates assemblages
In early August 2019, benthic macroinvertebrates individuals were sampled along the banks of the Reque river to evaluate the diversity of aquatic biota. There were 16 sample stations from the Reque town (upstream, yellow symbols) to the river mouth (downstream, orange symbols) as illustrated in Figure 1 items F and G. Based on the relative location of the flood protection structure, stations were defined to evaluate the influence of hydraulic infrastructure in aquatic biota.
At each sampling station, in an area of 1 m 2 , bottom substrate was removed and filtered using a kick net of a mesh diameter of 300 μm, and then filtered into sieves of an aperture of less than 425 μm to pass through the coarse sediment. Each macroinvertebrate individual collected was fixed in situ and preserved with a 3:1 concentration of water and 96% ethanol in plastic bottles of 30 ml (Samanez et al., 2014). In the laboratory, all specimens were identified at the family's level, using available taxonomic keys (Domínguez and Fernández 2009) and reference collections (Loayza-Muro and La Matta 2015; Roldán 1996).
In addition, taxa richness S and abundance N indices were calculated to analyze the distribution and richness of macroinvertebrate communities in the wetland. Then, the Shannon-Wiener Index (H') and the improved Shannon-Wiener Index (B) were calculated to assess families' abundance and richness (Pielou 1996;Wang et al. 2009). In addition, equitability J' (evenness) of families was estimated from Shannon index H to evaluate the distribution of families in a community (Il et al. 2010).

DPSIR (Drivers, Pressures, States, Impacts, Responses) Assessment
Developed by the European Environment Agency, the DPSIR is a conceptual model that describes the relationships between human beings and the environment under a cause-effect framework. The DPSIR model consists in a chain of causal links starting with a "driver" (a demand that influences the ecosystem, e.g., agriculture) that can cause "pressures" (when demand can modify the ecosystem, e.g., changes in land use) on the "state" of the environment (physical, chemical, and biological) thus generating "impacts" on ecosystems and human health. Finally, the understanding of the resulting impacts leads to "responses" that can be displayed in terms of decisionmaking, target settings, prioritization projects, and policies for risk mitigation and conservation purposes (Kristensen 2004). This study considered an improved DPSIR model developed by Müller and Burkhard (2012), which is based on the "ecosystem service cascade" diagram, to evaluate the state of biophysical processes and ecosystem functions of the Eten wetland and how changes in the environment can impact ecosystem services and human benefits.

Results
(1) Hydro-geomorphic classification and identification of expected ecosystem services of the Eten wetland Based on the hydrogeomorphic classification, the Eten wetland corresponds to an Estuarine system, in which the dominant water source is the overbank flow from the Reque river, and the hydrodynamic flow structure is dominated by bidirectional and horizontal flow. The Reque river describes a delta with small streams flowing into it, forming an estuary that encloses the Eten wetland and describes a wedge-shaped section that is seasonally connected to the ocean ( Figure 3). Moreover, the geomorphic setting is described as a blind estuary with a characteristic funnel shape in the river mouth and foredunes formed along the coastline (Fairbridge 1980). In estuarine environments, such as the Eten wetland, it is expected to have ecosystem services such as groundwater recharge, sediment retention, biological diversity, flood control, carbon sequestration, nutrient cycling, storm attenuation, sediment retention, and water regulation (Barbier 2019;Camacho-Valdez et al. 2013;Costanza et al. 1997;Dugan 1992;Perillo et al. 2009). These ecosystem services were then categorized as provisioning, regulating, cultural, and supporting services, following the Millennium Ecosystem Assessment (MEA 2005). Table 1 presents the results of this classification as an input in the DPSIR model.

Characterization of environmental, biological, and biological parameters
As illustrated in Figure 3, two types of hydraulic structures along the river have been identified: two bridges (the Reque and Eten bridges), and the 3.5-km length dike built along the river. Additionally, intensive agricultural activities demanded channeling riverbanks to derive water inputs for croplands. The lack of waste disposal facilities also resulted in the accumulation of liquid and solid waste within the wetland.

River morphology
The morphology of the Reque river has changed over the years (Figure 4). A comparison between 1949 and 2019 indicated that the Reque river was used to create deltas generating a variety of landforms and ridges of water and sediment, and average sinuosity values indicated a moderate meandering form (1.20-1.38) (Abad, Paredes, and Montoro 2010;Odgaard and Abad 2008), especially from the Reque to Ciudad Eten town. In this section, the averaged channel width is 120 m (100 m) at the upstream (downstream) end, which indicates a channel constraint due to flooding protection structures and floodplain occupation for agricultural purposes.
Aerial photographs from August 2019 illustrate the current state of the Reque river morphology, starting from the Reque bridge until the river reaches the ocean. Figure 1 shows the data obtained from the photogrammetry and the sections to be analyzed: Figure 1B, a section without river embankment, the Reque river presents a wider floodplain upstream in the side in which there are less fluvial constraints (e.g., agricultural lands). Hillshades (azimuth 45°) created from the Digital Elevation Model (DEM) of aerial photographs shows a naturally meandering river morphology based on the observed paleochannel patterns (dotted yellow and red lines). Paleochannels evidenced the natural extension of the river floodplain, which was covered by croplands in different seasons. Vegetation along the stream is represented by agricultural lands, shrubs, and small trees. • In Figure 1C, a section located near Cascajales and a one-sided dike, aerial photographs illustrate the types of sand bars along the Reque stream. The presence of point and channel bars evidenced that the river still maintains its meandering dynamics. However, the dike on the right side constrains the sinuosity of the river, thus reducing channel migration. • In Figure 1D, a section with both-sides river embankment, a bridge and channelization for agriculture, the river morphology presents a different dynamic. Here, the river channel presents less curvature and more mid-channel bars created by the new equilibrium condition along the stream. In this section, there are more croplands that channelize the stream to abstract water Spatial analysis along the Reque river indicated that embankment for flood protection changed river morphology from a meandering to almost a braided stream, in which the channel is constrained to follow only the left margin of its natural  course. Braided rivers are characterized by a wider channel prone to bank that can also increase turbidity (Ashmore 2013). The accumulated plant and waste debris also trap the sediment load and, consequently, increase the flow velocity. It is evidenced by the low vegetation cover in Figure 1D. On the one hand, river constraining can directly impact the natural flow and vegetation cover of the Reque river.
On the other hand, bank erosion might also increase the maintenance cost of river infrastructure because it reduces the effectivity of dikes and flood protection.
The Reque river represents a crucial component of the Eten wetland because it represents the main source of freshwater input, and its morphology defines the wetland structure. However, lateral connectivity has been reduced due to flood protecting structures that constrain the river channel and agricultural land occupation. Channelization also disturbed the riverbanks, thus increasing bank erosion.

Land Cover and Land Use analysis
The qualitative analysis of temporal changes in land cover and land use shows a rapid increase in agricultural areas from 1949 to 2014 (Figure 4). Moreover, aerial photographs suggest soil salinization of agricultural lands. Figure 3 shows saline water springs in crop fields and the presence of salinity tolerant plants. In estuarine environments, river receives both saline and freshwater inputs from water, depending on the season (Phillips 2018). Figure 5 summarizes the particle size distribution sampled at different locations. Bedforms (ripples) of 2-5 cm, typically characterized by fine-to-medium sand mainly carried by wind and low-velocity water flow, had a mean size grain (D50) was 247.60 μm. Sediment texture in croplands was composed of fine sand. The D50 found in crops was 189.1 μm which might be caused by tillage activity (removal and preparation of arable soils) and soil erosion due to river channelization. Moreover, finest material was found on the right side of the Reque river. Therefore, it is suggested that crop areas are an important source of finer sediment material in the Reque river, which is transported by irrigation-water flow.

Analysis of sediment grain size and composition
The substrate in marshes was mainly medium and very fine sand and the D50 was like croplands (189.9 μm). It might be because sampling sites in marshes were located between the river and crops. Riverbanks presented greater medium coarse and fine sand content. The mean grain size D50 is about 217.8 μm and slightly tended to become finer near the river mouth in the estuary. Finally, when the Reque river reaches the sea, it creates a barrier where sediments are deposited. Sediment texture, composed of a medium-coarse to fine sand, is uniformly distributed along the estuary mouth. The average D50 value is 203.4 μm along the shoreline, not so far from the riverbank composition.

Physical-chemical parameter of water quality
As indicated in Table 2, the pH values of all sites were greater than 8 during the dry season, indicating a low alkaline environment. Previous studies indicated that a decline of pH at some levels can disturb biological responses in the aquatic environment (Ringwood and Keppler, 2002;Wright and Worrall, 2001); therefore, pH values are important to monitor the biological diversity of the Eten wetland.
The average value of dissolved oxygen (DO) was 8.975. Differences between oxbow lakes (Figure 1.E) and the river values were very small and there was no evidence of hypoxia in water bodies (<2-3 mg/l O 2 ). Conductivity decreased from the mouth to upstream and increased from the mouth to the oxbow lakes. Salinity and conductivity measures are directly related because dissolved ions increase both values, which suggests that lakes receive more brackish than freshwater inputs with increasing chloride values. The biochemical oxygen demand (BOD) values along the Reque river were less than 2 mg/L. However, BOD increased from the mouth to the oxbow lakes. Values in farther oxbow lakes (W1 and W2) were more than 400 mg/L. Water temperatures have an average of 22°C. Warmer temperatures have been registered within the boundaries of the wetland, mainly in the oxbow lakes far from the river mouth. The chemical oxygen demand (COD) along the Reque river was less than 6 mg/L, while values in farther oxbow lakes (W1 and W2) were more than 1,000 mg/ L. This can be because oxbow lakes are located near the urban areas, thus receiving more domestic wastewater inputs (Figure 3), which also increased the concentration of sulfate and bicarbonate in water.
Water pollution and salinization are also evidenced by the results of the water quality evaluation. For instance, even though higher values of chloride are a consequence of brackish water inputs during coastal flooding and animals and plants are commonly adapted, concentration values could be increased by sewage and wastewater contamination, as illustrated in Figure 2.

Aquatic biota: estuarine macroinvertebrates assemblages
A total of 309 samples of macroinvertebrates were collected from the 16 sampling sites. They belong to 17 families distributed heterogeneously along the Reque river. Sampling sites have been characterized in two sections: from Reque to Cascajales (upstream) ( Figure 1F) and Ciudad Eten to the ocean (downstream) ( Figure 1G) Based on their relative location from the flood protection structure, this division was considered to compare the influence of hydraulic infrastructure in aquatic biota.
Results indicated that taxa richness and Shannon-Wiener index values upstream (from Reque to Ciudad Eten) were greater than downstream (from Ciudad Eten until the river reaches the ocean). It might be because the upstream zone presents a lower disturbance in the river due to channel embankment. Moreover, evenness values suggest a constant population size for families collected downstream. Nevertheless, upstream registered a disparity between the number of individuals within each family. The results indicated that the downstream area harbors fewer families than the area above the tidal limit (upstream), but the abundance of individuals was higher (mainly tolerant families). According to Zhao et al. (2017), the heterogeneity in families' distribution might be due to several causes including: 1) lower temperatures upstream, 2) type of sediment material in the substrate, and 3) few riparian vegetations along the riverbanks in the wetland region. In the Eten wetland, the loss of vegetation due to bank erosion influenced the aquatic biota. Figure 6 shows the relative abundance of macroinvertebrate families. Chironomidae and Baetidae represent the most abundant families, followed by the Palaemonidae. Both Chironomidae and Baetidae families are considered detritus feeders and algae, and during their adult phase as flies they become vectors for allergies (Hall and Gerhardt 2002) (Hall and Gerhardt 2002). Chironomidae individuals are tolerant to environmental stress and their presence also suggests a high organic matter (Roldán 1996), elevated salinity, and heavy metal concentrations in the stream (Loayza-Muro et al. 2010;Mendelssohn et al. 2014). In case of the Palaemonidae, also known as freshwater shrimp or prawn, they play an important role in fishery and aquaculture, and are sensitive to environmental stress (Roldán 1996). Moreover, Ephemeroptera, Plecoptera, and Trichoptera (EPT) orders also positively affect primary production positively and are considered vital role as food for stream fishes. However, they are very sensitive to sedimentation. Upstream, individuals sampled included Odontoceridae taxa in the order Trichoptera upstream, and Baetidae taxa in the order Ephemeroptera downstream.
Results in this section indicate that the diversity of families and abundance of individuals are influenced by changes in the environment. Modifications along the river channels, such as embankment and river channelization disturb the grain-size of sediments (increase in finer sediment concentration), and physical-chemical parameters, such as temperature, DO, COD, and BOD which result in alterations in aquatic biota.

DPSIR (Drivers, Pressures, States, Impacts, Responses) Assessment
The model uses both hydrogeomorphic classification and environmental characterization as inputs to define the drivers (threats), pressures, states, and services of the Eten wetland. Based on the hydro-geomorphic classification (section 4.1), 13 ecosystem services that the Eten wetland possibly provides were identified. The environmental characterization (section 4.2) identified 3 main threats (agricultural, livestock, and farming activities; dikes and hydraulic infrastructure for flood protection; and population growth and rapid urbanization) that influence wetland status, generating changes in river morphology, land use, sediment grain size, water quality, and aquatic biota.
As a first step in the assessment, each wetland condition was linked to one or more threats. For instance, agricultural activities modify the composition of sediment transported by increasing fine material concentration due to river channeling and uncontrolled water harvesting. Then, for each threat, we established a linkage to pressures and impacts generated based on the information collected. Following the previous example, increasing the number of finer sediments in the stream might lead to an increase in turbidity (pressure) and a reduction in aquatic habitats (impact). The last step was to correlate how the impacts affect the ecosystem services and human benefits identified. For example, the reduction in aquatic habitats can modify ecosystem services and leads, for example, to a reduction in wildlife and biological diversity, or the disturbance of nutrient cycling. Table 3 shows a summary of the causal links identified among the components of the DPSIR model for the Eten wetland.

Threat 1: Agricultural, livestock, and farming activities
Pressures: Uncontrolled water harvesting, cultivation of fields and croplands, river channelization and water withdrawal State: The results of the environmental characterization indicated that coastal dunes and the river floodplain have been heavily reduced due to changes in land use by agricultural activities. Consequently, this modified the sediment load and river morphology. Moreover, uncontrolled water harvesting demanded river channelization and the modification of riverbanks.
Impacts: Water exchange between the river flow and croplands also increased ion concentration; thus, it generated water pollution and fostering an abundance of tolerant biota. These pressures resulted in the alteration of water balance, soil salinization, and disturbance of the hydrological connectivity and nutrient cycling. Affected ecosystem services included provisioning, regulating, and supporting. Based on the hydrogeomorphic classification, it is expected that these impacts could affect groundwater recharge, sediment retention, waste processing, and the biological diversity of the Eten wetland.

Threat 2: Dikes and hydraulic infrastructure for flood protection
Pressures: Construction of levees and dikes, and channel embankment State: Results of remote sensing analysis indicated that river channeling and embankment constrained the lateral dynamics of the Reque river. When the dynamics of river morphology is reduced, it causes the loss of lateral connectivity between the river and the wetland.
Impacts: Impacts on ecosystem services resulted in the alteration of water balance and the disturbance of flow regime and hydrological connectivity (possible impacts). Moreover, channel embankment can cause the capability of storm protection, thus leading to human and economic losses during flooding and extreme events. Modification of river morphology also increased the supply of very fine sediment and the velocity of river flow, thus generating bank erosion that reduced plant communities and disturbed biota and primary production population in the riverbanks. Results indicated that ecosystem values affected would be water regulation, flood control, biological diversity, and storm protection.

Threat 3: Population growth and rapid urbanization
Pressures: Waste disposal, bulldozing of marshes, land claim, and uncontrolled groundwater pumping State: The evaluation of changes in land use and water quality indicated that the urban sprawl and rapid urbanization reduced the Eten wetland area and polluted soils and water bodies. The contamination of the Reque river and lakes within the wetland is worsening due to the lack of sewage and sanitary infrastructure in the area. Uncontrolled groundwater pumping could also foster the infiltration of pollutants into the aquifer, but more measurements should be done.
Impacts: Pressures impacted the Eten wetland by reducing wetland area, plant communities, habitats, and esthetic value. Therefore, the ecosystem services that could be affected include biological diversity; the cultural, educational, and research value; and the capability of waste processing and biomass export.

Geomorphic and environmental features of the Eten wetland
Results of both the hydro-geomorphic and environmental characterizations indicated that the Reque river plays an important role in defining the landscape of the Eten wetland. Aerial photographs also illustrated the dynamism of the Reque river and its natural floodplain extension. The meandering nature of the Reque river determined the location of agricultural lands because it transported not only water but also nutrients that fertilize soils during flooding events. In the last 70 years, the agricultural land area increased in 125% which fostered the demand for irrigation water. Because the river represents the main source of freshwater, channelization is a common practice for local farmers to redirect water. These modifications disturbed not only channel morphology but also habitats in overflow areas. Besides, the environmental characterization indicates that urban sprawl has caused the contamination and desiccation of water bodies and soil in areas near the urban areas.
The analysis using photogrammetry and GIS-tools indicated that embankment and river channelization constrained the lateral connectivity of the river, thus modifying floodplain morphology by increasing flow velocity and reducing sinuosity. While riverbanks from Reque to Ciudad Eten are prone to erosion, the region between the Ciudad Eten and the river mouth presents bank failures that led to the loss of riparian vegetation and aquatic communities. Therefore, bank erosion generated more sediment supply (mainly silt and clay) that increased turbidity.

Impacts of human activities on ecosystem services
The presence of hydraulic structures for flood protection in the Eten wetland caused the reduction of floodplain ecosystem functions, such as water quality and biological diversity. Losing the capacity for lateral migration reduces the formation of geomorphic structures, thus affecting flow intrusion into the floodplain and leading the soil to desiccation. Moreover, the hydrologic connectivity within the floodplain allows the streamflow to trap sediment-associated nutrients, but in crop areas trapped material also includes agriculture-derived pollutants (Hupp, Pierce, and Noe 2009).
Benthic macroinvertebrates play an important role in supporting the food web by cycling nutrients and contributing to the productivity of higher trophic levels (Zhao et al. 2017). Riverbanks provide macroinvertebrates of ecotone habitats, but they are highly sensitive to any variability of the flow, sediment concentration, and nutrients transported. Changes in water quality were evaluated based on the type of families found and how they were distributed along the Reque river. This evaluation indicated that the presence of hydraulic structures (dikes) and agricultural activities reduced the aquatic habitats and increased the diversity of tolerant families, thus disturbing the biological diversity. Moreover, urbanization in wetlands tends to generate negative impacts of aquatic biota, as indicated in Novoa et al. (2020) and Vizcardo and Gil-Kodaka (2015).
Results also evidenced that changes in land use disturbed the natural landscape and geomorphology, thus affecting support and regulation services, such as groundwater recharge and water regulation. Additionally, agricultural activities caused salinization of soils and desiccation of lakes due to the disconnection of the river with the floodplain. Based on municipality records, these problems are currently exacerbated by the lack of regulations, sanitation facilities, and culture of water.

Applicability for sustainable ecosystem management and conservation plans
Nowadays, most Peruvian coastal wetlands face similar issues (agricultural and urban development); however, the ecosystem response in each environment may be different. Therefore, understanding the natural processes and functions of a coastal wetland allows the valuation of ecosystem services, and its vulnerability and capability to respond. Conservation strategies and sustainable management can then be implemented based on a better understanding of the functioning of the environment and the pressures of human activities.
First, it is important to characterize the type of coastal wetland because it determines the geomorphic settings that describe the landscape. The hydro-geomorphic characterization demonstrated the role of the Reque river in determining riverine communities along the alluvial corridor and fertilizing the floodplain due to the lateral connectivity with the Eten wetland. For instance, other types of coastal wetlands are depressional (Pantanos de Villa) and lacustrine (Albufera de Medio Mundo) wetlands which might be influenced by other components, such as precipitation and groundwater; thus, the geomorphic features and ecosystem functions will be different. Once the type and features are described, potential ecosystem services of a coastal wetland can be defined. For example, depressional wetlands such as swamps and peatlands provide the service of carbon sequestration and storage. Second, an environmental characterization assesses the status of the components of a coastal wetland. By evaluating the changes in the geomorphic features (defined in the previous step); the quality of water, soil, air, plant, and animal communities; and landscape over time, it is possible to define the main threats and pressures through a cause-and-effect relationship. Third, the DSPIR framework provides a wide overview of how threats and pressures can disturb the ecosystem services and the impacts generated. Based on this information, decision-makers can respond by prioritizing conservation projects, establishing regulations to reduce the pressures, and other actions to mitigate impacts. Moreover, the framework can be applied as many times as required to evaluate the implementation of actions and responses.
The proposed methodology aims to be a tool to identify the most effective responses to reduce environmental impacts in ecosystem services, with a scope that this is a cyclic process that needs to be reinforced over time. Applied to a Peruvian coastal wetland with little field data, this study provides a rapid diagnosis of ecosystem status and potential impacts on services. Moreover, the environmental characterization applied in the methodology is not limited to the areas described above, but also nutrient analysis, soil quality, water balance, sediment transport, among other disciplines. This versatile methodology, however, can be improved by increasing the time, area of study, and field data collection. On the other hand, results are also presented in plain language for non-technical decision-makers, which helps them to easily evaluate responses and actions.

Conclusions
This study proposes a methodology to assess the impacts on ecosystem services based on the understanding of geomorphic features and status of a Peruvian coastal wetland. By combining remote sensing applications and field surveys, this study presents an approach to improve coastal wetland management through the identification of the most affected functions in terms of geomorphological changes and environmental parameters.
Considering that the study area had very limited local data, it was challenging to characterize the landscape features using only available literature data and open-source GIS tools and products. The study can be improved by acquiring high-resolution multitemporal imagery and field surveys to map the full extension of the wetland. Furthermore, the environmental characterization consisted of spatial analysis that identified the main threats and pressures in the Eten wetland. However, the accuracy of the results should be complemented with more field evidence and spatiotemporal monitoring. This study addresses the need to enhance decision-making in terms of conservation and sustainability. Based on the information derived from this methodology, it is possible to establish actions to mitigate impacts.
The results of this study indicate that agricultural and urban development generates critical pressures on the ecosystem services of the Eten wetland. Uncontrolled water abstractions and overbank modifications can lead to the loss of the wetland and its ecological value. Indeed, it is required to both prioritize conservation projects and to foster more research, monitoring, and educational activities to promote the value of the Eten wetland in the region. The methodology also aims to be versatile and provide a rapid diagnosis. Researchers can add more environmental parameters and enhance hydrogeomorphic characterization with more field data. The scope of the study can be easily adapted to inland wetlands as well.