Open Journal of Modern Hydrology
Vol.3 No.4(2013), Article ID:38886,15 pages DOI:10.4236/ojmh.2013.34027

Spatial and Temporal Variation of Stable Isotopes in Precipitation across Costa Rica: An Analysis of Historic GNIP Records

Ricardo Sánchez-Murillo1,2*, Germain Esquivel-Hernández1,3, Kristen Welsh2,4, Erin S. Brooks5, Jan Boll2,5, Rosa Alfaro-Solís1,3, Juan Valdés-González1,3

1Escuela de Química, Universidad Nacional, Heredia, Costa Rica; 2Waters of the West-Water Resources Program, University of Idaho, Moscow, USA; 3Laboratorio de Química de la Atmósfera, Escuela de Química, Universidad Nacional, Heredia, Costa Rica; 4Division of Research and Development, Centro Agronómico Tropical de Investigación y Enseñanza (CATIE), Turrialba, Costa Rica; 5Department of Biological and Agricultural Engineering, University of Idaho, Moscow, USA.

Email: *

Copyright © 2013 Ricardo Sánchez-Murillo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received September 1st, 2013; revised October 1st, 2013; accepted October 8th, 2013

Keywords: Stable Isotopes; Costa Rica; Precipitation; GNIP; HYSPLIT Model


The location of Costa Rica on the Central American Isthmus creates unique microclimate systems that receive moisture inputs directly from the Caribbean Sea and the Pacific Ocean. In Costa Rica, stable isotope monitoring was conducted by the International Atomic Energy Agency and the World Meteorological Association as part of the worldwide effort entitled Global Network of Isotopes in Precipitation. Sampling campaigns were mainly comprised of monthly-integrated samples during intermittent years from 1990 to 2005. The main goal of this study was to determine spatial and temporal isotopic variations of meteoric waters in Costa Rica using historic records. Samples were grouped in four main regions: Nicoya Peninsula (d2H = 6.65d18O − 0.13; r2 = 0.86); Pacific Coast (d2H = 7.60d18O + 7.95; r2 = 0.99); Caribbean Slope (d2H = 6.97d18O + 4.97; r2 = 0.97); and Central Valley (d2H = 7.94d18O + 10.38; r2 = 0.98). The water meteoric line for Costa Rica can be defined as d2H = 7.61d18O + 7.40 (r2 = 0.98). The regression of precipitation amount and annual arithmetic means yields a slope of −1.6‰ d18O per 100 mm of rain (r2 = 0.57) which corresponds with a temperature effect of −0.37‰ d18O/˚C. A strong correlation (r2 = 0.77) of −2.0‰ d18O per km of elevation was found. Samples within the Nicoya Peninsula and Caribbean lowlands appear to be dominated by evaporation enrichment as shown in d-excess interpolation, especially during the dry months, likely resulting from small precipitation amounts. In the inter-mountainous region of the Central Valley and Pacific slope, complex moisture recycling processes may dominate isotopic variations. Generally, isotopic values tend to be more depleted as the rainy season progresses over the year. Air parcel back trajectories indicate that enriched isotopic compositions both in Turrialba and Monteverde are related to central Caribbean parental moisture and low rainfall intensities. Depleted events appear to be related to high rainfall amounts despite the parental origin of the moisture.

1. Introduction

Tropical regions cover approximately 50% of Earth’s landmass and are home to three-quarters of the human population. During the last decade, the scientific community, environmental institutions, governments, and communities have increased their awareness of tropical climate variability based on the premise that changes in regional and global circulation trends may lead to intensification (i.e., greater rainfall intensity) or weakening (i.e., prolonged drought) of the hydrological cycle in particular regions [1]. For example, several studies have shown the potential high vulnerability of tropical ecosystems to increasing mean surface temperatures [2-4]. In addition, extreme climate events could be even more severe than those experienced to date [5]. Climate anomalies will increase the vulnerability of certain regions and communities to changes in magnitude, timing, and duration of hydrological responses.

Understanding precipitation patterns, particularly how they differ spatially and temporally, will provide further information on regional and global hydrologic processes and enable better planning and preparation for future potential changes in climate. A better understanding of the factors that control precipitation patterns will inform stakeholders and environmental agencies so that they can prioritize efforts and resources in ungauged basins where potential future droughts or floods may drastically affect ecological assemblages and socio-economic activities. In Central America, many socio-economic activities, especially agriculture, tourism, and hydropower generation, are dependent on seasonal cycles of precipitation. Therefore, in order to determine the interannual variability in precipitation, it is necessary to describe the regional annual cycle precisely and examine the mechanisms that control it [6].

