Evapotranspiration on Natural and Reclaimed Coral Islands in the South China Sea

remote sensing

Article Evapotranspiration on Natural and Reclaimed Coral Islands in the South China Sea

Shengsheng Han 1,2 , Suxia Liu 1,2,*, Shi Hu 1, Xianfang Song 1,2 and Xingguo Mo 1,2

1 Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China; [email protected] (S.H.); [email protected] (S.H.); [email protected] (X.S.); [email protected] (X.M.) 2 Sino-Danish Center, College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China * Correspondence: [email protected]

Abstract: Studies of evapotranspiration on remote tropical coral islands are important to explore and sustain scarce freshwater resources. However, there is a significant knowledge gap between research to evaluate evapotranspiration based on remote sensing methods and the influences of different land use types on water dynamics on reclaimed coral islands. This study applied the remote-sensing- based Vegetation Interfaces Processes (VIP-RS) model to estimate actual evapotranspiration (ETa) on Zhaoshu Island, Yongxing Island, and Yongshu Island in the South China Sea from 2016 to 2019. The results showed that the average annual ETa of Zhaoshu Island, Yongxing Island, and Yongshu Island was 685 mm, 530 mm, and 210 mm, respectively. Annual transpiration (Ec) and soil evaporation (Es) exhibited similar patterns on the natural islands; however, Es controlled the water consumption on the reclaimed islands. Water dynamics exhibited seasonal fluctuations due to the uneven distribution of precipitation (PRP). However, ETa of the natural islands was higher than PRP in the dry season,

indicating vegetation has to absorb water from the groundwater to sustain growth. The results also agreed with the analysis of dominant driving factors based on partial correlation analysis, which

