Next Article in Journal
Fuzzy AHP Integrated with GIS Analyses for Drought Risk Assessment: A Case Study from Upper Phetchaburi River Basin, Thailand
Previous Article in Journal
Reclaiming Indigenous Planning as a Pathway to Local Water Security
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Integrated Approach for Evaluating Water Quality between 2007–2015 in Santa Cruz Island in the Galapagos Archipelago

by
Cristina Mateus
1,2,*,
Christian A. Guerrero
1,
Galo Quezada
3,
Daniel Lara
3 and
Valeria Ochoa-Herrera
1,2,4
1
Instituto Biósfera, Colegio de Ciencias e Ingeniería, Universidad San Francisco de Quito, Diego de Robles y Vía Interoceánica, Quito 17-1200-841, Ecuador
2
Galapagos Science Center, Universidad San Francisco de Quito and University of North Carolina at Chapel Hill, Puerto Baquerizo Moreno, San Cristobal, Galapagos 200150, Ecuador
3
Dirección del Parque Nacional Galápagos, Ministerio del Ambiente, Av. Charles Darwin S/N, Puerto Ayora, Santa Cruz, Galápagos 200350, Ecuador
4
Department of Environmental Sciences and Engineering, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
*
Author to whom correspondence should be addressed.
Water 2019, 11(5), 937; https://doi.org/10.3390/w11050937
Submission received: 21 February 2019 / Revised: 10 April 2019 / Accepted: 23 April 2019 / Published: 3 May 2019
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Water quality in Galápagos has been deteriorating by increased human impacts over the past few decades. Water quality is a key environmental component and issue in need to be evaluated in the Pelican Bay Watershed, the biggest urban and economic development of Santa Cruz Island, for better management and regulation of water resources. This study assesses coastal and ground water bodies of Pelican Bay by employing a 9-year dataset obtained during a local water quality monitoring program conducted by the Galápagos National Park. Physical-chemical and microbial parameters were evaluated with respect to national and international water quality standards. A statistical integrated approach was performed to calculate environmental background levels of water quality parameters and to explore their seasonal and spatial variation. In addition, a sensitivity analysis was conducted to evaluate the impact of changes in tourism and residents in San Cruz Island in the degradation of water sources. Results highlighted are: (a) water is not suitable for drinking and domestic use at some inland sites; (b) saline water is used for irrigation in the highlands; (c) the presence of parameters of concern at coastal sites represent a risk for human and ecosystem health; (d) background levels may serve for defining site-specific limits to control water quality, and; (e) the influence of population change on water quality conditions varied at each site with a higher effect at coastal sites relatively to inland sites. This study provided valuable information of the water quality status in Santa Cruz Island and can serve as a baseline for effective water management and control of pollution.

Graphical Abstract

1. Introduction

Growing water quality impairment represents a global scale problem threatening human development and ecosystem integrity. Alterations of water can pose a threat to public health and affect human access to safe water for drinking and domestic purposes, food production, and recreation [1]. In aquatic ecosystems, poor water quality can lead to detrimental effects such as the increase in algal blooms, aquatic species mortality, increase in fish disease, among others [2,3]. Within the Galapagos Islands, human activities have contributed to major alterations in coastal and ground water quality [4,5,6], resulting in public health problems such as a high incidence of respiratory and gastrointestinal parasitic infection in residents and tourists [4]. In addition, as observed in other parts of the world, poor water quality can be impaired further by natural-climate driven processes [7,8,9,10]. This study seeks to evaluate coastal and ground water quality on the Pelican Bay watershed in Santa Cruz Island to gain a better understanding of the current status of water quality along with seasonal and spatial variations.
Increasing permanent human population and tourism in the Galapagos exerts a significant pressure over the fragile ecosystem of the islands and enhances water degradation [5,11]. Main groundwater sources for Santa Cruz Island are contaminated due to several factors such as the location of the basal aquifer beneath dense urban settlements, the lack of effective wastewater treatment plants, and seawater intrusion [12]. In fact, seawater may be migrating inland due to overexploitation of the aquifer [5,13,14]. In particular, malfunctioning septic tanks are responsible for most of the fecal contamination in the island [4]. Main sources of fecal contamination for coastal marine waters in Puerto Ayora, the main urban area on Santa Cruz Island, include: land runoff of sewage waters during wet weather, submarine and boat discharges, and overland pipe sewage effluents [5]. Additionally, high levels of E. coli in the coastal waters seem to be related to groundwater leachate flowing seaward during low tide [15]. Yet, more specific monitoring is needed to understand the complete picture of potential water impairment in the Galapagos Islands and to identify the implications of basin sensitivity on water quality to increases in tourism and residents.
Besides the anthropogenic factors, water quality can be directly affected by several climate-related mechanisms in both the short and long term. For example, water quality variation in the Huai River Basin in China is significantly correlated with precipitation and strongly affected by air temperature [7]. Other studies [8,9] found that intense rainfalls and runoff make pathogens more susceptible to rise in concentration in superficial waters. Moreover, water temperature increases reduce gases solubility and favors many physical-chemical processes such as dissolution, solubilization, complexation, nitrification, degradation, and evaporation occurring in water bodies, leading to elevated levels of dissolved substances found in temperate and tropical regions [10]. Additionally, mineralization increases leads to a greater release of nitrogen, phosphorus, and carbon from soil that is then transported by runoff [16]. Warmer temperatures may also impact water bodies with long residence times, which would imply eutrophication problems in the future [16]. Regarding precipitation, literature studies have demonstrated that pollutant transport is enhanced by heavy rainfall events, which increase runoff and erosion, and in many cases cause higher levels of turbidity and organic matter [10]. Therefore, similar climate-related processes as those described herein may affect water quality in the Galapagos Islands. Although, other factors, such as the direction of the groundwater flow, the distance from the sea, and changes in altitude, have been observed to influence changes in water temperature, pH, conductivity, and dissolved oxygen (DO) in Santa Cruz Island [13], there is still the need to investigate water quality variability due to seasonal and natural processes in the Galapagos.
Determining environmental background levels constitutes a helpful comparative tool that can be utilized to define thresholds values for specific water bodies [17,18] and defines if water characteristics are influenced by either natural processes or human activities. Background concentration is defined as “a concentration of a substance in a particular environment that is indicative of minimal influence by human (anthropogenic) sources” [19]. Therefore, the environmental background levels can reveal natural conditions unaffected by anthropogenic influence [17]; thus, it enables us to distinguish values of a particular parameter that are product of natural processes from anthropogenic impacts [20,21,22]. Environmental background levels and the upper limit background levels (or thresholds) have been established worldwide for many waterbodies [17,22]. However, the application of large-scale background values to evaluate local status of water quality can lead to misunderstandings because local background is constrained by local conditions [17,22] such as water–rock interactions, chemical and biological processes, residence time of water, recharge inputs and transport processes [18,21]. For that reason, statistical methods are preferred for the calculation of local background levels.
This study evaluates coastal and ground water status in Pelican Bay, a human-influenced watershed in Santa Cruz Island, using a 9-year (2007–2015) dataset obtained during a local water quality monitoring program conducted by the Galápagos National Park. Physical-chemical and microbiological parameters were used to assess water quality and its suitability for designated purposes at each study site, such as potable water, irrigation and/or recreation based on national and international guidelines. Additionally, environmental background levels were calculated to reveal natural conditions unaffected by anthropogenic contributions within the 9-year study period. To evaluate whether or not changes in tourism and residents are responsible for degradation of the water sources between 2007 and 2015, a sensitivity analysis was conducted. Finally parameters were analyzed to describe their seasonal and spatial variations.

2. Materials and Methods

2.1. Study Area

Water quality parameters at seven sites located in Santa Cruz Island in the Galápagos Archipelago at the equator in the Eastern Pacific Ocean, about 1000 km off the coast of Ecuador, were analyzed in this study. The sites were classified into two groups (Figure 1): (a) Inland sites that measured groundwater in the aquifer: Pozo Profundo (I1), INGALA (I2), Misión San Francisco (I3); and (b) coastal sites: Las Ninfas Lagoon (C1), Academy Bay in (C2), Academy Bay out (C3), and Muelle de los Pescadores (C4). All sites are located within Pelican Bay watershed in Santa Cruz Island. This watershed occupies an area of 43 km2 including the concentration of the main population and economic development of the island [23]. Water quantity is an issue in this watershed due to the demand for an extensive agricultural area (irrigation and cattle) that occupies 22.5 km2, and for municipal and domestic supply for the two main towns, Puerto Ayora and Bellavista [23]. In addition, water quality at the study sites, as noted above, is subject to perturbations associated with human activities.
Inland sites (Figure 1) provide access to the basal aquifer water of Santa Cruz, which consists of brackish water resulting from seawater laterally intruding the highly pervious fractured basalt, at least 9 km inland [14]. Water is extracted for different purposes from this aquifer. For example, site I1 is a well built on the area of agricultural lands at 4.7 km from the sea and 157 m.a.s.l. The well has been operated since 2002 by Bellavista municipality, supplying water to the local population of approximately 1721 inhabitants [13]. Water from this source is considered safe for domestic use since no coliform contamination has been detected yet; however, high contamination risk is present due to its proximity to the main road, urban settlements and saline intrusion [23]. Site I2 corresponds to a crevice located in Puerto Ayora’s periphery at 1.1 km from the sea and 23 m.a.s.l., and it is surrounded by urban cover and mixed vegetation. It was tapped in the 1980s and provided water for drinking and domestic use to Puerto Ayora until the detection of E. coli and other coliforms determined in 2011 [4,14]. Currently, its use is limited to supply water for irrigation and cleaning purposes [13] while it is facing risk of an increase of fecal coliform contamination [24]. Finally, I3 is a crevice within Puerto Ayora urban area, at 0.5 km from the sea and 20 m.a.s.l. Its water was distributed through the municipal water supply and designated for private desalinization until its closure in 2011 due to the presence of extremely high levels of E. coli and total coliforms [4,12,13]. Despite that inland sites are located at different elevations, extracted and sampled water belong to the same basal aquifer whose water level is virtually the same as sea level [25].
Coastal sites are located within Puerto Ayora’s coastline (Figure 1). Site C1 consists in a small semi-enclosed basin with brackish waters that, due to its salinity and tidal influence it approximates to an estuary. This place was a popular swimming spot for tourists and residents until signs of human fecal contamination were detected [26] which in part was due to an inadequate management of stormwater runoff in the surrounding area [5]. Nowadays swimming and other recreational activities are prohibited since they represent a high potential risk of waterborne disease incidence [26]. Moving seaward from C1 are sites corresponding to Academy bay, a narrow and long indentation of Puerto Ayora’ shoreline, where C2 covers the enclosed water and C3 the outlet. On the other hand, C4 site is in the southeast bay and represents the last point of freshwater outflow to the sea. Both, C1 and C4, host surrounding mangroves, though their ecological integrity was disturbed due to urban development.