Across the globe, spatial variation of stable isotopes in the water cycle occurs due to the movement of air masses and fractionation during phase changes, such as condensation and evaporation, creating a directional spatial pattern across regions and continents [7]. Light stable isotope compositions of tropical meteoric waters have proven to be an important indicator of current climate variability [8-16]. In particular, d18O values have provided novel insights on El Niño/Southern Oscillation dynamics [12,17-19]. In temperate regions, isotopic variations in precipitation frequently have been correlated with mean surface air temperatures [20-25], whereas, in the tropics, several authors have reported “the amount effect” as the main controlling factor [16,24,26,27]. These consistent processes create relatively strong spatial patterns of isotopic compositions across landscapes [28]. Large-scale regional and global patterns can be depicted through isoscapes that illustrate the spatial and temporal variation of isotopes over a geographic region.

Despite novel advances to understand stable isotope precipitation dynamics in the tropics, there is consensus among scientists in the published literature [13,16,27] on the urgent need for long-term sampling networks that may lead to a better spatial and temporal resolution, particularly in mountainous regions such as the Central American Continental Divide (CACD). In the CACD orographic effects, moisture recycling and canopy interception, intense evapotranspiration, and microclimates are relevant factors in isotopic variation. Furthermore, recent studies have demonstrated that variations in the stable isotope composition of precipitation occur even between storm events [29-32], emphasizing the importance of event-based sampling.

Stable isotope research in Costa Rica has been limited to date with no studies conducted on a regional scale examining trends in precipitation. Lachniet and Patterson (2002) [27] provided a preliminary evaluation of d18O and d2H in Costa Rican surface waters (n = 63 from rivers and lakes). Their study found regional trends in the isotopic composition of surface waters inversely correlated to altitude ranges (−1.4 d18O ‰/km) and evaporation enrichment, particularly in the northwestern region. Rhodes et al. (2006) studied the seasonal variation of isotopic composition in tropical montane forests of Monteverde’s National Park located on the continental divide. Based on seasonal d-excess variations, the authors concluded that water evaporated from land is an important component of the water budget to the region during the transitional and dry seasons. Reynolds and Fraile (2009) [33] conducted a groundwater recharge study on one of the main aquifers of the Central Valley of Costa Rica that provides water to approximately 30% of the population. Based upon precipitation, well, and spring records they were able to identify local recharge areas mainly dominated by Pacific-originated storms. Greater understanding of spatial and temporal isotopic variations in tropical regions with similar topographic, vegetation, and climatic attributes can enhance hydrological studies of this nature.

Costa Rica’s climate is mainly controlled by four phenomena: northeast trades-winds (i.e., alisios), the shifts of the Intertropical Convergence Zone (ITCZ), cold continental outbreaks, and indirect influence of Caribbean cyclones [34]. These circulation processes produce three distinct seasons. The wet season (May-October) corresponds to the time when the ITCZ travels over Costa Rica, and precipitation is characterized by heavy convective rain storms. The transitional (NovemberJanuary) and dry (February-April) seasons comprise the months when the ITCZ is located to the south of Costa Rica [35]. During the wet season, the air masses arriving in Costa Rica can be classified as continental winds, reaching Costa Rica’s Central Valley from the Pacific Ocean. Most of these air masses move over industrial and urban areas of the country. In the dry season, trade winds bring air masses to Costa Rica from the Caribbean Sea, and they move over non-industrial and less populated areas. Air masses arriving over Costa Rica during the transitional season can be transported by continental or trade-winds, and a weakening in the trade winds reaching Costa Rica is usually observed during the afternoon. Due to the transport pattern of the air masses observed in Costa Rica throughout the year, the isotopic composition of precipitation in Costa Rica is expected to vary depending on its origin from either the Caribbean Sea or the Pacific Ocean and also on the air trajectory across the country (i.e. industrial-urban areas versus remote-forested areas). Annual precipitation varies from ~1500 mm in the northwestern region and Nicoya Peninsula to ~7000 mm on the Caribbean side of the Talamanca range. Temperature seasonality is low. The mean annual temperatures vary from around 27˚C on the coastal lowlands, 20˚C in the Central Valley, and below 10˚C at the summits of the highest mountain range (~4000 m.a.s.l.).

The Global Network for Isotope in Precipitation (GNIP) [36] data set provides isotopic values of precipitation from 46 stations across Costa Rica between 1990 and 2005 (intermittent collection over this period). Most of the reported stations were located within the Central Valley (greatest populated area) and the North Pacific region. Therefore, a considerable gap exists along the Pacific and Caribbean shores, continental divide, and the northern and southern regions. Herein, we present an analysis of the spatial and temporal variation of stable isotopes in meteoric waters in Costa Rica and determine which broad climatic or topographic factors may influence isotopic variations across the region. Results of this analysis provide baseline historical spatial-temporal characterization of meteoric waters in Costa Rica that has been lacking to date.

2. Methods and Materials

2.1. Study Area