Citation: Han, S.; Liu, S.; Hu, S.; demonstrated that the Normalized Difference Vegetation Index (NDVI) is the most important factor Song, X.; Mo, X. Evapotranspiration that influences ETa, while relative humidity (RH) controlled the bare land or sparsely vegetated on Natural and Reclaimed Coral areas on the reclaimed islands. The setting of different land use types showed that vegetation and Islands in the South China Sea. built-up or hardened roads took control of evapotranspiration and rainwater collection, respectively, Remote Sens. 2021, 13, 1110. which play important roles in water dynamics on corals islands. The evaluation of ETa based on https://doi.org/10.3390/rs13061110 a remote-sensing-based model overcame the difficulty in fieldwork observation, which improves the certainty and accuracy at a spatial scale. In addition, it gave us a new reference to protect and Received: 9 February 2021 manage scarce freshwater resources properly. Accepted: 7 March 2021 Published: 15 March 2021 Keywords: evapotranspiration; coral islands; water balance; VIP-RS model

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional afﬁl- 1. Introduction iations. Coral islands are widely distributed in the Paciﬁc, Indian, and Atlantic Oceans between 30◦N and 30◦S (Figure1a). Most of them have ecosystems vulnerable to both natural processes and anthropogenic activities [1,2]. Fresh water is one of the most important resources for the natural ecosystem and the people inhabiting the coral islands. In addition, Copyright: © 2021 by the authors. it also attracts much attention for its sensitivity to climate change [3,4]. Freshwater resources Licensee MDPI, Basel, Switzerland. on coral islands differ from those on the continent: (1) Freshwater resources on coral islands This article is an open access article distributed under the terms and mostly come from precipitation. However, it is easy for precipitation to penetrate the thin, conditions of the Creative Commons permeable unsaturated zone. There is nearly no ponding water on coral islands [5,6]. Attribution (CC BY) license (https:// (2) The fresh water is mainly stored in the form of groundwater, whose availability is creativecommons.org/licenses/by/ highly dependent on its resilience in drought periods. When droughts occur, the water of 4.0/). the coral islands gotten from precipitation decreases, and the vegetation may increase the

Remote Sens. 2021, 13, 1110. https://doi.org/10.3390/rs13061110 https://www.mdpi.com/journal/remotesensing Remote Sens. 2021, 13, 1110 2 of 18

consumption of water from groundwater [7,8]. (3) The low land surface elevation of atoll islands also makes fresh water vulnerable to inundation by seawater, such as storm surges, wave setup, extremely high tides, and tsunamis [2,9,10]. (4) The high population densities of some coral islands have led to unsustainable rates of extraction. Changes in land use types and contamination from human settlements can also signiﬁcantly affect freshwater resources [3,11].

Figure 1. (a) Distribution of coral reefs in the world, (b) location of coral islands in the South China Sea, and land use maps of (c) Zhaoshu Island, (d) Yongxing Island, and (e) Yongshu Island in 2019. Note: The distribution of coral reefs in (a) is made using Natural Earth (free vector and raster map data @ naturalearthdata.com (accessed on 15 December 2020)).

Actual evapotranspiration (ETa), which represents the sum of water evaporation and transpiration from a surface area to the atmosphere, is the least quantified freshwater component of atoll islands’ water budget analysis because of the inherent difficulties in measuring it [1,5]. ETa could be estimated from potential evapotranspiration (ETp) by means of a water balance, taking account of evapotranspiration losses from interception and soil water stores [5], where ETp has been investigated based on the Priestley–Taylor formula, the Penman formula, and pan evaporation [12–14]. The ETa of coral islands could also be estimated by field observations such as eddy covariance, heat-dissipative sap flow, and monitor wells [15–19]. However, these studies only concentrated on a few coral islands, and it was hard to go further in such studies due to the difficulty in data collection and equipment maintenance. Satellite remote sensing data, which consider the vegetation and land use type changes at a long temporal scale [20,21], could be convenient to study hydrological processes on remote islands; however, it has rarely been used in ETa evaluation on coral islands. In this study, the remote-sensing-based Vegetation Interfaces Processes (VIP-RS) model [22–24] was used to evaluate the ETa of coral islands in the South China Sea from 2016 to 2019. The VIP-RS model takes account of the detailed physical mechanism of the exchanges of energy and water between the land surface and atmosphere, assimilated with remote sensing information. We aimed to (1) explore the ETa and water dynamics on natural and reclaimed coral islands in the South China Sea and dominant driving factors Remote Sens. 2021, 13, 1110 3 of 18

for ETa using partial correlation analysis and (2) compare the differences in ETa between natural and reclaimed coral islands and analyze the inﬂuences of ETa on different land use types according to the scenario analysis of land use types. This method considers the spatial changes in vegetation and land use types and reduces the hard work of ﬁeld measurements, helping us evaluate the characteristics of freshwater resources accurately and improve water resource management.

2. Materials and Methods 2.1. Study Area The South China Sea (Figure1b), situated between the Karimata Strait (~3 ◦S) and the middle of the Taiwan Strait (~23.5◦N), has numerous coral islands. These coral islands develop a wide variety of scenery and wildlife and provide convenience to people with fish, tourism, and energy. To study the water cycle in natural and reclaimed coral islands in the South China Sea, three typical islands were selected. Zhaoshu Island (latitude 16◦580N, longitude 112◦160E) and Yongxing Island (latitude 16◦500N, longitude 112◦200E) were chosen as typical natural islands with different island sizes. Yongshu Island (latitude 9◦370N, longitude 112◦580E), reclaimed in 2015, was chosen as a case of reclaimed islands. These three islands could represent most coral islands in the South China Sea in island size and type. Zhaoshu Island (Figure1c) has a natural area of 0.22 km 2, with an area of about 0.07 km2 reclaimed in the west. Lush vegetation is distributed in a ring and band in the natural parts, and the reclaimed area is occupied with greens, a built-up area, and roads. The elevation is approximately 4.4 m, with the core of the island being somewhat low. The soil type is calcareous sand, which is marine sediment originating from coral and other creatures’ debris, with high permeability. It has a tropical monsoon climate that is warm all year, with dry and wet seasons, and the average annual temperature is above 26 °C, while the daily temperature variation in a year is only about 6 °C. The rainfall is mainly controlled by typhoons and convective rains, with an average annual rainfall of about 1500 mm, 85% of which is concentrated in June–November [25]. The water table is between 2.3 and 3.0 m and is influenced by the tides. Yongxing Island (Figure1d) is the largest coral island of the Xisha Archipelago, which has an area of 3.58 km2, with the reclaimed area in the northeast. It is 18 km away from Zhaoshu Island, with approximately 3000 people living there. It has similar climate conditions, morphology, and depositional environments as Zhaoshu Island. However, its altitude is as high as 8.2 m. The vegetation is more diverse than that on Zhaoshu Island, but it has suffered some destruction because of human activities. Yongshu Island (Figure1e) is a reclaimed island that is located in Yongshu Atoll. It was reclaimed from August 2014 to September 2015, with the land area increased from less than 0.2 km2 to 2.23 km2. The length is 3.71 km, and the maximum width is 0.99 km [26]. The topography is flat, with an elevation normally 3 to 4 m above the mean sea level. It was reclaimed stratified, which is similar to the formation process of natural coral islands who were deposited by coral debris under the influence of storms or tidal waves [14,16], but the reclaimed island has accelerated the process. Therefore, the coral sand properties of the reclaimed islands are similar to those of natural calcareous sands, but the vegetation is sparse and needs time to recover under natural and artificial processes [27].

2.2. Model Description The VIP-RS model (Mo et al., 2015) is based on the VIP-distributed ecohydrological dynamic model (Mo and Liu, 2001, Mo et al., 2004). The model inputs include land use types, atmospheric forcing, Digital Elevation Model (DEM), and the Normalized Difference Vegetation Index (NDVI) data series. The model was run from 01 August 2015, to 31 December 2019, at a daily scale, with a spatial resolution of 10 m. The simulated period was chosen for two reasons: (1) the availability of image data from Sentinel-2 is from July 2015, and (2) the reclaimed construction of Yongshu Island was completed in July 2015. Remote Sens. 2021, 13, 1110 4 of 18

In the VIP-RS model, actual ETa consists of vegetation transpiration (Ec), soil evaporation (Es), and canopy interception evaporation (Ei). Ec and Es based on the energy budgets of canopy and soil surface calculated separately with the Penman–Monteith equation. Ec is estimated from potential transpiration (Ecp), mediated by plant physiological characteristics and environmental stresses, and expressed as follows:

Ec = Ecp fw ft (1)

where fw and ft are stress functions of the atmospheric water vapor pressure deﬁcit and air temperature, respectively, and they are calculated using algorithms proposed by Mu et al. [28] and Zhang et al. [29]. Ecp is based on the Penman–Monteith equation:

1 E = ∆R + f ρc D/r /(∆ + γη) (2) cp λ nc c p a

−1 where Rnc is the net radiation absorbed by the canopy (MJ d ), fc is the fractional cover of −1 ◦ −1 −3 vegetation, cp is the specific heat capacity of air (J kg C ), ρ is the air density (kg m ), −1 λ is the latent heat of vaporization of water (J kg ), ra is the aerodynamic resistance between the canopy and the reference height (s m−1), and D is the saturated water vapor pressure deficit of air (hPa). ∆ is the slope of the curve relating saturated water vapor pressure to temperature (hPa ◦C−1), and γ is the psychrometric constant (hPa ◦C−1). η is the ratio of the minimum stomatal resistance (s m−1) of a natural plant functional type to that of a reference crop. It was proposed by Mo et al. (2015) to avoid bias from considering the crop reference transpiration for other vegetation types, as done by most remote sensing evapotranspiration models. The minimum stomatal resistance was adopted from Leuning et al. (2008) and Bastiaanssen et al. (2008) [30,31]. Soil evaporation (Es) is restricted by surface potential evaporation (Esp) and soil moisture exfiltration (Eex), Es = min Esp, Eex (3) 1 E = ∆(R − G) + (1 − f )ρc D/r /(∆ + γ) (4) sp λ ns c p as 0.5 0.5 Eex = S t − (t − 1) (5)

where Rns is the net radiation absorbed by the soil surface, G is the soil heat ﬂux, ras is the aerodynamic resistance between reference height and soil surface, S is the soil-controlled exﬁltration volume and usually falls by 3–5 mm d−0.5 (in this study, a value of 4 mm d−0.5 is set), and t is the time (days) elapsed since the day following rainfall [32]. Net radiation Rn is partitioned into energy for the soil and the canopy (Rns for the soil and Rnc for the canopy): Rnc = fcRn (6)

Rns = Rn − Rnc (7)

The vegetation fraction ( fc) is estimated using an empirical formula based on the remote sensing vegetation index [33]:

β NDVImax − NDVI fc = 1 − ( ) (8) NDVImax − NDVImin

where β is an empirical constant (from 0.6 to 1.2), which was 0.9 in this study; NDVImax represents NDVI values of the full vegetation cover, and NDVImin represents NDVI values of bare soil.

2.3. Data The land use classification was partitioned through artificial interpretation according to field investigation and images of Google Earth. First, we investigated different plant Remote Sens. 2021, 13, 1110 5 of 18

communities around the whole islands and created the land use type classification according to plant growth characteristics and stand structure. The basis of the land use type classification and field investigation numbers is shown in Table1. Second, we confirmed different plant species’ distributions and their land use types and drew the maps of land use types from 2016 to 2019, manually combing images of Google Earth. The map was made at an annual scale for rare changes in land use types within a year. Figure1c–e showed the map of the land use types on Zhaoshu Island, Yongxing Island, and Yongshu Island, respectively, in 2019.

Table 1. Basis of land use types classiﬁcation and ﬁeld investigation numbers.

Land Use Types Basis of Classiﬁcation Sample Numbers Evergreen forest Mainly the community of Cocos nucifera on Yongxing Island 15 Evergreen shrubs Mainly the community of Scaevola sericea on Zhaoshu Island 12 Mainly a mixture of different vegetation species (common vegetation species, such as Calophylluminophyllum, Mixed forests Terminalia catappa, Cocos nucifera, Clerodendrum inerme, 31 Scaevola sericea, Pandanus tectorius, Ceodes grandis, Guettarda speciosa, Canavalia maritima, etc.) Grassland The lawn beside the airport, roads, or built-up area 12 Built-up or hardened roads Construction land or concrete roads 12 Barren or sparsely vegetated area Bare land or sparsely vegetated areas 15

Meteorological data were obtained from the China Meteorological Administration at a daily scale at the Yongxing Island meteorological station (ID: 59981) and the Yongshu Island meteorological station (ID: 59985). Atmospheric forcing variables included air temperature (TA), maximum air temperature (TMX), minimum air temperature (TMN), relative humidity (RH), wind speed (WIN), sunshine duration (SHIN), precipitation (PRP), and atmospheric pressure (AP). The missing SHIN data at the Yongxing Island meteorological station from 2018 to 2019 were replaced by the data at the Qionghai Island meteorological station (ID: 59985; 19◦140N, 110◦280E). SHIN is used to calculate solar radiation (Ra) in the VIP-RS model, so the missing SHIN data at the Yongshu Island meteorological station were replaced by Ra, which was obtained from the modeling results of Hargreaves’ radiation formula [34]. According to the station meteorological data from 2010 to 2019 (Figure2), precipitation and temperatures are similar at both meteorological sites. The annual TA ranges from 25.8 ◦C to 28 ◦C in different years, while the average monthly TA ranges from 24 ◦C to 30 ◦C, with its maximum value in June and minimum value in January. The annual PRP varies from 900 mm to 1900 mm, with a large ﬂuctuation in different years. Similar to TA, the average monthly PRP has seasonal differences, with the maximum value of 270 mm in July and the minimum value of less than 10 mm in February.

Figure 2. (a) Average monthly and (b) annual precipitation (PRP) and temperature (TA) values during 2010–2019 from Yongxing Island and Yongshu Island meteorological stations. Remote Sens. 2021, 13, 1110 6 of 18

The NDVI was extracted using Sentinel-2 images by Google Earth Engine. The data spanned from July 2015 to December 2019, with a monthly temporal composite at 10 m spatial resolution. The temporal changes in the NDVI in Figure3 show that Zhaoshu Island has the largest spatial average NDVI for its large area of vegetation, followed by Yongxing Island and Yongshu Island. The NDVI in these islands exhibits a seasonal change and a slow ﬂuctuated increase tendency during the simulated period. We noticed that a rapid decline in the NDVI occurred between September 2018 and May 2019 in Zhaoshu Island, which may be related to drought and vegetation destruction caused by insufﬁcient precipitation. Figure4 demonstrates that an NDVI above 0.4 has a larger ratio to the entire island area and a lower degree of fragmentation in Zhaoshu Island than in Yongxing Island. The reason for the differences between the two natural islands is the high-intensity human activities in Yongxing Island. The NDVI in reclaimed areas ranges from 0.0 to 0.2, except the grassland beside the airport and scattered green areas ranging between 0.2 and 0.4.

Figure 3. Spatial average monthly Normalized Difference Vegetation Index (NDVI) on Zhaoshu Island, Yongxing Island, and Yongshu Island from July 2015 to December 2019.

Figure 4. Spatial distribution of the average NDVI on Zhaoshu Island, Yongxing Island, and Yongshu Island from July 2015 to December 2019. The altitude was set as a 2 m height for the whole area of islands because the small altitude differences on coral islands have a rare inﬂuence on simulated results.

2.4. Model Validation The model validation could be divided into four parts: (1) The potential evaporation was validated by pan evaporation measured with pan Φ20 at a monthly scale at the Yongxing Island meteorological station from 2008 to 2013 because pan evaporation data at the Yongxing Island meteorological station after 2014 are missing. Potential evapotranspiration was only driven by meteorological data, and the changes in vegetation and land use types have no influences on potential evapotranspiration. The validation of the simulated potential evapotranspiration from 2008 to 2013 confirmed the applicability of the model under tropical coral island meteorological conditions, although the validated periods had different time spans than our study periods. (2) The daily simulated ETa of grasses and Es of bare land in July 2018 and August 2018 were validated by field-collected ETa using a specially constructed micro-lysimeter in July 2018 and August 2018 [35]. The micro-lysimeter was constructed from a polyvinyl chloride (PVC) pipe with a height of 150 mm and a diameter of 101.6 mm, and the sampling processes could refer to Plauborg, (1995) [36]. The Es of bare land was measured from 22 Remote Sens. 2021, 13, 1110 7 of 18

June 2018 to 6 July 2018 during the periods with no precipitation, with the ETa of grasses being measured between 22 July 2018 and 28 July 2018 [35]. (3) The daily simulated Ec of evergreen forests, evergreen shrubs, and mixed forests was validated with the Ec measured by Jie et al. (2015) [37], who measured the leaf Ec of Calophylluminophyllum, Terminalia catappa, Cocos nucifera, Clerodendrum inerme, Scaevola sericea, Pandanus tectorius, and Canavalia maritima with the Li-Cor 6400 Portable Photo- synthesis System on Shanhu Island of the Xisha Archipelago. In this study, evergreen forests were mainly the Cocos nucifera community on Yongxing Island, evergreen shrubs represented the community of Scaevola sericea on Zhaoshu Island, and mixed forests were mixtures of different vegetation species (Table1). We compared the simulated average Ec of evergreen forests, evergreen shrubs, and mixed forests from 31 August 2015 to 2 September 2015 with the Ec of Cocos nucifera and Scaevola sericea and the average Ec of all vegetation species, as measured by Jie et al. (2015) [37] on 1 September 2015. The measured branch Ec at a daily scale was transferred from the leaf Ec using the leaf area index (LAI) [38] and the sinusoidal distributional shape of ETp [38,39]. (4) To validate the model performance on coral islands at a temporal scale, we simulated the monthly ETa and Ec of the Cocos nucifera plantation on Vanuatu in the South Pacific to compare with the works of Roupsard et al. (2007) [19], who monitored the monthly ETa by eddy covariance and the monthly Ec by heat-dissipative sap flow, soil heat flux and, stand heat storage from 2001 to 2004. Vanuatu (15◦260N, 167◦110E) is a coral island having a similar distance to the equator as Yongxing Island (16◦500N, 112◦200E), and the climate of Vanuatu is tropical monsoons with seasonally fluctuating precipitation, which is similar to the climate in the South China Sea [19]. Although the meteorology of Vanuatu has some differences in temporal distribution, the validation of evapotranspiration in Vanuatu could also illustrate the applicability of the VIP-RS model in tropical coral islands. The meteorological data for simulating ETa in Vanuatu were obtained from the website of the National Oceanic and Atmospheric Administration (ID: 915510-99999/915540-99999), and the LAI was set as 3.0, as given by Roupsard et al. (2007) [19].

2.5. Dominant Driving Factor Analysis Dominant driving factor analysis can explore the main driving factors that control ETa. The analysis process includes simple linear regression, the first difference de-trending method, normalization, and partial correlation analysis, which are all performed in each grid over the whole region. The dominant driving factor is set as the factor that has the highest partial correlation coefficient. The specific methods used in this study are as follows: First, we used the first-difference de-trending method (i.e., the difference in values from one month to the previous month) to get rid of the trend influences on evapotranspiration [40,41]. Second, we normalized each variable to facilitate the comparison between variables of different units and orders of magnitude:

xi − xmin Xs = (9) xmax − xmin Third, since climatic factors have correlations with each other, the relationships between the de-trending climatic factors and the de-trending ETa are quantiﬁed by partial correlation analysis. Partial correlation analysis explores the relationship between two variables independent of the inﬂuences of other factors [40,42], which is shown as follows:

Ys = a1X1s + a2X2s + a3X3s + ······ (10)

The ETa of the coral islands was simulated from 2016 to 2019 with a 10 m spatial resolution and a daily time step. Partial correlation analysis was used to investigate the relationships between ETa, the NDVI, and climatic factors at the monthly scale. The Remote Sens. 2021, 13, 1110 8 of 18

factor that had the highest partial correlation coefﬁcient was identiﬁed as the dominant climatic factor.

3. Results 3.1. Evaluation of the Model

The simulated water ETp and pan evaporation (pan coefficient of 0.7) agreed well as they demonstrated similar inter-annual and seasonal variations with a maximum value of about 170 mm/month in July and a minimum value of about 80 mm/month in January (Figure5a). The comparison (Figure5b) between pan evaporation and simulated ETp exhibited a determination coefficient of 0.82 within the significance of 0.01, with the root- mean-square error (RMSE) being 10.54. The results indicated that the model is sensitive to meteorological conditions and can simulate the potential evaporation well.

Figure 5. (a) Temporal pattern and (b) comparison of simulated potential evaporation and pan evaporation at the Yongxing Island meteorological station.

Figure6 exhibits the evapotranspiration comparison of different land use types. The average daily Ec of evergreen forests, evergreen shrubs, and mixed forests was about 1.9 mm/day, 2.6 mm/day, and 1.51 mm/day, respectively. The average daily ETa of grassland was 2.1 mm/day, with the Es of bare land being approximately 0.9 mm/day. Evergreen shrubs showed the largest Ec because the dominant species of Scaevola sericea on Zhaoshu Island had a high transpiration capacity, as shown by Jie et al. (2015) [37], and its distribution is on a natural area with almost no human activities. Error bars showed that simulated values had a smaller ﬂuctuation than measured values, the reason is that the measured evapotranspiration may show some differences at the inner-daily scale, versus the simulated values showed stable values. The comparison results demonstrated that the model can simulate evapotranspiration of different land use types well.

Figure 6. Comparisons of simulated and measured evapotranspiration of ﬁve types of land use (BS(Es): evaporation of bare land; EF(Ec): transpiration of evergreen forests; ES(Ec): transpiration of evergreen shrubs; GL(ETa): evapotranspiration of grassland; MF(Ec): transpiration of mixed forests). Error bars are standard deviations.

To validate the performance of the model at the temporal scale, we compared ETa and Ec between simulated and measured values at the Cocos nucifera plantation on Vanuatu Remote Sens. 2021, 13, 1110 9 of 18

from 2001 to 2004 (Figure7). The measured ETa of Cocos nucifera had an average value of 2.6 mm/day, with the highest value of approximately 3.4 mm/day in December and the lowest value of 1.8 mm/day in June. The measured Ec of Cocos nucifera had an average value of 1.7 mm/day, with the largest value of approximately 2.3 mm/day in November versus the lowest value of 1.3 mm/day in June. The simulated ETa and Ec showed a similar ﬂuctuation with measured values in different months, with the relative errors of ETa and Ec being less than 0.1 mm/day. The results demonstrated that the model can simulate ETa and Ec at a monthly scale well.

Figure 7. (a) Comparison of simulated and measured monthly average actual evapotranspiration (ETa) and (b) actual transpiration (Ec) of Cocos nucifera on Vanuatu. Error bars are standard deviations. 3.2. Evapotranspiration Characteristics on Different Islands

Figure8 shows the spatial patterns of the average annual ETa from 2016 to 2019, which was similar to the distribution of the NDVI. ETa in the natural area showed higher values than that in reclaimed areas. Zhaoshu Island had the largest average annual ETa of 685 mm for its largest proportion of the natural vegetation area. The average annual ETa of Yongxing Island was 530 mm. Although Yongxing Island also had a large area of natural vegetation, the vegetation was somewhat destroyed because of the inhabitants’ activities. Yongshu Island was a reclaimed island that had the lowest ETa of 210 mm because of the high proportion of built-up or hardened roads. The average annual ETa of different land use types is shown in the bar ﬁgures in Figure8 . We found that evergreen forests on Yongxing Island and evergreen shrubs on Zhaoshu Island had the largest average annual ETa, ranging from 1000 mm/year to 1200 mm/year. The ETa of grassland and mixed forests ranged from 600 mm/year to 800 mm/year versus the ETa of bare land or sparsely vegetated areas, which ranging from 400 mm/year to 600 mm/year on Yongxing Island and Zhaoshu Island. The ETa of grassland, mixed forests, and bare land or sparsely vegetated land ranged from 300 mm/year to 600 mm/year on Yongshu Island, which was smaller than that on Yongxing Island and Zhaoshu Island. The differences may be caused by the somewhat small vegetation cover in the reclaimed island. The ETa of bare land or sparsely vegetated land was similar to the ETa of grassland and mixed forests for high values of Es in bare land. The ETa of built-up or hardened roads was less than 50 mm/year on all three islands. Remote Sens. 2021, 13, 1110 10 of 18

Figure 8. Spatial patterns of annual average evapotranspiration and evapotranspiration of different land use types on (a) Zhaoshu Island, (b) Yongxing Island, and (c) Yongshu Island from 2016 to 2019. EF: evergreen forests; BL: built-up or hardened roads; ES: evergreen shrubs; MF: mixed forests; GL: grassland; BS: barren or sparse vegetation.

Figure9 shows the simulated daily ETa of different land use types on the three islands. The daily ETa exhibited seasonal fluctuations due to the conditions of energies and water. ETa was high from July to September, while it was low from January to March, which was caused by the change in meteorological conditions (Figure2). At a daily scale, ETa rose rapidly after rainfall, especially for the bare land or sparsely vegetated areas, due to increases in Es and Ei. After a period of no rainfall, the ETa of bare land or sparsely vegetated areas dropped to approximately 1 mm/day, while evergreen forests and evergreen shrubs had a value of at least 2 mm/day, even in the dry season, and the ETa of grassland or mixed forests fell in between. The ETa of built-up or hardened roads exhibited small values during a few days after rainfall. ETa was influenced by PRP on the reclaimed island more than that on the natural island because Es and Ei are sensitive in a short time, which also indicated that the water deficit should be considered for ETa of the reclaimed island. We partially enlarged the detail of ETa from 20 March 2019 to 30 April 2019 when these three islands experienced two periods of rainfall after more than 30 days of no rainfall. One rainfall level is 37 mm on 31 March 2019, and the other is 13 mm on 16 April 2019. When it rained, we could see that ETa rose rapidly for all land use types. Es and Ei displayed an obvious increase with rainfall; however, Ec somewhat decreased because of stomatal closure. After rainfall, the ETa of the land use types with large vegetation cover had a high value, and that of bare land or sparsely vegetated areas dropped to a low level rapidly.

3.3. Dominant Driving Factors for ETa

ETa is mainly influenced by meteorological factors and NDVI conditions [43–45]. Gao et al. [46] found that PRP controlled ETa changes in most parts of Chinese areas, especially in arid and semi-arid regions. Coral islands in the South China Sea have a tropical oceanic monsoon climate, where vegetation is somewhat similar to the tropical dry evergreen forests in Southeast Asia [47]. Pour et al. [48] demonstrated that the minimum temperature (31.5–48.2%) is the most influencing factor, followed by wind speed (15.1– 32.8%), in defining ETp in peninsular Malaysia. To identify the dominant driving factors for ETa of the coral islands in the South China Sea, partial correlation analysis was used based on the detrending data series for each grid. Remote Sens. 2021, 13, 1110 11 of 18

Figure 9. Evapotranspiration of different land use types at a daily scale on Zhaoshu Island, Yongxing Island, and Yong- shu Island.

Table2 shows the spatial distribution of dominant driving factors within the significance of 0.01. The NDVI was the main influencing factor for ETa of Zhaoshu Island and Yongxing Island. On Yongshu Island, RH controlled the ETa of mixed forests and bare land or sparsely vegetated areas, while the NDVI controlled the ETa of grassland. The results revealed that the NDVI is the dominant driving factor in the natural islands. However, for reclaimed islands, RH was another dominant driving factor, except for the NDVI. The reason for the difference is that RH controlled the diffusion of soil water to air in bare land and sparsely vegetated areas. Salvucci and Gentine [49] also demonstrated that ETa is controlled by the gradient of humidity near the land surface, which reflects a balance of radiative heating and cooling, conversion to latent heat, canopy and soil heating, and the turbulent transport of sensible heat into the atmosphere.

Table 2. Dominant impact variables of evapotranspiration of different land use types.

Land Use Types Yongxing Island Yongshu Island Zhaoshu Island EF/ES NDVI - NDVI MF NDVI RH NDVI GL NDVI NDVI - BS NDVI RH NDVI Note: EF: evergreen forests; ES: evergreen shrubs; MF: mixed forests; GL: grassland; BS: barren or sparse vegetation.

3.4. Evapotranspiration Distribution on the Three Islands

PRP, the only source of fresh water on the coral islands, is mainly consumed by ETa, inﬁltration, and rainwater collection or water ﬂowing into the sea. Storage changes within Remote Sens. 2021, 13, 1110 12 of 18

the unsaturated zone and surface runoff can be neglected [1,5] on coral islands due to the low water-binding capacity of gravelly and sandy soil texture, both for natural and for reclaimed islands [26,50]. In this study, the water budget analysis could be written as ∆V = PRP − Eta, where PRP is rainfall (mm); ETa is the actual evapotranspiration (mm), which contains Ec, Es, and Ei; and ∆V is the total amount of infiltration, rainwater collection, and water flowing into the sea. Figure 10 shows the annual water balance from 2016 to 2019; the average PRP was about 1443 mm/year on the three islands. The average annual Ec on Zhaoshu Island and Yongxing Island was 319 mm/year and 246 mm/year, respectively, while the Es on Zhaoshu Island and Yongxing Island was 293 mm/year and 306 mm/year, respectively. The Ec on Yongshu Island was 48 mm/year, which was far smaller than the Es of 174 mm/year. The annual Ei on the three islands ranged from 20 mm to 70 mm, which is a low value. Figure 11 shows the ratio of Ec, Es, Ei, and ∆V to precipitation on Zhaoshu Island, Yongxing Island, and Yongshu Island from 2016 to 2019. Ec and Es occupied more than 20% of PRP on Zhaoshu Island and Yongxing Island. However, the ratio of Ec and Es to PRP was 3% and 13%, respectively. PRP had a somewhat large fluctuation in different years, and ETa had small fluctuations compared to PRP, indicating that ETa is restricted more by vegetation or energy than by PRP. The results were consistent with the analysis of Table2.

Figure 10. Water distribution analysis at the annual scale from 2016 to 2019 on Zhaoshu Island, Yongxing Island, and Yongshu Island.

Figure 11. Average annual ratio of Ec, Es, Ei, and ∆V to precipitation on Zhaoshu Island, Yongxing Island, and Yongshu Island from 2016 to 2019.

The average monthly water balance analysis demonstrated a seasonal distribution over the simulated period (Figure 12). The wet season had PRP between 120 mm/month and 220 mm/month from June to November on Zhaoshu Island and Yongxing Island. In comparison, the dry season had PRP between 2 mm/month and 90 mm/month from December to May. Yongshu Island had a similar seasonal distribution, but PRP had a larger ﬂuctuation than that on Zhaoshu Island and Yongxing Island. In the dry season, PRP on Zhaoshu Island and Yongxing Island was 10 mm/month to 30 mm/month less than ETa. Vegetation had to absorb water from groundwater to sustain growth, which has been found in many studies on coral islands [4,5]. However, on the reclaimed island, the groundwater Remote Sens. 2021, 13, 1110 13 of 18

could be recharged by PRP through the entire year. Comparing the water balance at the annual scale, we found that a relatively large fluctuation of PRP mainly occurs in the wet season. In the wet season, PRP was sufficient for ETa. In addition, ETa could be limited by PPR in the dry season. However, the changes in ETa due to fluctuations of rainfall in the dry season were small for the low values of ETa. The results further proved that ETa is controlled more by energy and the NDVI than by PRP.

Figure 12. Water distribution analysis at the monthly scale from 2016 to 2019 on Zhaoshu Island, Yongxing Island, and Yongshu Island.

4. Discussion 4.1. Evapotranspiration Differences in Natural and Reclaimed Areas of Coral Islands Reclaimed coral islands have somewhat large areas of hardening and barren land or sparsely vegetated areas compared to natural islands, which indicates that reclaimed islands consume small volumes of water due to the low Ec caused by the lack of evergreen forests or shrubs. Figure 13 compares the differences in ETa between natural areas and reclaimed areas. ETa had similar seasonal fluctuations at temporal series on different islands. However, ETa in natural areas was 25–50 mm/month higher than that in reclaimed areas. Natural and reclaimed coral islands face different water resource problems. The most important freshwater problem on natural islands is water contamination and sustainable exploitation of water resources [1]. However, for newly reclaimed coral islands that have no fresh groundwater, it will take 20 years to form a stable freshwater lens [26]. The most important thing is to evaluate the spatial and temporal distribution of the fresh water of coral islands and the available water resources [51]. Falkland and Woodroffe [4] predicted an approximately 50% increase in groundwater recharge if all coconut trees were removed from Tarawa Atoll and Kiritimati. However, it could bring new problems if we simply decrease the vegetation cover. For the people on coral islands, however, a low vegetation cover could increase the groundwater recharge, which have a positive influence on the freshwater storage [52,53]. It could also provide shelters to seabirds, which could promote soil-forming processes, and decrease the surface temperature to make the islands livable [54–56]. However, the proper way to manage the limited fresh water and how the different land use types influence the water dynamics on the reclaimed islands are still not clear. Remote Sens. 2021, 13, 1110 14 of 18

Figure 13. Comparison of evapotranspiration in natural and reclaimed areas of Zhaoshu Island, Yongxing Island, and Yongshu Island from 2016 to 2019.

4.2. Influences of Land Use Types on Water Dynamics in Reclaimed Islands To study the influences of land use types on water dynamics in reclaimed coral islands, we attempted to set different land uses on Yongshu Island to study its possible water transport process. The land use types on Yongshu Island in 2019 were set as a reference. The ratio of trees and shrubs, grassland, barren or sparsely vegetated areas, and built-up or hardened roads to the whole island area was 0.08, 0.19, 0.22, and 0.51, respectively, with Ec, Es, Ei, and ∆V in 2019 of 98 mm, 161 mm, 24 mm, and 1160 mm, respectively. The water distribution changes are shown in Table3 when we covered one of these land use types over the entire island. If we cover the whole islands with trees and shrubs, Ec, Es, and Ei increase to 520 mm, 253 mm, and 118 mm, respectively, with their values increasing from 57% to 430%. In such a setting, ETa is larger than ∆V obviously and vegetation consumes a relatively large proportion of precipitation. If we cover the entire coral island with grassland, Ec, Es, and Ei increase to 193 mm, 307 mm, and 46 mm, respectively, with their values increasing from 89% to 93%, and ∆V decreases to 897 mm. In this situation, Es becomes the largest ETa component, and ∆V is the largest consumption of precipitation. If we cover the entire coral island with barren or sparsely vegetated areas, the overall ETa still exhibits a 57% increase, although the area of vegetation decreases. The reason for the increase in ETa can be attributed to the reduction in the area of built-up or hardened roads. The results showed that ETa of Yongshu Island in 2019 was still at a low value due to the restriction of nearly 51% built-up or hardened roads. Built-up or hardened roads influence the water cycle in another way. They cut off the interchange of water between air and the land surface, which can benefit rainwater collection. For reclaimed coral islands, it will take a long time to form fresh groundwater, so a rainwater collection system is an important measure to collect the limited freshwater for the people living there [57]. The water of the rainwater collection system mainly comes from built-up or hardened roads. Setting Yongshu Island in 2019 as a reference, if 50% of the precipitation on built-up or hardened roads were collected by the rainwater collection system, the rainwater collection system could collect 332 mm of fresh water in a year, which is about 0.74 million squares of water. The effect of built-up or hardened roads on the rainwater collection system needs further research for reclaimed islands.

Table 3. Inﬂuence of changes in land use types on water dynamics of coral islands.

Land Use Types Ec (mm)/Ratio Es (mm)/Ratio Ei (mm)/Ratio ∆V (mm)/Ratio Trees and shrubs 520/430% 253/57% 118/389% 552/−50% Grassland 193/97% 307/90% 46/89% 897/−23% Barren or sparsely 89/−9% 319/98% 27/13% 1007/−13% vegetated areas Built-up or hardened roads 0/−100% 24/−85% 0/−100% 1419/22% Note: The ratio means the relative value of the difference between the evapotranspiration of different land use settings and evapotranspiration in 2019 to evapotranspiration in 2019. Remote Sens. 2021, 13, 1110 15 of 18

4.3. Implications and Limitations

To evaluate ETa of coral islands conveniently, free and relatively high-resolution satellite products are the first reference. We compared the NDVI products that extracted Landsat series images (30 m spatial resolution) and Sentinel-2 images (10 m spatial resolution). However, the products of Landsat-5 (1985 to 2011) had a wrong geospatial calibration, which causes deviations in the coordinates of images at different time, sand products of Landsat-7 (1999 to now) had some Scan Line Corrector failure. Landsat-8 (2014 to now) can simulate ETa well, but it has a similar on-orbit time with Sentinel-2 (2015 to now), which has a high temporal-spatial resolution. Another limitation is the shortage of data in the wet season due to the influence of cloud cover. In this study, we use 31 days of temporal-scale data to run the model, although the products can provide a sight image of approximately a week. Land use types were obtained using artificial interpretation methods because of the difficulty in distinguishing different land use types from the satellite remote sensing image due to limits of spatial resolution. We tried to classify land use types using supervised classification and unsupervised classification methods, but the results were unsatisfactory, especially for differences in build-ups and bare lands in reclaimed areas, so an effective algorithm should be developed for coral islands. The simulated results were validated with data of field samplings and other literature; however, uncertainties still exist due to the limited hydrological observations and intrinsic uncertainties of the hydrological model. Some plants close pores to defend against drought on coral islands in the South China Sea [37], which is similar to the study of Pitman [47], who found that the stomatal conductance of planta may reduce by 50% when the soil water potential reaches −1.2 (±0.07) MPa. Stomatal conductance is an important parameter during simulated processes. However, the model does not differentiate stomatal conductance changes between the dry season and the wet season, which may lead to some errors for some vegetation species. So, it is necessary to adjust sensitive parameters for coral islands referring to dominant driving factors. Even so, the evaluation of ETa based on the remote-sensing-based model makes up for the difficulty in fieldwork, which improves the certainty and accuracy at a spatial scale. In addition, it gives us a new reference to protect and manage scarce freshwater resources properly.

5. Conclusions

This study evaluated ETa and water dynamics on natural and reclaimed coral islands in the South China Sea from 2016 to 2019 based on a simplified remote-sensing-based VIP model. The simulated results showed that the water dynamics on coral islands have a seasonal distribution. The average annual PRP was about 1443 mm/year on the three islands, with 1088 mm in the wet season and 355 mm in the dry season. The average annual ETa of Zhaoshu Island, Yongxing Island, and Yongshu Island was 685 mm, 530 mm, 210 mm, respectively. Ec and Es had a similar value, which ranged from 246 mm to 319 mm (approximately 20% of PRP) on Zhaoshu Island and Yongxing Island. However, Es controls water consumption on Yongshu Island. ETa is restricted more by vegetation or energy than by PRP, which is consistent with the partial correlation analysis result that the NDVI is the dominant factor for the whole area of the natural islands and grassland of the reclaimed island, RH is the dominant driving factor for mixed forests and barren land or sparsely vegetated areas on Yongshu Island. Approximately half of PRP is consumed by evapotranspiration on the natural islands, while approximately 20% of PRP can be consumed by evapotranspiration on the reclaimed island. Under the background of global climate changes, more and more islands are constructed to defend against the possible impact of extreme climate events. Through setting different land uses on the entire islands, we quantitatively evaluated the effects of vegetation and built-up or hardened roads to water cycles on coral islands. Although further studies need to be carried out using the VIP-RS model to optimize parameters due to the specific condition of coral islands, the evaluation of ETa based on a remote-sensing-based model makes up for the difficulty in fieldwork, improving the certainty and accuracy Remote Sens. 2021, 13, 1110 16 of 18

at a spatial scale. In addition, it gives us a new reference to protect and manage scarce freshwater resources properly.

Author Contributions: Conceptualization, S.L., X.M., and S.H. (Shengsheng Han); methodology, X.M. and S.H. (Shi Hu); software, X.M.; validation, X.M. and S.H. (Shengsheng Han); writing—original draft preparation, S.H. (Shengsheng Han); writing—review and editing, S.L., S.H. (Shengsheng Han), X.M., and X.S.; visualization, S.H. (Shengsheng Han) and S.L.; supervision, S.L. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDA13010303). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The data presented in this study are available on request from the corresponding author. The data are not publicly available due to restrictions of privacy. Acknowledgments: We thank Meng Bai, Xuejuan Chen, and Yi Liu in Institute of Geographic Sciences and Natural Resources Research, CAS for helping the first author with model simulation and analysis during preliminary article writing. We also thank all the team members of the project for pertinent comments in multiple group meetings. Conflicts of Interest: The authors declare no conflict of interest.