Climatic Conditions

Due to its location, Santa Cruz Island is characterized by cool and dry climatic conditions influenced by extreme inter-annual variability and spatial variability between coast and highlands [14]. Two seasons in the island, as well as in all the Archipelago, are clearly defined: the cool garúa season, from June to December, and the hot invierno, from January to May [27]. During the garúa season, average temperatures are 22.6 ± 3.4 °C and the cool air brings a moisture-laden mist above approximately 250 m, whereas below this elevation the lowland areas remain dry [28]. During the invierno season, rainfall is convective; the amount of precipitation is positively correlated with sea surface temperatures [27], the climate is warm (26.0 ± 2.0 °C) with occasional heavy rain showers, and dry spells result in strong evapotranspiration [14]. Spatial variability has also been identified: annual rainfall ranges from 502 mm on the coast to 3068 mm in the highlands [14] whereas the temperature gradient is of −0.8 °C per 100 m of elevation [29]. El Niño Southern Oscillation (ENSO) also induces extreme inter-annual climatic variations.

2.2. Study Approach

2.2.1. Water Quality Parameters

Analyses employed a 9-year (2007–2015) dataset obtained from the routine water quality monitoring program carried out by the Galapagos National Park in Santa Cruz Island to monitor the overall status of inland and coastal waters [30]. Data consisted in a total of 11 water quality parameters presented in Table A1: electrical conductivity (EC), salinity (Sal), dissolved oxygen (DO), water temperature (WT), pH, turbidity (Turb), fecal coliforms (FC), total coliforms (TC), nitrite (NO2-N), nitrate (NO3-N), and total phosphorus (TP).
Over the 2007–2015 sampling period, monitoring based on field measurements and sampling was done by the Galapagos National Park once a month or a few times in the year between 9:00 and 11:30 am at all study sites [30]. Briefly, field measurements of DO, pH, EC, salinity, temperature, and turbidity were conducted in situ using field portable instruments in triplicates (Table A1). According to Galapagos National Park, nitrate, nitrite, and total phosphorus were analyzed using a spectrophotometer, and fecal and total coliforms were detected by using a portable microbiology laboratory [30]. Overall, sampling, preservation, transportation and analysis of the water samples were as per standard methods [31] as described by Galapagos National Park [30].
The main focus of this study was to complete a cursory analysis of available, but limited, data to evaluate water quality in Santa Cruz Island. We conducted a statistical analysis based on data collected by a third party in order to calculate environmental background levels of water quality variables and to explore seasonal and spatial variations. The analytical techniques applied for data collection and data monitoring were not evaluated in this study. While monitoring was carried out monthly, some parameters, such as pH and DO, may present changes in narrower time-scales like hours or minutes. Therefore, the frequency of the data analyzed herein may not be enough to represent short-term variability. Additionally, data may be impaired by the presence of outliers, missing values and possible measurement limitations that were not appropriately verified when they were reported by the Galapagos National Park [30].

2.2.2. Meteorological Variables

Basic meteorological records from 2007–2015 were obtained from existing long-term manual stations operated by the Charles Darwin Foundation, which are located at 2 m.a.s.l. (CDF1 Charles Darwin) and 180 m.a.s.l. (CDF2 Bellavista), respectively (Figure 1) [25]. Assessed variables include: mid-day air and sea temperatures, and total daily rainfall. Air temperature was recorded daily in the shade 2 m above the ground surface, sea temperature was recorded at the closest coast of CFD1 in a bucket pulled from the sea, and total daily-rainfall was the sum of measurements taken at 06h00, 12h00 and 18h00 from a Hellmann rain gauge 1.5 m above the ground [32].

2.2.3. Evaluating Water Suitability for Designated Water Uses

Suitability of ground and coastal water for different designated uses was evaluated for each parameter considered in this study (Table A1). These parameters were compared to standards of three different guidelines on water quality: (a) Unified Text about Secondary Environmental Legislation of Ecuador 2015 [33], (b) WHO guidelines for drinking-water quality [34] and (c) Hawaii Administrative Rules Title 11, Chapter 54 [35]. Thus, for inland sites, parameters were evaluated comparing Ecuadorian water quality criteria for drinking and domestic uses and criteria for protection of agriculture [33], and WHO drinking water guidelines [34]. For coastal sites, parameters were evaluated with specific criteria for estuaries and embayments from water quality standards of Hawaii Administrative Rules Title 11, Chapter 54 [35]. These standards were chosen due to hydrogeological similarities between Hawaii and Galapagos established in previous studies [14,36]. Ecuadorian guidelines for preservation of wild and aquatic life in marine waters and estuaries, and for primary and secondary contact recreational waters [33] were also used. Guideline values and criteria considered for the analysis on both inland and coastal sites were those that matched the available data. In order to assess DO data of coastal sites, DO was converted to percent saturation units through Benson and Krause equations [37].

2.2.4. Statistical Techniques for Environmental Background Level Determination

Several different approaches exist to determine environmental background levels of water quality variables in water bodies. For this study, we used two statistical techniques presented by Matschullat et al. [38], using two model-based objective methods to define background levels of water quality: (a) the iterative 2σ technique, and; (b) the distribution function (df) technique. The data was statistically processed using an Excel application developed by the University of Zagreb that uses an algorithm designed to estimate background levels by means of the interative 2σ and the distribution function [20]. Both techniques aimed to define the background by approaching a normal range [20,21,38] so then, anomalous and background subpopulations could be differentiated. For most water quality parameters, background levels are those with lowest values and most frequent occurrence, whereas anomalies due to anthropogenic influence lie separated and are less numerous. After calculating background levels, the Lilliefors test, an adaptation of the Kolmogorov–Smirnov test, was applied for each data series to evaluate the goodness of fit to a normal distribution by generating the t-statistic as suggested by Nakić et al. and Urresti-Estala et al. [20,21]. If the t-statistic value of a parameter is below the critical t, background values are considered to fit a normal distribution with 95% confidence level, suggesting they have been estimated correctly and the utilized technique was appropriate for such parameter [21].
The iterative 2σ and the distribution function (df) techniques estimate the threshold between background and anomalies. The iterative 2σ calculates the threshold as the upper and lower values of the background variation since both high and low values may define anomalies in certain parameters [21]. This technique constructs a normal distributions around the mode of the original dataset and calculates the mean and standard deviation of this normal distribution, then the 2σ range is stablished, corresponding to the background levels, and the values outside this range are removed [38]. The distribution function (df) defines the upper limit of background levels by taking the positive asymmetry of a normal curve as anomalies that may be the result of anthropogenic influence [38]. Specifically, this technique calculates the median of the data and eliminates the values above it, what is left is reflected towards values above the median to stablished a range from which a median and a standard deviation is calculated [38]. Advantages toward utilizing these techniques have been addressed in previous studies. For example, the targeted dataset does not require a specific prior distribution, such as to be normal or lognormal distributed [21]. However, they have been showed to be more successfully applied in unimodal and skewed distributions [20]. Another advantage is that both techniques can be applied to relatively small data populations (n > 30) in contrast to other techniques that require hundreds or thousands of data to be successfully applied.
For practical application in water quality control, background levels were calculated for all water quality parameters except for microbial parameters (TC and FC) because the utilized techniques have only been applied to geochemical parameters [20,21,38].

2.2.5. Estimating Sensitivities to Changes Population

To evaluate whether or not changes in population are responsible for the degradation of the water sources between 2007 and 2015, we a conducted sensitivity analysis (S) that takes into account changes in population and changes in parameter concentration over the 9-year period of study. We calculated S using a simplified approach with linear regression analysis applied to data at the annual time scale to develop a general understanding of the relative sensitivity of water quality parameters to changes in population at each study site. Regression analysis methods to evaluate the significant influencing factors of ambient water quality conditions and estimate the magnitude of their impacts have been applied in other areas [39,40,41,42,43,44]. However, in reporting the results, we acknowledge that their interpretation is greatly limited by assuming that any interactions across space and time are not fully represented in this analysis.
Changes in population was calculated as the number of permanent residents in Santa Cruz Island and the number of tourists that stay in the island per day. The number of permanent residents was obtained from the last three population census done in Galápagos in 2001, 2010, and 2015 [45]. We first estimated the growth rate between census and then used a regression analysis to estimate the slope of the curve between 2001 and 2010, and 2010 and 2015. Thus, the number of permanent residents between 2007–2015 was obtained. The number of tourists staying on the island was obtained from the Galapagos Tourism Observatory [46]. Tourists landing in Baltra Airport in Santa Cruz Island but not staying in Puerto Ayora overnight, were not considered in our study.
The sensitivity analysis evaluates the strength of relationship between the annual percent change in population and annual percent change in water quality parameters based on the R2 and the size of the effect of population change on water quality based on the ß1 (slope). Therefore, the sensitivity analysis will also determine whether or not sources of contamination during the 9-year period are due to anthropogenic influences (strong relationship and large effect: high R2 and ß1), natural variability (weak relationship and low effect: low R2 and ß1) or both (weak relationship and large effect: low R2 and high ß1). No relationship (R2 < 0.2) means that any change in parameter concentrations of water quality between 2007 and 2015 is not related with changes in population.

3. Results

3.1. Parameters Exceeding Water Quality Criteria and Guidelines

