Simulation and Analysis of Conjunctive Use with MODFLOW's Farm Process

The extension of MODFLOW onto the landscape with the Farm Process (MF-FMP) facilitates fully coupled simulation of the use and movement of water from precipitation, streamflow and runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. This allows for more complete analysis of conjunctive use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within “water-balance subregions” comprised of one or more model cells that can represent a single farm, a group of farms, or other hydrologic or geopolitical entities. Simulation of micro-agriculture in the Pajaro Valley and macro-agriculture in the Central Valley are used to demonstrate the utility of MF-FMP. For Pajaro Valley, the simulation of an aquifer storage and recovery system and related coastal water distribution system to supplant coastal pumpage was analyzed subject to climate variations and additional supplemental sources such as local runoff. For the Central Valley, analysis of conjunctive use from different hydrologic settings of northern and southern subregions shows how and when precipitation, surface water, and groundwater are important to conjunctive use. The examples show that through MF-FMP's ability to simulate natural and anthropogenic components of the hydrologic cycle, the distribution and dynamics of supply and demand can be analyzed, understood, and managed. This analysis of conjunctive use would be difficult without embedding them in the simulation and are difficult to estimate a priori.     

A conceptual–numerical model to simulate hydraulic head in aquifers that are hydraulically connected to surface water bodies

In this paper, we present a conceptual‐numerical model that can be deduced from a calibrated finite difference groundwater‐flow model, which provides a parsimonious approach to simulate and analyze hydraulic heads and surface water body–aquifer interaction for linear aquifers (linear response of head to stresses). The solution of linear groundwater‐flow problems using eigenvalue techniques can be formulated with a simple explicit state equation whose structure shows that the surface water body–aquifer interaction phenomenon can be approached as the drainage of a number of independent linear reservoirs. The hydraulic head field could be also approached by the summation of the head fields, estimated for those reservoirs, defined over the same domain set by the aquifer limits, where the hydraulic head field in each reservoir is proportional to a specific surface (an eigenfunction of an eigenproblem, or an eigenvector in discrete cases). All the parameters and initial conditions of each linear reservoir can be mathematically defined in a univocal way from the calibrated finite difference model, preserving its characteristics (geometry, boundary conditions, hydrodynamic parameters (heterogeneity), and spatial distribution of the stresses). We also demonstrated that, in practical cases, an accurate solution can be obtained with a reduced number of linear reservoirs. The reduced computational cost of these solutions can help to integrate the groundwater component within conjunctive use management models. Conceptual approximation also facilitates understanding of the physical phenomenon and analysis of the factors that influence it. A simple synthetic aquifer has been employed to show how the conceptual model can be built for different spatial discretizations, the parameters required, and their influence on the simulation of hydraulic head fields and stream–aquifer flow exchange variables. A real‐world case was also solved to test the accuracy of the proposed approaches, by comparing its solution with that obtained using finite‐difference MODFLOW code. Copyright © 2011 John Wiley & Sons, Ltd.     

Simple water balance modelling of surface reservoir systems in a large data-scarce semiarid region / Modélisation simple du bilan hydrologique de systèmes de réservoirs de surface dans une grande région semi-aride pauvre en données

Abstract

Abstract Water resources in dryland areas are often provided by numerous surface reservoirs. As a basis for securing future water supply, the dynamics of reservoir systems need to be simulated for large river basins, accounting for environmental change and an increasing water demand. For the State of Ceará in semiarid Northeast Brazil, with several thousands of reservoirs, a simple deterministic water balance model is presented. Within a cascade-type approach, the reservoirs are grouped into six classes according to storage capacity, rules for flow routing between reservoirs of different size are defined, and water withdrawal and return flow due to human water use is accounted for. While large uncertainties in model applications exist, particularly in terms of reservoir operation rules, model validation against observed reservoir storage volumes shows that the approach is a reasonable simplification to assess surface water availability in large river basins. The results demonstrate the large impact of reservoir storage on downstream flow and stress the need for a coupled simulation of runoff generation, network redistribution and water use.     

Role of reservoir regulation and groundwater feedback in a simulated ground-soil-vegetation continuum: A long-term regional scale analysis

The regional terrestrial water cycle is strongly altered by human activities. Among them, reservoir regulation is a way to spatially and temporally allocate water resources in a basin for multi-purposes. However, it is still not sufficiently understood how reservoir regulation modifies the regional terrestrial- and subsequently, the atmospheric water cycle. To address this question, the representation of reservoir regulation into the terrestrial component of fully coupled regional Earth system models is required. In this study, an existing process-based reservoir network module is implemented into NOAH-HMS, that is, the terrestrial component of an atmospheric–hydrologic modelling system, namely, the WRF-HMS. It allows to quantitatively differentiate role of reservoir regulation and of groundwater feedback in a simulated ground-soil-vegetation continuum. Our study focuses on the Poyang Lake basin, where the largest freshwater lake of China and reservoirs of different sizes are located. As compared to streamflow observations, the newly extended NOAH-HMS slightly improves the streamflow and streamflow duration curves simulation for the Poyang Lake basin for the period 1979–1986. The inclusion of reservoir regulation leads to major changes in the simulated groundwater recharges and evaporation from reservoirs at local scale, but has minor effects on the simulated soil moisture and surface runoff at basin scale. The performed groundwater feedback sensitivity analysis shows that the strength of the groundwater feedback is not altered by the consideration of reservoir regulation. Furthermore, both reservoir regulation and groundwater feedback modify the partitioning of the simulated evapotranspiration, thus affecting the atmospheric water cycle in the Poyang Lake region. This finding motivates future research with our extended fully coupled atmospheric–hydrologic modelling system by the community.     

A stochastic model for optimal operation of inter-basin water allocation systems: a case study

This paper presents a new methodology for optimal operation of inter-basin water transfer systems by conjunctive use of surface water resources in water donor basin and groundwater resources in water receiving basin. To incorporate the streamflow uncertainty, an integrated stochastic dynamic programming (ISDP) model is developed. In the ISDP, the monthly inflow to the reservoir in the water donor basin, the water storage of the reservoir, and the water storage of the aquifer in the water receiving basin are considered as state variables. A water allocation optimization model is embedded in the main structure of ISDP and a new ensemble streamflow prediction model based on K-nearest-neighbourhood algorithm is also developed and linked to the ISDP. By using a new reoptimization process, the ISDP model provides monthly policies for water allocation to users in water donor and receiving basins. As water users can form a coalition to increase their benefits, several solution concepts in cooperative game theory, namely Nash–Harsanyi, Shapley, Nucleolus, Weak Nucleolus, Proportional Nucleolus, Separable Costs Remaining Benefits (SCRBs) and Minimum Costs Remaining Savings are utilized to determine the profit of each water user. In the last step, stakeholders make negotiation over these solution concepts using the Fallback bargaining theory to reach a unanimous agreement on the final distribution of the total benefit. The methodology is applied to an inter-basin water transfer project and the results show that the Shapley and SCRB solutions concepts can provide better distributions for the total benefit and the total benefit of water users is increased by a factor of 1.6 when they participate in a grand coalition.     

Evaluation of Bias Associated with Capture Maps Derived from Nonlinear Groundwater Flow Models

The impact of groundwater withdrawal on surface water is a concern of water users and water managers, particularly in the arid western United States. Capture maps are useful tools to spatially assess the impact of groundwater pumping on water sources (e.g., streamflow depletion) and are being used more frequently for conjunctive management of surface water and groundwater. Capture maps have been derived using linear groundwater flow models and rely on the principle of superposition to demonstrate the effects of pumping in various locations on resources of interest. However, nonlinear models are often necessary to simulate head‐dependent boundary conditions and unconfined aquifers. Capture maps developed using nonlinear models with the principle of superposition may over‐ or underestimate capture magnitude and spatial extent. This paper presents new methods for generating capture difference maps, which assess spatial effects of model nonlinearity on capture fraction sensitivity to pumping rate, and for calculating the bias associated with capture maps. The sensitivity of capture map bias to selected parameters related to model design and conceptualization for the arid western United States is explored. This study finds that the simulation of stream continuity, pumping rates, stream incision, well proximity to capture sources, aquifer hydraulic conductivity, and groundwater evapotranspiration extinction depth substantially affect capture map bias. Capture difference maps demonstrate that regions with large capture fraction differences are indicative of greater potential capture map bias. Understanding both spatial and temporal bias in capture maps derived from nonlinear groundwater flow models improves their utility and defensibility as conjunctive‐use management tools.     

Assessing the restoration time of surface water and groundwater systems under groundwater pumping

