Assessing hydro-morphological changes in Mediterranean stream using curvilinear grid modeling approach - climate change impacts

The objective of this work was the estimation of time-space hydraulic (water depth, flow velocity) and morphological (sediment transport and bank erosion) characteristics in the downstream part of a Mediterranean stream under current and future climatic conditions. The two-dimensional hydraulic model MIKE 21C was used, which has been developed specifically to simulate 2D flow and morphological changes in rivers. The model is based on an orthogonal curvilinear grid and comprises two parts: (a) the hydrodynamic part and (b) the morphological changes part. The curvilinear grid and the bathymetry file were generated using a very high-resolution DEM (1 m × 1 m). Time series discharge data from a hydrometric station introduced in the hydrodynamic part of the model. Regarding the morphological part of the model, field measurements of suspended sediment concentration and of bank erosion were used. The model was calibrated and verified using field data that were collected during high and low flow discharges. Model simulation was in good agreement with field observations as indicated by a variety of statistical measures. Next, for predicting the riverbank change, future meteorological data and river flow data for the next 10 years (2017–2027) were employed. These data series were created according to a lower and a higher emission climate change scenario. Based on the results, an increase in rainfall intensity may cause significant changes in river banks after 10 years (more than 5 m of soil loss in river meanders). Using the obtained simulation results, extreme hydrological events such as floods transporting large sediment loads and changes in river morphology can be monitored. The proposed methodology was applied to the downstream part of the Koiliaris River Basin in Crete, Greece.     

Effect of irrigation pumpage during drought on karst aquifer systems in highly agricultural watersheds: example of the Apalachicola-Chattahoochee-Flint river basin,southeastern USA

In the Apalachicola-Chattahoochee-Flint (ACF) river basin in Alabama, Georgia, and Florida (USA), population growth in the city of Atlanta and increased groundwater withdrawal for irrigation in southwest Georgia are greatly affecting the supply of freshwater to downstream regions. This study was conducted to understand and quantify the effect of irrigation pumpage on the karst Upper Floridan Aquifer and river–aquifer interactions in the lower ACF river basin in southwest Georgia. The groundwater MODular Finite-Element model (MODFE) was used for this study. The effect of two drought years, a moderate and a severe drought year, were simulated. Comparison of the results of the irrigated and non-irrigated scenarios showed that groundwater discharge to streams is a major outflow from the aquifer, and irrigation can cause as much as 10 % change in river–aquifer flux. The results also show that during months with high irrigation (e.g., June 2011), storage loss (34 %), the recharge and discharge from the upper semi-confining unit (30 %), and the river–aquifer flux (31 %) are the major water components contributing towards the impact of irrigation pumpage in the study area. A similar scenario plays out in many river basins throughout the world, especially in basins in which underlying karst aquifers are directly connected to a nearby stream. The study suggests that improved groundwater withdrawal strategies using climate forecasts needs to be developed in such a way that excessive withdrawals during droughts can be reduced to protect streams and river flows.     

Influences of the unsaturated,saturated, and riparian zones on the transport of nitrate near the Merced River,California, USA

Transport and transformation of nitrate was evaluated along a 1-km groundwater transect from an almond orchard to the Merced River, California, USA, within an irrigated agricultural setting. As indicated by measurements of pore-water nitrate and modeling using the root zone water quality model, about 63% of the applied nitrogen was transported through a 6.5-m unsaturated zone. Transport times from recharge locations to the edge of a riparian zone ranged from approximately 6 months to greater than 100 years. This allowed for partial denitrification in horizons having mildly reducing conditions, and essentially no denitrification in horizons with oxidizing conditions. Transport times across a 50–100-m-wide riparian zone of less than a year to over 6 years and more strongly reducing conditions resulted in greater rates of denitrification. Isotopic measurements and concentrations of excess N₂ in water were indicative of denitrification with the highest rates below the Merced River. Discharge of water and nitrate into the river was dependent on gradients driven by irrigation or river stage. The results suggest that the assimilative capacity for nitrate of the groundwater system, and particularly the riverbed, is limiting the nitrate load to the Merced River in the study area.     

Modeling groundwater/surface-water interactions in an Alpine valley (the Aosta Plain,NW Italy): the effect of groundwater abstraction on surface-water resources