Parameters at inland sites were evaluated according to values exceeding the WHO guidelines and Ecuadorian criteria of water quality (Table 1). Fecal coliform levels are alarmingly high at I3, and therefore not recommended for drinking nor domestic use. Results suggested that water at inland sites may be not suitable for drinking and domestic use due to its Sal and Turb levels. Salinity levels at I2 and I3 (90% and 100% of values exceeding) exceeded acceptable levels in the WHO water quality guidelines [34]. Site I1 can also be considered unpalatable but in less extent (26%). Turb levels at all inland sites exceeded guidelines, particularly I3 (18% of values exceeding). Furthermore, some values of EC and Sal exceeded irrigation water limits, making its water not recommended for neither drinking nor irrigation. Similar to Sal, site I2 and I3 have more EC values within severe limits relative to I1. Inland sites levels of NO3-N, NO2-N, are within safe limits of WHO guidelines and Ecuadorian criteria for drinking water, and irrigation water for agriculture and water for livestock, respectively. Therefore, they are not presented in Table 1.
Parameters at coastal sites were evaluated according to values exceeding Hawaiian rules and Ecuadorian criteria (Table 1). Coastal site C1, frequently exceeded the guidelines for: DO, EC, FC, the sum of nitrogen [NO3+NO2]-N and TP. Levels of DO exceeded under the three criteria considered though exceedance was higher at C1 when compared to Ecuadorian criteria of recreational waters (77% of values exceeding). For this criteria, high levels of fecal coliforms at C1 (88% of values exceeding) were also reported. Nutrient levels as [NO3+NO2]-N and TP often exceeded the Hawaiian rules (99–100% of the data) at both C1 and C4. However, the concentrations of NO3-N and NO2-N comply with the Ecuadorian criteria (13 and 0.2 mg/L, respectively) for the preservation of wild and aquatic life in marine waters and estuaries. Levels of Turb at C4 were also high for an embayment (99% of values exceeding). Finally, levels of TC fell below Ecuadorian criteria for secondary contact waters (2000 MPN/100 mL) in C1.

3.2. Environmental Background Levels

Environmental background levels of the physical-chemical parameters were determined using two techniques: (a) the iterative 2σ and (b) the distribution function (df) as described in Section 2.2.5. The outer thresholds of background values for the parameters considered are reported in Table 2. Best estimations are indicated in bold and shaded in grey according to the nature and distribution of original data, and recommendations proposed in previous studies [20,21,38]. Variables that passed the Lilliefors statistic test are also indicated (*) which means that the data fits a normal distribution and that the techniques applied herein are appropriate for such parameters. Overall, results demonstrated that: (1) background levels calculation fit better for coastal sites than for inland sites; (2) the techniques applied herein performed different for each parameter and study site, and; (3) while some parameters exceeded water quality criteria, background levels suggested that those parameters have been varying naturally since 2007.
As demonstrated by coastal sites (C1–C4), background ranges yielded a greater number of parameters that passed the Lilliefors test with at least one of the techniques considered suggesting a better fit for coastal sites relative to inland sites.
In most cases, both the iterative 2σ and distribution function (df) techniques produced similar results, however, differences in their performance was observed. In general, results suggested that: (a) the iterative 2σ technique usually produced the best estimations for EC, Turb and Sal; (b) both techniques were equally suitable for TP, NO2-N, and NO3-N, and; (c) none of the techniques are a good fit for WT, pH, and DO. Performance differences depend primarily on the frequency distribution of the evaluated datasets. These differences might be provided by the degree of human disturbance, along with the hydrodynamic and hydrochemical characteristics of the water body [21]. Parameters on which the iterative 2σ worked well (EC, Turb, and Sal) presented unimodal or skewed distributions in which background values are calculated around this unique mode and they follow a normal distribution (Figure A1a), thus making them suitable for the application of this technique [20]. In contrast, background values for EC, Turb, and Sal calculated using the df technique fails Lilliefors test because positive anomalies in data cannot be distinguished and thus upper thresholds are overestimated, making backgrounds to not follow a normal distribution (Figure A1b). On the other hand, parameters that showed overlapping polymodal distributions, which may include right skewness, with prevailing lower values representative of background levels (NO3-N, NO2-N, and TP), were suitable to apply the distribution function technique (Figure A2b) [20,38]. When the distribution of the original dataset is dispersed and has larger number of small values that may represent background levels, the df technique seems to work remarkably good relative to the iterative 2σ technique (Figure A2). The parameters WT, DO, and pH did not present any of the mentioned distributions in most cases, therefore most of the calculated background values across study sites did not pass the Lilliefors test. In the case of WT and DO, distributions different from normal and skewed were documented on these parameters previously [47].
Differences in the range of the background values were also observed between study sites. Background values for parameters such as Turb, Sal, and EC were higher in coastal sites, which may be associated with its proximity to the sea. In contrast, TP background is higher among inland sites which could be a result of land human activities. For NO3-N and NO2-N ranges of background values were very similar.
Levels for Sal, EC, Turb, [NO3 + NO2]-N, and TP exceeded water quality criteria at some study sites. However, levels of these parameters are considered to be varying naturally according to background values calculated for the 9-year time period studied herein. For example, EC and Sal levels make water at inland sites unsuitable for irrigation, however, those levels fall mostly within calculated background levels (Table 2) suggesting that have been varying naturally since 2007, reflecting natural brackish conditions. High levels of Turb at coastal sites largely exceeded Hawaiian criteria. However, according to the background levels, is subject to natural increases that surpasses Hawaiian specific criteria for estuaries (3 NTU). Finally, while [NO3 + NO2]-N and TP largely exceeded Hawaiian criteria for coastal waters, background levels of these parameters show higher natural-occurring values than the limit set for coastal waters in Hawaii. This suggests that the use of Hawaii’s water quality rules for evaluating our coastal study sites may lead to misunderstandings about the degree of anthropogenic water quality alteration.

3.3. Relative Influence of Population Change on Water Quality Parameters

The sensitivity analysis provided useful information about the significant influence of population change on water quality conditions over the 9-year study period. Results highlighted that (a) the strength of the relationship and the influence that changes in population has over water quality parameters varied at each site; (b) the strength of the relationship was stronger at coastal sites relatively to inland sites, and; (c) anthropogenic influence on nutrients (TP, NO3-N and NO2-N) seems to be higher the closer they get to the coast (C2, C3, and C4).
Higher sensitivity to anthropogenic activities was observed at coastal sites relative to inland sites between 2007 and 2015. Changes in water quality parameters at inland sites exhibited relative weak relationship with percent change in population, based on the low R2 (Table 3). For example, a large change in population size at I3 resulted in a small change (low ß1, small effect) in FC concentrations. Thus, results suggested that changes in the number of habitants at inland sites has a small impact on water quality parameters. On the other hand, based on the large number of high R2 found at coastal sites (Table 3), results suggested a strong relationship between changes in population and water quality parameters at these sites. However, the size of the effect of anthropogenic activities varied between parameters, based on the high or low ß1 coefficient (slope). For example, at C4, a small change in population (high R2) resulted in a big change in Turb values (high ß1, large effect) and a small change in NO3-N concentrations (low ß1, small effect). Results also highlighted that the strength of the relationship between population changes and nutrients ([NO3 + NO2]-N and TP) increases the closer the site is to coastal waters.

3.4. Seasonal Variations

Data variability within invierno and garúa seasons differed for each parameter suggesting that some of them may be influenced by seasons. While WT, EC, Sal, Turb, NO3-N, TC, and FC present seasonal variations (Figure 2), DO, pH, TP, and NO2-N did not show any important differences between seasons (Figure A3). In general, higher variation happened during the garúa season relative to the invierno season with the exception of TC and FC in which higher variability happened during invierno. Results for WT, Turb, NO3-N, and TC are shown in Figure 2, and results for EC, Sal, and FC are presented in Figure A4.
Variations observed in pH are within normal pH levels of coastal water and coastal aquifers and may not be a problem for ecosystems health and recreational activities [48]. However, a 10% exceedance was found at C1 (Figure A3) which still falls within estuarine pH levels that generally average from 7.0 to 8.6 [48]. Literature suggests that pH anomalies may be caused by bacterial activity, chemical constituents in runoff and sewage overflows which can greatly stress most aquatic plants and animals [48]. However, further analysis needs to be done to better understand the possible causes of the anomalies found in our results.
Seasonal analysis shows that coastal WT variability is consistent with seasonal changes, where higher WT concurred with the hot invierno season and the lower temperatures fell during the garúa season (Figure 3). However, differences between WT observations and measurements of air and sea surface temperatures from the Charles Darwin Foundation (CDF) stations were observed within study sites. Water temperature at coastal sites (Figure 3a–d) behaved similar to air and sea surface temperature between January to May. However, it appears to be one month lag between high peaks for air temperature and coastal water temperature. High peak for air temperature is in March whereas high peak for coastal water is in April. From September through December, a steep rise of WT would be expected according to the seasonal variation seen for air and sea surface temperature; however, the high temperature in September compared to the neighboring months breaks this pattern since unusual high water temperature were reported for this month at many years of the study period. Water temperature regime at inland sites (Figure 3e–g) presented a weaker seasonal variation (also seen in Figure 2b) and only I3 which is Misión San Francisco (Figure 3e) seems to receive seasonal influence that affects its temperature. Periodic annual fluctuation of these inland sites has a noticeable amplitude less than 1 °C.

3.5. Spatial Variations

