基于SAR的拉合尔地面沉降监测与沉降因素评价（英文）

合理开发和管理自然资源对于实现可持续发展至关重要,地面沉降作为目前影响城市可持续发展的最重要地质灾害,其与人类活动有着紧密的关系,如地下水开采,基础设施开发导致的施工超载等。巴基斯坦第二大城市拉合尔出现了显著的地面沉降现象,与此同时,拉合尔也正在经历城市化过程中的人口密度急剧增加、大规模城市建设和地下水开采,通过分析拉合尔地面沉降的不同空间模式与人类活动的对应关系,将有助于合理规划拉合尔城市的发展。本文利用永久散射体合成孔径雷达干涉测量(PSI)技术对拉合尔2018年8月-2019年8月的Sentinel-1数据进行In SAR处理,获取了该城市的地表形变速率为-114-15毫米/年,并结合拉合尔的土地覆盖图、交通网络和水道、土壤类型、人口密度和现场点数据,分析各种因素在地面抬升或沉降发生中的作用。结果显示,沉降与各种参数（如地下水抽取和地下水位降低、土壤类型变化、土地覆盖变化、地表水通道和人口密度）之间建立了很强的相关性,其中拉合尔中部人口密集,地下水开采严重,使其成为该地区沉降最为严重的区域,而在拉维河附近人口较少的农村地区,由于通过运河系统,地下水得到河流补给,观察到了地表抬升。     

Effect of the overconsolidation ratio of soils in surface settlements due to tunneling

Construction of urban tunnels requires the control of surface subsidence to minimize any disturbance to nearby buildings and services. Past study of surface subsidence has been limited to mainly empiri...     

合成孔径雷达干涉测量（InSAR）技术在地表沉降监测中的应用

近年来,合成孔径雷达干涉测量（InSAR）技术应用于地表沉降监测研究领域的方法与实例在国内外工程实践和科研活动中经常出现。近几十年来由于地下水的长期超量开采,华北平原已成为世界上超采地下水最严重、地下水降落漏斗面积最大、地面沉降面积最大、沉降类型最复杂的地区,地表沉降监测周期长、成本高、数据量大,通过In SAR技术的应用可有效提高地表沉降观测的精确性,为预防和减少经济损失发挥积极重要作用。     

用多角度遥感数据反演混合像元组分温度的可行性分析

利用多通道信息反演陆面温度的传统方法,由于通道间信息的高度相关使反演结果的精度难以提高,并且无法得到混合像元中的组分温度。文中以非同温混合像元热红外辐射方向性规律为基础,建立了连续植被类型非同温混合像元热红外辐射模型,用蒙特卡罗方法模拟分析了连续植被的组分有效比辐射率与变量的关系。结果表明:连续植被热辐射亮度值是组分温度T_v、T_s,叶面积指数LAI,叶倾角分布LAD,以及单叶面比辐射率ε_v(θ)和土壤表面比辐射率ε_v(θ)的函数。利用先验知识对变量进行分析后表明,6个变量中只有T_v、T_s、LAI、ε_v(0)需要反演。由于多角度间组分辐射亮度的相关性低,从理论上讲只需要4个角度的辐射亮度观测值就可以解出4个未知量而达到反演组分温度的目的。这4个角度数据除选择垂直方向上辐射亮度数据外,应在30°～50°视角范围内选择另外3个热辐射亮度数据。     

LAND SUBSIDENCE IN THE LOS BANOS–KETTLEMAN CITY AREA,CALIFORNIA: PAST AND FUTURE OCCURRENCE

Land subsidence in the western San Joaquin Valley has been studied and monitored for almost six decades. Especially in the 1950s and 1960s, land subsidence had tremendous negative environmental and socioeconomic implications. The rate of land subsidence in the Los Banos–Kettleman City area, however, has decreased during the last three decades as a result of a decrease in groundwater withdrawal, following the importation of surface water to the Valley in late 1960s. The land subsidence in the Valley potentially could be a serious issue if the existing surface water supplies cannot meet increasing water demands during future critical drought periods. In this respect, this paper proposes a preliminary analysis to predict magnitudes of subsidence for the period 2000 to 2040, based on historical occurrences. In this study, extrapolation from earlier subsidence is based on: (1) knowledge of recent post-drought records at extensometers, (2) assumption of a single time porosity scenario, and (3) the premise that pumpage from the confined aquifer gradually will decrease as a result of progressive water quality degradation. Predicted magnitudes of land subsidence in the Los Banos–Kettleman City area are expected to be in the range 0.5 to 4.0 m for the next 4 to 5 decades under the scenario considered in this paper. [Key words: subsidence, groundwater extraction, aquitard, inelastic compaction, active porosity, Los Banos-Kettleman City, California.]     