References 1. Werner, A.D.; Sharp, H.K.; Galvis, S.C.; Post, V.E.A.; Sinclair, P. Hydrogeology and management of freshwater lenses on atoll islands: Review of current knowledge and research needs. J. Hydrol. 2017, 551, 819–844. [CrossRef] 2. Storlazzi, C.D.; Gingerich, S.B.; van Dongeren, A.; Cheriton, O.M.; Swarzenski, P.W.; Quataert, E.; Voss, C.I.; Field, D.W.; Annamalai, H.; Piniak, G.A.; et al. Most atolls will be uninhabitable by the mid-21st century because of sea-level rise exacerbating wave-driven flooding. Sci. Adv. 2018, 4, eaap9741. [CrossRef] 3. White, I.; Falkland, T. Management of freshwater lenses on small Pacific islands. Hydrogeol. J. 2010, 18, 227–246. [CrossRef] 4. Falkland, A.; Woodroffe, C. Geology and hydrogeology of Tarawa and Christmas Island, Kiribati. In Developments in Sedimentology; Elsevier: Amsterdam, The Netherlands, 2004; Volume 54, pp. 577–610. 5. Falkland, A.; Custodio, E. Hydrology and Water Resources of Small Islands: A Practical Guide; Unesco: Paris, France, 1991. 6. Robins, N. A review of small island hydrogeology: Progress (and setbacks) during the recent past. Q. J. Eng. Geol. Hydrogeol. 2013, 46, 157–165. [CrossRef] 7. White, I.; Falkland, T.; Scott, D. Droughts in Small Coral Islands: Case Study, South Tarawa, Kiribati; Unesco: Paris, France, 1999. 8. Van der Velde, M.; Green, S.R.; Vanclooster, M.; Clothier, B.E. Sustainable development in small island developing states: Agricultural intensification, economic development, and freshwater resources management on the coral atoll of Tongatapu. Ecol. Econ. 2007, 61, 456–468. [CrossRef] 9. Post, V.E.A.; Bosserelle, A.L.; Galvis, S.C.; Sinclair, P.J.; Werner, A.D. On the resilience of small-island freshwater lenses: Evidence of the long-term impacts of groundwater abstraction on Bonriki Island, Kiribati. J. Hydrol. 2018, 564, 133–148. [CrossRef] 10. Swanepoel, C.R. Effect of Coral Reefs on Wave Attenuation and Erosion: Mnemba Island, Zanzibar; University of Kwazulu-Natal: Durban, South Africa, 2017. 11. Hejazian, M.; Gurdak, J.J.; Swarzenski, P.; Odigie, K.O.; Storlazzi, C.D. Land-use change and managed aquifer recharge effects on the hydrogeochemistry of two contrasting atoll island aquifers, Roi-Namur Island, Republic of the Marshall Islands. Appl. Geochem. 2017, 80, 58–71. [CrossRef] 12. Nullet, D. Water balance of pacific atolls 1. JAWRA J. Am. Water Resour. Assoc. 1987, 23, 1125–1132. [CrossRef] 13. Jocson, J.M.U.; Jenson, J.W.; Contractor, D.N. Recharge and aquifer response: Northern Guam Lens Aquifer, Guam, Mariana Islands. J. Hydrol. 2002, 260, 231–254. [CrossRef] 14. Woodroffe, C.D.; Falkland, A.C. Geology and hydrogeolgy of the cocos(Keeling) islands. Dev. Sedimentol. 1997.[CrossRef] 15. Falkland, A. Climate, Hydrology, and Water Resources of the Cocos (Keeling) Islands; Atoll Research Bulletin: Washington, DC, USA, 1994. [CrossRef] 16. Falkland, A.C. Hydrology and Water Management on Small Tropical Islands. Hydrol. Warm Humid Reg. 1993, 216, 263–303. 17. Krauss, K.W.; Duberstein, J.A.; Cormier, N.; Young, H.S.; Hathaway, S.A. Proximity to encroaching coconut palm limits native forest water use and persistence on a Pacific atoll. Ecohydrology 2015, 8, 1514–1524. [CrossRef] 18. Comte, J.C.; Join, J.L.; Banton, O.; Nicolini, E. Modelling the response of fresh groundwater to climate and vegetation changes in coral islands. Hydrogeol. J. 2014, 22, 1905–1920. [CrossRef] Remote Sens. 2021, 13, 1110 17 of 18