Physical and chemical parameters were analyzed to explore patterns of spatial variations within Pelican Bay water quality. For WT, EC, Sal, Turb, and TP, important differences between coastal and inland sites were observed (Figure 4) and higher variability was observed in coastal sites relative to inland sites, with the exception of TP. For DO, pH, NO3-N, NO2-N, TC, and FC no clear spatial pattern was found (Figure A5).
Results demonstrate that coastal sites are subject to a higher variability of WT than inland sites with seasonal and spatial components (Figure 3a and Figure 4a, respectively). These differences in WT are likely to be attributed to the distance since water is at 20, 28, and 180 m below ground surface at I3, I2, and I1, respectively. While coastal WT variations happen because they correspond to surface water that receives higher heat fluxes driven by seasonal fluctuation [49], inland WT variation may be associated to its distance below the ground level rather than seasonality, which concurs with Kurylyk et al. study [50] that reported that low groundwater temperatures are related to higher topographic areas and higher groundwater temperatures are derived from high air temperatures in sites near coastal areas. However, our results at I2 differ with Szczucińska and Wasielewski [51] which demonstrated that seasonal fluctuations in temperature can affect a depth of 20–30 m in soils depending on the heat transfer ability of the upper soils or rocks.
As expected, higher variability was also observed for EC, Sal, and Turb at coastal sites relatively to inland sites. For EC (Figure 4b), a relatively big difference with distance from the sea was seen with decreases in the mean of 52.9 mS/cm and 26.9 mS/cm from C4 to C1 at coastal sites, and 3.49 mS/cm and 2.74 mS/cm from I3 to I1 at inland sites. Something similar occurred for Sal (Figure 4c), with decreases in the mean of 32.5 g/L and 19.5 g/L from C4 to C1, and 2.4 g/L and 1.0 g/L from I3 to I1. Turb levels (Figure 4d) also showed to be significantly higher and more dispersed among coastal sites, being C4 the one with the highest turbidity (IQR: 3.2–6.4 NTU). While Turb levels seem to decrease with proximity to the sea at coastal sites, they are likely to increase at inland sites.
In contrast, TP presented higher levels and wider variation in inland sites relative to coastal sites, being I2 the one with the highest TP (IQR: 0.7–1.1 mg/L) (Figure 4e). However, those levels may not represent any harm for human health in waters intended for drinking and domestic purposes. Higher inland TP levels relative to coastal TP levels suggests that inland might be a source of phosphorus to coastal waters. This was suggested in previous studies [52] in which TP might migrate from the aquifer to the coasts by processes involved in inland water discharges into the sea. Nonetheless, this cannot be confirmed by this study and further research should be made.

4. Discussion

Water quality monitoring is an essential part of effective water management at local and regional level. Our results demonstrated that data from the 9-year local monitoring program of Santa Cruz Island provided valuable information that allowed the characterization of physical-chemical and microbial properties of study sites and the analysis of water suitability for designated water uses as described in detail in the sections below. To the best of our knowledge, this is the first study that has investigated water quality variability due to seasonal, anthropogenic and natural processes in the Galapagos Islands.

4.1. Suitability Analysis for Drinking and Domestic Use

Results suggested that water at inland sites may be not suitable for drinking and domestic use due to its Sal, EC, Turb, and FC levels (Figure 5).
Sal and EC levels at inland sites reflect seawater intrusion in the basal aquifer of the island. Seawater intrusion has been reported previously [14] with linkages to human over-pumping of the aquifer. This relationship is not uncommon and has been documented in other places [21,53,54]. However, while it is true that human overexploitation of water induces seawater intrusion, there are other factors involved such as dispersive mixing, tidal effects, surface hydrology, and geological characteristics [53] that needs to be evaluated to elucidate seawater intrusion processes in the island. According to the sensitivity analysis, anthropogenic activities at these sites have limited effect on these parameters between 2007 and 2015 (Table 3). Results also suggested that the variability of these parameters may be related to seasons. Both parameters showed higher variabilities during the garúa relative to invierno season at most sites, and variability in garúa surpassed upper thresholds of background levels. These results suggest that garúa may be influencing EC and Sal parameters and could be a direct consequence of Humboldt current influence, producing higher concentrations than those considered to be natural.
According to WHO guidelines for drinking-water quality [34], increases in Turb levels may affect negatively effective disinfection processes, representing a critical issue if these waters come to get polluted by pathogenic bacteria. Fortunately, at I1 (Pozo Profundo) in which Turb values exceed the 1 NTU limit of WHO guidelines during the invierno season (Figure 2), no coliform contamination has been detected yet. While high contamination risk at I1 is present due to its proximity to the main road, urban settlements and saline intrusion [23], results from the sensitivity analysis suggested that anthropogenic activities has no effect on Turb levels at this site (Table 3). Heavy rainfall events can lead to higher Turb levels due to increased runoff and erosion [10]. This might be the case for I1 since this site is prone to receive the highest infiltration of water originated from rainfalls [13].
FC levels are alarmingly high at I3 (Misión San Francisco), which concurs with previous research [4,12] that indicated that I3 has the highest load on fecal indicator bacteria contamination among all crevices operated by Puerto Ayora’s Municipality for water extraction [4]. This site was originally distributed through the municipal water supply and designated for private desalinization [4,13,26]. Nowadays, I3 is closed for extraction [4]. It was suggested that this contamination is related to its proximity to urban areas and an inefficient sanitary infrastructure that generates leakages from the septic tanks to the basal aquifer [4]. However, our results suggested that population increases in this area between 2007–2015 has a relatively small influence on FC concentrations (Table 3). Nevertheless, effective wastewater management practices for Puerto Ayora are still needed.

4.2. Suitability Analysis for Irrigation

EC and Sal levels (Figure 5) make inland sites unsuitable for irrigation and they seemed to be enhanced during the garúa season. However, these levels are considered natural according to background values calculated for 2007–2015 period (Table 2) with no significant influence from changes in anthropogenic activities (Table 3). High salinity water is responsible for altering soil properties with negative impacts on crops [53]. Since saline water has been used for irrigation at these sites for at least a decade, long term use of saline water for irrigation needs to be addressed. Using a proper irrigation schedule with saline water to ensure sustainability and adoption of salt-tolerant crops (such as quinoa) has been recommended [47]. However, while water management strategies would further facilitate the use of saline waters for irrigation [55], there is the need to understand how exactly salinity affects soils, and thus crops while constantly monitoring salinity levels to ensure they stay within the acceptable range.

4.3. Suitability Analysis for Wild and Aquatic Life in Coastal Waters and Estuaries

In coastal sites, DO, FC, [NO3 + NO2]-N, TP, and Turb are parameters of concern for its different usage purposes (Figure 5). While [NO3 + NO2]-N, TP, and Turb largely exceeded Hawaiian criteria for coastal waters, background levels of these parameters between 2007 and 2015 show higher natural-occurring values than the limit set for coastal waters in Hawaii. This indicates that the use of Hawaii’s water quality rules for evaluating coastal study sites in Galápagos may lead to misunderstandings about the degree of anthropogenic water quality alteration. However, results from the sensitivity analysis suggested that increases in anthropogenic activities strongly influenced changes in the concentration of these parameters in coastal areas and, the closer it gets to the coast, the higher the effect is (Table 3). The influence of human activities on the mentioned parameters is possible since diffuse pollution sources such as agricultural runoffs may be causing some enrichment of nutrients in water and FC. Furthermore, higher concentrations of nutrients (NO3, NO2, and TP) and FC in coastal areas may be related with feces contamination due to untreated wastewater discharges and/or the presence of animal feces. Therefore, results from this analysis suggest that outer thresholds of background levels along with the sensitivity analysis of parameters to changes in population, may serve for defining site-specific limits to control water quality.
Concentrations of DO maintain similar levels across study sites, and those concentrations also reflect prevailing aerobic conditions in the water. However, C1 (Las Ninfas Lagoon) is an exception to this because DO concentrations at this site were lower and highly variable compared to the rest of the study sites (Figure A3a). Low levels of DO may represent a problem for ecosystem preservation, threatening flora and fauna of the estuary, as well as recreational activities [33]. In estuaries, particularly, they lead to stressful conditions for benthic macrofaunal [56]. Additionally, many outliers were observed for DO, which could be a result of simultaneously rapid DO decreases and increases processes derived from bacterial degradation and oxygen release from phytoplankton, respectively [57].
Evidence of high fecal contamination was also found at C1 through the presence of fecal coliform bacteria (Figure A3f). This site was a popular swimming spot for tourists and residents until recreational activities were prohibited due to the presence of human fecal contamination [2,9]. Changes in FC concentrations between 2007–2015 are related with increases in anthropogenic activities at this site (Table 3). It was also observed that higher variability of FC happened during invierno relative to garúa (Figure A4e). This concurs with previous studies that showed evidence of fecal contamination at this site originated from discharges of sewage and agricultural runoff [2,9]. Higher concentration of coliform bacteria proceeding from heavy rainfall events have been reported in other studies [8,9] suggesting that stormwater runoff can be a primary source of coliform bacteria. Furthermore, FC concentrations at C1 largely exceeded Ecuadorian criteria for recreational waters (88% of values exceeding), representing high potential risks of waterborne illness [58]. Additional management practices should be adopted to attenuate fecal contamination at recreational coastal sites.
At coastal sites, NO3-N have higher variability and reach higher values during the cool and dry garúa season (Figure 3) suggesting that rainfall runoff during invierno may tentatively be producing a dilution effect of NO3-N. Increased NO3-N levels were reported previously at Academy bay waters (C2 and C3 area) and Cartago Bay in San Cristobal Island [59]. Human activities such as agriculture, combustion, and sewage discharges may be causing higher-than-natural enrichment of NO3-N, especially at coastal sites. For example, Howarth et al. [52] reported high levels of NO3-N in the U.S. coastal waters resulting from natural origin, and other sources such as that from native fauna, fossil fuel combustion, and N-compounds from fertilizers. Bacteria action on nitrogenous matter concentrated in sewage can also increase concentrations of NO3-N [57]. While our results suggested a strong relationship between anthropogenic activities and NO3-N concentrations during the 9-year period (Table 3), further research needs to be done to be able to explain the processes causing NO3-N enrichment at these sites.

4.4. Environmental Background Levels

Background levels are a useful tool that serves as a reference to assess the level of water affectation since values outside the range of the background can occur due to atypical events. Our results showed that while some parameters (Sal, EC, Turb, [NO3 + NO2]-N, and TP) exceeded water quality criteria, those levels were considered to be varying naturally since 2007 according to the calculated background values. However, our data only captures the 2007–2015 period. Therefore, to reveal real natural conditions unaffected by anthropogenic contributions, data should be compared with historical water sources data and data from pristine sites. Since data for that period is not available, these findings provide valuable information on the water quality status in Santa Cruz Island over the 9-year period analyzed herein and can be used as a baseline for future analysis. Thus, the importance of maintaining long term water quality monitoring to keep track of water status and assess whether alterations are due to impacts of human-induced activities or natural climate-driven processes is clearly highlighted.
Additionally, our results suggested that background levels calculation worked better with parameters at coastal sites than inland sites suggesting that processes underlying coastal sites provide more suitable conditions for the application of the techniques applied herein. While previous research addressed the use of the iterative 2σ and the distribution function techniques in groundwater bodies only, our results suggested that these techniques can also be considered for evaluating coastal water bodies.
Finally, results highlighted the importance of calculating background values for specific sites since the conditions and distribution of the parameters may be different and could lead to misinterpretations when assessing the degree of water alteration [17,22]. This was the case for coastal TP levels that exceed almost 100% with respect to the limits of Hawaii while background calculations suggested that levels were subject to natural increases that surpasses Hawaiian specific criteria for estuaries (3 NTU). Thus, the application of large-scale background values to evaluate local status of water quality is not recommended since local background is constrained by local conditions [17,18,21,22].