东江流域地表水功能区水质对土地利用的响应

以广东省东江流域2013年地表水功能区水质监测数据和Landsat8 OLI遥感影像为数据源,通过遥感信息提取、GIS空间分析和数理统计方法,对不同空间尺度河岸缓冲带土地利用结构与水体中的高锰酸盐指数（CODMn）、氨氮（NH3-N）和总磷（TP）进行相关性分析。结果表明：1）在水功能区中,开发利用区河岸带的城市化程度较高,水质状况相对较差;保护区、保留区、缓冲区河岸带用地类型以耕地和林地为主,水质状况明显优于开发利用区。2）耕地、林地和草地对水质具有正效应,水体和建设用地对水质具有负效应;同时,建设用地对水质影响作用最强,其次是水体和耕地。3）CODMn与耕地、林地、草地、建设用地的相关性远高于NH3-N和TP,3种水质指标与水体相关性都较大。4）耕地和草地与水质在小尺度范围内相关性最为显著;林地和建设用地与水质相关性随河岸缓冲带尺度增大而增大,水体与水质相关性随河岸缓冲带尺度不同呈波动性变化。     

新疆天山北坡平原区地下水可开采量确定方法浅析

在天山北坡以选取的典型区昌吉市和阜康市为例,利用动态相关法、典型流域开采模型法、可开采系数法分析确定地下水可开采量。对这3种方法分别分析了其计算过程中的制约因素,并将计算的结果和以往成果进行了比较。结果表明：动态相关法计算结果准确度较高,成果可靠性高,且在以往的新疆地下水资源评价成果中,至今还未采用过,故推荐动态相关法为新疆天山北坡平原区地下水可开采量确定方法的首选。     

香港现代海平面变化及其影响

潮汐站资料分析表明香港地区海平面每年以１ｍｍ的速度上升。然而，在沿海填海地区，由于受长期地面沉陷的影响，下沉速度超过上升速度，导致海平面明显相对上升。受大规模地面沉陷的影响。填海区在台风暴潮期间极易受海水浸淹。为了有助于未来海岸的开发，建议利用现有技术水平的测量方法，包括卫星测高和激光测距方法等，对填海区进行长期连续的地面沉陷监测。对下沉问题提出预警，以便采取措施，如修筑海提，建抽水站等，以抵御海     

黑河流域土地利用变化对地下水资源的影响

Land use and land cover changes have a great impact on the regional hydrological process. Based on three periods of remote sensing data from the 1960s and the long-term observed data of groundwater from the 1980s, the impacts of land use changes on the groundwater system in the middle reach of Heihe River Basin in recent three decades are analyzed by the perspective of groundwater recharge and discharge system. The results indicate that with the different intensities of land use changes, the impacts on the groundwater recharge were 2.602 × 10^8 m^3/a in the former 15 years （1969-1985） and 0.218 × 10^8 m^3/a in the latter 15 years （1986-2000）, and the impacts on the groundwater discharge were 2.035 × 10^8 m^3/a and 4.91 × 10^8 m^3/a respectively. When the groundwater exploitation was in a reasonable range less than 3.0 × 10^8 m^3/a, the land use changes could control the changes of regional groundwater resources. Influenced by the land use changes and the large-scale exploitation in the recent decade, the groundwater resources present apparently regional differences in Zhangye region. Realizing the impact of land use changes on groundwater system and the characteristics of spatial-temporal variations of regional groundwater resources would be very important for reasonably utilizing and managing water and soil resources.     

基于FEFLOW和GIS技术的矿区地下水动态模拟及预测