19. Roupsard, O.; Bonnefond, J.M.; Irvine, M.; Berbigier, P.; Nouvellon, Y.; Dauzat, J.; Taga, S.; Hamel, O.; Jourdan, C.; Saint-Andre, L.; et al. Partitioning energy and evapo-transpiration above and below a tropical palm canopy. Agric. For. Meteorol. 2006, 139, 252–268. [CrossRef] 20. Veettil, B.K.; Van, D.D.; Quang, N.X. Biodiversity and nature conservation in island ecosystems: Spatiotemporal changes in Socotra archipelago (Yemen). J. Coast. Conserv. 2020, 24, 1–9. [CrossRef] 21. Somching, N.; Wongsai, S.; Wongsai, N.; Koedsin, W. Using machine learning algorithm and landsat time series to identify establishment year of para rubber plantations: A case study in Thalang district, Phuket Island, Thailand. Int. J. Remote Sens. 2020, 41, 9075–9100. [CrossRef] 22. Mo, X.; Liu, S.; Lin, Z.; Wang, S.; Hu, S. Trends in land surface evapotranspiration across China with remotely sensed NDVI and climatological data for 1981–2010. Hydrol. Sci. J. 2015, 60, 2163–2177. [CrossRef] 23. Mo, X.G.; Liu, S.X. Simulating evapotranspiration and photosynthesis of winter wheat over the growing season. Agric. For. Meteorol. 2001, 109, 203–222. [CrossRef] 24. Mo, X.G.; Beven, K. Multi-objective parameter conditioning of a three-source wheat canopy model. Agric. For. Meteorol. 2004, 122, 39–63. [CrossRef] 25. Zuo, X.; Su, F.; Wu, W.; Chen, Z.; Shi, W. Spatial and temporal variability of thermal stress to China’s coral reefs in South China Sea. Chin. Geogr. Sci. 2015, 25, 159–173. [CrossRef] 26. Sheng, C.; Han, D.M.; Xu, H.H.; Li, F.C.; Zhang, Y.F.; Shen, Y.Q. Evaluating dynamic mechanisms and formation process of freshwater lenses on reclaimed atoll islands in the South China Sea. J. Hydrol. 2020, 584, 124641. [CrossRef] 27. Zhu, C.Q.; Qin, Y.; Meng, Q.S.; Wang, X.Z.; Wang, R. Formation and sedimentary evolution characteristics of Yongshu Atoll in the South China Sea Islands. Ocean Eng. 2014, 84, 61–66. [CrossRef] 28. Mu, Q.; Heinsch, F.A.; Zhao, M.; Running, S.W. Development of a global evapotranspiration algorithm based on MODIS and global meteorology data. Remote Sens. Environ. 2007, 111, 519–536. [CrossRef] 29. Zhang, K.; Kimball, J.S.; Nemani, R.R.; Running, S.W. A continuous satellite-derived global record of land surface evapotranspiration from 1983 to 2006. Water Resour. Res. 2010, 46.[CrossRef] 30. Leuning, R.; Zhang, Y.Q.; Rajaud, A.; Cleugh, H.; Tu, K. A simple surface conductance model to estimate regional evaporation using MODIS leaf area index and the Penman-Monteith equation. Water Resour. Res. 2008, 44.[CrossRef] 31. Bastiaanssen, W.G.M.; Cheema, M.J.M.; Immerzeel, W.W.; Miltenburg, I.J.; Pelgrum, H. Surface energy balance and actual evapotranspiration of the transboundary Indus Basin estimated from satellite measurements and the ETLook model. Water Resour. Res. 2012, 48.[CrossRef] 32. Choudhury, B.J.; DiGirolamo, N.E. A biophysical process-based estimate of global land surface evaporation using satellite and ancillary data—I. Model description and comparison with observations. J. Hydrol. 1998, 205, 164–185. [CrossRef] 33. Li, F.; Kustas, W.P.; Prueger, J.H.; Neale, C.M.; Jackson, T.J. Utility of remote sensing–based two-source energy balance model under low-and high-vegetation cover conditions. J. Hydrometeorol. 2005, 6, 878–891. [CrossRef] 34. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. 35. Han, S.; Liu, S.; Song, X.; Mo, X.; Yang, L. Field evapotranspiration experiment in Zhaoshu Island of Xisha Islands, South China sea. Geogr. Res. 2021, 40, 174–186, (in Chinese with abstract in English). [CrossRef] 36. Plauborg, F. Evaporation from bare soil in a temperate humid climate—Measurement using micro-lysimeters and time domain reﬂectometry. Agric. For. Meteorol. 1995, 76, 1–17. [CrossRef] 37. Li, J.; Liu, N.; Ren, H.; Shen, W.; Jian, S. Ecological Adaptability of Seven Plant Species to Tropical Coral Island Habitat. Ecol. Environ. Sci. 2016, 25, 790–794, (in Chinese with abstract in English). [CrossRef] 38. Šim ˚unek,J.; Šejna, M.; Saito, H.; Sakai, M.; Genuchten, M.T.V. The Hydrus-1D Software Package for Simulating the Movement of Water, Heat, and Multiple Solutes in Variably Saturated Media; Version 4.17; HYDRUS Software Series 3; Department of Environmental Sciences, University of California Riverside: Riverside, CA, USA, 2013; p. 315. 39. Fayer, M.J. UNSAT-H Version 3.0: Unsaturated Soil Water and Heat Flow Model, Theory, User Manual, and Examples; Report PNNL- 13249; Paciﬁc Northwest National Laboratory: Richland, WA, USA, 2000. 40. Nicholls, N. Increased Australian wheat yield due to recent climate trends. Nature 1997, 387, 484–485. [CrossRef] 41. Lobell, D.B.; Field, C.B. Global scale climate–crop yield relationships and the impacts of recent warming. Environ. Res. Lett. 2007, 2.[CrossRef] 42. Bai, M.; Mo, X.; Liu, S.; Hu, S. Contributions of climate change and vegetation greening to evapotranspiration trend in a typical hilly-gully basin on the Loess Plateau, China. Sci. Total. Environ. 2019, 657, 325–339. [CrossRef] 43. Pei, T.T.; Wu, X.C.; Li, X.Y.; Zhang, Y.; Shi, F.Z.; Ma, Y.J.; Wang, P.; Zhang, C.C. Seasonal divergence in the sensitivity of evapotranspiration to climate and vegetation growth in the Yellow River Basin, China. J. Geophys. Res. Biogeosci. 2017, 122, 103–118. [CrossRef] 44. She, D.; Xia, J.; Zhang, Y. Changes in reference evapotranspiration and its driving factors in the middle reaches of Yellow River Basin, China. Sci. Total Environ. 2017, 607–608, 1151–1162. [CrossRef] 45. Chen, X.J.; Mo, X.G.; Hu, S.; Liu, S.X. Contributions of climate change and human activities to ET and GPP trends over North China Plain from 2000 to 2014. J. Geogr. Sci. 2017, 27, 661–680. [CrossRef] Remote Sens. 2021, 13, 1110 18 of 18