4.5. Relative Influence of Population Change on Water Quality Paramertes

Results suggested that the regression analysis method applied for this study could potentially serve as an effective tool to evaluate the significant influence of population change on water quality conditions and estimate the magnitude of its impacts. However, the temporal and spatial correlations are not fully represented in this analysis. To identify the sources of pollutant loads from point and non-point pollution sources (e.g., agriculture, rural domestic, industries, municipal waste, or animal waste), spatial regression methods have been proposed as possible effective methods [60] and should be considered for future analysis. Nevertheless, results from the sensitivity analysis not only supports information obtained from background levels analysis, but provides valuable information that may serve as a scientific basis for water quality management, future research, and decision making [39,40,43] for Santa Cruz Island.

5. Conclusions

A nine year dataset was successfully evaluated providing valuable information of the water quality status in Santa Cruz Island that could serve as a baseline for effective water management in the Pelican Bay watershed. The influence of human activities such as discharges of untreated sewage, combustion, stormwater and agricultural runoffs could cause some enrichment of parameters such as nutrients and coliforms posing a risk for human and ecosystem health. Additional management practices should be considered for recreational coastal sites with high concentrations of coliforms as well as drinking and domestic water with high concentrations of Sal, EC, Turb and FC. Long term use of saline water for irrigation also needs to be addressed. Our results have elucidated that outer thresholds of background levels along with the sensitivity analysis of parameters to changes in population, may serve for defining site-specific limits to control water quality and to distinguish values that are a product of natural processes from anthropogenic activities. However, more specific monitoring is needed to understand the potential water impairment in the Galapagos Islands. Finally, the establishment of long-term water quality monitoring programs in the Galapagos island is a milestone contribution to the preservation and conservation of water resources in the Archipelago.

Author Contributions

Conceptualization, C.M., C.A.G., V.O.-H.; Methodology, C.M., C.A.G.; Formal Analysis, C.M., G.Q., D.L., V.O.-H.; Investigation, C.M., C.A.G., V.O.-H.; Resources, G.Q., D.L.; Data Curation, G.Q., D.L.; Writing—Original Draft Preparation, C.A.G., C.M.; Writing—Review and Editing, C.M., C.A.G., V.O.-H.; Supervision, C.M.; Project Administration, C.M., V.O.-H.; Funding Acquisition, C.M., G.Q., D.L.

Funding

This research project was funded by the 2017 and 2018 Collaboration Grants-USFQ.

Acknowledgments

We would like to express our sincere gratitude to the Galapagos National Park for providing the monitoring dataset that was used in this study. We would especially like to thank Javier López from the Galapagos National Park who performed the sampling and data collection. The authors gratefully acknowledge Kara DiFrancesco for her constructive comments on a previous version of this paper. Publication of this article was funded by the Universidad San Francisco de Quito USFQ Research Publication Fund

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Appendix A

Table A1. Water quality parameters measured by Galapagos National Park during 2007–2015 at all study sites, including their abbreviations, units, analytical methods, and data availability 1.
Table A1. Water quality parameters measured by Galapagos National Park during 2007–2015 at all study sites, including their abbreviations, units, analytical methods, and data availability 1.
ParametersAbbreviationsUnitsTime Coverage
I1I2I3C1C2C3C4
Water temperatureWT°C
TurbidityTurbNTU
Total ColiformsTCMPN 100 mL−1n.d.n.d.n.d.
Fecal ColiformsFCMPN 100 mL−1n.d.n.d.n.d.
SalinitySalg L−1
Dissolved oxygenDOmg L−1
Electrical ConductivityECmS cm−1
pHpHpH unit
NitrateNO3-Nmg L−1*******
NitriteNO2-Nμg L−1*******
Total PhosphorusTPmg L−1*******
Note: 1 Retrieved from the original planning document [30]. ✓ 91% time coverage: No data for March 2008 and February, March, May, June, July, September, October, and November 2009 (n = 99). * 30%—time coverage: Available data for January, April, August, and December of each year (n = 36). n.d.: no data.
Figure A1. Background concentrations for EC at I1 using: (a) the iterative 2σ technique (background levels: 2.1–3.2 mS/cm; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.094, T = 0.058), and; (b) the distribution function (background levels: 2–3.3 mS/cm; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.085, T = 0.130).
Figure A1. Background concentrations for EC at I1 using: (a) the iterative 2σ technique (background levels: 2.1–3.2 mS/cm; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.094, T = 0.058), and; (b) the distribution function (background levels: 2–3.3 mS/cm; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.085, T = 0.130).
Water 11 00937 g0a1
Figure A2. Background concentrations for NO3-N at C1 using: (a) the iterative 2σ technique (background levels: 0–0.1 mg/L; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.152, T = 0.155), and (b) the distribution function (background levels: 0–0.1 mg/L; Lilliefors test statistic [α = 0.5]: Tcrit. = 0.137, T = 0.086).
Figure A2. Background concentrations for NO3-N at C1 using: (a) the iterative 2σ technique (background levels: 0–0.1 mg/L; Lilliefors test statistic [α = 0.05]: Tcrit. = 0.152, T = 0.155), and (b) the distribution function (background levels: 0–0.1 mg/L; Lilliefors test statistic [α = 0.5]: Tcrit. = 0.137, T = 0.086).
Water 11 00937 g0a2
Figure A3. Seasonal variation: Invierno (January–May) and garúa (June–December) variations for (a) DO coastal, (b) DO inland, (c) pH coastal, (d) pH inland, (e) TP coastal, (f) TP inland, (g) NO2-N inland, and (h) NO2-N coastal. The scale for coastal TO (C1, C2, C3, and C4) differs from the scale for inland TP (I1, I2, and I3).
Figure A3. Seasonal variation: Invierno (January–May) and garúa (June–December) variations for (a) DO coastal, (b) DO inland, (c) pH coastal, (d) pH inland, (e) TP coastal, (f) TP inland, (g) NO2-N inland, and (h) NO2-N coastal. The scale for coastal TO (C1, C2, C3, and C4) differs from the scale for inland TP (I1, I2, and I3).
Water 11 00937 g0a3
Figure A4. Seasonal variation. Invierno (January–May) and garúa (June–December) variations for (a) EC coastal, (b) EC inland, (c) Sal coastal, (d) Sal inland, and (e) FC coastal (left) TC inland (right), and (f) FC inland. The scale at coastal sites differs from the scale at inland sites.
Figure A4. Seasonal variation. Invierno (January–May) and garúa (June–December) variations for (a) EC coastal, (b) EC inland, (c) Sal coastal, (d) Sal inland, and (e) FC coastal (left) TC inland (right), and (f) FC inland. The scale at coastal sites differs from the scale at inland sites.
Water 11 00937 g0a4
Figure A5. Spatial variation for (a) DO, (b) pH, (c) NO3-N, (d) NO2-N, (e) TC, and (f) FC. Study sites are ordered by proximity to the sea, being C2 the closest to the sea and I2 the furthest. The upper and lower quartiles are represented by the ends of the box, the median is the horizontal line between box ends, the mean is marked by “x”, the whiskers indicate variability outside the upper and lower quartiles and, finally, points outside whiskers are considered as outliers. TC and FC data only exists for C1 at the coastal sites.
Figure A5. Spatial variation for (a) DO, (b) pH, (c) NO3-N, (d) NO2-N, (e) TC, and (f) FC. Study sites are ordered by proximity to the sea, being C2 the closest to the sea and I2 the furthest. The upper and lower quartiles are represented by the ends of the box, the median is the horizontal line between box ends, the mean is marked by “x”, the whiskers indicate variability outside the upper and lower quartiles and, finally, points outside whiskers are considered as outliers. TC and FC data only exists for C1 at the coastal sites.
Water 11 00937 g0a5