A groundwater flow model of the Alpine valley aquifer in the Aosta Plain (NW Italy) showed that well pumping can induce river streamflow depletions as a function of well location. Analysis of the water budget showed that ～80% of the water pumped during 2 years by a selected well in the downstream area comes from the baseflow of the main river discharge. Alluvial aquifers hosted in Alpine valleys fall within a particular hydrogeological context where groundwater/surface-water relationships change from upstream to downstream as well as seasonally. A transient groundwater model using MODFLOW2005 and the Streamflow-Routing (SFR2) Package is here presented, aimed at investigating water exchanges between the main regional river (Dora Baltea River, a left-hand tributary of the Po River), its tributaries and the underlying shallow aquifer, which is affected by seasonal oscillations. The three-dimensional distribution of the hydraulic conductivity of the aquifer was obtained by means of a specific coding system within the database TANGRAM. Both head and flux targets were used to perform the model calibration using PEST. Results showed that the fluctuations of the water table play an important role in groundwater/surface-water interconnections. In upstream areas, groundwater is recharged by water leaking through the riverbed and the well abstraction component of the water budget changes as a function of the hydraulic conditions of the aquifer. In downstream areas, groundwater is drained by the river and most of the water pumped by wells comes from the base flow component of the river discharge.     

Particle tracking analysis of river–aquifer interaction via bank infiltration techniques

Induced bank infiltration (BI) is commonly implemented in other countries, but remains new and unexplored in Malaysia. Increasing river pollution could affect drinking water resources. Given the threat of pollution to raw water sources, applying induced BI to sustain water management is essential. This paper presents a case study of the BI method, which evaluates the effects of groundwater pumping and BI operation on the installation of wells as well as determines the effect of pumping rate on flow paths, travel time, the size of the pumping and capture zone delineation, and groundwater mixing in a pumping well in Jenderam Hilir, Malaysia. The proposed method performs infiltration safely and achieves the ideal pumping rate. Numerical modeling packages, MODFLOW and MODPATH (particle tracking) were used. Results indicate that the migration of river water into the aquifer is generally slow and depends on the pumping rate and distance from well to the river. Most water arrives at the well by the end of a pumping period of 1–5 days at 3,072 m³/day for test wells DW1 and DW2, and during simultaneous pumping for DW2 and PW1 for a well located 36 and 18 m, respectively, from the river. During the 9.7-day pumping period, 33 % of the water pumped from the DW1 well was river water, and 38 % from DW2 throughout 4.6 days was river water. The models provide necessary information for water operators in the design and construction of pumping and sampling schedules of BI practices.     

Groundwater–river interaction and management in the context of inter-basin transfers

This study is focused on the management of coupled groundwater–river interaction (GRI) in the context of inter-basin transfers. The semi-coupled models, HEC-RAS and MODFLOW, were utilized in this study to replicate the River Tees experiment that had been conducted in 1976. Then the ‘surrogate’ models were applied to two “what if” scenarios in order to understand and manage GRI behavior and the nature of the reservoir control system. The results on the River Tees showed that the total temporary bankside storage volume was up to 17.5 % of the total reservoir release, given about 22 h of travel time. In addition, zones 4 and 5 of the river model with significant alluvium deposits were the best storage sites with values of 4.5 and 6.7 % of the total reservoir release, respectively, estimated as temporary losses. The final direct runoff received downstream at broken scar was about 95 % of the Cow Green total reservoir release. Furthermore, the closure rate proved to be vital to compute the effective return flow in respect of this study area. Finally, with careful reservoir and groundwater management, the River Tees proved to be an ideal natural aqueduct for inter-basin transfers when used as part of a major regional water resources scheme.     

Effect of the addition of granitic powder to an acidic soil from Galicia (NW Spain) in comparison with lime