46. Gao, G.; Chen, D.; Xu, C.-y.; Simelton, E. Trend of estimated actual evapotranspiration over China during 1960–2002. J. Geophys. Res. 2007, 112.[CrossRef] 47. Pitman, J.I. Ecophysiology of tropical dry evergreen forest, Thailand: Measured and modelled stomatal conductance of Hopea ferrea, a dominant canopy emergent. J. Appl. Ecol. 1996, 33, 1366–1378. [CrossRef] 48. Pour, S.H.; Abd Wahab, A.K.; Shahid, S.; Bin Ismail, Z. Changes in reference evapotranspiration and its driving factors in peninsular Malaysia. Atmos. Res. 2020, 246.[CrossRef] 49. Salvucci, G.D.; Gentine, P. Emergent relation between surface vapor conductance and relative humidity profiles yields evaporation rates from weather data. Proc. Natl. Acad. Sci. USA 2013, 110, 6287–6291. [CrossRef] 50. Bailey, R.T.; Jenson, J.W.; Olsen, A.E. Estimating the Ground Water Resources of Atoll Islands. Water 2010, 2, 1. [CrossRef] 51. Han, D.; Cao, G.; Song, X. Formation processes and influencing factors of freshwater lens in artificial island of coral reef in South China Sea. Acta Geogr. Sin. 2020, 75, 1053–1064, (in Chinese with abstract in English). [CrossRef] 52. Wallace, C.D.; Bailey, R.T.; Arabi, M. Rainwater catchment system design using simulated future climate data. J. Hydrol. 2015, 529, 1798–1809. [CrossRef] 53. White, I.; Falkland, T.; Perez, P.; Dray, A.; Metuter, T.; Metai, E.; Overmars, M. Challenges in freshwater management in low coral atolls. J. Clean. Prod. 2007, 15, 1522–1528. [CrossRef] 54. Morrison, R.J.; Woodroffe, C.D. The Soils of Kiritimati (Christmas) island, Kiribati, central Pacific: New information and comparison with previous studies. Report 2009, 63, 397–411. [CrossRef] 55. Wu, L.; Liu, X.; Xu, L. Change of Organic δ13C in Ornithogenic Sediments of the Xisha Archipelago, South China Sea and its Implication for Ecological Development. Acta Geol. Sin. Engl. Ed. 2017, 91, 1109–1119. [CrossRef] 56. Graham, N.A.J.; Wilson, S.K.; Carr, P.; Hoey, A.S.; Jennings, S.; MacNeil, M.A. Seabirds enhance coral reef productivity and functioning in the absence of invasive rats. Nature 2018, 559, 250–253. [CrossRef][PubMed] 57. Wallace, C.D. Atoll Island Freshwater Resources: Modeling, Analysis, and Optimization; Colorado State University: Fort Collins, OR, USA, 2015.