References

  1. Michalak, A.M. Study role of climate change in extreme threats to water quality. Nature 2016, 535, 349–350. [Google Scholar] [CrossRef]
  2. Whitehead, P.G.; Wilby, R.L.; Battarbee, R.W.; Kernan, M.; Wade, A.J. A review of the potential impacts of climate change on surface water quality. Hydrol. Sci. J. 2009, 54, 101–123. [Google Scholar] [CrossRef] [Green Version]
  3. Perkins, D.L.; Kann, J.; Scoppettone, G.G. The Role of Poor Water Quality and Fish Kills in the Decline of Endangered Lost River and Shortnose Suckers in Upper Klamath Lake; U.S. Geological Survey: Klamath Falls, OR, USA, 2000; p. 40. [Google Scholar]
  4. Liu, J.; d’Ozouville, N. Water Contamination in Puerto Ayora: Applied Interdisciplinary Research Using Escherichia Coli as an Indicator Bacteria; Galapagos Report 2011–2012; GNPD, GCREG, CDF and GC: Puerto Ayora, Galápagos, Ecuador, 2011; pp. 76–83. [Google Scholar]
  5. Walsh, S.J.; McCleary, A.L.; Heumann, B.W.; Brewington, L.; Raczkowski, E.J.; Mena, C.F. Community Expansion and Infrastructure Development: Implications for Human Health and Environmental Quality in the Galápagos Islands of Ecuador. J. Lat. Am. Geogr. 2010, 9, 137–159. [Google Scholar] [CrossRef]
  6. Overbey, K.N.; Hatcher, S.M.; Stewart, J.R. Water quality and antibiotic resistance at beaches of the Galápagos Islands. Front. Environ. Sci. 2015, 3. [Google Scholar] [CrossRef]
  7. Jun, X.; Shubo, C.; Xiuping, H.; Rui, X.; Xiaojie, L. Potential Impacts and Challenges of Climate Change on Water Quality and Ecosystem: Case Studies in Representative Rivers in China. J. Resour. Ecol. 2010, 1, 31–35. [Google Scholar]
  8. Kistemann, T.; Classen, T.; Koch, C.; Dangendorf, F.; Fischeder, R.; Gebel, J.; Vacata, V.; Exner, M. Microbial Load of Drinking Water Reservoir Tributaries during Extreme Rainfall and Runoff. Appl. Environ. Microbiol. 2002, 68, 2188–2197. [Google Scholar] [CrossRef] [Green Version]
  9. Pednekar, A.M.; Grant, S.B.; Jeong, Y.; Poon, Y.; Oancea, C. Influence of Climate Change, Tidal Mixing, and Watershed Urbanization on Historical Water Quality in Newport Bay, a Saltwater Wetland and Tidal Embayment in Southern California. Environ. Sci. Technol. 2005, 39, 9071–9082. [Google Scholar] [CrossRef] [PubMed]
  10. Delpla, I.; Jung, A.-V.; Baures, E.; Clement, M.; Thomas, O. Impacts of climate change on surface water quality in relation to drinking water production. Environ. Int. 2009, 35, 1225–1233. [Google Scholar] [CrossRef]
  11. Alava, J.J.; Palomera, C.; Bendell, L.; Ross, P.S. Pollution as an Emerging Threat for the Conservation of the Galapagos Marine Reserve: Environmental Impacts and Management Perspectives. In The Galapagos Marine Reserve; Denkinger, J., Vinueza, L., Eds.; Springer International Publishing: Cham, Switzerland, 2014; pp. 247–283. ISBN 978-3-319-02768-5. [Google Scholar]
  12. López, J.; Rueda, D. Water Quality Monitoring System in Santa Cruz, San Cristóbal, and Isabela; Galapagos Report 2009–2010; CDF, GNP & CGG: Puerto Ayora, Galapagos, Ecuador, 2010; pp. 103–107. [Google Scholar]
  13. CISPDR. Plan Reginal de los Recursos Hídricos de las Islas Galápagos; Changjiang Institute of Survey Planning Design and Research (CISPDR): London, UK, 2015. [Google Scholar]
  14. Violette, S.; d’Ozouville, N.; Pryet, A.; Deffontaines, B.; Fortin, J.; Adelinet, M. Hydrogeology of the Galápagos Archipelago: An Integrated and Comparative Approach Between Islands. Galápagos A Nat. Lab. Earth Sci. Geophys. Monogr. 2014, 204, 167–185. [Google Scholar]
  15. Palmer, M.; Armijos, E.; Moir, F.C. Selection of an ocean outfall site for the town of Puerto Ayora in the Galapagos Islands. In Proceedings of the XXVIII Congreso Interamericano de Ingeniería Sanitaria y Ambiental, Cancún, México, 27–31 October 2002. [Google Scholar]
  16. Malmaeus, J.M.; Blenckner, T.; Markensten, H.; Persson, I. Lake phosphorus dynamics and climate warming: A mechanistic model approach. Ecol. Model. 2006, 190, 1–14. [Google Scholar] [CrossRef]
  17. Sun, L.; Gui, H.; Peng, W. Geostatistical analyses of iron in shallow groundwater and their application for establishing an environmental baseline: A case study. Water Pract. Technol. 2014, 9, 450. [Google Scholar] [CrossRef]
  18. Gałuszka, A.; Migaszewski, Z.M.; Dołęgowska, S.; Michalik, A.; Duczmal-Czernikiewicz, A. Geochemical background of potentially toxic trace elements in soils of the historic copper mining area: A case study from Miedzianka Mt., Holy Cross Mountains, south-central Poland. Environ. Earth Sci. 2015, 74, 4589–4605. [Google Scholar] [CrossRef]
  19. Water Resources Glossaries. Available online: https://water.usgs.gov/water-basics_glossary.html#B (accessed on 7 April 2019).
  20. Nakić, Z.; Posavec, K.; Bačani, A. A Visual Basic Spreadsheet Macro for Geochemical Background Analysis. Ground Water 2007, 45, 642–647. [Google Scholar] [CrossRef] [PubMed]
  21. Urresti-Estala, B.; Carrasco-Cantos, F.; Vadillo-Pérez, I.; Jiménez-Gavilán, P. Determination of background levels on water quality of groundwater bodies: A methodological proposal applied to a Mediterranean River basin (Guadalhorce River, Málaga, southern Spain). J. Environ. Manag. 2013, 117, 121–130. [Google Scholar] [CrossRef] [PubMed]
  22. Zgłobicki, W.; Lata, L.; Plak, A.; Reszka, M. Geochemical and statistical approach to evaluate background concentrations of Cd, Cu, Pb and Zn (case study: Eastern Poland). Environ. Earth Sci. 2011, 62, 347–355. [Google Scholar] [CrossRef]
  23. d’Ozouville, N.; Deffontaines, B.; Benveniste, J.; Wegmüller, U.; Violette, S.; de Marsily, G. DEM generation using ASAR (ENVISAT) for addressing the lack of freshwater ecosystems management, Santa Cruz Island, Galapagos. Remote Sens. Environ. 2008, 112, 4131–4147. [Google Scholar] [CrossRef]
  24. d’Ozouville, N. Manejo de Recursos Hídricos: Caso de la Cuenca de Pelican Bay; Informe Galápagos 2007–2008; DPNG, CGREG, FCD y GC: Puerto Ayora, Ecuador, 2007. [Google Scholar]
  25. Pryet, A.; Domínguez, C.; Tomai, P.F.; Chaumont, C.; d’Ozouville, N.; Villacís, M.; Violette, S. Quantification of cloud water interception along the windward slope of Santa Cruz Island, Galapagos (Ecuador). Agric. For. Meteorol. 2012, 161, 94–106. [Google Scholar] [CrossRef] [Green Version]
  26. Stumpf, C.H.; González, R.A.; Noble, R.T. Chapter 10: Investigating the Coastal Water Quality of the Galapagos Islands, Ecuador. In Science and Conservation in the Galapagos Islands: Frameworks & Perspectives; Walsh, S.J., Mena, C.F., Eds.; Springer: New York, NY, USA, 2013; pp. 173–184. ISBN 2195-1055. [Google Scholar]
  27. Trueman, M.; d’Ozouville, N. Characterizing the Galapagos terrestrial climate in the face of global climate change. Galapagos Res. 2010, 67, 12. [Google Scholar]
  28. Adelinet, M.; Fortin, J.; d’Ozouville, N.; Violette, S. The relationship between hydrodynamic properties and weathering of soils derived from volcanic rocks, Galapagos Islands (Ecuador). Environ. Geol. 2008, 56, 45–58. [Google Scholar] [CrossRef]
  29. d’Ozouville, N. Etude du fonctionnement hydrologique dans les îles: Caractérisation d’un milieu volcani ue insulaire et préalable à la gestion de la ressource. Ph.D. Thesis, Université Pierre et Marie Curie, Paris, France, 2007. [Google Scholar]
  30. Zapata, F. Planificación del Monitoreo de la Contaminación Costera y Calidad del Agua en Puerto Ayora, Galápagos; Monitoreo de Ambientes Marinos; Parque Nacional Galápagos-JICA: Santa Cruz, CA, USA, 2004. [Google Scholar]
  31. Standard Methods for the Examination of Water and Wastewater, 22nd ed.; Rice, E.W.; American Public Health Association (Eds.) American Public Health Association: Washington, DC, USA, 2012; ISBN 978-0-87553-013-0. [Google Scholar]
  32. Snell, H.; Rea, S. The 1997–1998 El Niño in Galapagos: Can 34 years of data estimate 120 years of pattern. In Climate Change Vulnerability Assessment of the Galápagos Islands; Larrea Oña, I., Di Carlo, G., Eds.; WWF: Grand, Switzerland; Conservation International: Arlington, VA, USA, 2011; p. 10. ISBN 978-9942-03-454-0. [Google Scholar]
  33. MAE. Acuerdo 097-A, Anexo 1 del Libro VI del Texto Unificado de Legislación Secundaria del Ministerio del Ambiente: Norma de Calidad Ambiental y de Descarga de Efluentes al Recurso Agua; MAE: Quito, Ecuador, 2015. [Google Scholar]
  34. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011.
  35. Water Quality Standards. In Hawaii Administrative Rules Title 11, Chapter 54—Water Quality Standards; 2015. Available online: epa.gov (accessed on 24 April 2019).
  36. MacDonald, G.; Abbott, A.; Peterson, F.L. Volcanoes in the Sea; The Geology of Hawaii; University of Hawaii Press: Honolulu, HI, USA, 1983; p. 517. [Google Scholar]
  37. Benson, B.B.; Krause, D. The Concentration and Isotopic Fractionation of Gases Dissolved in Freshwater in Equilibrium with the Atmosphere. 1. Oxygen. Limnol. Oceanogr. 1984, 25, 662–671. [Google Scholar] [CrossRef]
  38. Matschullat, J.; Ottenstein, R.; Reimann, C. Geochemical background—Can we calculate it? Environ. Geol. 2000, 39, 990–1000. [Google Scholar] [CrossRef]
  39. Bu, H.; Meng, W.; Zhang, Y.; Wan, J. Relationships between land use patterns and water quality in the Taizi River basin, China. Ecol. Indic. 2014, 41, 187–197. [Google Scholar] [CrossRef]
  40. Shi, W.; Xia, J.; Zhang, X. Influences of anthropogenic activities and topography on water quality in the highly regulated Huai River basin, China. Environ. Sci. Pollut. Res. 2016, 23, 21460–21474. [Google Scholar] [CrossRef] [PubMed]
  41. Halstead, J.A.; Kliman, S.; Berheide, C.W.; Chaucer, A.; Cock-Esteb, A. Urban stream syndrome in a small, lightly developed watershed: A statistical analysis of water chemistry parameters, land use patterns, and natural sources. Environ. Monit. Assess. 2014, 186, 3391–3414. [Google Scholar] [CrossRef] [PubMed]
  42. Ahamad, K.U.; Raj, P.; Barbhuiya, N.H.; Deep, A. Surface Water Quality Modeling by Regression Analysis and Artificial Neural Network. In Advances in Waste Management; Kalamdhad, A., Singh, J., Dhamodharan, K., Eds.; Springer: Singapore, 2019; pp. 215–230. [Google Scholar]
  43. Tasker, G.D.; Driver, N.E. Nationwide Regression Models for Predicting Urban Runoff Water Quality at Unmonitored Sites. JAWRA J. Am. Water Res. Assoc. 1988, 24, 1091–1101. [Google Scholar] [CrossRef]
  44. Christensen, V.; Xiaodong, J.; Ziegler, A. Regression Analysis and Real-Time Water-Quality Monitoring to Estimate Constituent Concentrations, Loads, and Yields in the Little Arkansas River, South-Central Kansas, 1995–1999; U.S. Geological Survey & U.S. Department of Interior: Washington, DC, USA, 2000. [Google Scholar]
  45. Instituto Nacional de Estatisticas y Censos (INEC). Censo de Población Y Vivienda 2011: Quito, Ecuador. Proyección a 2018. Available online: http://www.ecuadorencifras.gob.ec/ (accessed on 8 July 2018).
  46. Ministerio de Turismo Observatorio de Turismo de Galápagos. Available online: https://www.observatoriogalapagos.gob.ec/ (accessed on 1 April 2019).
  47. Kovač, Z.; Nakić, Z.; Pavlić, K. Influence of groundwater quality indicators on nitrate concentrations in the Zagreb aquifer system. Geol. Croat. 2017, 70, 93–103. [Google Scholar] [CrossRef]
  48. New Zealand’s Environmental Reporting Series: Environment Aotearoa 2015; Ministry for the Environment & Statistics New Zealand: Wellington, New Zealand, 2015.
  49. Kaplan, D.M.; Largier, J.L.; Navarrete, S.; Guiñez, R.; Castilla, J.C. Large diurnal temperature fluctuations in the nearshore water column. Estuar. Coast. Shelf Sci. 2003, 57, 385–398. [Google Scholar] [CrossRef]
  50. Kurylyk, B.L.; Bourque, C.P.-A.; MacQuarrie, K.T.B. Potential surface temperature and shallow groundwater temperature response to climate change: An example from a small forested catchment in east-central New Brunswick (Canada). Hydrol. Earth Syst. Sci. 2013, 17, 2701–2716. [Google Scholar] [CrossRef]
  51. Szczucińska, A.M.; Wasielewski, H. Seasonal Water Temperature Variability of Springs from Porous Sediments in Gryżynka Valley, Western Poland. Quaest. Geogr. 2013, 32. [Google Scholar] [CrossRef]
  52. Howarth, R.W.; Sharpley, A.; Walker, D. Sources of nutrient pollution to coastal waters in the United States: Implications for achieving coastal water quality goals. Estuaries 2002, 25, 656–676. [Google Scholar] [CrossRef]
  53. Werner, A.D.; Bakker, M.; Post, V.E.A.; Vandenbohede, A.; Lu, C.; Ataie-Ashtiani, B.; Simmons, C.T.; Barry, D.A. Seawater intrusion processes, investigation and management: Recent advances and future challenges. Adv. Water Res. 2013, 51, 3–26. [Google Scholar] [CrossRef] [Green Version]
  54. Molinari, A.; Chidichimo, F.; Straface, S.; Guadagnini, A. Assessment of natural background levels in potentially contaminated coastal aquifers. Sci. Total Environ. 2014, 476–477, 38–48. [Google Scholar] [CrossRef]
  55. Wang, Q.; Huo, Z.; Zhang, L.; Wang, J.; Zhao, Y. Impact of saline water irrigation on water use efficiency and soil salt accumulation for spring maize in arid regions of China. Agric. Water Manag. 2016, 163, 125–138. [Google Scholar] [CrossRef]
  56. Diaz, R.J.; Rosenberg, R. Marine benthic hypoxia: A review of its ecological effects and the behavioral responses of benthic macrofauna. Oceanogr. Morine Biol. An Annu. Rev. 1995, 33, 245–303. [Google Scholar]
  57. Parameters of Water Quality: Interpretation and Standards; Environmental Protection Agency: Dublin, Ireland, 2001; p. 133.
  58. Colford, J.M.; Wade, T.J.; Schiff, K.C.; Wright, C.C.; Griffith, J.F.; Sandhu, S.K.; Burns, S.; Sobsey, M.; Lovelace, G.; Weisberg, S.B. Water Quality Indicators and the Risk of Illness at Beaches with Nonpoint Sources of Fecal Contamination. Epidemiology 2007, 18, 27–35. [Google Scholar] [CrossRef] [PubMed]
  59. Werdeman, J. Effects of Populated Towns on Water Quality in Neighboring Galapagos Bays. Available online: http://courses.washington.edu/ocean443/2005_6/DraftProposals/Joni.pdf (accessed on 24 April 2019).
  60. Yang, X.; Liu, Q.; Luo, X.; Zheng, Z. Spatial Regression and Prediction of Water Quality in a Watershed with Complex Pollution Sources. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Pelican Bay Watershed: Urban, Rural and Agricultural Areas, Two Weather Stations (purple triangle), Four Coastal Sites (blue circle) and Three Inland Sites (red square). Coastal sites are: Las Ninfas Lagoon (C1), Academy Bay in (C2), Academy Bay out (C3) and Muelle de los Pescadores (C4). Inland sites are: Pozo Profundo (I1), INGALA (I2), Misión San Francisco (I3). CDF1 and CDF2 represent the Charles Darwin Foundation weather stations located in the Charles Darwin Research Station and Bellavista, respectively. Top inset: Galapagos Archipelago. Bottom inset: Santa Cruz Island.