The Central Valley is a heavily populated area containing four major cities of Costa Rica located within its boundaries (~60% population) along with significant industrial activity. Weather conditions in the valley are mainly influenced by easterly Caribbean trade-winds during the dry season (February-April) and by frequent westerly continental air masses during the wet season (May-October). About ~50% of GNIP were located within the valley limits. In this region, mean annual precipitation ranges from ~3500 mm in the highlands to ~2000 mm within the foothills (Figure 2), and the area experiences 3 - 4 dry months over the summer. This inter-mountainous valley region has an elevation between 1000 and 2500 meters m.a.s.l, with average temperatures that range from 14˚C to 18˚C. The northern Pacific region of the Nicoya Peninsula receives on average 1500 mm of precipitation and suffers extensive dry periods that can reach up to 6 months. The Caribbean slope is characterized by a short dry season (1 - 2 months); mean annual precipitation in this region ranges from 4500 to 6000 mm.

2.2. Historic GNIP Records Analysis

We analyzed stable isotope records (n = 679) of 46 monitoring stations across Costa Rica obtained from the GNIP data base (Table 1, Figure 1). Station altitude ranges from 0 - 3000 m.a.s.l with a mean altitude of 930 m.a.s.l. Sampling frequency was typically monthly except for one station (Monteverde) where event-based sampling was conducted by Rhodes et al. (2006). The entire data set is comprised of three relative short-term sampling campaigns: 1990-1992, 1997-1998, and 2002- 2005. Weighted d18O and d2H means were calculated for 42 stations (87%). Linear regression analysis of elevation and arithmetic d18O values was conducted considering 18 stations which are unlikely to be influenced by evaporation enrichment: La Selva, Turrialba, Santa María, San Pablo, Santa Lucía, Birrí, Monteverde, Zurquí, Heinz Hoffman, Monte de la Cruz, Pacayas, Fraijanes, Vara Blanca, La Giralda, Sacramento, Paso Llano, Poás, and Irazú. Likewise, linear regression analysis of precipitation amount and d18O values was determined using data from 11 stations which are representative of all four regions: Turrialba, Puriscal, Fraijanes, Orotina, Santa Barbara, Rancho Redondo, Puntarenas, Hacienda Tempisque, and Sacramento.  

Spatial analysis was conducted using the Atlas of Costa Rica 2008 database [37]. The natural neighbor procedure in ArcGIS 9.3 was applied to calculate deuterium excess, d18O, and d2H interpolations. High station densities within the Central Valley and Nicoya Peninsula may overemphasize the isotopic signature of these regions. Data limitations constrained this and additional analyses. First, the majority of the stations (61%) generally reported isotopic values over a 1 year period with monthly values ranging from n = 4 to 28 with a mean of n = 15. Second, most of the reported stations were located within the Central Valley (the most populated area) and the North Pacific region (Figure 1). Therefore, a considerable gap exists along the Pacific and Caribbean shores, continental divide, and the northern and southern regions (Figure 1). Third, there was an omission of temperature values which prevented calculation of d18O and temperature correlation.

2.3. HYSPLIT Model Simulations

To complement the historic GNIP station data analysis, air parcel back trajectories were calculated for particular locations using the hybrid single particle Lagrangian integrated trajectory model (HYSPLIT) [38]. Isotopic values in precipitation from two distinct study locations, Monteverde (continental divide) and Turrialba (Caribbean Slope), were analyzed. Both a depleted and enriched isotopic event during the 2011 and 2012 water years were chosen from these two locations. Single (24 hours) five days back trajectories were calculated using the vertical velocity model calculation method. Trajectory ensembles were done using the GDAS1 meteorological database. Back trajectories were calculated starting at 2300 UTC. Rainfall rate, downward radiation,

Figure 1 . (Top) Digital elevation map of Costa Rica and location of historic GNIP stations; (bottom) Inset showing Central Valley stations. DEM source: Atlas of Costa Rica 2008 (ITEC); 100 m resolution.

Table 1. Summary of GNIP stations operated in Costa Rica from 1990 to 2005.

relative humidity, and elevation of the mixing layer depth were also computed.

2.4. Stable Isotope Analysis

Precipitation was collected on event-basis. Monteverde sampling was carried out during July, August, and September of 2011. Turrialba’s sampling period used in this study was from September to December, 2012. The collector for individual events consisted of a 10.16 cm plastic funnel coupled with a metal filter mesh (i.e., to prevent external contamination). The funnel was connected to 4 L high density polyethylene (HDPE) container. An approximately 3 cm mineral oil layer was added to prevent fractionation according with standard sampling protocols. Later, mineral oil was separated using a 250 - 500 mL separatory funnel. Samples were stored upside down in a HDPE bottles with conic and polyseal inserts and parafilm seals until analysis. Stable isotope analyses were conducted at the Chemistry School of the National University (Heredia, Costa Rica) and the Idaho Stable Isotope Laboratory (Moscow, Idaho) using a Cavity Ring Down Spectroscopy (CRDS) water isotope analyzer L2120-I (Picarro, CA) and L1120-i (Picarro, CA), respectively. Ratios of d18O/d16O and d2H/d1H are expressed in delta units (‰, parts per mil) relative to Vienna Standard Mean Ocean Water (V-SMOW). The current analytical analyzers precision is 0.1‰ d18O/d16O and 0.5‰ d2H/d1H.

3. Results and Discussion

3.1. Precipitation

Precipitation variability in Costa Rica reflects regional and local complex dynamics. Figure 2 shows mean annual precipitation throughout Costa Rica and precipitation amounts recorded during isotopic sampling campaigns. GNIP precipitation records are consistent with the Atlas of Costa Rica reported values and ranged from 5074 mm (Zurquí) to 897 mm (Hacienda Tempisque) with a mean of 2479 mm (Table 1, Figure 2). The lowlands of the Nicoya Peninsula is the region with the lowest precipitation regime (<1500 mm) while the Talamanca Range and northeastern Caribbean lowlands received the greatest amount, between 4500 - 8000 mm

Figure 2. Mean annual precipitation throughout Costa Rica with GNIP stations showing precipitation recorded during sampling periods.

annually. These precipitation records provide a reliable basis for further stable isotope spatial analysis.

3.2. (18O and (2H Isotopic Variations

Stable isotope signatures in precipitation across Costa Rica were variable during the sampling campaigns. d18O composition ranged from −1.8‰ to −10.1‰ with an arithmetic mean of −6.7 ± 1.4 (‰) (Figure 3, Table 1). Weighted d18O values range from −3.2‰ to −11.1‰ with a mean of −7.4 ± 1.8 (‰) (Table 1). Deuterium composition varied from −17.3‰ to −70.4‰ with an arithmetic mean of −43.2 ± 10.3 (‰) (Figure 3, Table 1). Weighted d2H values ranged from −19.7‰ to −76.8‰ with a mean of −49.0 ± 13.7 (‰) (Table 1). Figure 4 shows the regional water meteoric lines for each of the four main regions for the GNIP Stations (Nicoya Peninsula, Pacific Coast, Caribbean Slope, and Central Valley). The best fit straight line to the data for the Nicoya Peninsula yields an arithmetic water line equation of d2H = 7.84d18O + 7.60 (Figure 4), whereas the weighted meteoric water line yields a distinct best fit straight line equation of d2H = 6.65d18O − 0.13 (Figure 4), likely resulting from small precipitation amounts. For the Pacific Coast, the best fit straight line to precipitation-weighted annual means is d2H = 7.60d18O + 7.95 (Figure 4), which is consistent with the Global Meteoric Water Line (GMWL) of d2H = 8d18O + 10 [39]. Interestingly, Pacific Coast stations are close to equilibrium conditions rather than exhibiting evaporation enrichment. The Caribbean Slope precipitation-weighted water line d2H = 6.98d18O + 4.97 (Figure 4) reflects orographic distillation of Caribbean-sourced moisture. Overall, the best fit line to the precipitation-weighted data for all regions combined yields a water line of d2H = 7.61d18O + 7.40 (r2 = 0.98) which corresponds closely to the GMWL (Figure 5). Enriched values above the annual mean of −6.7‰ d18O were often observed during January to March. The greater variation occurred in May and June when the ITCZ starts its migration over Costa Rica. Towards the middle of the wet season (August, September, and October) d18O variation appears to decrease due to increase in precipitation amounts. The intensification of northeast trade winds in December produced a decrease in precipitation, and consequently, enrichment in isotopic values (Figure 5).

3.3. The Altitude Effect

The variability in isotopic composition of precipitation with elevation is not consistent across the GNIP stations indicating a complex system. Figure 6 shows three lon-

Figure 3. Spatial variations of d18O (top) and d2H (bottom) arithmetic means.

Figure 4. Meteoric water lines for four major geographic regions of Costa Rica.

gitudinal transects. For the Caribbean slope, there is an evident orographic distillation effect from Puerto Limón (0 m.a.s.l.) to Taus (900 m.a.s.l.) of −3.5‰ d18O/km. The ~75 km Pacific slope transect from Herradura (3 m.a.s.l.) to Paso Llano (2397 m.a.s.l.) represents a more complex system. For instance, a rainout effect is observed up to the Santa Bárbara station (1102 m.a.s.l.) (Figure 6). However, at the two higher elevation stations (Birrí and Heinz Hoffman) isotopic compositions are more enriched, suggesting another moisture source. A similar phenomena is also observed along the Nicoya Peninsula transect of ~66 km from the coast to the Rincón de la Vieja volcano (Figure 6).

Considering the 18 stations mainly from the Central Valley cluster, the correlation of elevation with arithmetic mean annual d18O is moderately strong (r2 = 0.77). The best fit straight line has a slope of −0.2‰ d18O per 100 m, which is close to the range of 0.1‰ and 0.36‰ reported by Leibungdgut et al. (2009) [40]. Since the altitude effect is the combined result of the temperature effect and moisture depletion by adiabatic cooling, a temperature effect can be calculated by knowing the atmospheric lapse rate. Lachniet and Paterson (2002) [27] reported a lapse rate of −5.4˚C/km for the San José area. Therefore, an altitude effect of −0.2‰/km corresponds to −0.37‰/˚C, in agreement with temperature effects reported in the pioneer work done by Dansgaard [23].

3.4. The Amount Effect and d-Excess Isoscape

One important factor controlling isotopic variations in the tropics is the “amount effect” [23,24,41]. This effect has been attributed to three main factors according to Scholl et al. (2009) [42]: 1) the isotopic composition of the condensation in a cloud decreases as cooling and “rainout” occurs; 2) smaller raindrops equilibrate to a larger degree with the water vapor and temperature conditions below the cloud; and 3) small raindrops evaporate more than larger raindrops on their way to the land