The impact of diffuse pollution, agricultural land use and climate change on the long-term response of subsurface–surface water quality is not well understood, but is a prerequisite for evaluation of water management options. The goal of this study is to model geochemical evolution of water chemistry from the infiltration through soil into the unsaturated zone, transport through bedrocks and granular aquifers to a river in order to identify zones of steep concentration gradients and high dynamics under transient flow conditions. A numerical model was constructed comprising a 2-D 1,500 m × 150 m vertical cross-section of typical sedimentary rock formations, a glacio-fluvial quaternary gravel aquifer in the valley and soil layers. The model coupled saturated/un-saturated flow and reactive transport under steady state and transient conditions. Geochemical interactions, include intra-aqueous kinetic reactions of oxygen with dissolved organic matter, as well as kinetics of carbonate dissolution/precipitation. This model section was chosen to provide insight in to the principal processes and time scales affecting water chemistry along different flow paths. The numerical simulator MIN3P was used, a finite volume program for variably saturated subsurface flow and multi-component reactive transport. The results show that subsurface water residence times range from approximately 2 to 2,000 years. Different zones are to be expected with respect to the development of mineral equilibria; namely, purely atmospherically influenced, as well as open and closed system carbonate dissolution. Short-term responses to daily averaged changes in precipitation, however, are only visible to some extent in the shallower and near-river parts of flow system and solute loads. This can most likely be explained by directional changes in flow paths, indicating that equilibrium geochemical condition predominate at the hillslope scale, i.e. water quality depends on transport pathways rather than on kinetic effects. The extent of reducing conditions is controlled by the presence of organic-rich layers (i.e. peat deposits), the dissolution kinetics of aquifer organic matter and the subsequent mixing with oxygenated water by hydrodynamic dispersion.     

Plio-Pleistocene landscape evolution in Northern Switzerland

Re-evaluation of the river history, palaeosurface levels and exhumation history in northern Switzerland for the last 10 million years reveals that distinct morphotectonic events about 4.2 and 2.8 million years ago (Ma) caused major reorganisation of river networks and morphosculpture. As a result of the earlier formation of the Swiss Jura, potential relief energy in the piggy-back North Alpine Foreland Basin (NAFB) of northern central Switzerland south of the Jura fold belt was built up after 11–10 Ma. It was suddenly released by river capture at about 4.2 Ma when the Aare-Danube was captured by a tributary of the Rhône-Doubs river system which rooted southeast of the Black forest. This event triggered rapid denudation of weakly consolidated Molasse sediments, in the order of about 1 km, as constrained by apatite fission track data from drillholes in the NAFB. Likely mechanisms of river capture are (a) headward erosion of Rhône-Doubs tributaries, (b) uplift and rapidly increasing erosion of the Swiss Alps after about 5.3 Ma, and (c) gravel aggradation at the eastern termination of the Jura fold belt in the course of eastward and northward tilt of the piggy-back NAFB. A morphotectonic event between 4.2 and 2.5 Ma, probably at about 2.8 Ma, caused a phase of planation, accompanied by local gravel aggradation and temporary storage of Alpine debris. Between 2.8 and 2.5 Ma, the Aare-Rhône river system is cannibalised by the modern Rhine River, the latter later connecting with the Alpine Rhine River.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Quantitative and qualitative simulation of groundwater by mathematical models in Rafsanjan aquifer using MODFLOW and MT3DMS

Agriculture sector by using 80% of freshwater is the greatest water consumer in Iran. Excessive use of agricultural fertilizers in last decade, caused accumulation of enormous amounts of salts and subsequence declined the physical properties of soil. In desert and dry regions such as Rafsanjan plain, use of the groundwater resources is more than the surface water resources. Therefore, information about the quality of these resources remains a necessary task for optimum management, protection of water resources, and stopping the future damages. In this study, the groundwater quantity and quality of Rafsanjan plain was investigated by MODFLOW and MT3DMS. The presented quantitative model for this aquifer was compared by observed data and calibrated. This model was used to predict a 10-year period. Results show that water elevation decreases approximately 15 m for 10 years to come in this plain. Qualitative model results show that most quality parameters will increase. Electrical conductivity will increase more than other parameters. As values of this parameter will reach 16,000 µs/l for next 5 years. Therefore, we suggest that exploitation of water from these resources should be reduced and discharge from some of agricultural wells stop; also we suggested that recharge to groundwater resources should be increased and agricultural activities should be limited or improved using of modern irrigation systems in this plain.     

Assessment of flood inundation mapping of Surat city by coupled 1D/2D hydrodynamic modeling: a case application of the new HEC-RAS 5

