Next Article in Journal
Performance of Conservation Techniques for Semiarid Environments: Field Observations with Caatinga, Mulch, and Cactus Forage Palma
Next Article in Special Issue
Stable Isotopes of Precipitation in China: A Consideration of Moisture Sources
Previous Article in Journal
Assessing Groundwater Mineralization Process, Quality, and Isotopic Recharge Origin in the Sahel Region in Africa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Local Meteoric Water Line of Northern Chile (18° S–30° S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation

1
Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze, 43124 Parma, Italy
2
Servicio Nacional de Geología y Minería, Avenida Santa María #0104, 7520405 Santiago, Chile
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 791; https://doi.org/10.3390/w11040791
Submission received: 13 March 2019 / Revised: 11 April 2019 / Accepted: 11 April 2019 / Published: 16 April 2019
(This article belongs to the Special Issue Use of Water Stable Isotopes in Hydrological Process)

Abstract

:
In this study, a revision of the previously published data on hydrogen (2H/1H) and oxygen (18O/16O) stable isotope ratio of precipitation in northern Chile is presented. Using the amount-weighted mean data and the combined standard deviation (related to both the weighted mean calculation and the spectrometric measurement), the equation of the local meteoric line calculated by error-in-variables regression is as follows: Northern Chile EIV-LMWL: δ2H = [(7.93 ± 0.15) δ18O] + [12.3 ± 2.1]. The slope is similar to that obtained by ordinary least square regression or other types of regression methods, whether weighted or not (e.g., reduced major axis or major axis) by the amount of precipitation. However, the error-in-variables regression is more accurate and suitable than ordinary least square regression (and other types of regression models) where statistical assumptions (i.e., no measurement errors in the x-axis) are violated. A generalized interval of δ2H = ±13.1‰ is also proposed to be used with the local meteoric line. This combines the confidence and prediction intervals around the regression line and appears to be a valid tool for distinguishing outliers or water samples with an isotope composition significantly different from local precipitation. The applicative examples for the Pampa del Tamarugal aquifer system, snow samples and the local geothermal waters are discussed.

1. Introduction

To fit a model to data pairs, the linear least squares regression is by far the most widely used modeling method (e.g., [1]). The ordinary least squares regression (OLSR) minimizes the sum of the squared vertical distances between the y data values and the corresponding y values on the fitted line (the predictions). The OLSR design assumes that there is no variation in the independent variable (x) because it is controlled by the researcher. This is also known as ‘Model I’ regression [2]. In contrast, in ‘Model II’ regression, both the response and the explanatory variables of the model are random (i.e., not controlled by the researcher), and there are errors associated with the measurements of both x and y [2]. Among the ‘Model II’ regression, major axis regression (MA) assumes equal variance on both variables and minimizes the orthogonal (perpendicular) distances from the data points to the fitted line. To the best of our knowledge, one of the first suggestions for applying orthogonal regression to the precipitation amount weighted stable isotope ratios of the oxygen (δ18O, the independent variable x) and hydrogen (δ2H, the dependent variable y) of water molecules in order to determine a better global meteoric water line (GMWL) dates back to the 1980s [3]. Ten years later, the International Atomic Energy Agency (IAEA) proposed the use of the reduced major axis (RMA) regression for calculation of the local meteoric water lines (LMWL) [4]. However, in this latter publication, RMA was wrongly defined as an orthogonal regression. In fact, RMA minimizes the sum of the areas (thus using both vertical and horizontal distances of the data points from the resulting line) rather than the least squares sum of the squared vertical distances, as in OLSR [5]. Meanwhile, the local meteoric water lines in the GNIP dataset of IAEA are calculated through precipitation amount-weighted least squares regression (PWLSR) [6,7]. This regression approach is considered to be the most suitable for representing the isotope composition of groundwater because these are recharged by important rainfall events [6,7]. Alternatively, the precipitation weighted RMA (PWRMA) appears to be most suitable for coastal, island, and Mediterranean sites [6]. Another proposed approach for the calculation of LMWL was the generalized least squares regression model (GENLS), an error-in-variables regression (EIV) that considers the combined standard deviation of the δ18O and δ2H values [8]. Within the statistical literature, synonyms of this approach are “Deming” [9] or “error-in-variables” [10] regression, the former often distinguished as simple (when the data errors are constant among all measurements for each of the two variables) and general (different data error at each observation) [11]. It should be noted that in the case where the variance ratio is equal to 1, Deming regression is equivalent to orthogonal regression. From a general statistical and chemometrics point of view, comparisons between the different approaches can be found in the existing literature [12]. In contrast, in the specific case of water isotope geochemistry and the related meteoric water line, a comparison between error-in-variables and other regression methods has never been presented or discussed. Moreover, despite the advantages identified by previous publications of the alternative regression methods [4,6,7,13], most of the studies on the stable isotope ratio of the waters still use OLSR to calculate the LMWLs. While this does not necessarily mean that the use of OLSR is wrong [7], it has been shown, especially in the hydrology study of arid regions, that the differences between the OLSR approach and other regression approaches can be significant, particularly in terms of the slope of the LMWL [7].
In this manuscript, we compare the regression lines obtained from OLSR, RMA, MA, PWLSR, PWRMA, and PWMA with that from the error-in-variables (EIV) approach, which is very similar to the method followed for GENLS [8]. Here, the existing data on δ18O and δ2H values of precipitation in northern Chile were chosen. This brief study does not attempt to be mathematically or statistically exhaustive and the reader is advised to refer to the referenced literature for more details. In fact, the main purpose here is to provide the results of alternative regression methods applied to stable isotope ratios of precipitation and to provide guidance and advice when the obtained local meteoric water lines are employed for the interpretation of the isotope data of groundwater.

2. Materials and Methods

Northern Chile extends southwards from approximately 18° S to 30° S and includes the two macro-regions of Norte Grande (Regions XV, I and II; Figure 1B) and part of Norte Chico (Regions III and IV; Figure 1B) [14,15,16].
This presents a narrow strip of land (~300 km) between the Andean ridge and the Pacific coast [17]. The latitudinal range of this northern part of Chile aligns with a major extension of the arid desert climate (BW climate zone [20,21]), followed by a tundra climate (ET climate zone [20,21]), and a semi-arid climate in a minor extension [14,21,22]. The Atacama Desert, the driest place on earth with mean precipitation as low as 10 mm per decade, and the higher elevations in the Andes are striking examples of the BW and the ET climates, respectively. The Amazonian influence in the South American monsoon system, with concentrated summer rainfall or dry winters, affects the climate in this area from the north boundary to Nevado Ojos del Salado [22], a stratovolcano with an altitude of 6893 m, the highest point in Chile and the second highest outside of Asia.
The previous studies on isotope composition of meteoric water revealed a general increase in δ18O on the Andean plateau (the “Altiplano”, Figure 1C) from north to south, concomitant with an increase in aridity and decrease in convective moistening (amount effect; [23]). Stable isotope and seasonal precipitation patterns suggest an eastern provenance of the vast majority of moisture that falls as precipitation across the Andean Plateau and Western Cordillera (Figure 1C), with Pacific-derived moisture contributing a minor amount at low elevations near the coast ([23,24]). However, over most regions, the δ18O signal of precipitation is influenced by a combination of factors ([24]). The δ18O -depleted values observed in the high altitude area (i.e., the Andean plateau) were related to processes that affect the air masses that (i) originated over the Atlantic Ocean, (ii) cross the Amazon Basin (continental effect), (iii) ascended the Andes (altitude effect) and (iv) precipitated (convective effect) in the Andean plateau ([25]). In particular, over the eastern Andes, precipitation at low elevations has δ18O from −2 to −8‰, but δ18O becomes more depleted toward the west as vapor is lifted across the Eastern Cordillera of the Andes ([26,27,28,29]). The dominant Altiplano summer rain has δ18O values from −8‰ to −15‰, with the expectation that precipitation in the driest places (e.g., Atamaca desert) should be moderately to strongly depleted at all elevations up to approximately −18‰ ([28,30]).
During the last forty years, several meteoric water lines calculated through OLSR regression were proposed to describe the best fit for the oxygen and hydrogen stable isotope ratios of the precipitation in northern Chile (Table 1, [31,32,33,34,35,36,37,38,39]).
However, not all the studies show the oxygen and hydrogen stable isotope ratios of the water molecule along with precipitation amount [31,32,35,40,41,42]. “La Serena” is the northern Chile station (Region IV, Figure 1B) monitored by the Global Network of Isotopes in Precipitation (GNIP dataset) from 1988 until 2015 [37].
For this study, a dataset of 32 stations has been constructed through collecting previously published data, considering only the isotope values of precipitation with the measured amount and recalculating the amount-weighted mean in order to have only one δ18O and δ2H data pair per station (Table 2; Supplementary File S1).
Assuming x i to be the measurement of the stable isotope ratios expressed as δ18O and δ2H in ‰ versus SMOW (standard mean ocean water)—determined on a specific precipitation amount p i (mm) collected over a specific time period—the standard deviation related to the amount weighted mean x ¯ w [4] of the meteoric stations with more than one accumulation periods of precipitation ( N a . p . > 1) was calculated as follows [44]:
σ w = i = 1 N a . p . p i ( x i - x ¯ w ) 2 ( N a . p . - 1 ) i = 1 N a . p . p i N a . p .
It should be noted that the accumulation periods, which coincides with the total number of the isotope data pairs and precipitation amounts available per each station, are not uniformly distributed through the year and not equally long for most of the stations considered (Supplementary File S1).
In the existing literature, there are several definitions of ‘combined standard deviation’ or ‘combined uncertainty.’ Here, we deal with two different ‘uncertainty’ sources [8,45], that is, a standard deviation related to the amount-weighted mean calculation σ w , as described above, and a standard deviation related to the mass spectrometric measurements σ ms , which differs depending on the publication (Supplementary File S1). One of the easiest ways to calculate the combined standard deviation σ csd involves foreseeing the combination of σ w and σ ms following the law of propagation of uncertainty, that is, the square root of the sum-of-the-squares [45,46,47]:
σ csd = ( σ w ) 2 + ( σ ms ) 2
The combined standard deviations, associated with the δ18O and δ2H amount-weighted mean data in 30 of the total 32 stations, were used to calculate the local meteoric water line through EIV regression. In two stations, Quisquiro [32] and La Serena [37], the calculation of the standard deviation of the weighted mean σ w was performed differently for two reasons. In the case of Quisquiro, the complete dataset is unavailable; therefore, the standard deviation was calculated among the amount-weighted data corresponding to the extreme seasons (winter and summer, [32]). Meanwhile, in the case of La Serena, the σ w was calculated among the years in which more than 70% of precipitation was analyzed for a given isotope composition [4]. Following this, the σ csd was calculated through using the σ ms declared in the analytical section of the publication dataset [32] or elsewhere [25]. With this, the σ csd on δ18O and δ2H is within the mean of the other stations.
The calculation of the EIV regression through different codes employ similar methods based on York and Williamson’s works [48,49,50,51]. Generally speaking, the process of searching for best straight-line parameters requires an iterative procedure in which each new slope estimate is used to modify a weight term, which combines the reciprocal variance of the x and y variables, in a series of convergent fitting cycles [11,52,53,54,55]. The results of this method are compared with those obtained by OLSR, RMA, MA, PWLSR, PWRMA, and PWMA [6]. In accordance with Crawford et al. [6], the Student’s t-test was used to check the difference between the two regression slopes. Furthermore, the obtained EIV meteoric water line was also compared with previously published isotope data on precipitation without amount data [28,30,36,56,57,58,59,60], snow [43,58,59,61,62,63,64,65,66,67,68], groundwater [69], and geothermal water [70] samples collected up to this point in northern Chile.

3. Results

The parameters of the northern Chile meteoric water lines obtained through the different software and regression methods are listed in Table 3.
The values of slope and intercept through EIV appear to be slightly higher than those obtained by other methods that do not take into account the combined standard deviation of the oxygen and hydrogen measurements. However, there is no difference with other regression methods, aside from the results from Cantrell’s bivariate worksheet (EIV-d, Table 3b), which was tested on atmospheric chemical data [54]. The regression results of EIV-d method include: (i) a lower standard error in comparison with those from other codes used for EIV calculation; and (ii) a higher slope when compared to PWLSR, but not significantly different (p > 0.05). The difference between the slopes are further smoothed when the sample #9 is removed (Na.p. = 1; Supplementary File S2). It should be noted that lower standard errors on slope and intercept similar to those of Cantrell’s bivariate worksheet are obtained in the other EIV models introducing a scale factor in the maximum likelihood method, which is estimated by the square root of the reduced chi-square statistic that results from the fit analysis [11,53]. Lower than half standard error are also obtained by “Deming regression” with a constant variance ratio of λ = δ18Ovar2Hvar = 0.019 (EIV-e, Table 3b), that is the mean ratio of the squared combined sigma (Table 2). Therefore, we consider as most representative of the samples the regression results with the standard errors on slope and intercept calculated without scale factor and with combined sigma for each station (EIV results from “a” to “c” in Table 3b):
Northern Chile EIV-LMWL: δ2H = [(7.93 ± 0.15) δ18O] + [12.3 ± 2.1]
The upper and lower confidence intervals at 95% probability of the slope and intercept in the EIV regressions (results from “a” to “c” in Table 3b) include all the previous published values obtained by OLSR (Table 1) and all the other slopes and intercepts of the meteoric water lines calculated by different regression methods (Table 3a). This is clear and straightforward, in particular looking to the standard errors of slope and intercept, which are higher in EIV regression results.

4. Discussion

In isotope hydrology applications, it is often desirable to know whether a water sample from the surface (river, lake, or lagoon) or from an aquifer is shifted in a statistically significant way in relation to the local meteoric water line. This is particularly true in the following: (i) arid zones, where fractionation due to evaporation produces isotopically enriched residual waters; and (ii) geothermal areas, where hot water samples could be isotopically enriched due to water-rock interaction. Northern Chile offers an ideal location for this kind of study because both of these conditions exist here [65,72].
The calculated confidence intervals obtained through EIV regression represent the scatter of the weighted mean data around the fitted line, disregarding the combined standard deviation (Table 3). In OLSR, prediction intervals for y at a specific x value estimate the bands around the fitted line in which the future observation will fall, with a certain probability (usually 95% as confidence interval), given what has already been observed with the sample dataset employed for the regression (t-distribution multiplied by the standard error of the fit). The OLSR prediction interval can be quickly calculated by using commercial or free codes (e.g., [71,73]). However, OLSR does not take into account the measurement error in x when predicting y, either in the combined standard deviations of the data in the fittings or the prediction interval calculations. Moreover, an officially recognized method for the calculation of prediction intervals related to EIV regression still does not exist.
An approximate prediction area around the EIV local meteoric water line could be traced through using an empirical-graphical method, joining the outer limits of the bars representing the combined standard deviations around the mean of the isotope data of precipitation. However, its irregular shape, related to the non-homogeneous values of the combined standard deviations of the regressed stations, necessitates an ad-hoc correction. Specifically speaking, the lower sector of the line (approximately, δ18O < −12‰ and δ2H < −80‰) shows numerous monitored stations with wider standard deviations—especially in terms of hydrogen composition δ2H—in comparison with the higher sector of the line (Figure 2A).
The wider combined standard deviation of the samples that fall within the lower sector is mainly due to the irregular/inhomogeneous frequency of the sampling and to the different accumulation periods of precipitation in the same station (mixing of monthly and seasonal data). In contrast, the stations on the higher sector of the line, which are not depleted by the “altitudinal effect” (lower isotope values at a higher elevation, [74,75]), show a minor combined standard deviation. This is probably due to the fact that only four stations for rainwater are located at an altitude of lower than 2500 m, and monitored by a more regular frequency of sampling (e.g., La Serena station).
An innovative method has recently been proposed for the calculation of a “generalized interval,” which combines the confidence and predictive concepts in EIV [76,77,78]. Using this approach, the obtained interval forcing a mean on 12 values per station (qx = qy = 12; i.e., as supposing a sampling of precipitation extended to 12 months) is 8.9 ± 4.2‰ up and down the fitted line, which is fairly similar to the mean combined standard deviation of δ2H (9.1 ± 4.8‰, excluding the station #9 with Na.p. = 1, Table 2; Supplementary File S3). In contrast, the prediction intervals calculated through the jackknife method [71]—a statistical approach which simulate the resampling of the collected data taking into account their variance—and through Deming regression [71]—which necessitate a constant variance ratio (constant λ [71], calculated on the variances ratio of the 32 samples)—give a prediction interval of ±1.53‰, which is narrower than the combined standard deviation (Supplementary File S3).
Considering that the OLSR and EIV meteoric water lines of northern Chile are not different in terms of slope, the stable isotope ratio data on 115 precipitation samples available from the existing literature (Supplementary File S4; [28,30,36,56,57,58,59,60,68,79]) could help with checking the effectiveness of prediction and generalized intervals. While all these isotope data do not report the amount of precipitation and are therefore not included in the EIV regression calculation, most of them are clustered between or extended over the enriched and depleted sector of the meteoric water line, thus completing the station’s gap. The upper and lower 95% prediction intervals on the δ2H values related to the OLSR regression of this latter dataset are shifted ±13.4‰ up and down respectively on the regressed line (Supplementary File S4). This latter value is not significantly different, neither from the upper/lower δ2H combined standard deviation values of the weighted-amount dataset (±13.8‰), nor from the upper/lower extremes of the generalized interval (±13.1‰). This is particularly true if we take into account that the unweighted rainfall dataset pertains to single rain events, meaning it is quite normal that their predicted δ2H values are shifted towards the higher value of the combined standard deviation of the weighted rainfall dataset.
To verify the effectiveness of this approach, we compared the EIV regression and the above described generalized interval with the isotope composition of the following samples from northern Chile: (i) the groundwater from Pampa del Tamarugal [69] (Figure 1B); (ii) the snow samples (Supplementary File S5; [43,58,59,61,62,63,64,65,66,67,68]); and (iii) the geothermal waters samples [70].
In terms of the first groundwater sample type (i), in previous studies, the clustering of the δ18O and δ2H values of the groundwater from Pampa del Tamarugal falling on the right side of the meteoric water line has led several authors to consider this aquifer system as isotopically different from local precipitation [35,80]. This has been attributed to the following: (i) a recharge in climatic conditions different from the current ones [35]; and (ii) to the evaporation that affects the precipitation in the unsaturated zone of the recharge area [80]. A recent re-evaluation of the isotope data published up to that point shows that the isotope composition of the groundwater from that area is effectively parallel to the meteoric water line, while it was noted that the previous hypothesis on climate variation must be reconsidered [69]. However, the “uncertainty wings” around the meteoric water line have not been traced in any of the above-described cases.
In Figure 2B it is shown how the recent calculated regression line (in orange) representing the groundwater in Pampa del Tamarugal [69] falls within the generalized intervals of the EIV meteoric water line. This demonstrates that the mean composition of the groundwater from this area is not significantly different from the precipitation. This does not exclude the fact that some groundwater samples could be significantly affected by evaporation during recharge. Indeed, the right side of the generalized interval of the EIV regression is the location of the enriched precipitation samples associated with the evaporation during rainfall. Furthermore, the regression line of the groundwater samples should have a proper prediction interval, which certainly does not cover exactly that of the rainwater on that side. However, it is beyond the scope of this study to show also the prediction interval related to groundwater regression.
In terms of the snow samples, all but two of the 61 snow samples from northern Chile fall within the area between the generalized intervals determined in this study, confirming that this area could be approximated as an uncertainty ribbon around the EIV regression (Figure 2B). It should be noted that this dataset includes both samples from precipitation samplers and the Andean snow cover. The two samples out of this area are one outlier (left side of the meteoric water line, Figure 2B) and one isotopically enriched sample from the so-called “penitents” of the Parinacota volcano ([64], Figure 1B—Region I, Supplementary File S5). The penitentes are an annual phenomenon of snow sublimation related to solar radiation and arid climates; consequently, the isotope composition of these samples could be enriched, in particular, if sampled from surface layers [64,67,81]. The volume weighted mean and the related combined standard deviation calculated on the isotope composition of the snow samples from the Cerro Tapado glacier ([67], Figure 1B—Region IV, Supplementary File S5) fall exactly on the EIV-LMWL’s extension beyond the most negative value and within the generalized interval (Figure 2B). As expected, the penitentes from the same area are enriched in comparison with the mean isotope composition of the local snow due to sublimation (Figure 2B, [67]). Moreover, although their mean values fall close to EIV-LMWL, the combined sigma of the penitentes can be higher than the generalized interval (up to ±16‰; Figure 2B, Supplementary File S5).
Finally, geothermal water samples from a recent study were plotted with the EIV regression line and the generalized interval, showing more clearly the distinction between the meteoric waters and the hydrothermal water samples that can be isotopically enriched as a result of high-temperature isotopic fractionations between water–rock, mixing with andesitic water or water–gas interactions [70,74] (Figure 2B). Specifically speaking, the Torta de Tocorpuri (TT in Figure 2B) and Pampa Lirima (PL in Figure 2B) samples are not different from meteoric waters. According to the Tassi et al. [70], the limited temperature and involvement of the so-called “andesitic water” does not produce a relevant shift from the meteoric composition in these two sample groups (Figure 2B). In particular, the PL sample group lies precisely between Pampa del Tamarugal line and the EIV local meteoric water line (Figure 2B).

5. Conclusions

The regression analysis based on the error-in-variables approach (EIV) produces a meteoric water line with a slope similar to those obtained with ordinary least squares regression or other regression approaches.
This most likely occurs because the variance ratio (i.e., the square of the combined standard deviations) is high (δ2Hvar18Ovar >> 3; [2]).
Further studies that focus on other areas and which use a greater number of amount-weighted data pairs (for example, the GNIP/WMO dataset [37]) could be useful in terms of verifying this.
The present study offers an EIV local meteoric water line that summarizes all the previously published isotope data of precipitation from northern Chile. Moreover, the use of the innovative generalized interval (locally, δ2H = ±13.1‰ up and down the fitted values on the line), which is a good compromise between the confidence and the prediction intervals, appears to be a useful tool for distinguishing significantly different data from precipitation.
The case of the groundwater system of Pampa del Tamarugal, snow samples and the local geothermal waters is emblematic. An efficient application of the generalized interval or prediction band could be extended to waters that have been subjected to evaporative fractionation (surface water or precipitation samples from desert climate areas), while it could be used to distinguish the outliers.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/11/4/791/s1: Supplementary File S1—Table 2 raw data. Dark blue-highlighted bold values are the data shown in Table 2. Light blue-highlighted values are the results of partial calculations or additional parameters. Other colors depict parameters from the previously published tables (e.g., yellow-highlighted annual weighted mean in stations 1–9, 10–11, and 12–24) or warnings (grey-highlighted discharged data in stations 29–31 and 32; orange-highlighted measurement error in all the sheets). Supplementary File S2—Regression results as in Table 3 but excluding the sample #9 (Na.p. = 1) from Table 2 dataset. Supplementary File S3—Intervals of the meteoric water line calculated using Table 2 dataset by BivRegBLS [77] and Real Statistics Using Excel© [71]. In the former, code admits only Na.p. > 1 stations, whereas the latter use constant variance ratio (λ = δ18Ovar2Hvar = 0.019). Supplementary File S4—δ18O and δ2H of single precipitation without amount data [28,30,36,56,57,58,59,60,68,79]; OLSR results and prediction band calculated on the isotope composition of the single precipitation dataset. Supplementary File S5—Sheet 1: δ18O and δ2H of single snow samples (collected from snow cover or by precipitation sampler; [43,58,59,61,62,63,64,65,66,68]). Sheet 2: volume weighted mean and combined standard deviation on snow and penitentes [67].

Author Contributions

Study design, data search, data and interpretation, writing the manuscript—original draft preparation, communicating with the journal, T.B.; data search, components for discussion on results and conclusions, writing—review and editing, J.C.; components for discussion on results and conclusions, writing—review and editing, P.I. and E.S.

Funding

This research received no external funding.

Acknowledgments

We are indebted to Bernard Francq, Université Catholique de Louvain—Belgium, who gave a great deal of assistance in the R code and for a critical revision of a first version of the manuscript. Special thanks to Javier Uribe and Jorge Gironás, Pontificia Universidad Católica de Chile, for the raw data on precipitation collected at Salar del Huasco basin. We extend acknowledgments and thanks to all the authors who published isotope data of precipitation in the studied area. Finally, we would like to thank the Editor, Polona Vreča, and two anonymous reviewers for their helpful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NIST/SEMATECH. NIST/SEMATECH e-Handbook of Statistical Methods. Available online: https://www.itl.nist.gov/div898/handbook/index.htm (accessed on 8 February 2019).
  2. Legendre, P.; Legendre, L. Numerical Ecology, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2012; Volume 24. [Google Scholar]
  3. Yurtsever, Y.; Gat, J.R. Atmospheric waters. In Stable Isotope Hydrology: Deuterium and Oxygen-18 in the Water Cycle - Technical Report Series N°210; Gonfiantini, R., Ed.; International Atomic Energy Agency: Vienna, Austria, 1981; pp. 103–142. [Google Scholar]
  4. IAEA. Statistical Treatment of Data on Environmental Isotopes in Precipitation—Technical Report Series n° 331; International Atomic Energy Agency—IAEA: Vienna, Austria, 1992; p. 781. [Google Scholar]
  5. Harper, W.V. Reduced Major Axis Regression: Teaching Alternatives to Least Squares. In Proceedings of the Ninth International Conference on Teaching Statistics (ICOTS9), Flagstaff, AZ, USA, 13–18 July 2014. [Google Scholar]
  6. Crawford, J.; Hughes, C.E.; Lykoudis, S. Alternative least squares methods for determining the meteoric water line, demonstrated using GNIP data. J. Hydrol. 2014, 519, 2331–2340. [Google Scholar] [CrossRef]
  7. Hughes, C.E.; Crawford, J. A new precipitation weighted method for determining the meteoric water line for hydrological applications demonstrated using Australian and global GNIP data. J. Hydrol. 2012, 464, 344–351. [Google Scholar] [CrossRef]
  8. Argiriou, A.A.; Lykoudis, S. Isotopic composition of precipitation in Greece. J. Hydrol. 2006, 327, 486–495. [Google Scholar] [CrossRef]
  9. Deming, W.E. Statistical Adjustment of Data; John Wiley & Sons, Inc.: New York, NY, USA, 1943. [Google Scholar]
  10. Gillard, J. An overview of linear structural models in errors in variables regression. REVSTAT–Stat. J. 2010, 8, 57–80. [Google Scholar]
  11. Maind. Sigmaplot di Systat Software. Available online: http://www.maind.it/document/SigmaplotStatistica.pdf (accessed on 27 December 2018).
  12. Francq, B.G.; Govaerts, B.B. Measurement methods comparison with errors-in-variables regressions. From horizontal to vertical OLS regression, review and new perspectives. Chemom. Intell. Lab. Syst. 2014, 134, 123–139. [Google Scholar] [CrossRef]
  13. Argiriou, A.A.; Salamalikis, V.; Dotsika, E. A Total Weighted Least Squares Method for the Determination of the Meteoric Water Line of Precipitation for Hydrological Purposes. In Perspectives on Atmospheric Sciences; Karacostas, T.S., Bais, A., Nastos, P.T., Eds.; Springer: Cham, Switzerland, 2017; pp. 233–238. [Google Scholar]
  14. Garreaud, R.D. The climate of northern Chile: Mean state, variability and trends. Rev. Mex. De Astron. Y Astrofísica - Ser. De Conf. 2011, 41, 5–11. [Google Scholar]
  15. Harrison, S.; Glasser, N.F. The Pleistocene glaciations of Chile. In Developments in Quaternary Sciences; Elhers, J., Gibbard, P.L., Hughes, C.E., Eds.; Elsevier: Amsterdam, The Netherlands, 2011; Volume 15, pp. 739–756. [Google Scholar]
  16. Yáñez, E.; Lagos, N.A.; Norambuena, R.; Silva, C.; Letelier, J.; Muck, K.P.; San Martin, G.; Benítez, S.; Broitman, B.R.; Contreras, H.; et al. Impacts of climate change on marine fisheries and aquaculture in Chile. In Climate Change Impacts on Fisheries and Aquaculture; Phillips, B.F., Pérez-Ramírez, M., Eds.; John Wiley & Sons: Chichester, UK, 2018; pp. 239–332. [Google Scholar]
  17. SERNAGEOMIN. Mapa geológico de Chile: Versión digital. Escala 1:1’000’000, 1.0 ed.; Servicio Nacional de Geología y Minería, Subdirección Nacional de Geología: Gobierno de Chile, Santiago, 2003; Publicación Geológica Digital No. 4 (CD-ROM, versión 1, 2003). [Google Scholar]
  18. DeCelles, P.G.; Carrapa, B.; Horton, B.K.; McNabb, J.; Gehrels, G.E.; Boyd, J. The Miocene Arizaro Basin, central Andean hinterland: Response to partial lithosphere removal. In Geodynamics of a Cordilleran Orogenic System: The Central Andes of Argentina and Northern Chile: Geological Society of America; DeCelles, P.G., Ducea, M.N., Carrapa, B., Kapp, P.A., Eds.; The Geological Society of America: Boulder, CO, USA, 2015; Volume 212, pp. 359–386. [Google Scholar]
  19. Schilling, F.R.; Trumbull, R.B.; Brasse, H.; Haberland, C.; Asch, G.; Bruhn, D.; Mai, K.; Haak, V.; Giese, P.; Muñoz, M.; et al. Partial melting in the Central Andean crust: A review of geophysical, petrophysical, and petrologic evidence. In The Andes—Active Subduction Orogeny; Oncken, O., Chong, G., Franz, G., Giese, P., Götze, H.J., Ramos, V.A., Strecker, M.R., Wigger, P., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 459–474. [Google Scholar]
  20. Geiger, R. Klassifikation der klimate nach W. Köppen. Landolt-Börnstein–Zahlenwerte Und Funkt. Aus Phys. Chem. Astron. Geophys. Und Tech. 1954, 3, 603–607. [Google Scholar]
  21. Beck, H.E.; Zimmermann, N.E.; McVicar, T.R.; Vergopolan, N.; Berg, A.; Wood, E.F. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 2018, 5. [Google Scholar] [CrossRef]
  22. Sarricolea, P.; Herrera-Ossandon, M.; Meseguer-Ruiz, Ó. Climatic regionalisation of continental Chile. J. Maps 2017, 13, 66–73. [Google Scholar] [CrossRef]
  23. Bershaw, J.; Saylor, J.E.; Garzione, C.N.; Leier, A.; Sundell, K.E. Stable isotope variations (δ18O and δD) in modern waters across the Andean Plateau. Geochim. Cosmochim. Acta 2016, 194, 310–324. [Google Scholar] [CrossRef]
  24. Vuille, M.; Bradley, R.S.; Werner, M.; Healy, R.; Keimig, F. Modeling δ18O in precipitation over the tropical Americas: 1. Interannual variability and climatic controls. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  25. Sánchez-Murillo, R.; Aguirre-Dueñas, E.; Gallardo-Amestica, M.; Moya-Vega, P.; Birkel, C.; Esquivel-Hernández, G.; Boll, J. Isotopic characterization of waters across Chile. In Andean Hydrology; Rivera, D.A., Godoy-Faundez, A., Lillo-Saavedra, M., Eds.; CRC Press: New York, NY, USA, 2018; pp. 205–230. [Google Scholar]
  26. Fiorella, R.P.; Poulsen, C.J.; Pillco-Zolá, R.S.; Barnes, J.B.; Tabor, C.R.; Ehlers, T.A. Spatiotemporal variability of modern precipitation δ18O in the central Andes and implications for paleoclimate and paleoaltimetry estimates. J. Geophys. Res. Atmos. 2015, 120, 4630–4656. [Google Scholar] [CrossRef]
  27. Gonfiantini, R.; Roche, M.A.; Olivry, J.C.; Fontes, J.C.; Zuppi, G.M. The altitude effect on the isotopic composition of tropical rains. Chem. Geol. 2001, 181, 147–167. [Google Scholar] [CrossRef]
  28. Jordan, T.E.; Herrera, L.C.; Godfrey, L.V.; Colucci, S.J.; Gamboa, P.C.; Urrutia, M.J.; González, L.G.; Jacob, F.P. Isotopic characteristics and paleoclimate implications of the extreme precipitation event of March 2015 in northern Chile. Andean Geol. 2019, 46, 1–31. [Google Scholar] [CrossRef]
  29. Rohrmann, A.; Strecker, M.R.; Bookhagen, B.; Mulch, A.; Sachse, D.; Pingel, H.; Alonso, R.N.; Schildgen, T.F.; Montero, C. Can stable isotopes ride out the storms? The role of convection for water isotopes in models, records, and paleoaltimetry studies in the central Andes. Earth Planet. Sci. Lett. 2014, 407, 187–195. [Google Scholar] [CrossRef]
  30. Burgener, L.; Huntington, K.W.; Hoke, G.D.; Schauer, A.; Ringham, M.C.; Latorre, C.; Díaz, F.P. Variations in soil carbonate formation and seasonal bias over >4 km of relief in the western Andes (30° S) revealed by clumped isotope thermometry. Earth Planet. Sci. Lett. 2016, 441, 188–199. [Google Scholar] [CrossRef]
  31. Aravena, R.; Suzuki, O.; Peña, H.; Pollastri, A.; Fuenzalida, H.; Grilli, A. Isotopic composition and origin of the precipitation in Northern Chile. Appl. Geochem. 1999, 14, 411–422. [Google Scholar] [CrossRef]
  32. Chaffaut, I. Précipitations d’altitude, eaux souterraines et changements climatiques de l’altiplano Nord-Chilien; Universite Paris Sud, U.F.R. Scientifique D’Orsay: Paris, France, 1998. [Google Scholar]
  33. Chaffaut, I.; Coudrain-Ribstein, A.; Michelot, J.L.; Pouyaud, B. Précipitation d’altitude du nord-Chili, origine des sources de vapeur et données isotopiques. Bull. De L’institut Français D’études Andin. 1998, 27, 367–384. [Google Scholar]
  34. DGA. Diagnóstico de disponibilidad hídrica en la cuenca del río Lauca, región de Arica y Parinacota; Ministerio de Obras Públicas, Dirección General de Aguas - Arica y Parinacota, XV Región: Arica, Chile, 2015. [Google Scholar]
  35. Fritz, P.; Suzuki, O.; Silva, C.; Salati, E. Isotope hydrology of groundwaters in the Pampa del Tamarugal, Chile. J. Hydrol. 1981, 53, 161–184. [Google Scholar] [CrossRef]
  36. Herrera, C.; Pueyo, J.J.; Sáez, A.; Valero-Garcés, B.L. Relación de aguas superficiales y subterráneas en el área del lago Chungará y lagunas de Cotacotani, norte de Chile: Un estudio isotópico. Rev. Geológica De Chile 2006, 33, 299–325. [Google Scholar] [CrossRef]
  37. IAEA/WMO. Global Network of Isotopes in Precipitation. The GNIP Database. Available online: https://nucleus.iaea.org/wiser (accessed on 15 January 2019).
  38. Squeo, F.A.; Aravena, R.; Aguirre, E.; Pollastri, A.; Jorquera, C.B.; Ehleringer, J.R. Groundwater dynamics in a coastal aquifer in north-central Chile: Implications for groundwater recharge in an arid ecosystem. J. Arid Environ. 2006, 67, 240–254. [Google Scholar] [CrossRef]
  39. Troncoso, R.; Castro, R.; Lorca, M.E.; Espinoza, C.; Pérez, Y. Análisis Preliminar de la Composición Isotópica Oxígeno 18 – Deuterio de las Aguas de la Cuenca del Río Copiapó, Región de Atacama: Una Contribución al Conocimiento del Sistema Hidrogeológico. In Proceedings of the XIII Congreso Geológico Chileno, Antofagasta, Chile, 5–9 August 2012; pp. 774–776. [Google Scholar]
  40. Uribe, J.; Muñoz, J.F.; Gironás, J.; Oyarzún, R.; Aguirre, E.; Aravena, R. Assessing groundwater recharge in an Andean closed basin using isotopic characterization and a rainfall-runoff model: Salar del Huasco basin, Chile. Hydrogeol. J. 2015, 23, 1535–1551. [Google Scholar] [CrossRef]
  41. Salazar, C.M.; Rojas, L.B.; Pollastri, A. Evaluación de recursos hídricos en el sector de Pica hoya de la Pampa del Tamarugal I region; Ministerio de Obras Públicas, Dirección General de Aguas, CCHEN: Santiago, Chile, 1998; p. 98. [Google Scholar]
  42. Aravena, R.; Peña, H.; Grilli, A.; Suzuki, O.; Mordeckai, M. Evolución isotópica de las lluvias y origen de las masas de aire en el Altiplano chileno. In Isotope Hydrology Investigations in Latin America; International Atomic Energy Agency—IAEA: Vienna, Austria, 1989; pp. 129–142. [Google Scholar]
  43. Herrera, C.; Custodio, E.; Chong, G.; Lambán, L.J.; Riquelme, R.; Wilke, H.; Jódar, J.; Urrutia, J.; Urqueta, H.; Sarmiento, A.; et al. Groundwater flow in a closed basin with a saline shallow lake in a volcanic area: Laguna Tuyajto, northern Chilean Altiplano of the Andes. Sci. Total Environ. 2016, 541, 303–318. [Google Scholar] [CrossRef]
  44. NIST. DATAPLOT Reference Manual - weighted standard deviation. Available online: https://www.itl.nist.gov/div898/software/dataplot/refman2/ch2/weightsd.pdf (accessed on 28 December 2018).
  45. Taylor, B.N.; Kuyatt, C.E. Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results; National Institute of Standards and Technology—NIST: Gaithersburg, MD, USA, 1994. [Google Scholar]
  46. Leito, I.; Jalukse, L.; Helm, I. Estimation of Measurement Uncertainty in Chemical Analysis (Analytical Chemistry) Course; University of Tartu: Tartu, Estonia, 2018. [Google Scholar]
  47. Bell, S. A Beginner’s Guide to Uncertainty of Measurement; National Physical Laboratory: Teddington, UK, 1999. [Google Scholar]
  48. Williamson, J.H. Least-squares fitting of a straight line. Can. J. Phys. 1968, 46, 1845–1847. [Google Scholar] [CrossRef]
  49. York, D. Least-squares fitting of a straight line. Can. J. Phys. 1966, 44, 1079–1086. [Google Scholar] [CrossRef]
  50. York, D. Unified equations for the slope, intercept, and standard error of the best straight line. Am. J. Phys. 2004, 72, 367–375. [Google Scholar] [CrossRef]
  51. York, D. Least squares fitting of a straight line with correlated errors. Earth Planet. Sci. Lett. 1968, 5, 320–324. [Google Scholar] [CrossRef]
  52. Originlab. Algorithms (Fit Linear with X Error). Available online: https://www.originlab.com/doc/Origin-Help/Ref-Linear-XErr#Fit_Parameters (accessed on 27 December 2018).
  53. Systat. Using SigmaStat Statistics in SigmaPlot; Systat Software: San Jose, CA, USA, 2013; p. 470. [Google Scholar]
  54. Cantrell, C.A. Technical Note: Review of methods for linear least-squares fitting of data and application to atmospheric chemistry problems. Atmos. Chem. Phys. 2008, 8, 5477–5487. [Google Scholar] [CrossRef]
  55. Sturm, P. bfsl: Best-Fit Straight Line, 0.1.0; CRAN.R-project.org. 2018. Available online: https://cran.r-project.org/web/packages/bfsl/index.html (accessed on 3 February 2019).
  56. Aravena, R.; Suzuki, O.; Pollastri, A. Coastal fog and its relation to groundwater in the IV region of northern Chile. Chem. Geol. 1989, 79, 83–91. [Google Scholar] [CrossRef]
  57. Cifuentes, J.L.; Cervetto, M.M.; López, L.A.; Fuentes, F.C.; Feuker, P.; Espinoza, C. Análisis preliminar de Isótopos Estables en aguas subterráneas, superficiales y lluvia de la Pampa del Tamarugal. In Proceedings of the XIV Congreso Geológico Chileno, La Serena, Chile, 4–8 October 2015; pp. 277–280. [Google Scholar]
  58. DGA. Levantamiento hidrogeológico para el desarrollo de nuevas fuentes de agua en áreas prioritarias de la zona norte de Chile, Regiones XV, I, II y III, Etapa 2; Ministerio de Obras Públicas, Dirección General de Aguas, Departamento de Estudios y Planificación; Pontificia Universidad Católica de Chile, Departamento de Ingeniería Hidráulica y Ambiental: Santiago, Chile, 2009. [Google Scholar]
  59. López, L.A.; Cifuentes, J.L.; Fuentes, F.C.; Neira, H.A.; Cervetto, M.M.; Troncoso, R.A.; Feuker, P. Hidrogeología de la Cuenca de la Pampa del Tamarugal, Región de Tarapacá; Gobierno de Chile, Servicio Nacional de Geología y Minería: Santiago, Chile, 2017; p. 186. [Google Scholar]
  60. Durán, L.V.L. Hidrogeoquímica de fuentes termales en ambientes salinos relacionados con salares en los Andes del Norte de Chile; Universidad de Chile: Santigo, Chile, 2016. [Google Scholar]
  61. Giggenbach, W. The isotopic composition of waters from the El Tatio geothermal field, Northern Chile. Geochim. Cosmochim. Acta 1978, 42, 979–988. [Google Scholar] [CrossRef]
  62. Godfrey, L.V.; Jordan, T.E.; Lowenstein, T.K.; Alonso, R.L. Stable isotope constraints on the transport of water to the Andes between 22° and 26°S during the last glacial cycle. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2003, 194, 299–317. [Google Scholar] [CrossRef]
  63. Munoz-Saez, C.; Manga, M.; Hurwitz, S. Hydrothermal discharge from the El Tatio basin, Atacama, Chile. J. Volcanol. Geotherm. Res. 2018, 361, 25–35. [Google Scholar] [CrossRef]
  64. Peña, H. Mediciones de 180 y 2H en “Penitentes” de Nieve. In Isotope Hydrology Investigations in Latin America; International Atomic Energy Agency—IAEA: Vienna, Austria, 1988; Volume IAEA-TECDOC-502, pp. 143–154. [Google Scholar]
  65. Cortecci, G.; Boschetti, T.; Mussi, M.; Herrera Lameli, C.; Mucchino, C.; Barbieri, M. New chemical and isotopic data on waters of El Tatio Geothermal Field, Northern Chile. Geochem. J. 2005, 39, 547–571. [Google Scholar] [CrossRef]
  66. Alpers, C.N.; Whittemore, D.O. Hydrogeochemistry and stable isotopes of ground and surface waters from two adjacent closed basins, Atacama Desert, northern Chile. Appl. Geochem. 1990, 5, 719–734. [Google Scholar] [CrossRef]
  67. Sinclair, K.E.; MacDonell, S. Seasonal evolution of penitente glaciochemistry at Tapado Glacier, Northern Chile. Hydrol. Process. 2016, 30, 176–186. [Google Scholar] [CrossRef]
  68. Cervetto, M.M. Caracterización hidrogeológica e hidrogeoquímica de las cuencas: Salar de Aguas calientes 2, Puntas negras, Laguna Tuyajto, Pampa Colorada, Pampa Las Tecas y Salar el Laco, II región de Chile; Universidad de Chile: Santiago, Chile, 2012. [Google Scholar]
  69. Scheihing, K.; Moya, C.; Struck, U.; Lictevout, E.; Tröger, U. Reassessing Hydrological Processes That Control Stable Isotope Tracers in Groundwater of the Atacama Desert (Northern Chile). Hydrology 2018, 5, 3. [Google Scholar] [CrossRef]
  70. Tassi, F.; Aguilera, F.; Darrah, T.; Vaselli, O.; Capaccioni, B.; Poreda, R.J.; Huertas, A.D. Fluid geochemistry of hydrothermal systems in the Arica-Parinacota, Tarapacá and Antofagasta regions (northern Chile). J. Volcanol. Geotherm. Res. 2010, 192, 1–15. [Google Scholar] [CrossRef]
  71. Zaiontz, C. Real Statistics Using Excel. Available online: www.real-statistics.com (accessed on 3 February 2019).
  72. Boschetti, T.; Cortecci, C.; Barbieri, M.; Mussi, M. New and past geochemical data on fresh to brine waters of the Salar de Atacama and Andean Altiplano, northern Chile. Geofluids 2007, 7, 35–50. [Google Scholar] [CrossRef]
  73. Verschuuren, G. Excel 2013 for Scientists; Holy Macro! Books: Chicago, IL, USA, 2014; p. 250. [Google Scholar]
  74. Clark, I. Groundwater Geochemistry and Isotopes; CRC Press, Taylor & Francis Group: Boca Raton, FL, USA, 2015. [Google Scholar]
  75. Gat, J.R. Isotope Hydrology: A Study of the Water Cycle; Imperial College Press: London, UK; Singapore, 2010. [Google Scholar]
  76. Francq, B.; Govaerts, B. How to regress and predict in a Bland–Altman plot? Review and contribution based on tolerance intervals and correlated-errors-in-variables models. Stat. Med. 2016, 35, 2328–2358. [Google Scholar] [CrossRef] [PubMed]
  77. Francq, B.G.; Berger, M. BivRegBLS: Tolerance Intervals and Errors-in-Variables Regressions in Method Comparison Studies, 1.0.0; CRAN.R-project.org. 2017. Available online: https://CRAN.R-project.org/package=BivRegBLS (accessed on 3 February 2019).
  78. Berger, M.; Francq, B. BivRegBLS: A new R package in method comparison studies with tolerance intervals and (correlated)-errors-in-variables regressions. In Proceedings of the Chimiométrie XVIII, Paris, France, 30 January–1 February 2017. [Google Scholar]
  79. Lorca, M.E. Hidrogeología e hidrogeoquímica de la cuenca de la Quebrada Paipote, Región de Atacama; Universidad de Chile: Santiago, Chile, 2011. [Google Scholar]
  80. Aravena, R. Isotope hydrology and geochemistry of northern Chile groundwaters. Bull. De L’institut Français D’études Andin. 1995, 24, 495–503. [Google Scholar]
  81. Stichler, W.; Schotterer, U.; Fröhlich, K.; Ginot, P.; Kull, C.; Gäggeler, H.; Pouyaud, B. Influence of sublimation on stable isotope records recovered from high-altitude glaciers in the tropical Andes. J. Geophys. Res. Atmos. 2001, 106, 22613–22620. [Google Scholar] [CrossRef]
Figure 1. (A) Localization of the area of this study. (B) Morphostructural domains (shades of brown) and administrative regions (modified from [17]). In blue: Pampa de Tamarugal aquifer. (C) Morpho-tectonic units of the central Andes (modified from [18]). SBS: Santa Barbara System; SP: Sierras Pampeanas. The “Modern Forearc” coincides with the Precordillera and the Coastal Cordillera [19].
Figure 1. (A) Localization of the area of this study. (B) Morphostructural domains (shades of brown) and administrative regions (modified from [17]). In blue: Pampa de Tamarugal aquifer. (C) Morpho-tectonic units of the central Andes (modified from [18]). SBS: Santa Barbara System; SP: Sierras Pampeanas. The “Modern Forearc” coincides with the Precordillera and the Coastal Cordillera [19].
Water 11 00791 g001
Figure 2. (A) δ2H and δ18O amount-weighted mean of the samples (open circles; Table 2) along with the combined standard deviation on hydrogen and oxygen isotope composition (y and x bars, respectively); error-in-variables local meteoric water line (EIV-LMWL solid line) and the generalized interval (dashed lines; δ2H = ±13.1‰ up and down the fitted values on the line). (B) Pampa de Tamarugal aquifer (thick-orange line, [69]), snow and penitentes (dark and light blue diamonds, respectively [41,43,59,61,63,64,65,66,67,68]) and hydrothermal groundwater (red triangles, [70]; PL: Pampa Lirima, TT: Torta de Tocorpuri) samples compared with the lines described in ‘A’. Arrows explain the isotope effects in hydrothermal groundwater [70,74]. y and x bars in snow and penitentes depict the combined sigma related to the volume weighted mean (samples from Cerro Tapado [67], Supplementary File S5). Light-blue diamonds without error bars are single samples of the penitentes from Parinacota volcano ([64], Supplementary File S5).
Figure 2. (A) δ2H and δ18O amount-weighted mean of the samples (open circles; Table 2) along with the combined standard deviation on hydrogen and oxygen isotope composition (y and x bars, respectively); error-in-variables local meteoric water line (EIV-LMWL solid line) and the generalized interval (dashed lines; δ2H = ±13.1‰ up and down the fitted values on the line). (B) Pampa de Tamarugal aquifer (thick-orange line, [69]), snow and penitentes (dark and light blue diamonds, respectively [41,43,59,61,63,64,65,66,67,68]) and hydrothermal groundwater (red triangles, [70]; PL: Pampa Lirima, TT: Torta de Tocorpuri) samples compared with the lines described in ‘A’. Arrows explain the isotope effects in hydrothermal groundwater [70,74]. y and x bars in snow and penitentes depict the combined sigma related to the volume weighted mean (samples from Cerro Tapado [67], Supplementary File S5). Light-blue diamonds without error bars are single samples of the penitentes from Parinacota volcano ([64], Supplementary File S5).
Water 11 00791 g002
Table 1. Slope and intercept values of the previously published northern Chile meteoric water lines [31,32,33,34,35,36,37,38,39]. All the shown data were calculated by OLSR. The IAEA/GNIP data pairs are referred to as the “La Serena” station (Figure 1B). For this station, other two equations are also proposed on the same database, δ2H = [(7.64 ± 0.39) δ18O] + [9.44 ± 2.32] and δ2H = [(8.04 ± 0.45) δ18O] + [9.98 ± 2.59] by PWLSR and RMA methods [6,7], respectively, both applied to all the available data (N = 40).
Table 1. Slope and intercept values of the previously published northern Chile meteoric water lines [31,32,33,34,35,36,37,38,39]. All the shown data were calculated by OLSR. The IAEA/GNIP data pairs are referred to as the “La Serena” station (Figure 1B). For this station, other two equations are also proposed on the same database, δ2H = [(7.64 ± 0.39) δ18O] + [9.44 ± 2.32] and δ2H = [(8.04 ± 0.45) δ18O] + [9.98 ± 2.59] by PWLSR and RMA methods [6,7], respectively, both applied to all the available data (N = 40).
Slope±s.e.Intercept±s.e.NR2ReferenceYearRef. [#]
7.8 -10.3 -29-Fritz et al.1981[35]
8.2±0.1716.6±2.2390.98Chaffaut et al.1998[32,33]
7.8 -9.7 -129-Aravena et al.1999[31]
7.9 -14 ---Herrera et al.2006[36]
7.7 -9.6 --0.91Squeo et al.2006[38]
7.95 -14.9 --0.99DGA2015[34]
8.07 -13.5 -230.98Troncoso et al.2012[39]
7.52±0.467.18±2.64400.87IAEA/WHO2015[37]
mean7.87 11.97
std.de.0.21 3.22
-: data not calculated or not declared; s.e.: standard error.
Table 2. Amount-weighted oxygen and hydrogen stable isotope ratios (δ18O and δ2H in permil versus the standard mean ocean water SMOW, respectively) of the previously published precipitation samples collected in northern Chile [31,32,35,37,40,41,42,43]. Calculations are described in the Supplementary File S1. Samples were mainly collected in Region I (from #1 to #27), Region II (from #28 to #31) and Region IV (#32) (Figure 1B).
Table 2. Amount-weighted oxygen and hydrogen stable isotope ratios (δ18O and δ2H in permil versus the standard mean ocean water SMOW, respectively) of the previously published precipitation samples collected in northern Chile [31,32,35,37,40,41,42,43]. Calculations are described in the Supplementary File S1. Samples were mainly collected in Region I (from #1 to #27), Region II (from #28 to #31) and Region IV (#32) (Figure 1B).
ID #LocationCoordinatesElevation Metersδ18O (a.w.m. ± c.s.)δ2H (a.w.m. ± c.s.)Na.p.Precipitation Amount (Mean Value) MillimetersReferencesRef. [#]
LatitudeLongitude‰ vs. SMOW‰ vs. SMOW
1Apacheta Tapa−19.5891−68.99404350−17.48 ± 1.14−122.3 ± 8.9473.73Fritz et al. 1981[35]
2Chusmiza−19.7060−69.19063360−7.51 ± 0.22−48.9 ± 1.03103.93Fritz et al. 1981[35]
3Alto Mocha−19.8181−69.29872590−7.13 ± 0.31−43.7 ± 1.5233.90Fritz et al. 1981[35]
4Mocha−19.8436−69.28402200−4.91 ± 0.90−29.2 ± 6.7239.65Fritz et al. 1981[35]
5Collacagua−20.0703−68.82593915−15.00 ± 1.22−108.1 ± 8.6782.23Fritz et al. 1981[35]
6Huasco−20.3144−68.80653800−17.13 ± 0.65−127.7 ± 9.3334.63Fritz et al. 1981[35]
7Indio Muerto−20.3610−68.96154135−18.27 ± 0.35−129.9 ± 3.4453.05Fritz et al. 1981[35]
8Tambillos−20.4672−69.14153300−8.02 ± 0.27−50.1 ± 2.1377.13Fritz et al. 1981[35]
9Apacheta Mama *−20.4734−69.15804115−19.90 ± 0.15−146.0 ± 1.0177.00Fritz et al. 1981[35]
10Pumire−19.0955−69.11084200−2.25 ± 0.78−11.5 ± 3.1310.00Salazar et al. 1998; Aravena et al. 1989[41,42]
11Coposa−20.7089−68.69423460−16.66 ± 2.01−119.8 ± 16.299.17Salazar et al. 1998; Aravena et al. 1989[41,42]
12Collaguasi−20.9833−68.70004250−16.96 ± 1.99−122.3 ± 16.0513.30Aravena et al. 1999[31]
13Ujina−20.9833−68.60004200−11.98 ± 1.57−83.3 ± 12.0818.79Aravena et al. 1999[31]
14Puchuldiza−19.4000−68.95004150−14.55 ± 0.93−98.4 ± 11.01511.88Aravena et al. 1999[31]
15Pampa Lirima−19.8200−68.90004100−16.97 ± 1.23−121.6 ± 9.7279.86Aravena et al. 1999[31]
16Colchane−19.7003−68.88333965−14.95 ± 2.28−106.3 ± 18.2451.90Aravena et al. 1999[31]
17Collacagua−20.0333−68.85003990−12.68 ± 1.34−88.3 ± 10.51415.04Aravena et al. 1999[31]
18Cancosa−19.8500−68.60003800−14.00 ± 2.45−96.4 ± 20.7528.50Aravena et al. 1999[31]
19Huaytane−19.5433−68.60423720−14.10 ± 1.48−100.9 ± 12.1532.20Aravena et al. 1999[31]
20Copaquire−20.9305−68.89223490−12.47 ± 0.98−84.6 ± 7.8165.75Aravena et al. 1999[31]
21Poroma−19.8667−69.18332880−5.00 ± 0.66−28.1 ± 6.098.01Aravena et al. 1999[31]
22Parca−20.0167−68.85002570−6.64 ± 0.90−40.9 ± 7.844.88Aravena et al. 1999[31]
23Huatacondo−20.9167−69.05002460−8.26 ± 1.26−53.8 ± 9.8126.21Aravena et al. 1999[31]
24Camina−19.3116−69.42992380−7.02 ± 1.00−42.8 ± 8.1411.50Aravena et al. 1999[31]
25Sillillica −20.1738−68.73774270−15.96 ± 1.07−110.9 ± 7.9395.30Uribe et al. 2015[40]
26Altos del Huasco −20.3221−68.90223784−14.82 ± 1.57−104.8 ± 12.3351.80Uribe et al. 2015[40]
27Diablo Marca −20.0505−68.99624603−18.43 ± 1.59−128.8 ± 11.8378.60Uribe et al. 2015[40]
28Quisquiro−23.2100−67.25004260−13.40 ± 2.12−91.1 ± 13.33614.30Chaffaut 1998[32]
29Tuyajto−23.9435−67.59214040−6.23 ± 0.50−32.4 ± 3.1213.30Herrera et al. 2016[43]
30Pampa Colorada−23.8588−67.47744426−8.80 ± 0.86−51.3 ± 7.6328.50Herrera et al. 2016[43]
31Aguas Calientes-3−23.9171−67.69713900−6.00 ± 0.51−31.9 ± 4.5321.67Herrera et al. 2016[43]
32La Serena−29.8981−71.2425142−5.78 ± 1.32−35.9 ± 9.93443.10IAEA/WHO 2015[37]
a.w.m. ± c.s.: amount-weighted mean ± combined sigma (see methods paragraph for details); Na.p.: Number of accumulation periods of precipitation. This value coincides with the total number of the available isotope data pairs per station (Supplementary File S1); * sample with only one accumulation period (Na.p. = 1). In this case, the combined sigma (c.s.) contain only the standard deviation related to spectrometric measurement (Supplementary File S1). The regression has been calculated both including (Table 3) and excluding (Supplementary File S2) this sample.
Table 3. Regression results obtained by different methods applied to the data of Table 2. (a) Crawford et al.’s non-weighted (OLSR, MA, RMA) and amount weighted (PWLSR, PWRMA, PWMA) regressions [6,7]. (b) Ordinary Least Squares Regression (OLSR) and error-in-variables regression (EIV) results by different codes [52,53,55,71]. (c) Comparison between regression’s slopes obtained by EIV and PWLR, PWRMA, and PWRMA.
(a) regression results according to Crawford et al. 2014
(a) regression results according to Crawford et al. 2014
ModelsNSlopeInterceptrmSSEavt-Valuep
Values.e.Values.e.
OLSR327.780.1011.31.21.0009--
RMA327.800.0911.51.21.00040.13180.896
MA327.820.1011.71.21.00090.25930.797
PWLSR327.740.0910.71.31.00540.30470.763
PWRMA327.760.0910.91.31.00210.18100.858
PWMA327.780.0911.11.31.00000.06080.952
Local Meteoric Water Line Freeware©—Australian Nuclear Science and Technology Organisation. OLSR: Ordinary Least Square Regression; RMA: Reduced Major Axis; MA: Major Axis; PWLSR: Precipitation Weighted Least Square Regression; PWRMA: Precipitation Weighted Reduced Major Axis. PWMA: Precipitation Weighted Major Axis. s.e.: standard error; rmSSEav: root mean sum of squared error average; Comparison between OLSR and the other regression models listed in (a): t-value: Student’s t-test value; p: probability distribution by Microsoft®Office Excel (d.f. = N − 2).
(b) EIV regression results compared with OLSR
(b) EIV regression results compared with OLSR
ModelsNSlopeIntercept95% C.I. Slope95% C.I. InterceptrGOFrmSSEt-Valuep
Values.e.Values.e.UpperLowerUpperLower
OLSR327.790.1011.31.27.987.5913.88.70.99777.17402.6784--
EIV-a327.930.1512.32.18.247.6316.78.00.99770.14570.38170.82920.414
EIV-b327.930.1512.32.18.247.6216.67.90.99770.14570.38170.81190.423
EIV-c327.930.1512.32.1-----0.14570.38170.82890.414
EIV-d327.930.0612.30.8-----0.14570.38171.31850.197
EIV-e327.810.1211.51.68.057.5614.88.2- -0.12560.901
OLSR: here calculated by Microsoft®Office Excel—2013 (GOF parameter from Origin®Pro). EIV-a: Origin®Pro; EIV-b: SigmaPlot©; EIV-c: BFSL; EIV-d: Cantrell 2008; EIV-e: Real Statistics Using Excel© with a mean variances ratio of λ = Xvar/Yvar = 0.01913. C.I.: Confidence Intervals; r: Pearson’s correlation coefficient; GOF: Goodness of Fit; rmSSE: root mean sum of squared error (this parameter is called “standard error” of regression in Microsoft®Excel’s OLSR; in other regression codes, it could be easily calculated by the square root of the GOF parameter).
(c) EIV regression results compared with PWLSR, PWRMA, and PWMA
(c) EIV regression results compared with PWLSR, PWRMA, and PWMA
Modelst-ValuepModelst-ValuepModelst-Valuep
PWLSR--PWRMA--PWMA--
EIV-a1.08430.287EIV-a0.98830.331EIV-a0.89460.378
EIV-b1.06620.295EIV-b0.97050.340EIV-b0.87700.387
EIV-c1.08400.287EIV-c0.98810.331EIV-c0.89440.378
EIV-d1.74810.091EIV-d1.59340.122EIV-d1.44100.160
EIV-e0.41410.682EIV-e0.30290.764EIV-e0.19490.847

Share and Cite

MDPI and ACS Style

Boschetti, T.; Cifuentes, J.; Iacumin, P.; Selmo, E. Local Meteoric Water Line of Northern Chile (18° S–30° S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation. Water 2019, 11, 791. https://doi.org/10.3390/w11040791

AMA Style

Boschetti T, Cifuentes J, Iacumin P, Selmo E. Local Meteoric Water Line of Northern Chile (18° S–30° S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation. Water. 2019; 11(4):791. https://doi.org/10.3390/w11040791

Chicago/Turabian Style

Boschetti, Tiziano, José Cifuentes, Paola Iacumin, and Enricomaria Selmo. 2019. "Local Meteoric Water Line of Northern Chile (18° S–30° S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation" Water 11, no. 4: 791. https://doi.org/10.3390/w11040791

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