Figure 1. Pelican Bay Watershed: Urban, Rural and Agricultural Areas, Two Weather Stations (purple triangle), Four Coastal Sites (blue circle) and Three Inland Sites (red square). Coastal sites are: Las Ninfas Lagoon (C1), Academy Bay in (C2), Academy Bay out (C3) and Muelle de los Pescadores (C4). Inland sites are: Pozo Profundo (I1), INGALA (I2), Misión San Francisco (I3). CDF1 and CDF2 represent the Charles Darwin Foundation weather stations located in the Charles Darwin Research Station and Bellavista, respectively. Top inset: Galapagos Archipelago. Bottom inset: Santa Cruz Island.
Water 11 00937 g001
Figure 2. Seasonal variation. Invierno (January–May) and garúa (June–December) variations for (a) WT coastal, (b) WT inland, (c) Turb coastal, (d) Turb inland, (e) NO3-N coastal, (f) NO3-N inland, (g) TC coastal (left) TC inland (right), and (h) TC inland. The scale at coastal sites for NO3-N and TC differ from the scale at inland sites.
Figure 2. Seasonal variation. Invierno (January–May) and garúa (June–December) variations for (a) WT coastal, (b) WT inland, (c) Turb coastal, (d) Turb inland, (e) NO3-N coastal, (f) NO3-N inland, (g) TC coastal (left) TC inland (right), and (h) TC inland. The scale at coastal sites for NO3-N and TC differ from the scale at inland sites.
Water 11 00937 g002
Figure 3. Mean water temperature (WT) from 2007 to 2015 (solid gray line), mean air temperature (dotted line), sea surface temperature (dashed line), and rainfall (bars) at coastal sites: (a) Muelle de los Pescadores (C4), (b) Academy Bay out (C3), (c) Academy Bay in (C2), and (d) Las Ninfas (C1); and inland sites: (e) Misión San Francisco (I3), (f) INGALA (I2), and (g) Pozo Profundo (I1).
Figure 3. Mean water temperature (WT) from 2007 to 2015 (solid gray line), mean air temperature (dotted line), sea surface temperature (dashed line), and rainfall (bars) at coastal sites: (a) Muelle de los Pescadores (C4), (b) Academy Bay out (C3), (c) Academy Bay in (C2), and (d) Las Ninfas (C1); and inland sites: (e) Misión San Francisco (I3), (f) INGALA (I2), and (g) Pozo Profundo (I1).
Water 11 00937 g003
Figure 4. Spatial variation for (a) WT, (b) EC, (c) Sal, (d) Turb, and (e) TP. Study sites are ordered by proximity to the sea, being C2 the closest to the sea and I2 the furthest. The upper and lower quartiles are represented by the ends of the box, the median is the horizontal line between box ends, the mean is marked by “x”, the whiskers indicate variability outside the upper and lower quartiles and, finally, points outside whiskers are considered as outliers.
Figure 4. Spatial variation for (a) WT, (b) EC, (c) Sal, (d) Turb, and (e) TP. Study sites are ordered by proximity to the sea, being C2 the closest to the sea and I2 the furthest. The upper and lower quartiles are represented by the ends of the box, the median is the horizontal line between box ends, the mean is marked by “x”, the whiskers indicate variability outside the upper and lower quartiles and, finally, points outside whiskers are considered as outliers.
Water 11 00937 g004
Figure 5. Parameters exceeding guidelines criteria at inland sites (I1, I2, I3) and coastal sites (C1–C4). Bottom left inset: Santa Cruz location in the Galápagos Archipelago (top) and study areas location in Santa Cruz Island (bottom). Bottom middle inset: Location for each parameter within the inset figures showing the results inside the map: Sal—Salinity in the top left side (circles in purple tones), Turb—Turbidity in the top right side (circles in grey tones), DO—Dissolved Oxygen in the bottom left side (circles with multiple colors), TC—Total Coliforms (brown rectangle) and FC—Fecal Coliforms (red rectangle) in the bottom right side. The range of values and color tones are described in the figure legend in the bottom right inset.
Figure 5. Parameters exceeding guidelines criteria at inland sites (I1, I2, I3) and coastal sites (C1–C4). Bottom left inset: Santa Cruz location in the Galápagos Archipelago (top) and study areas location in Santa Cruz Island (bottom). Bottom middle inset: Location for each parameter within the inset figures showing the results inside the map: Sal—Salinity in the top left side (circles in purple tones), Turb—Turbidity in the top right side (circles in grey tones), DO—Dissolved Oxygen in the bottom left side (circles with multiple colors), TC—Total Coliforms (brown rectangle) and FC—Fecal Coliforms (red rectangle) in the bottom right side. The range of values and color tones are described in the figure legend in the bottom right inset.
Water 11 00937 g005
Table 1. Percentage of parameter measurements from study sites that exceed water quality guidelines. Percentage (%) of exceedance represents the fraction of a certain parameter that exceed a specified water quality criteria value. FC: fecal coliforms; Turb: turbidity; Sal: salinity; EC: electrical conductivity; DO: dissolved oxygen; TP: total phosphorus; WT: water temperature.
Table 1. Percentage of parameter measurements from study sites that exceed water quality guidelines. Percentage (%) of exceedance represents the fraction of a certain parameter that exceed a specified water quality criteria value. FC: fecal coliforms; Turb: turbidity; Sal: salinity; EC: electrical conductivity; DO: dissolved oxygen; TP: total phosphorus; WT: water temperature.
ParametersGuideline/Criteria Values% of Exceedance
I1I2I3
(A) Drinking and domestic use
FC a1000 MPN/100 mL--53
pH a6–9--2
Sal bUnpalatable: >1 g/L2690100
Turb b<1 NTU1864
(B) Irrigation water c
pH6–9--2
FC1000 MPN/100 mL--53
ECModerate: 0.7–3 mS/cm 835133
Severe: >3 mS/cm174967
SalModerate: 0.45–2 g/L99766
Severe: >2 g/L02494
ParametersGuideline/Criteria Values% of Exceedance
C1C2C3C4
(A) Preservation of wild and aquatic life in marine waters and estuaries d
DO>60% saturation15331
pH6.5–9.52-11
TurbShall not deviate more than 5% from background levels416714
(B) Primary and Secondary contact recreational water e
FC200 MPN/100 mL88n.d.n.d.n.d.
TC2000 MPN/100 mL4n.d.n.d.n.d.
DO>80% saturation7727158
pH6.5–8.3 (primary)3231
6–9 (secondary)2--1
(C) Specific criteria for estuaries and embayment f
[NO3+NO2]-N25 μg/L (estuaries)10010099-
20 μg/L (embayments)---99
TP50 μg/L971008699
Turb3 NTU8204099
pH7-8.6 and shall not deviate more than 0.5 units from background levels10---
shall not deviate more than 0.5 units from 8.1-322018
DO>75% sat.6914106
WTshall not deviate more than 1 °C from background levels5211
Salshall not deviate more than 10% C from background levels131599
Note: a Acuerdo 097-A, Anexo 1 del Libro VI del Texto Unificado de Legislación Secundaria del Ministerio de Ambiente, Tabla 1, 2015; b Guidelines for Drinking-water Quality; 4th ed.; World Health Organization, 2011; c Acuerdo 097-A, Anexo 1 del Libro VI del Texto Unificado de Legislación Secundaria del Ministerio de Ambiente, Tabla 4, 2015; d Acuerdo 097-A, Anexo 1 del Libro VI del Texto Unificado de Legislación Secundaria del Ministerio de Ambiente, Tabla 2, 2015; e Acuerdo 097-A, Anexo 1 del Libro VI del Texto Unificado de Legislación Secundaria del Ministerio de Ambiente, Tabla 6, 2015; f Specific criteria for estuaries and embayment from water quality standards of Hawaii Administrative Rules Title 11, Chapter 54.
Table 2. Background levels at Inland Sites (I1, I2, I3) and Coastal Sites (C1, C2, C3, C4) using the iterative 2σ technique (2σ) and the distribution function (d.f.) technique.
Table 2. Background levels at Inland Sites (I1, I2, I3) and Coastal Sites (C1, C2, C3, C4) using the iterative 2σ technique (2σ) and the distribution function (d.f.) technique.
(a) Inland Sites
ParametersI1I2I3
d.f.d.f.d.f.
WT (°C)22.4–24.318.8–28.023.4–24.223.0–24.623.8–24.422.9–25.5
Turb (NTU)0.2–0.70.2–0.80.1–0.6 *0.1–0.60.1–0.6 *0.1–0.7 *
Sal (g/L)0.7–1.2 *0.7–1.11.1–1.80.8–2.02.1–2.62.0–2.6
DO (mg/L)7.2–8.36.6–9.07.6–8.37.2–8.87.6–8.3 *7.2–8.8
EC (mS/cm)2.1–3.2 *2.0–3.32.2–3.7 *2.2–3.7 *2.0–4.7 *1.9–5.1 *
pH7.0–8.36.8–8.6 *7.4–8.27.2–8.47.5–8.27.5–8.2 *
NO3-N (mg/L)0.0–0.1 *0.0–0.1 *0.0–0.1 *0.0–0.1 *0.1–0.4 *0.1–0.4 *
NO2-N (μg/L)2.6–8.13.2–6.82.6–8.02.9–7.12.9–7.4 *2.7–7.3
TP (mg/L)0.0–1.3 *0.0–1.4 *0.0–1.6 *0.0–1.7 *0.0–1.1 *0.0–1.3 *
(b) Coastal Sites
ParametersC1C2C3C4
d.f.d.f.d.f.d.f.
WT (°C)23.1–26.521.7–27.721.8–27.2 *21.5–27.9 *21.7–27.3 *21.3–27.9 *21.9–27.6 *21.3–28.3 *
Turb (NTU)0.0–2.3 *0.1–2.3 *0.3–2.9 *0.4–3.2 *0.0–4.70.0–5.21.2–7.5 *1.2–8.4 *
Sal (g/L)16.4–21.6 *13.3–25.116.8–24.1 *16.7–25.9 *32.4–33.2 *28.0–37.632.4–33.2 *28.0–37.6
DO (mg/L)4.7–7.44.5–7.96.5–7.85.5–8.57.5–8.37.5–8.3 *7.5–8.26.1–9.5
EC (mS/cm)21.0–29.4 *20.2–30.326.4–31.7 *25.6–33.5 *47.1–53.9 *39.5–61.146.8–55.0 *39.5–61.1
pH7.0–8.1 *6.4–8.67.4–8.27.1–8.57.5–8.2 *7.2–8.47.5–8.27.1–8.5
NO3-N (mg/L)0.0–0.10.0–0.1 *0.0–0.10.0–0.1 *0.0–0.1 *0.0–0.1 *0.0–0.1 *0.0–0.1 *
NO2-N (μg/L)3.9–9.1 *3.7–10.3 *3.5–8.8 *3.9–8.1 *3.3–9.5 *3.4–10.6 *3.4–10.0 *3.7–10.3 *
TP (mg/L)0.0–0.2 *0.0–0.2 *0.0–0.2 *0.0–0.3 *0.0–0.5 *0.0–0.6 *0.0–0.5 *0.0–0.6 *
Note: The * indicates variables that passed the Lilliefors test. Values shaded in grey represent the best estimations according to the nature and distribution of data which means that the data fits a normal distribution and that the techniques applied herein are appropriate for such parameters.
Table 3. Sensitivity Analysis for Inland and Coastal Sites: Represented by the Strength of the Relationship (R2) and the Size of the Effect (slope, ß1) between water quality parameters and forcing variables for number of habitants in the island taking into account both, number of tourists and local population.
Table 3. Sensitivity Analysis for Inland and Coastal Sites: Represented by the Strength of the Relationship (R2) and the Size of the Effect (slope, ß1) between water quality parameters and forcing variables for number of habitants in the island taking into account both, number of tourists and local population.
Inland SitesCoastal Sites
SiteParameterß1ßoR2SiteParameterß1ßoR2
I1TEMP−0.920.040.29C1TEMP0.140.030.01
EC0.100.030.16EC0.050.030.07
PH0.200.030.02PH0.580.030.58
TURB0.000.320.00TURB0.020.030.19
SAL0.070.040.06SAL0.190.030.24
DO−0.470.040.25DO0.100.030.13
NO2-N0.020.030.00NO2-N0.100.040.17
NO3-N−0.010.030.01NO3-N0.000.030.00
TP−0.030.020.07TP0.010.030.00
TC0.000.030.03TC0.010.030.03
FC0.000.030.03FC0.060.030.29
I2TEMP0.040.030.00C2TEMP0.080.030.00
EC0.150.040.30EC0.190.040.72
PH0.700.030.21PH0.840.030.43
TURB0.360.030.14TURB0.080.040.38
SAL0.000.030.00SAL0.090.030.05
DO0.870.040.47DO0.200.030.46
NO2-N0.010.030.01NO2-N0.080.030.08
NO3-N0.010.030.00NO3-N0.040.030.38
TP0.020.030.02TP0.020.040.72
TC0.010.030.04TCNANANA
FC0.000.320.00FCNANANA
I3TEMP0.600.030.02C3TEMP0.280.030.09
EC0.090.040.11EC0.020.030.10
PH0.450.030.03PH1.580.030.21
TURB0.040.040.30TURB0.030.030.02
SAL0.110.030.09SAL3.280.040.79
DO0.420.030.04DO0.230.040.22
NO2-N0.090.030.12NO2-N0.050.030.30
NO3-N0.020.030.29NO3-N0.070.040.52
TP0.170.040.26TP0.080.030.42
TC0.050.040.21TCNANANA
FC0.100.040.34FCNANANA
C4TEMP0.350.040.14
EC0.020.030.09
PH0.570.030.15
TURB0.160.050.46
SAL0.470.030.05
DO0.500.040.22
NO2-N0.040.030.13
NO3-N0.080.040.49
TP0.020.030.02
TCNANANA
FCNANANA
Note: Values with the strongest relationships are highlighted in dark gray, values with weak relationships are highlighted in light gray.

Share and Cite

MDPI and ACS Style

Mateus, C.; Guerrero, C.A.; Quezada, G.; Lara, D.; Ochoa-Herrera, V. An Integrated Approach for Evaluating Water Quality between 2007–2015 in Santa Cruz Island in the Galapagos Archipelago. Water 2019, 11, 937. https://doi.org/10.3390/w11050937

AMA Style

Mateus C, Guerrero CA, Quezada G, Lara D, Ochoa-Herrera V. An Integrated Approach for Evaluating Water Quality between 2007–2015 in Santa Cruz Island in the Galapagos Archipelago. Water. 2019; 11(5):937. https://doi.org/10.3390/w11050937

Chicago/Turabian Style

Mateus, Cristina, Christian A. Guerrero, Galo Quezada, Daniel Lara, and Valeria Ochoa-Herrera. 2019. "An Integrated Approach for Evaluating Water Quality between 2007–2015 in Santa Cruz Island in the Galapagos Archipelago" Water 11, no. 5: 937. https://doi.org/10.3390/w11050937

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop