Data Descriptor: Daily observations of stable isotope ratios of rainfall in the tropics

We present precipitation isotope data (δ2H and δ18O values) from 19 stations across the tropics collected from 2012 to 2017 under the Coordinated Research Project F31004 sponsored by the International Atomic Energy Agency. Rainfall samples were collected daily and analysed for stable isotopic ratios of oxygen and hydrogen by participating laboratories following a common analytical framework. We also calculated daily mean stratiform rainfall area fractions around each station over an area of 5° x 5° longitude/latitude based on TRMM/GPM satellite data. Isotope time series, along with information on rainfall amount and stratiform/convective proportions provide a valuable tool for rainfall characterisation and to improve the ability of isotope-enabled Global Circulation Models to predict variability and availability of inputs to fresh water resources across the tropics.

across the tropics. Several sampling stations with no prior rainfall isotope observations were established. It is noteworthy that this data set is the first systematic effort to capture ground-based daily isotope record of tropical rainfall in all phases (i.e. neutral, warm, and cool) of the El Niño-Southern Oscillation (ENSO). Furthermore, the 2015-2016 ENSO event was one of the strongest on record 3 .
The availability of Global Circulation Models (GCMs) incorporating isotope physics provided a promising framework to study isotopic variability in precipitation [4][5][6][7][8] . To evaluate the stable isotope fields simulated by such GCMs, the IAEA's Global Network of Isotopes in Precipitation (GNIP) database 9,10 has recurrently been used as the major comprehensive source of data. While isotope-enabled GCMs simulate modern isotopic patterns in global precipitation on monthly and inter-annual time scales reasonably well over mid and high latitude regions, simulations of tropical rainfall need further development 6,8 . A major reason for this problem is that the physical mechanisms controlling the spatial and temporal isotopic variations in tropical precipitation are still not fully understood. In addition, the poor spatial and temporal coverage of precipitation isotope data in tropical regions poses a challenge to explore the primary drivers of isotope variability. Therefore, an expanded monitoring network in the tropics is required, not only to identify the key controls of isotopic variability, but to improve the reproducibility of GCMs simulations. This improvement may help elucidate the mechanisms (e.g. sub-cloud evaporation, moisture convergence and entrainment) that control isotopic changes in precipitation and enhance the capabilities of climate models to predict variability and availability of fresh water resources. In addition, better understanding of the controls on precipitation isotope variability will lead to improvements in our interpretation of isotope-based proxies in terrestrial and maritime paleoarchives (e.g., caves, corals, lake sediments, and potentially tree-rings) 11 .
In the tropics, precipitation isotopes mainly correlate negatively with precipitation amounts on a monthly scale. This empirical low-latitude inverse correlation between water isotope ratios and the amount of rainfall, known as the 'amount effect' 12 , has been used as a rationale to infer wet and dry paleo-hydroclimate periods based on available maritime and terrestrial paleo proxies across the tropics 13 . However, the amount effect is not universal over the tropical regions. For example, in SE Asia the effect is relatively strong at two marine island stations (Palau and Bali) but relatively weak at continental coastal stations (e.g., Bangkok and Da Nang) 14 . Furthermore, the correlation between precipitation amount and isotopic composition is commonly weaker or non-existent on a daily basis compared to monthly time scales 2,14,15 . This is because the linear regression approach neglects other processes such as moisture convergence and entrainment, resulting in weak correlations across the tropics with >80% of the variance unexplained when using daily data, whereas stronger correlations (30-70% variance explained) are reported when computing monthly means 2 . These complexities indicate that the isotopic variations in precipitation, even at tropical islands where a pronounced amount-effect is observed, are not directly controlled by rainfall amount but is rather influenced by the other convection related processes (e.g. cloud microphysics, cloud type, moisture transport). This conclusion is supported by several recent studies which have demonstrated that the isotopic variability is associated with regional, rather than local, convective activity (e.g. in North Africa 16 , Tibet 17 , Borneo 1 , Pacific Ocean 18,19 , India 20,21 , Australia 22 , Brazil 23 and Costa Rica 15 ). In Costa Rica, isotopic variations also exhibited more complex interactions between variable moisture sources, humidity and orographic lifting 15 . A gradual decrease in δ 18 O values of precipitation with cumulative rainfall along upstream air mass trajectories over several days was reported in Tibet 17 , Northern Australia 22 , and Southeast Asia 24 . On the other hand, in some amount-effect-dominated regions, changes in moisture source and air mass trajectories largely influenced sub-monthly or seasonal isotopic variations (e.g. East Africa 25,26 , Japan 27 , Namibia 28 ).
It is well known that intra-seasonal isotopic variations are clearly seen in many tropical regions (e.g. Borneo 1 , Western Pacific 19 , Northern Australia 22 ). These variations typically manifest as negative excursions in δ 18 O and δ 2 H values and typically, but not exclusively, correspond to the wet phases of the Madden Julian Oscillation (MJO). Cyclonic lows (i.e., hurricanes, typhoons) may also produce extreme rainfall and negative δ 18 O and δ 2 H anomalies. Within the MJO wet phase, large organised convective cloud systems, referred to as mesoscale convective systems (MCSs), account for a large portion of tropical rainfall. Stratiform rainfall associated with MCSs has been shown to be mainly associated with large negative excursions of δ 18 O and δ 2 H values in tropical rainfall 19,29 . Negative isotope anomalies in stratiform rainfall have been linked to deposition of 18 O-and 2 H-depleted water vapour onto ice particles at altitude which, as they fall, aggregate and melt at mid tropospheric levels 29,30 . Most recently, the relative influence of bulk precipitation microphysics, cloud type, and surface moisture transport on precipitation amounts and 18 O/ 16 O ratios was assessed in the tropics 2 . This analysis showed that bulk precipitation microphysics and cloud type (i.e., stratiform rain fraction) exert comparable influences on the isotopic composition of precipitation, whereas moisture transport plays an important secondary role in regions of deep atmospheric convection. Reduced Outgoing Longwave Radiation (OLR) values are indicative of stronger MCS activity and often correlate with low δ 18 O values in rainfall 20,22 . Based on these results, we can hypothesise that stratiform rainfall is a major driver of isotopic variability over the amount-effect dominated region. However, currently both spatial and temporal coverage of precipitation isotope data is insufficient to adequately test the hypothesis.
Here we present stable isotope data for daily precipitation collected at nineteen stations at both maritime and continental locations within the tropics; ten stations are (near-) coastal and nine are located from ≈ 80 to 600 km inland. We also calculated daily mean stratiform rainfall area fractions at these stations using TRMM and Global Precipitation Mapping (GPM) satellite observations.

Methods
Sampling and analysis. Details of the sampling stations including location, updated Köppen-Geiger climate zone classification 31 , mean annual precipitation and temperature, sampling period and number of samples collected are provided in Table 1 and Fig. 1. Rainfall samples were usually collected at 9am local time but variations of up to a few hours occurred on some occasions for practical reasons. calculation of Stratiform Rainfall fraction. The daily mean stratiform rainfall area fraction (F st ) was calculated to examine the influence of stratiform rainfall on the daily isotopic variability. F st is defined as the average percent of rainfall area covered by stratiform rainfall over the 5° x 5° longitude/latitude box centered over each isotope monitoring station. We used the Ku-band Precipitation Radar (KuPR) convective/stratiform classification data from version 5, level 2 product of GPM (Global Precipitation Measurement) Core Observatory (https://pmm.nasa.gov/GPM), which is a successor of the TRMM Precipitation Radar. The GPM satellite flies at an altitude of 407 km in a non-sun-synchronous orbit and completes roughly 16 orbits per day between 65° N and 65° S. The KuPR's horizontal footprint size along a track (swath width) is 245 km. Based on this swath width, 164 orbits are required to cover the whole equator (roughly 40,000 km) corresponding to around 10 days. This    Table 1.
www.nature.com/scientificreports www.nature.com/scientificreports/ means that only a few GPM orbits pass through the domain (5° x 5° longitude/latitude box) within 10 days and so we cannot obtain daily convective/stratiform classification data from the satellite. As an alternative approach, we used area averaged (5° x 5° box) daily precipitation to estimate the daily F st . In the tropics, stratiform rainfall area is significantly larger than the convective rainfall area, and the increase in rain area is more likely to result in larger total rain amounts 33,34 . We found statistically significant relationships of area-averaged rainfall amount (P area ) to F st over the study domains at fourteen stations (i.e., Cairns, Darwin, 28 Millas, Heredia, Barisal, Cox´s Bazar, Port Blair, Nanyang Tech., Mulu, Ho Chi Minh City (HCM City), Kuala Terengganu., Kuching, Abetifi, and Amedzofe) from 2014 to 2017 (p < 0.05, R2 range = 0.25-0.72, mean R 2 = 0.53). As shown in Fig. 2, combining data from the 14 stations show a strong correlation between P area and F st (R 2 = 0.87, p < 0.01). However, at five inland or slightly extra-tropical stations (i.e., Rio Claro, Addis Ababa, Windhoek, Nagoya and SP Reyes) the relationships were weak, and we excluded these stations from further analysis. We then applied the correlation shown in Fig. 2 to a satellite-based daily precipitation dataset to estimate the daily F st over the study domain at 14 stations. A Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis product, namely TMPA 3B42, was used to estimate area-averaged daily precipitation at each station. The TMPA 3B42 product has a 3-hourly temporal resolution and a 0.25° spatial resolution 35 . Original data are available online at http://disc. gsfc.nasa.gov/datacollection/TRMM_3B42_V7.shtml. Figure 3 shows the correlation between the 10-day moving average rainfall isotope data and stratiform rainfall area fractions >0.01 for 14 tropical stations. Since there is a time lag for organized convective cloud systems (MCSs) in  Table 2. Investigator, sampling and analytical information. www.nature.com/scientificreports www.nature.com/scientificreports/ the upwind region to arrive at rainfall sampling stations, the moving average improves correlation compared to using the daily data at most stations. The stations Heredia and 28 Millas in Costa Rica, Darwin and Cairns in Australia, Amedzofe in Ghana, Nanyang Tech in Singapore, Mulu in Malaysia and HCM City in Vietnam had the strongest correlations (R 2 = 0.28-0.58) while the remaining stations had weaker or insignificant correlations (p > 0.05) ( Table 3). Figure 4 shows a time series for the 10-day moving averages of rainfall isotope data and stratiform rainfall area fraction at Cairns, where a relatively strong anti-phase variation of these two parameters (R 2 = 0.43) is observed.

Technical Validation
All laboratory analyses were carried out using infrared laser-absorption spectrometry either by Off-Axis Integrated Cavity Output Spectroscopy (Los Gatos Research, San Jose, CA, USA) or by Cavity Ring Down Spectroscopy (Picarro, Santa Clara, CA, USA). These spectrometers can produce accurate and precise results provided that volatile organic compounds do not cause spectral interferences. However, such interferences are generally absent from rainfall samples and can easily be monitored using instrument software.
The datasets presented here were quality controlled by the individual investigators and laboratories. During the CRP participants were provided with guidance on the production, use and calibration of in-house (secondary) standard waters traceable to the International Measurement Standards VSMOW and SLAP (IAEA 2019) to  Table 1 for sampling period for each station and Table 3 for linear coefficients and correlation coefficient (R 2 ).  www.nature.com/scientificreports www.nature.com/scientificreports/ ensure analytical accuracy and optimise quality control procedures. Most of the CRP's laboratories participated in the IAEA Water Isotope Inter-Comparison WICO2016 36 as well as an unofficial CRP-wide inter-comparison.

Data Records
Data sets are available at figshare.com (https://figshare.com/s/fdfabb43a844cad530a5). The file 'CRP isotopes' contains stable isotope data in daily rainfall at 19 stations sampled within the period 2012 to 2017. They are recorded in the following order: Local sampling start date/time (YYYY-MM-DDTHH:MM), local sampling end date/time (YYYY-MM-DDTHH:MM), precipitation δ 18 O (‰ VSMOW), precipitation δ 2 H (‰ VSMOW), precipitation d-excess (‰ VSMOW), precipitation amount (mm). Note that at some stations only the sampling end date was recorded (sampling of 24-hr rainfall occurred at 9:00 local time). Empty cells indicate that no data was obtained.
The file 'CRP stratiform P' contains calculated stratiform precipitation area fraction for the subset of 14 stations for which this fraction could be calculated (five stations were excluded, see above). They are recorded in the following order: Local sampling start date/time (YYYY-MM-DDTHH:MM), local sampling end date/time (YYYY-MM-DDTHH:MM), sampled precipitation amount (mm), precipitation δ 2 H (‰ VSMOW), precipitation δ 18 O (‰ VSMOW), precipitation d-excess (‰ VSMOW), satellite observation date (YYYY-MM-DDT), daily mean area-averaged precipitation amount (mm), daily mean stratiform precipitation area-fraction, moving 10-day average precipitation δ 18 O (‰ VSMOW), moving 10-day average stratiform precipitation area-fraction. Stratiform rainfall area fractions <0.01 were disregarded due to their high uncertainty. Empty cells indicate that no data was obtained.

Usage notes
We encourage interested parties to contact the site investigators to explore possible collaboration opportunities based on these data. It is noted that some data have been published in peer-reviewed journals.