Surat city of India, situated 100 km downstream of Ukai Dam and 19.4 km upstream from the mouth of River Tapi, has experienced the largest flood in 2006. The peak discharge of about 25,770 m³ s^?1 released from the Ukai Dam was responsible for a disaster. To assess the flood and find inundation in low-lying areas, simulation work is carried out under the 1D/2D couple hydrodynamic modeling. Two hundred ninety-nine cross sections, two hydraulic structures and five major bridges across the river are considered for 1D modeling, whereas a topographic map at 0.5 m contour interval was used to produce a 5 m grid, and SRTM (30 and 90 m) grid has been considered for Surat and the Lower Tapi Basin. The tidal level at the river mouth and the release from the Ukai Dam during 2006 flood are considered as the downstream and upstream boundaries, respectively. The model is simulated under the unsteady flow condition and validated for the year 2006. The simulated result shows that 9th August was the worst day in terms of flooding for Surat city and a maximum 75–77% area was under inundation. Out of seven zones, the west zone had the deepest flood and inundated under 4–5 m. Furthermore, inundation is simulated under the bank protection work (i.e., levees, retaining wall) constructed after the 2006 flood. The simulated results show that the major zones are safe against the inundation under 14,430 m³ s^?1 water releases from Ukai Dam except for the west zone. The study shows the 2D capability of new HEC-RAS 5 for flood inundation mapping and management studies.     

Numerical modeling of the Rideau Valley Watershed

Using the Mike11 modeling system by the Danish Hydraulic Institute, a detailed model of the Rideau Valley Watershed was constructed. It includes 532 km of rivers and lakes, 106 basins, 122 bridges and culverts, and 20 water control structures. The model was calibrated using measured streamflow data for a time period of 5 years; additional 5 years of data was used for validation. Various methods, both qualitative and quantitative, were used to evaluate the model performance. It was found that the model can simulate the hydrological response with a reasonable to high degree of accuracy. This model is now being used for various watershed management purposes, including flood forecasting, dam safety assessment, quantification of wetland functions, and derivation of design flows.     

Chloride contamination in construction aggregates due to periodic saline water intrusion: a case study in the Kaluganga River Estuary, Sri Lanka

The Kaluganga River Estuary is one of the main sources of construction sand in Sri Lanka. Salt water intrusion along this estuary due to extensive sand mining has increased over the years. Thus, the focus of the current research is to understand the relationship between river sand mining, salt water intrusion, and the resultant effects on construction sand. Two surveys were conducted along the Kaluganga Estuary along an 11 km stretch from the river mouth at predetermined intervals to measure depth water quality profiles, and to collect sediment samples. These surveys were carried out during maximum spring tide; first in a dry period and then in a wet period, to understand hydrographic effects on the quality of river sands. Sand samples were analysed for absolute chloride content and grain size distribution. Results showed significant salt water intrusion during the dry period, averaging 2,307 μS cm^?1 in surface waters throughout the surveyed 11 km stretch along with 3,818 μS cm^?1 (average) in bottom waters up to 5.6 km upstream from the river mouth causing above normal chloride content in the bottom sandy sediments. The high chloride content in bottom sands was recorded up to 5.5 km from the river mouth making them unsuitable for construction purposes. However, during wet period, salt water intrusion levels in the bottom waters were insignificant (average 61 μS cm^?1) and the chloride content in bottom sediments was very low. This study highlighted the requirement for regulations on river estuary sandmining for construction purposes.     

A GIS-based DRASTIC model for assessing groundwater vulnerability in the Ordos Plateau, China

Groundwater plays a key role in arid regions as the majority of water is supplied by it. Groundwater pollution is a major issue, because it is susceptible to contamination from land use and other anthropogenic impacts. A study was carried out to build a vulnerability map for the Ordos Plateau using the DRASTIC model in a GIS environment. The map was designed to show the areas of the highest potential for groundwater pollution based on hydrogeological conditions. Seven environmental parameters, such as depth to water table, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and hydraulic conductivity of the aquifer, were incorporated into the DRASTIC model and GIS was used to create a groundwater vulnerability map by overlaying the available data. The results of this study show that 24.8 % of the study area has high pollution potential, 24.2 % has moderate pollution potential, 19.7 % has low pollution potential, and the remaining 31.3 % of the area has no risk of groundwater pollution. The regional distribution of nitrate is well correlated with the DRASTIC vulnerability index. In contrast to this, although the DRASTIC model indicated that the western part had no risk, nitrate concentrations were higher in some of these areas. In particular, higher nitrate concentrations were recorded along river valleys and around lakes, such as the Mulin River valley. This is mainly caused by the intensive agricultural development and favorable conditions for recharge along river valleys.     

Assessment of the effects of climate and land cover changes on river discharge and sediment yield,and an adaptive spatial planning in the Jakarta region

In Jakarta, climate change has been detected through rising air temperatures, increased intensity of rainfall in the wet season, and sea level rise. The coupling of such changes with local anthropogenic driven modifications in the environmental setting could contribute to an increased probability of flooding, due to increase in both extreme river discharge and sedimentation (as a result of erosion in the watersheds above Jakarta and as indicated by sediment yield in the downstream area). In order to respond to the observed and projected changes in river discharge and sediment yield, and their secondary impacts, adaptation strategies are required. A possible adaptation strategy is through policy making in the field of spatial planning. For example, in Indonesia, presidential regulation number 54 year 2008 (Peraturan Presiden Nomor 54 Tahun 2008—Perpres 54/2008) was issued as a reference for the implementation of water and soil conservation. This paper assesses the impact of climate and land cover change on river discharge and sediment yield, as well as the effects of Perpres 54/2008 on that river discharge and sediment yield. The spatial water balance model Spatial Tools for River Basins and Environmental and Analysis of Management Option was used for the runoff computations, whilst the Spatial Decision Assistance of Watershed Sedimentation model was used to simulate erosion, Sediment Delivery Ratio, and sediment yield. The computation period is from January 1901 to December 2005, at the scale of the following watersheds: Ciujung, Cisadane, Ciliwung, and Citarum. During the twentieth century, computed average discharge in the downstream area (near Jakarta) increased between 2.5 and 35 m³/s/month, and sediment yield increased between 1 × 10³ and 42 × 10³ tons/year. These changes were caused by changes in both land cover and climate, with the former playing a stronger role. Based on a computation under a theoretical full implementation of the spatial plan proposed by Perpres 54/2008, river discharge would decrease by up to 5 % in the Ciliwung watershed and 26 % in the Cisadane watershed. The implementation of Perpres 54/2008 could also decrease the sediment yield, by up to 61 and 22 % in the Ciliwung and Cisadane watersheds, respectively. These findings show that the implementation of the spatial plan of Perpres 54/2008 could significantly improve watershed response to runoff and erosion. This study may serve as a tool for assessing the reduction in climate change impacts and evaluating the role of spatial planning for adaptation strategies.     

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin,northwest China

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?10⁴ m³/year in 2017 to 11,875?×?10⁴ m³/year in 2021, and to 17,039?×?10⁴ m³/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?10⁴ m³/year in 2017 to 1857?×?10⁴ m³/year in 2021, and to 510?×?10⁴ m³/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.     

Applying image recording and identification for measuring water stages to prevent flood hazards

Information on river stages can be transmitted to relevant management offices over a network by using an automatic stage gauge so that management offices can effectively evaluate whether a river stage is exceeding the warning line and take action if necessary. However, current river stage measurement methods can present this information only in data form because the methods cannot simultaneously obtain images such as the rising or overflow of the river stage. In addition, the stage gauge can fail or be washed away when the river stage is high. To solve these problems, this research evaluates a novel measurement method involving a video surveillance system that exhibits features such as easy installation, low maintenance cost, and low failure possibility. Through on-site image recording, this measurement method involves using image identification technology to read water level figures automatically. This method offers instant river stage figures and on-site video so that disaster prevention measures can be implemented accordingly. The results of a dynamic water flow test conducted in an indoor experimental channel indicated that all of the average absolute error levels of river stage identification were less than ±1.2 %, meaning that the image identification technology could achieve identification results at any flow level. By contrast, the findings of a rainfall simulation experiment suggest that the average absolute error of river stage identification was less than ±2.5 %, meaning that the measurement technology and method used in this research are useful and feasible at various rainfall intensities.     

Runoff and water budget of the Liudaogou Catchment at the wind–water erosion crisscross region on the Loess Plateau of China