Figure 5. (Top) Meteoric water line for Costa Rica. Contributions from all four regions are color coded. Inset: distribution of d18O values; (Bottom) Temporal variation of d18O values.

surface. The slope of the best fit straight line to precipitation amount and d18O from 11 stations, which are representative of the four main regions described previously, was found to be −1.6‰ per 100 mm decrease in rainfall (r2 = 0.57). This effect is clearly present during the times of lowest precipitation (January to April).

The “amount effect” was also observed in the d-excess isoscape presented in Figure 7. D-excess composition ranged from 2.9‰ to 14.4‰ with an arithmetic mean of 10.8 ± 2.0 (‰) (Table 1, Figure 7). Weighted d-excess values range from 3.5‰ to 17.4‰ with a mean of 10.3 ± 2.1 (‰) (Table 1, Figure 7). Low d-excess (~3‰) is predominantly present within the Nicoya Peninsula and Caribbean lowlands and may also be a common phenomena across the northern lowlands and Nicaragua trough. High d-excess (>17‰) was found in the dense forested areas east of the Central Valley (Zurquí and Irazú mountain range) and the Pacific slope of the Monteverde tropical cloud forest (Figure 7), which may indicate successive evaporation and precipitation cycles. Based on d-excess values, Rhodes et al. (2006) [43] concluded that moisture flux from evaporation is an important input to the region during the transitional and dry seasons when trade winds from the Caribbean slope dominate. Therefore, it appears that intense moisture recycling processes govern isotopic variations within the

Figure 6. Longitudinal elevation and isotopic transects of Nicoya Peninsula, Pacific and Caribbean slopes.

inter-mountainous regions of Costa Rica.

3.5. Air Mass Back Trajectories and Precipitation Events

Figure 8 shows daily variation in isotopic composition in precipitation at Monteverde and Turrialba. For Monteverde, a remarkable depletion in δ18O was observed on 7/18/2011 and an enriched event was sampled on 7/26/ 2011 (Figure 8). For Turrialba, sampling coincided with the wet to dry season transition, but sampling events for the purposes of the air mass trajectories were selected

Figure 7. Preliminary d-excess isoscape for Costa Rica.

Figure 8. Daily d18O variations in precipitation at Monteverde and Turrialba stations.

during the wet and transitional seasons (September-December, 2011). Considerably depleted and enriched δ18O values were recorded on 10/17/2011 and 10/27/2011, respectively (Figure 8). Figure 9 shows air mass back trajectories of five days prior to the depleted or enriched precipitation event coupled with rainfall intensities. HYSPLIT simulated trajectories indicate that enriched isotopic compositions both in Turrialba and Monteverde correspond to central Caribbean-sourced moisture and small rainfall intensities (Figures 9(a) and (b)). Depleted events appear to correspond to high rainfall amounts despite the parental origin of the moisture. In the case of Turrialba, the depleted δ18O event was related to continental air masses arriving over Costa Rica from the southwest, producing convective precipitation as they travel along and over the Pacific Coast. Air masses intrusions from the Pacific Ocean are typically observed during the wet season and are less frequent when the wet to dry season transition begins in November. The depleted precipitation sampled in Monteverde was related to greater rainfall rates and parental moisture from the southeastern Caribbean basin.

4. Conclusions

Stable isotope signatures in precipitation and water vapor are useful indicators of climate variability and provide critical information regarding regional and global water processes. In temperate regions, isotopic sampling has been conducted for decades, providing novel insights on water cycle dynamics. However, isotopic records from the humid tropics, such as Central America, have been quite limited to date. Long-term d18O and d2H records


Figure 9. HYSPLIT model simulations including rainfall intensity in mm/hr. Depleted events: (a) Monteverde 7/18/11 and (b) Turrialba 10/17/11. Enriched events: (c) Monteverde 7/26/11 and (d) Turrialba 10/27/11. Five-day air mass back trajectories are color coded.

offer a unique opportunity to study precipitation and ground and surface water processes at a relatively low cost compared to continuous streamflow or climate monitoring, which is often not feasible in tropical developing countries due to high costs and limited resources. Based on historic GNIP records, we conducted a spatial and temporal characterization of stable isotopes in four main geographical regions of Costa Rica: Nicoya Peninsula, Pacific Coast, Caribbean slope, and Central Valley. Despite data limitations several correlations were found. The regression analysis of precipitation amounts with annual arimethic means in samples from all four regions yields a slope of −1.6‰ d18O per 100 mm of rain (r2 = 0.57), which corresponds with a temperature effect of −0.37‰ d18O/˚C. A strong correlation (r2 = 0.77) of −2.0‰ d18O per km of elevation was also found in these data.

D-excess has been identified as a potential indicator of recycled moisture. Therefore, meteoric d-excess isoscapes coupled with isotopic monitoring in surface waters may provide quantitative estimates of the contribution of recycled moisture, which would greatly enhance our understanding of water budgets in Costa Rica. D-excess interpolation from 46 GNIP stations suggests that evaporation enrichment often occurs over the Nicoya Peninsula, especially during periods of lowest precipitation (January-April). The Nicoya Peninsula, which receives less than 1500 mm annually, has been recently identified [44] as a hydrologically vulnerable region due to the increase of water demand and susceptibility to experience prolonged drought periods. Greater d-excess values were observed in the inter-mountainous region of the Central Valley and the Pacific slope. It is still unknown if moisture recycling is occurring due to evaporation within dense canopy covers (e.g., Zurquí-Irazú or Monteverde areas) or if instead moisture is available from multiple sources such as flood plains, dams or reservoirs, or soil moisture. However, this moisture input may play a significant role in late afternoon and overnight convective storms along the Pacific slope. 

The illustrative example using the HYSPLIT model indicates that enriched isotopic compositions both in Turrialba and Monteverde are related to central Caribbean parental moisture and small rainfall intensities. Depleted events appear to be related to high rainfall amounts despite the parental moisture origin. Further event-based sampling in strategic locations across the country coupled with air mass trajectory analysis is needed to evaluate the effects of parental moisture origin and isotope signatures. Overall, we believe that this preliminary stable isotope analysis across Costa Rica may provide a useful baseline for future isotopic studies in the region as well as enhance our understanding of moisture recycling processes within tropical inter-mountainous systems. Further research should evaluate spatial and temporal isotopic variations in surface and groundwater systems to ultimately determine critical recharge zones across the country.

5. Acknowledgements

We thank the International Atomic Energy Agency for granting access to historic GNIP records. A National Science Foundation-IGERT scholarship (Grant No. 0903479) US Borlaug Fellowship for Food Security supported the isotopic fieldwork in Turrialba, Costa Rica. Funding provided by the Chemistry School at the National University of Costa Rica and the National University Council of Costa Rica was essential to acquire a laser water isotope analyzer. The authors are also grateful for access to the field sites and logistical support by the Tropical Science Center and the Monteverde Cloud Forest Reserve.


  1. P. Kabat, R. E. Schulze, M. E. Hellmuth and J. A. Veraart, “Coping with Impacts of Climate Variability and Climate Change in Water Management: A Scoping Paper,” DWCReport No. DWCSSO-01 International Secretariat of the Dialogue on Water and Climate, Wageningen, 2003.
  2. R. K. Colwell, G. Brehm, C. L. Cardelús, A. C. Gilman and J. T. Longino, “Global Warming, Elevational Range Shifts, and Lowland Biotic Attrition in the Wet Tropics,” Science, Vol. 322, No. 5899, 2008, pp. 258-261.
  3. S. R. Loarie, P. P. Duffy, H. Hamilton, G. P. Asner, C. B. Field and D. D. Ackerly, “The Velocity of Climate Change,” Nature, Vol. 462, No. 7276, 2009, pp. 1052-1055.
  4. J. J. Tewksbury, R. B. Huey and C. A. Deutsch, “Putting the Heat on Tropical Animals,” Science, Vol. 320, No. 5881, 2008, pp. 1296-1297.
  5. Y. P. Zhou, K. M. Xu, Y. C. Sud and A. K. Betts, “Recent Trends of the Tropical Hydrological Cycle Inferred from Global Precipitation Climatology Project and International Satellite Cloud Climatology Project Data,” Journal of Geophysical Research, Vol. 116, No. D6, 2011, Article ID: D09101.
  6. V. Magaña, J. A. Amador and S. Medina, “The Midsummer Drought over Mexico and Central America,” Journal of Climate, Vol. 12, No. 6, 1999, pp. 1577-1588.<1577:TMDOMA>2.0.CO;2
  7. G. Bowen, “Statistical and Geostatistical Mapping of Precipitation Water Isotope Ratios,” In: J. B. West, et al., Eds., Isoscapes: Understanding Movement, Pattern, and Process on Earth through Isotope Mapping, Springer Science, New York, 2010, pp. 139-160.
  8. L. J. Araguás-Araguás, K. O. Froehlich and K. Rozanski, “Stable Isotope Composition of Precipitation over Southeast Asia,” Journal of Geophysical Research, Vol. 103, No. D22, 1998, pp. 28721-28742.
  9. K. M. Cobb, J. F. Adkins, J. W. Partin and B. Clark, “Regional-Scale Climate Influences on Temporal Variation of Rainwater and Cave Dripwater Oxygen Isotopes in Northern Borneo,” Earth and Planetary Science Letters, Vol. 263, No. 3-4, 2007, pp. 207-220.
  10. Y. Ishizaki, K. Yoshimura, S. Kanae, M. Kimoto, N. Kurita and T. Oki, “Interannual Variability of H218O in Precipitation over the Asian Monsoon Region,” Journal of Geophysical Research, Vol. 117, No. D16, 2012, Article ID: D16308.
  11. K. R. Johnson and B. L. Ingram, “Spatial and Temporal Variability in the Stable Isotope Systematics of Modern Precipitation in China: Implications for Paleoclimate Reconstructions,” Earth and Planetary Science Letters, Vol. 220, No. 3-4, 2004, pp. 365-377.
  12. M. Lachniet, “Sea Surface Temperature Control on the Stable Isotopic Composition of Rainfall in Panama,” Geophysical Research Letters, Vol. 36, No. 3, 2009, Article ID: L03701.
  13. M. Lachniet and W. Patterson, “Oxygen Isotope Values of Precipitation and Surface Waters in Northern Central America (Belize and Guatemala) Are Dominated by Temperature and Amount Effects,” Earth and Planetary Science Letters, Vol. 284, No. 3, 2009, pp. 435-446.
  14. M. Vuille, R. S. Bradley and F. Keimig, “Interannual Climate Variability in the Central Andes and Its Relation to Tropical Pacific and Atlantic Forcing,” Water Resources Research, Vol. 105, No. D10, 2000, pp. 12447-12460.
  15. M. Vuille, R. S. Bradley and F. Keimig, “Climatic Variability in the Andes of Ecuador and Its Relation to Tropical Pacific and Atlantic Sea Surface Temperature Anomalies,” Journal of Climate, Vol. 13, No. 14, 2000, pp. 2520-2535.<2520:CVITAO>2.0.CO;2
  16. M. Vuille, R. S. Bradley, R. Healy, M. Werner, D. R. Hardy, L. G. Thompson and F. Keimig, “Modeling d18O in Precipitation over the Tropical Americas 2: Simulation of the Stable Isotope Signal in Andean Ice Cores,” Journal of Geophysical Research, Vol. 108, No. D6, 2003, p. 4175.
  17. K. Ichiyanagi and M. D. Yamanaka, “Interannual Variation of Stable Isotopes in Precipitation at Bangkok in Response to El Niño Southern Oscillation,” Hydrological Processes, Vol. 19, No. 17, 2005, pp. 3413-3423.
  18. H. O. Panarello and C. Dapeña, “Large Scale Meteorological Phenomena, ENSO and ITCZ, Define the Paraná River Isotope Composition,” Journal of Hydrology, Vol. 365, No. 1-2, 2009, pp. 105-112.
  19. M. Vuille and M. Werner, “Stable Isotopes in Precipitation Recording South American Summer Monsoon and ENSO Variability: Observations and Model Results,” Climate Dynamics, Vol. 25, No. 4, 2005, pp. 401-413.
  20. P. K. Aggarwal, O. A. Alduchov, K. O. Froehlich, L. J. Araguas-Araguas, N. C. Sturchio and N. Kurita, “Stable Isotopes in Global Precipitation: A Unified Interpretation Based on Atmospheric Moisture Residence Time,” Geophysical Research Letters, Vol. 39, No. 11, 2012, Article ID: L11705.
  21. L. J. Araguás-Araguás, K. O. Froehlich and K. Rozanski, “Deuterium and Oxygen-18 Isotope Composition of Precipitation and Atmospheric Moisture,” Hydrological Processes, Vol. 14, No. 8, 2000, pp. 1341-1355.<1341::AID-HYP983>3.0.CO;2-Z
  22. G. Bowen, “Spatial Analysis of the Intra-Annual Variation of Precipitation Isotope Ratios and Its Climatological Corollaries,” Journal of Geophysical Research, Vol. 113, No. D5, 2008, Article ID: D05113.
  23. W. Dansgaard, “Stable Isotopes in Precipitation,” Tellus, Vol. 16, No. 4, 1964, pp. 436-468.
  24. K. Rozanski, L. J. Araguas-Araguas and R. Gonantini, “Isotopic Patterns in Modern Global Precipitation,” In: P. K. Swart, K. C. Lohmann, J. McKenzie and S. Savin, Eds., Climate Change in Continental Isotopic Records, American Geophysical Union, Washington DC, 1993, pp. 1-36.
  25. L. I. Wassenaar, P. Athanasopoulos and M. J. Hendry, “Isotope Hydrology of Precipitation, Surface and Ground Waters in the Okanagan Valley, British Columbia, Canada,” Hydrology, Vol. 411, No. 1, 2011, pp. 37-48.
  26. G. Bowen and J. Revenaugh, “Interpolating the Isotopic Composition of Modern Meteoric Precipitation,” Water Resources Research, Vol. 39, No. 10, 2003, p. 1299.
  27. M. Lachniet and W. Patterson, “Stable Isotope Values of Costa Rican Surface Waters,” Journal of Hydrology, Vol. 260, No. 1-4, 2002, pp. 135-150.
  28. G. J. Bowen, J. B. West and J. Hoogewerff, “Isoscapes: Isotope Mapping and Its Applications,” Geochemical Exploration, Vol. 102, No. 3, 2009, pp. 5-7.
  29. V. Barras and I. Simmonds, “Observation and Modeling of Stable Water Isotopes as Diagnostics of Rainfall Dynamics over Southeastern Australia,” Journal of Geophysical Research, Vol. 114, No. D23, 2009, pp. 1-17.
  30. H. Celle-Jeanton, R. Gonfiantini, Y. Travi and B. Sol, “Oxygen-18 Variations of Rainwater during Precipitation: Application of a Rayleigh Model to Selected Rainfalls in Southern France,” Journal of Hydrology, Vol. 289, No. 1-4, 2004, pp. 165-177.
  31. T. B. Coplen, P. J. Neiman, A. B. White, J. M. Landwehr, F. M. Ralph and M. D. Dettinger, “Extreme Changes in Stable Hydrogen Isotopes and Precipitation Characteristics in a Landfalling Pacific Storm,” Geophysical Research Letters, Vol. 35, No. 21, 2008, Article ID: L21808.
  32. N. C. Munksgaard, C. M. Wurster, A. Bass and M. I. Bird, “Extreme Short-Term Stable Isotope Variability Revealed by Continuous Rainwater Analysis,” Hydrological Processes, Vol. 26, No. 23, 2012, pp. 3630-3634.
  33. J. Reynolds-Vargas and J. Fraile, “Use of Stable Isotopes in Precipitation to Determine Recharge Areas within the Barva Aquifer, Costa Rica,” Latinoamerican Stable Isotopes Studies, International Atomic Energy Agency, Viena, 2006, pp. 83-96.
  34. P. Waylen, C. N. Caviedes and M. F. Quesada, “Interannual Variability of Monthly Precipitation in Costa Rica,” Journal of Climate, Vol. 9, No. 10, 1996, pp. 2506-2613.<2606:IVOMPI>2.0.CO;2
  35. A. J. Guswa, A. Rhodes and S. E. Newell, “Importance of Orographic Precipitation to the Water Resources of Monteverde, Costa Rica,” Advances in Water Resources, Vol. 30, No. 10, 2007, pp. 2098-2112.
  36. GNIP, “Global Network of Isotopes in Precipitation,” 2012.
  37. Instituto Tecnológico de Costa Rica (ITEC), “Atlas de Costa Rica. Escuela de Ingeniería Forestal. Cartago, Costa Rica,” Digital Atlas, ITEC, Cartago, Costa Rica, 2008.
  38. R. R. Draxler and G. D. Rolph, “HYSPLIT-Hybrid Single Particle Lagrangian Integrated Trajectory Model,” 2013.
  39. H. Craig, “Isotopic Variations in Meteoric Waters,” Science, Vol. 133, No. 3465, 1961, pp. 1702-1703.
  40. C. Leibundgut, P. Maloszewski and C. Kulls, “Tracers in hydrology,” Wiley-Blackwell, Oxford, 2009.
  41. R. Gonfiantini, M. A. Roche, J. C. Olivry, J. C. Fontes and G. M. Zuppi, “The Altitude Effect on the Isotopic Composition of Tropical Rains,” Chemical Geology, Vol. 181, No. 1, 2001, pp. 147-167.
  42. M. A. Scholl, J. B. Shanley, J. P. Zegarra and T. B. Coplen, “The Stable Isotope Amount Effect: New Insights from NEXRAD Echo Tops, Luquillo Mountains, Puerto Rico,” Water Resources Research, Vol. 45, No. 12, 2009, Article ID: W12407.
  43. A. L. Rhodes, A. J. Guswa and S. E. Newell, “Seasonal Variation in the Stable Isotopic Composition of Precipitation in the Tropical Montane Forests of Monteverde, Costa Rica,” Water Resources Research, Vol. 42, 2006, Article ID: W11402.
  44. MINAET-IMN, “Mejoramiento de las Capacidades Nacionales para la Evaluación de la Vulnerabilidad y Adaptación del Sistema híDrico al Cambio climÁtico en Costa Rica, Como Mecanismo para Disminuir el Riesgo al Cambio climÁtico y Aumentar el Índice de Desarrollo Humano,” In: J. A. Retana, Eds., Ministerio de Ambiente, Energía y Telecomunicaciones, Instituto Meteorológico Nacional San José, Costa Rica, 2012, 46 p.


*Corresponding author.