矿区地下水动态研究是山西亟待解决的重大需求问题。针对山西矿区水资源短缺与水环境恶化等现状,通过分析气象、水文、地质地貌、开采现状等监测数据资料,以FEFLOW模型和GIS技术为平台,构建山西矿区三维地下水数值模型,进行地下水动态研究,模拟并预测了4种不同情境下矿区地下水动态过程,从而定量分析采矿活动对地下水动态的影响,揭示矿区开采对地下水系统的作用机制。结果表明：当矿区开采强度提高10%、30%、50%时,地下水整体流场和运动趋势没有大的变化,但在矿区南部地下水流场发生突变,等值线形变,且在东南部形成一个水位变化剧烈的低水位带,迫使地下水流向发生偏转。水位分析表明采矿活动会造成地下水位整体下降,下降幅度与开采强度呈正相关关系。地下水系统均衡分析表明,当开采强度保持现状或提高10%,地下水系统仍然处于正均衡,补给量大于排泄量,当开采强度提高30%、50%时,系统转为负均衡状态,补给量小于排泄量,地下水水量大幅度减少。研究成果可为矿区制定合理的开采方案,保护矿区地下水资源提供参考依据,为有效遏制矿区水环境恶化、确保矿区水安全提供科技支撑。     

干旱内陆河流域人居环境适宜性评价——以石羊河流域为例(英文)

The study employs slope,aspect,relief degree of land surface,land use,vegetation index,hydrology and climate,as evaluation indexes to set up the Human Settlements Environmental Index(HEI) model to evaluate the environmental suitability for human settlements in the Shiyang River Basin.By using GIS spatial analysis technology,such as spatial overlay analysis,buffer analysis and density analysis,the environmental suitability of the human settlement spatial situation and spatial pattern are established to analyze their spatial distribution.The results show that the index of suitability for human settlements in the Shiyang River Basin is between 17.13 and 84.32.In general,suitability for human settlements decreases from the southwest to the northeast.Seen from an area pattern,the suitable region is mainly distributed in the Minqin oasis,Wuwei oasis and Changning basin,which are about 1080.01 km 2 and account for 2.59% of the total area.Rather and comparatively suitable region is mainly distributed around the counties of Gulang,Yongchang and north of Tianzhu,which is about 1100.30 km 2.The common suitable region is mainly distributed outside the counties of Yongchang,Jinchuan and most parts of Minqin County,which are about 23328.04 km 2,accounting for 56.08% of the total area.The unsuitable region is mainly distributed upstream and to the north of the river,which is about 9937.60 km 2,accounting for 23.89% of the total area.Meanwhile,the least suitable region is distributed around the Qilian Mountains,which are covered by snow and cold desert and lie in the intersecting area between the Tengger Desert and Badain Jaran Desert.The total area is about 6154.05 km 2,accounting for 14.79% of the total area.Suitable regions for human habitation are mainly distributed around rivers in the form of ribbons and batches,while others are scattered.The distribution pattern is identical to the residential spatial pattern.In addition,the relationships between HEI and other factors have been analyzed.There is a clear logarithmic correlation between the residential environment and population,that is,the correlation coefficient between the evaluation value and population density reaches 0.851.There is also a positive correlation between the residential environment and economy,which reaches an evaluation value of 0.845 between the residential environment and GDP.Results also show that the environment is out of bearing with the existing population in Shiyang River Basin.Spatial distribution of population is profoundly affected by severe environmental problems,such as the expanded deserts,the hilly terrain and the changing climate.Surface water shortage and slow economic growth are bottlenecks for suitable human settlement in the Shiyang River Basin.Combining these problems with planning for construction of new country and the exploitation of local land,some residential areas should be relocated to improve the residential environment.     

Quantitative estimation and spatial reconstruction of urban and rural construction land in Jiangsu Province, 1820–1985

Land cover is the most evident landscape signal to characterize the influence of human activities on terrestrial ecosystems. Since the industrial revolution, the expansion of construction land has profoundly changed the status of land use coverage and changes. This study is proposed to reconstruct the spatial pattern of construction land (urban construction land and rural settlement land) for five historical periods over the past 200 years in Jiangsu Province with 200 m × 200 m grids on the basis of quantitative estimation. Urban construction land is estimated based on data about city walls, four gates along walls, and other socio-economic factors. Rural settlement land is calculated based on the rural population and per capita housing allowance. The spatial pattern of historical construction land is simulated based on the distribution of modern construction land in 1985 with a quantitative-boundarysuitability control method and thorough consideration over connectivity of different land use types. The study concludes that: (1) the amount of construction land in Jiangsu Province is estimated at 963.46 km² in 1820, 1043.46 km² in 1911, 1672.40 km² in 1936, 1980.34 km² in 1952 and 10,687.20 km² in 1985; and (2) the spatial distribution of construction land features the great proclivity to water bodies and main roads and the strong polarization of existent residence. The results are verified directly and indirectly by applying the trend verification of construction land changes and patterns, the correlation analysis between rural settlement land and local arable land, and quantitative accuracy test of the reconstructed construction land to actual historical survey maps covering four sample regions in 1936.     

Application of Time-Series Model to Predict Groundwater Dynamic in Sanjiang Plain,Northeast China

To study the groundwater dynamic in the typical region of Sanjiang Plain, long-term groundwater level observation data in the Honghe State Farm were collected and analyzed in this paper. The seasonal and long-term groundwater dynamic was explored. From 1996 to 2008, groundwater level kept declining due to intensive exploitation of groundwater resources for rice irrigation. A decline of nearly 5 m was found for almost all the monitoring wells. A time-series method was established to model the groundwater dynamic. Modeled results by time-series model showed that the groundwater level in this region would keep declining according to the current exploitation intensity. A total dropdown of 1.07 m would occur from 2009 to 2012. Time-series model can be used to model and forecast the groundwater dynamic with high accuracy. Measures including control on groundwater exploitation amount and application of water saving irrigation technique should be taken to prevent the continuing declining of groundwater in the Sanjiang Plain.     

Collapsing loess-like soil in western Saudi Arabia

Earth fissures have developed at Wadi Al-Yutamah, western Saudi Arabia. The fissues are associated with land subsidence which is considered to be due to both rapid draw down of the ground-water level and hydrocompaction of the wadi soil after flooding. This phenomenon is relatively recent in the area.The wadi soil was investigated and classified in the field, and disturbed and undisturbed samples were collected for laboratory testing and analysis. The engineering properties of the wadi soils were determined, including in situ field density, specific gravity, liquid limit, plastic limit, plasticity index, shrinkage limit, and consolidation characteristics such as total settlement and coefficient of subsidence (collapse).Four types of silty soil of different colours were identified, in addition to yellowish-brown sand (representing the dunes) which accumulated on the surface of the study area. The silty soils include yellowish-brown silty sand, yellowish-gray silt, pale brown silt and yellowish-brown clayey silt. The results of X-ray diffraction analysis indicated that the clay is mainly kaolinite and illite with minor smectite. The dominant soil type in the study area is silt of low plasticity, high void ratio and low density which decreased with depth. This soil was classified as loess-like materials.The studied soils are of a collapsing type, and settlement is greatly increased by excessive wetting under constant pressure. The calculated coefficients of subsidence (collapse) of the soils at different depths generally increased with decreasing soil density and ranged between 3·1% and 10·8%. The wadi soils are considered to pose a moderate problem when wetted.     

Generic 3D interpolation of Holocene base‐level rise and provision of accommodation space,developed for the Netherlands coastal plain and infilled palaeovalleys

We present an interpolation model that describes Holocene groundwater level rise and the creation of accommodation space in 3D in the Rhine‐Meuse delta – the Netherlands. The model area (ca. 12 400 km²) covers two palaeovalleys of Late Pleistocene age (each 30 km wide) and the overlying Holocene deposits of the Rhine‐Meuse delta, the Holland coastal plain, and the Zuiderzee former lagoon. Water table rise is modelled from 10 800 to 1000 cal. BP, making use of age‐depth relations based on 384 basal peat index points, and producing output in the form of stacked palaeo groundwater surfaces, groundwater age‐depth curves, and voxel sets. These products allow to resolve (i) regional change and variations of inland water table slopes, (ii) spatial differences in the timing and pacing of transgression, and (iii) analysis of interplay of coastal, fluvial and subsidence controls on the provision of accommodation space. The interpolation model is a multi‐parameter trend function, to which a 3D‐kriging procedure of the residuals is added. This split design deploys a generic approach for modelling provision of accommodation space in deltas and coastal lowlands, aiming to work both in areas of intermediate data availability and in the most data‐rich environments. Major provision of accommodation space occurred from 8500 cal BP onwards, but a different evolution occurred in each of the two palaeovalleys. In the northern valley, creation of accommodation space began to stall at 7500 cal BP, while in the southern valley provision of new accommodation space in considerable quantities continued longer. The latter is due to the floodplain gradient that was maintained by the Rhine, which distinguishes the fluvial deltaic environment from the rest of the back‐barrier coastal plain. The interpolation results allow advanced mapping and investigation of apparent spatial differences in Holocene aggradation in larger coastal sedimentary systems. Furthermore, they provide a means to generate first‐order age information with centennial precision for 3D geological subsurface models of Holocene deltas and valley fills. As such, the interpolation is of use in studies into past and present land subsidence and into low land sedimentation.     

Models of the rapid post‐rift subsidence in the eastern Qiongdongnan Basin,South China Sea: implications for the development of the deep thermal anomaly

The Qiongdongnan Basin is one of the largest Cenozoic rifted basins on the northern passive margin of the South China Sea. It is well known that since the Late Miocene, approximately 10 Ma after the end of the syn‐rift phase, this basin has exhibited rapid thermal subsidence. However, detailed analysis reveals a two‐stage anomalous subsidence feature of the syn‐rift subsidence deficit and the well‐known rapid post‐rift subsidence after 10.5 Ma. Heat‐flow data show that heat flow in the central depression zone is 70–105 mW m^?2, considerably higher than the heat flow (<70 mW m^?2) on the northern shelf. In particular, there is a NE‐trending high heat‐flow zone of >85 mW m^?2 in the eastern basin. We used a numerical model of coupled geothermal processes, lithosphere thinning and depositional processes to analyse the origin of the anomalous subsidence pattern. Numerical analysis of different cases shows that the stretching factor β_s based on syn‐rift sequences is less than the observed crustal stretching factor β_c, and if the lithosphere is thinned with β_c during the syn‐rift phase (before 21 Ma), the present basement depth can be predicted fairly accurately. Further analysis does not support crustal thinning after 21 Ma, which indicates that the syn‐rift subsidence is in deficit compared with the predicted subsidence with the crustal stretching factor β_c. The observed high heat flow in the central depression zone is caused by the heating of magmatic injection equivalently at approximately 3–5 Ma, which affected the eastern basin more than the western basin, and the Neogene magmatism might be fed by the deep thermal anomaly. Our results suggest that the causes of the syn‐rift subsidence deficit and rapid post‐rift subsidence might be related. The syn‐rift subsidence deficit might be caused by the dynamic support of the influx of warmer asthenosphere material and a small‐scale thermal upwelling beneath the study area, which might have been persisting for about 10 Ma during the early post‐rift phase, and the post‐rift rapid subsidence might be the result of losing the dynamic support with the decaying or moving away of the deep thermal source, and the rapid cooling of the asthenosphere. We concluded that the excess post‐rift subsidence occurs to compensate for the syn‐rift subsidence deficit, and the deep thermal anomaly might have affected the eastern Qiongdongnan Basin since the Late Oligocene.     

Modeling Short-Term Spatial and Temporal Variability of Groundwater Level Using Geostatistics and GIS

Continuous depletion of groundwater levels from deliberate and uncontrolled exploitation of groundwater resources lead to the severe problems in arid and semi-arid hard-rock regions of the world. Geostatistics and geographic information system (GIS) have been proved as successful tools for efficient planning and management of the groundwater resources. The present study demonstrated applicability of geostatistics and GIS to understand spatial and temporal behavior of groundwater levels in a semi-arid hard-rock aquifer of Western India. Monthly groundwater levels of 50 sites in the study area for 36-month period (May 2006 to June 2009; excluding 3 months) were analyzed to find spatial autocorrelation and variances in the groundwater levels. Experimental variogram of the observed groundwater levels was computed at 750-m lag distance interval and the four most-widely used geostatistical models were fitted to the experimental variogram. The best-fit geostatistical model was selected by using two goodness-of-fit criteria, i.e., root mean square error (RMSE) and correlation coefficient (r). Then spatial maps of the groundwater levels were prepared through kriging technique by means of the best-fit geostatistical model. Results of two spatial statistics (Geary’s C and Moran’s I) indicated a strong positive autocorrelation in the groundwater levels within 3-km lag distance. It is emphasized that the spatial statistics are promising tools for geostatistical modeling, which help choose appropriate values of model parameters. Nugget-sill ratio (<0.25) revealed that the groundwater levels have strong spatial dependence in the area. The statistical indicators (RMSE and r) suggested that any of the three geostatistical models, i.e., spherical, circular, and exponential, can be selected as the best-fit model for reliable and accurate spatial interpolation. However, exponential model is used as the best-fit model in the present study. Selection of the exponential model as the best-fit was further supported by very high values of coefficient of determination (r ² ranging from 0.927 to 0.994). Spatial distribution maps of groundwater levels indicated that the groundwater levels are strongly affected by surface topography and the presence of surface water bodies in the study area. Temporal pattern of the groundwater levels is mainly controlled by the rainy-season recharge and amount of groundwater extraction. Furthermore, it was found that the kriging technique is helpful in identifying critical locations over the study area where water saving and groundwater augmentation techniques need to be implemented to protect depleting groundwater resources.     

Hydrogeomorphic control on tree growth responses in the Elton area of the Cheshire Saltfield, UK

Increment cores were sampled from oak (Quercus robur) and ash (Fraxinus excelsior) growing at Elton, an area of the Cheshire Saltfield that has experienced significant subsidence and damage to the natural and built environments in the latter part of the twentieth century. Ring-width measurements for Elton trees permitted the construction of one main site chronology (Elton) and four sub chronologies (Elton A, Elton B, Elton C, Elton ASH). Ring-width difference between these and a control chronology identified periods of sustained growth reduction in oak trees commencing in AD 1859/1861, 1886 and 1934. Growth reductions after 1934 are related to watertable draw down caused by brine pumping from a concentration of nine boreholes at Elton, up to 2 km from tree sampling locations. Growth reductions in 1859/1861 and 1886 are likely to be the result of earlier phases of brine pumping in the Wheelock Valley, up to 5 km to the east of Elton, and these reductions correlate well with historic records of subsidence and pumping activity. Cessation of pumping in 1977 led to a lagged growth recovery in oak trees between 1981 and 1986, indicating that an artificial drought had been imposed on the Elton area for a period in excess off 100-y. This research demonstrates a hydrological separation of surface water and groundwater in an area where salt beds are overlain by till and that ring-width records of Q. robur can be used to reconstruct watertable variability and also the spatial impact of solution mining.     

精明收缩视角下乡村人居空间变化特征及类型划分——以江苏省常熟市为例

为应对快速城镇化带来的乡村人口收缩的现实,苏南地区提出了一系列引导乡村空间收缩发展的政策。以常熟市为例,基于2005年和2017年乡村人口统计数据和农村居民点矢量数据,定量分析县、镇、村不同尺度下乡村人居空间变迁特征及差异并进行了类型划分。研究发现,常熟市乡村人居空间利用集而不约;乡村人居活动和乡村人居空间的变化特征存在不同尺度的差异性;镇域层面人居空间变化存在城镇化主导型、产业主导型、交通主导型、产业-交通复合主导型、政策主导型等五种类型;城镇化、产业、交通和政策是引起乡村人居空间变化的主要驱动因素。研究结论可以从精明收缩视角为苏南地区乡村人居空间的转型发展提供理论支撑。     

Characterization of Mechanisms and Processes Controlling Groundwater Recharge and its Quality in Drought-Prone Region of Central India (Buldhana,Maharashtra) Using Isotope Hydrochemical and End-Member Mixing Modeling

A thorough study on understanding of groundwater recharge sources and mechanisms was attempted by integrating the hydrogeological, geochemical and isotopic information along with groundwater dating and end-member mixing analysis (EMMA). This study was necessitated due to prolonged dryness and unavailability of freshwater in semi arid Deccan trap regions of Central India. In addition, groundwater resources are not characterized well in terms of their geochemical nature and recharge sources. The hydrogeochemical inferences suggest that aquifer I consists of recently recharged water dominated by Ca–Mg–HCO₃ facies, while groundwater in aquifer II shows water–rock interaction and ion exchange processes. Presence of agricultural contaminant, nitrate, in both aquifers infers limited hydraulic interconnection, which is supported by unconfined to semi-confined nature of aquifers. Groundwater in both aquifers is unsaturated with respect to carbonate and sulfate minerals indicating lesser water–rock interaction and shorter residence time. This inference is corroborated by tritium age of groundwater (aquifer I: 0.7–2 years old and aquifer II: 2–4.2 years old). Stable water isotopes (δ²H, δ¹⁸O) suggest that groundwater is a mixture of rainwater and evaporated water (surface water and irrigation return flow). EMMA analysis indicates three groundwater recharge sources with irrigation return flow being the dominant source compared to others (rainwater and surface waters). A conceptual model depicting groundwater chemistry, recharge and dynamics is prepared based on the inferences.