Since surface water and groundwater systems are fully coupled and integrated, increased groundwater withdrawal during drought may reduce groundwater discharges into the stream, thereby prolonging both systems’ recovery from drought. To analyze watershed response to basin-level groundwater pumping, we propose a modelling framework to understand the resiliency of surface water and groundwater systems using an integrated hydrologic model under transient pumping. The proposed framework incorporates uncertainties in initial conditions to develop robust estimates of restoration times of both surface water and groundwater and quantifies how pumping impacts state variables such as soil moisture. Groundwater pumping impacts over a watershed were also analyzed under different pumping volumes and different potential climate scenarios. Our analyses show that groundwater restoration time is more sensitive to variability in climate forcings as opposed to changes in pumping volumes. After the cessation of pumping, streamflow recovers quickly in comparison to groundwater, which has higher persistence. Pumping impacts on various hydrologic variables were also discussed. Potential for developing optimal conjunctive management plans using seasonal-to-interannual climate forecasts is also discussed.     

中国东南丘陵山区水质良好水库现状与天目湖保护实践

东南丘陵山区是我国水库分布最为集中的区域之一,这些水库在保障区域供水安全方面具有极其重要的作用,应该优先保护.然而,水库水环境保护正面临丘陵山区开发强度持续增加,开发方式和空间布局不合理,氮、磷污染及富营养化趋势严峻,缺乏完善监测和管理体系等众多问题.本文以2000年以来天目湖水库保护实践过程为例,从库体水环境治理、流域污染物削减和综合管理3个方面介绍天目湖沙河水库保护的措施和成效,在此基础上提出良好湖库优先保护的建议:建立具有部门协调能力的水库管理机构,实行基于湖库水生态目标的水质目标管理,治理丘陵山区茶果园的面源污染,注重流域生态系统整体的恢复,禁止上游水源涵养区和临湖地带的开发,划定生态保护红线,明确禁止和限制的开发类型与规模,加强湖库及流域的监测和预警,开展必要的水库水体治理工程,并针对性地制定湖库管理条例.天目湖十几年的保护实践中有效地解决了流域开发与水库水质保护之间的矛盾,使天目湖水质由快速恶化转为稳步好转,为东南丘陵山区经济发展过程中水库的环境保护探索了一条道路.     

Characteristics of chemistry and stable isotopes in groundwater of Chaobai and Yongding River basin,North China Plain

To identify the groundwater flow system in the North China Plain, the chemical and stable isotopes of the groundwater and surface water were analysed along the Chaobai River and Yongding River basin. According to the field survey, the study area in the North China Plain was classified hydrogeologically into three parts: mountain, piedmont alluvial fan and lowland areas. The change of electrical conductance and pH values coincided with groundwater flow from mountain to lowland areas. The following groundwater types are recognized: Ca? HCO₃ and Ca? Mg? HCO₃ in mountain areas, Ca? Mg? HCO₃ and Na? K? HCO₃ in piedmont alluvial fan areas, and HCO₃? Na in lowland areas. The stable isotope distribution of groundwater in the study area also has a good corresponding relation with other chemical characteristics. Stable isotope signatures reveal a major recharge from precipitation and surface water in the mountain areas. Chemical and stable isotope analysis data suggest that mountain and piedmont alluvial fan areas were the major recharge zones and the lowland areas belong to the main discharge zone. Precipitation and surface water were the major sources for groundwater in the North China Plain. Stable isotopic enrichment of groundwater near the dam area in front of the piedmont alluvial fan areas shows that the dam water infiltrated to the ground after evaporation. As a result, from the stable isotope analysis, isotope value of groundwater tends to deplete from sea level (horizontal ground surface) to both top of the mountain and the bottom of the lowland areas in symmetrically. This suggests that groundwater in the study area is controlled by the altitude effect. Shallow groundwater in the study area belongs to the local flow system and deep groundwater part of the regional flow system. Copyright © 2007 John Wiley & Sons, Ltd.     

Evapotranspiration and its main controlling mechanism over the desert riparian forests in the lower Tarim River Basin

Evapotranspiration(ET) and its controlling mechanism over the desert riparian forests in arid regions are the important scientific basis for the water resources managements of the lower reaches of the inland rivers of China. Nearly three years of continuous measurements of surface ET, soil water content at different depths and groundwater table over a typical Tamarix spp. stand and a typical Populus euphratica stand were conducted in the lower reach of the Tarim River. The ET seasonal trends in the growing season were controlled by plant phenology, and ET in non-growing season was weak. The diurnal variations of ET resulting from the comprehensive effects of all atmospheric factors were significantly related with reference ET. The spatial pattern of ET was determined by vegetation LAI, more vegetation coverage, more ET amount. Groundwater is the water source of surface ET, and the soil water in shallow layers hardly took part in the water exchange in the groundwatersoil-plant-air system. The temporal processes of ET over the Tamarix stand and the Populus stand were similar, but the water consumption of the well-grown Populus euphratica was higher than that of the well-grown Tamarix spp. Further analysis indicates that plant transpiration accounts for most of the surface ET, with soil evaporation weak and negligible; groundwater table is a crucial factor influencing ET over the desert riparian forests, groundwater influences the processes and amounts of ET by controlling the growth and spatial distribution of desert riparian forests; quantifying the water stress of desert riparian forests using groundwater table is more appropriate, rather than soil water content. Based on the understanding of ET and water movements in the groundwater-soil-plant-air system, a generalized framework expressing the water cycling and its key controlling mechanism in the lower reaches of the inland rivers of China is described, and a simple model to estimate water requirements of the desert riparian forests is presented.     

Quantifying Modern Recharge and Depletion Rates of the Nubian Aquifer in Egypt

Egypt is currently seeking additional freshwater resources to support national reclamation projects based mainly on the Nubian aquifer groundwater resources. In this study, temporal (April 2002 to June 2016) Gravity Recovery and Climate Experiment (GRACE)-derived terrestrial water storage (TWS_GRACE) along with other relevant datasets was used to monitor and quantify modern recharge and depletion rates of the Nubian aquifer in Egypt (NAE) and investigate the interaction of the NAE with artificial lakes. Results indicate: (1) the NAE is receiving a total recharge of 20.27 ± 1.95 km³ during 4/2002?2/2006 and 4/2008–6/2016 periods, (2) recharge events occur only under excessive precipitation conditions over the Nubian recharge domains and/or under a significant rise in Lake Nasser levels, (3) the NAE is witnessing a groundwater depletion of ? 13.45 ± 0.82 km³/year during 3/2006–3/2008 period, (4) the observed groundwater depletion is largely related to exceptional drought conditions and/or normal baseflow recession, and (5) a conjunctive surface water and groundwater management plan needs to be adapted to develop sustainable water resources management in the NAE. Findings demonstrate the use of global monthly TWS_GRACE solutions as a practical, informative, and cost-effective approach for monitoring aquifer systems across the globe.     

Sustainable conjunctive use of groundwater for additional irrigation

The key to ‘sustainable conjunctive use of groundwater for additional irrigation’ is the salt balance of groundwater below an irrigated field. This paper aims to develop a mathematical tool to study the accumulation of salt in the groundwater below an irrigated field as caused by irrigation recirculation. This study derives a salt balance of groundwater to ensure that the additional irrigation from groundwater remains possible in the future. The water and salt budgets by themselves do neither provide information concerning farmers' options nor on the limits of the individual terms in the budget equations. It is presumed that farmers will intuitively aim for (1) an optimal value of the actual evapotranspiration, and (2) a return flow as a feasible low fraction of the available water. We, therefore, derive the irrigation from groundwater Q as a consequence of the predefined farmers' aims to achieve a high actual evapotranspiration in combination with a given optimally used irrigation system. Our model concludes that the required amount of drainage is only dependent on the ratio of the salinity in the surface irrigation water and the acceptable salinity of the groundwater. The final salinity in the saturated zone only depends on salt‐carrying inflows and outflows. From the aforesaid model, it is further concluded that sustainable conjunctive use of groundwater for additional irrigation requires long‐term salt management, which should be founded on the essential controlling factors as derived in this paper. Copyright © 2013 John Wiley & Sons, Ltd.     

Isolating the impacts of anthropogenic water use within the hydrological regime of north India

The effects of anthropogenic water use play a significant role in determining the hydrological cycle of north India. This paper explores anthropogenic impacts within the region's hydrological regime by explicitly including observed human water use behaviour, irrigation infrastructure and the natural environment in the CHANSE (Coupled Human And Natural Systems Environment) socio-hydrological modelling framework. The model is constrained by observed qualitative and quantitative information collected in the study area, along with climate and socio-economic variables from additional sources. Four separate scenarios, including business as usual (BAU, representing observed irrigation practices), groundwater irrigation only (where the influence of the canal network is removed), canal irrigation only (where all irrigation water is supplied by diverted surface water) and rainfed only (where all human interventions are removed) are used. Under BAU conditions the modelling framework closely matched observed groundwater levels. Following the removal of the canal network, which forces farmers to rely completely on groundwater for irrigation, water levels decrease, while under a canal-only scenario flooding occurs. Under the rainfed-only scenario, groundwater levels similar to current business-as-usual conditions are observed, despite much larger volumes of recharge and discharge entering and leaving the system under BAU practices. While groundwater abstraction alone may lead to aquifer depletion, the conjunctive use of surface and groundwater resources, which includes unintended contributions of canal leakage, create conditions similar to those where no human interventions are present. Here, the importance of suitable water management practices, in maintaining sustainable water resources, is shown. This may include augmenting groundwater resources through managed aquifer recharge and reducing the impacts on aquifer resources through occasional canal water use where possible. The importance of optimal water management practices that highlight trade-offs between environmental impact and human wellbeing are shown, providing useful information for policy makers, water managers and users. © 2019 John Wiley & Sons, Ltd.     

Assessment of Managed Aquifer Recharge Site Suitability Using a GIS and Modeling

We completed a two‐step regional analysis of a coastal groundwater basin to (1) assess regional suitability for managed aquifer recharge (MAR), and (2) quantify the relative impact of MAR activities on groundwater levels and sea water intrusion. The first step comprised an analysis of surface and subsurface hydrologic properties and conditions, using a geographic information system (GIS). Surface and subsurface data coverages were compiled, georeferenced, reclassified, and integrated (including novel approaches for combining related datasets) to derive a spatial distribution of MAR suitability values. In the second step, results from the GIS analysis were used with a regional groundwater model to assess the hydrologic impact of potential MAR placement and operating scenarios. For the region evaluated in this study, the Pajaro Valley Groundwater Basin, California, GIS results suggest that about 7% (15 km²) of the basin may be highly suitable for MAR. Modeling suggests that simulated MAR projects placed near the coast help to reduce sea water intrusion more rapidly, but these projects also result in increased groundwater flows to the ocean. In contrast, projects placed farther inland result in more long‐term reduction in sea water intrusion and less groundwater flowing to the ocean. This work shows how combined GIS analysis and modeling can assist with regional water supply planning, including evaluation of options for enhancing groundwater resources.     

油田放射性异常成因研究

为探索地表放射性异常与地下油气藏的关系，在４种类型的９个油气藏上收集了３１０口井的测井数据，同时在５个油气藏上进行了地面放射性测量．上述资料证实，异常从油层到地表具有延续性，在同一油层上方放射性异常随测点与油层距离加大而减弱，复合式油气藏各油层异常独立存在，不同类型油藏异常总体形态相似，异常幅值与油藏圈闭条件有关．这些为油气物质垂向运移及放射性异常深部形成机制提供了直接证据．在此基础上，综合前人研究成果，提出了油田放射性异常物质来源主要是生油岩的观点和油气盆地内流体微运移循环的设想，建立了油田放射性异常综合成因模式．     

Seasonal expansion and contraction of stream networks in shallow groundwater systems

Surface water and groundwater are normally closely connected in areas with shallow aquifer systems. Stream systems can thus be considered as the outcrops of associated groundwater flows in areas with a shallow groundwater table and a previous subsurface. This situation prevails in sandy lowland areas where almost all rainfall percolates into the subsurface so that the surplus over evapotranspiration becomes part of a groundwater drainage system before it reappears at the surface in a stream. The stream network, being the interface with the groundwater system, must have the capacity to release the seasonally dependent precipitation surplus through the continuum of ground and surface waters. A river network therefore consists of a hierarchical system of different order and incision depth, of which the discharge-contributing component contracts and expands with the seasonal fluctuation in recharge and water table depth.

Coupling the mathematical expressions for groundwater drainage and stream flow enables development of a conjunctive model which relates the properties of a seasonally contracting and expanding stream network and related groundwater level fluctuation to the seasonal rainfall character for given geological and geomorphological conditions. This model further allows for assessment of drainage network response to a changing environment.     

Multi-objective optimization for sustainable groundwater resource management in a semiarid catchment

Abstract

Groundwater is an important water resource and its management is vital for integrated water resources development in semiarid catchments. The River Shiyang catchment in the semiarid area of northwestern China was studied to determine a sustainable multi-objective management plan of water resources. A multi-objective optimization model was developed which incorporated water supplies, groundwater quality, ecology, environment and economics on spatial and temporal scales under various detailed constraints. A calibrated groundwater flow model was supplemented by grey simulation of groundwater quality, thus providing two lines of evidence to use in the multi-objective water management. The response matrix method was used to link the groundwater simulation models and the optimization model. Multi-phase linear programming was used to minimize and compromise the objectives for the multi-period, conjunctive water use optimization model. Based on current water demands, this water use optimization management plan was able to meet ecological, environmental and economic objectives, but did not find a final solution to reduce the overall water deficit within the catchment.     

Evaluation of hydrologic data obtained from a local groundwater monitoring network in a metropolitan city,Korea

In the late 1980s, dramatic increases in water use caused over‐exploitation of groundwater resources and deterioration of water quality in Seoul metropolitan city. To monitor changes in quantity of groundwater resources and their quality, the metropolitan government established a local groundwater monitoring network in 1997 consisting of 119 monitoring wells. Groundwater resources in the urban area were affected by various human activities, including underground construction such as subways, pumping for public or private water use, leaky sewer systems and pavements. The variation patterns of the groundwater levels were mainly classified into four types, reflecting natural recharge due to rainfall events during the wet season, artificial recharge from leaky sewer or water supply systems, and heavy groundwater pumping for drainage or flood control purposes at underground construction sites. Significantly decreasing trends of groundwater levels in the suburbs of Seoul indicate groundwater use for various agricultural activities. Subway construction lowered the water level by an average of 25 m. Electrical conductivity values showed a wide range, from 100 to 1800 µS/cm (mean 470 µS/cm). Groundwater temperature generally showed a stable pattern, except for some sensitive increases at relatively shallow monitoring wells. Detailed analysis of the monitored groundwater data would provide some helpful implications for optimal and efficient management of groundwater resources in this metropolitan city. Copyright © 2005 John Wiley & Sons, Ltd.     

福建水库空间分布特征:沿海密度高水量少、内陆密度低水量多

水库是福建重要的水资源.通过2013-2015年遥感影像结合Google Earth和天地图提取福建水库3353座,分布在81个县区,总面积647.51 km2,约占全省土地面积的0.5%;其中面积≤1 km2水库3248座,总面积197.16 km2,面积1 km2水库105座,总面积450.35 km2.基于经验公式估算总蓄水量188.18亿m3,其中小型水库3078座(91.80%),蓄水总量37.06亿m3(19.69%),大中型水库275座(8.20%),蓄水总量151.12亿m3(80.31%).从空间分布格局来看,福建水库水资源空间分布不均,沿海六市水库密度大于三个内陆市,大中型水库主要分布于福建西北部,蓄水量呈现西北多、东南少的特点.单位陆地面积水库数量沿海城市县区均多于内陆,而单位人口水库数量则相反;单位面积水库蓄水量沿海与内陆差异不大,而人均蓄水量则沿海大部分县区远小于内陆.仅以水库作为供水水源,不能满足沿海地区用水,但内陆地区供水充足.水库蓄水对河流水体的平均滞留时间为0.053~0.341 a,除晋江流域受水库蓄水强烈影响外,其他流域受水库中等程度影响.     

Development and application of a stochastic optimization model for groundwater management: crop pattern and conjunctive use consideration

The need for irrigation water in arid and semi-arid regions is mostly supplied by groundwater. Furthermore, the agricultural development in these areas is not generally based on a comprehensive plan, which can cause aquifers depletion. On the other hand, to properly manage an aquifer and to have an optimal crop plan, the stochastic nature of the different parameters of a groundwater system such as groundwater recharge and water demands should be taken into consideration. In this paper, we develop an explicit stochastic optimization model for Firouzabad aquifer in Iran. This formulation is based on the first and second moment analysis for groundwater head which has been initially proposed for surface water resources management by Fletcher and Ponnambalam. We extend the model to create a new random withdrawal policy for conjunctive use setting in which the randomness in available precipitation is taken into account. The interesting point is that the model provides the respective probabilities of shortage and surplus without imposing the extra decision variables into the optimization model. A genetic-based algorithm is used to solve the stochastic nonlinear and non-convex formulation. The outcome results indicate that the current crop pattern should be changed, that is, the allocated areas of some crops have to be meaningfully reduced. Finally, to validate our model efficiency, we demonstrate that how much close the statistical characteristics obtained from the optimization model are to those estimated from the Monte Carlo simulation. Furthermore, the optimal benefits obtained using the proposed optimization model are as suitable as the benefits achieved using the corresponding Monte Carlo-based optimization model.