The objectives of this study were to examine the runoff characteristics and to estimate water budget at the wind–water erosion crisscross region on the Loess Plateau of China. A small catchment known as Liudaogou that has representative meteorological and hydrological conditions of the wind–water erosion crisscross region was chosen as the study location. A numerical model for rainfall-runoff was developed and verified; rainfall-runoff calculation for 5 years (2005–2009) was performed. The observed data and numerical result of the surface runoff were used for evaluating runoff characteristics and estimating the annual water budget. Runoff rate was proportional to average intensity of rain. Even though rainfall duration was for few minutes, surface runoff was generated by intensity of more than 2.6 mm × 5 min^?1, when rainfall duration exceeded 10 h; surface runoff was generated by an intensity of 0.6 mm × 5 min^?1, while annual runoff rate was 10–15 %. The unit area of 1 km² was adopted as the index area for estimating annual water budget. Runoff, evapotranspiration, variation of water storage, and habitant water consumption accounted for 20.4, 75.6, 0, and 4 % of the total annual precipitation, respectively. Results of this study provide the basis for further research on hydrology, water resources, and sustainable water development and utilization at the wind–water erosion crisscross region on the northern Loess Plateau where annual water resources are relatively deficient.     

Development of Tide–Surge Interaction Model for the Coastal Region of Bangladesh

In this paper, a two-dimensional, vertically integrated hydrodynamic model is developed taking into account entrained air bubbles during storm surges as well as incorporating inverted barometer, and river and land dynamics. The model is specifically designed for the coastal region of Bangladesh. A nested scheme method with a fine mesh scheme (FMS), capable of incorporating the complex coastline and all major offshore islands accurately, nested into a coarse mesh scheme (CMS) covering up to 15° N latitude in the Bay of Bengal is used. To incorporate the small and big offshore islands in the Meghna estuarine region with its complex coastline accurately, a very fine mesh scheme (VFMS) is again nested into the FMS. Along the northeast corner of the VFMS, the Meghna river discharge is taken into account. The coastal and island boundaries are approximated through proper stair steps. The model equations are solved by a semi-implicit finite difference technique using a staggered C-grid. A stable appropriate tidal condition over the model domain is generated by applying tidal forcing with the four major tidal constituents M ₂, S ₂, K ₁, and O ₁ along the southern open boundary of the CMS. This tidal regime is introduced as the initial state of the sea for nonlinear interaction of tide and surge. The model is applied to simulate water levels due to the interaction of tide and surge associated with the cyclones April 1991 and Aila at different coastal and island locations along the coast of Bangladesh. The results are found to be quite satisfactory with root mean square error of ～0.50 m as calculated for both the storm events. Tests of sensitivities on water levels are carried out for air bubbles, offshore islands, river discharge, inverse barometer, and grid resolution. The presence of air bubbles increases simulated water levels a little bit in our model, and the contribution of air bubbles in increasing water level is found around 2 %. Further, water levels are found to be influenced by offshore islands, river discharge, inverse barometer as well as grid resolution.     

Aquifer system incremental linear property and its application on groundwater management

Generally, an aquifer system coupled into a groundwater management model was regarded as a linear system. However, in terms of systems analysis, the aquifer system can be proven to be an incremental linear system rather than a linear system. For example, a confined aquifer system can be decomposed into two parts, one of which is a linear time invariability sub-system and another is a zero-input response. This system does not meet the additivity property of linear system, but satisfies the incremental linear system characteristics. In order to better understand, a case study of water resources management of Huaibei city within semiarid region, north Anhui province of China, is cited. Taking into account the water demand for satisfying the urban development in the next 15 years, three planning target years of water resources are preset as the present (2005), the short term (2010) and the long term (2020), respectively, and four hydrological years (e.g., wet year, mean year, dry year, and extremely dry year) are also defined by the rainfall data of many years. A groundwater management model based on linear programming is established. This model can deal with 12 possible scenarios (3 target years × 4 hydrological years), optimize the strategies of water resources development, integrate various kinds of water sources (e.g., groundwater, surface–water and additional water) and meet the water demand for the urban development of Huaibei city. In accordance with the groundwater management model solutions, the problem of groundwater drawdown funnels (groundwater overextraction funnels) which formed within the Huaibei downtown area for many years and led to some environmental and social issues will be solved over the whole planning period. More importantly, through statistically analyzing the model solutions, the relationships between the groundwater pumping (input signals) and groundwater level recovery (output signals) show up the characteristics of the incremental linear system.