A LANDSCAPE‐BASED APPROACH TO ESTIMATE RIPARIAN HYDROLOGICAL AND NITRATE REMOVAL FUNCTIONS1

ABSTRACT: This study evaluates a conceptual model developed for riparian zones in Ontario, Canada, that links landscape hydrogeological characteristics to riparian ground water hydrology and nitrate removal efficiency. Data from a range of riparian sites in the United States and Europe suggest that the riparian zone types identified in the model are consistent with patterns of riparian hydrology and nitrate flux and removal in many humid temperate landscapes. These data also support the view that a riparian width of less than 20 m is often sufficient for effective nitrate removal unless riparian sediments are coarse grained or nitrate transport occurs mainly in surface‐fed ground water seeps. This study assesses the possibility of using topographic, soil, surficial geology, and vegetation maps to determine landscape attributes linked by the model to riparian zone hydrological functioning and nitrate removal efficiency. Although mappable data can help in determining broad classes of riparian zones, field visits are necessary to determine non‐mappable riparian attributes such as seeps, organic horizons, and permeable sediment depth in the riparian zone. This research suggests that the conceptual model could be used for landscape management purposes in most temperate landscapes with minor modifications and that the hydrological component of the model could be adapted for contaminants other than nitrate.     

Identifying riparian sinks for watershed nitrate using soil surveys

The capacity of riparian zones to serve as critical control locations for watershed nitrogen flux varies with site characteristics. Without a means to stratify riparian zones into different levels of ground water nitrate removal capacity, this variability will confound spatially explicit source-sink models of watershed nitrate flux and limit efforts to target riparian restoration and management. We examined the capability of SSURGO (1:15 840 Soil Survey Geographic database) map classifications (slope class, geomorphology, and/or hydric soil designation) to identify riparian sites with high capacity for ground water nitrate removal. The study focused on 100 randomly selected riparian locations in a variety of forested and glaciated settings within Rhode Island. Geomorphic settings included till, outwash, and organic/alluvial deposits. We defined riparian zones with "high ground water nitrate removal capacity" as field sites possessing both >10 m of hydric soil width and an absence of ground water surface seeps. SSURGO classification based on a combination of geomorphology and hydric soil status created two functionally distinct sets of riparian sites. More than 75% of riparian sites classified by SSURGO as organic/alluviumhydric or as outwash-hydric had field attributes that suggest a high capacity for ground water nitrate removal. In contrast, >85% of all till sites and nonhydric outwash sites had field characteristics that minimize the capacity for ground water nitrate removal. Comparing the STATSGO and SSURGO databases for a 64000-ha watershed, STATSGO grossly under-represented critical riparian features. We conclude that the SSURGO database can provide modelers and managers with important insights into riparian zone nitrogen removal potential.     

LANDSCAPE ATTRIBUTES AS CONTROLS ON GROITHD WATER NITRATE REMOVAL CAPACITY OF RIPARIAN ZONES1

ABSTRACT: Inherent site factors can generate substantial variation in the ground water nitrate removal capacity of riparian zones. This paper examines research in the glaciated Northeast to relate variability in ground water nitrate removal to site attributes depicted in readily available spatial databases, such as SSUIRGO. Linking site‐specific studies of riparian ground water nitrate removal to spatial data can help target high‐value riparian locations for restoration or protection and improve the modeling of watershed nitrogen flux. Site attributes, such as hydric soil status (soil wetness) and geomorphology, affect the interaction of nitrate‐enriched ground water with portions of the soil ecosystem possessing elevated biogeochemical transformation rates (i.e., biologically active zones). At our riparian sites, high ground water nitrate‐N removal rates were restricted to hydric soils. Geomorphology provided insights into ground water flowpaths. Riparian sites located on outwash and organic/alluvial deposits have high potential for nitrate‐enriched ground water to interact with biologically active zones. In till deposits, ground water nitrate removal capacity may be limited by the high occurrence of surface seeps that markedly reduce the time available for biological transformations to occur within the riparian zone. To fully realize the value of riparian zones for nitrate retention, landscape controls of riparian nitrate removal in different climatic and physiographic regions must be determined and translated into available spatial databases.     

CHALLENGES IN MODELING HYDROLOGIC AND WATER QUALITY PROCESSES IN RIPARIAN ZONES1

This paper presents key challenges in modeling water quality processes of riparian ecosystems: How can the spatial and temporal extent of water and solute mixing in the riparian zone be modeled? What level of model complexity is justified? How can processes at the riparian scale be quantified? How can the impact of riparian ecosystems be determined at the watershed scale? Flexible models need to be introduced that can simulate varying levels of hillslope‐riparian mixing dictated by topography, upland and riparian depths, and moisture conditions. Model simulations need to account for storm event peak flow conditions when upland solute loadings may either bypass or overwhelm the riparian zone. Model complexity should be dictated by the level of detail in measured data. Model algorithms need to be developed using new macro‐scale and meso‐scale experiments that capture process dynamics at the hillslope or landscape scales. Monte Carlo simulations should be an integral part of model simulations and rigorous tests that go beyond simple time series, and point‐output comparisons need to be introduced. The impact of riparian zones on watershed‐scale water quality can be assessed by performing simulations for representative hillsloperiparian scenarios.     

Denitrification potential in relation to lithology in five headwater riparian zones

The influence of riparian zone lithology on nitrate dynamics is poorly understood. We investigated vertical variations in potential denitrification activity in relation to the lithology and stratigraphy of five headwater riparian zones on glacial till and outwash landscapes in southern Ontario, Canada. Conductive coarse sand and gravel layers occurred in four of the five riparian areas. These layers were thin and did not extend to the field-riparian perimeter in some riparian zones, which limited their role as conduits for ground water flow. We found widespread organic-rich layers at depths ranging from 40 to 300 cm that resulted from natural floodplain processes and the burial of surface soils by rapid valley-bottom sedimentation after European settlement. The organic matter content of these layers varied considerably from 2 to 5% (relic channel deposit) to 5 to 21% (buried soils) and 30 to 62% (buried peat). Denitrification potential (DNP) was measured by the acetylene block method in sediment slurries amended with nitrate. The highest DNP rates were usually found in the top 0- to 15-cm surface soil layer in all riparian zones. However, a steep decline in DNP with depth was often absent and high DNP activity occurred in the deep organic-rich layers. Water table variations in 2000-2002 indicated that ground water only interacted frequently with riparian surface soils between late March and May, whereas subsurface organic layers that sustain considerable DNP were below the water table for most of the year. These results suggest that riparian zones with organic deposits at depth may effectively remove nitrate from ground water even when the water table does not interact with organic-rich surface soil horizons.     

Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management1

Vidon, Philippe, Craig Allan, Douglas Burns, Tim P. Duval, Noel Gurwick, Shreeram Inamdar, Richard Lowrance, Judy Okay, Durelle Scott, and Steve Sebestyen, 2010. Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management. Journal of the American Water Resources Association (JAWRA) 46(2):278-298. DOI: 10.1111/j.1752-1688.2010.00420.x Abstract: Biogeochemical and hydrological processes in riparian zones regulate contaminant movement to receiving waters and often mitigate the impact of upland sources of contaminants on water quality. These heterogeneous processes have recently been conceptualized as “hot spots and moments” of retention, degradation, or production. Nevertheless, studies investigating the importance of hot phenomena (spots and moments) in riparian zones have thus far largely focused on nitrogen (N) despite compelling evidence that a variety of elements, chemicals, and particulate contaminant cycles are subject to the influence of both biogeochemical and transport hot spots and moments. In addition to N, this review summarizes current knowledge for phosphorus, organic matter, pesticides, and mercury across riparian zones, identifies variables controlling the occurrence and magnitude of hot phenomena in riparian zones for these contaminants, and discusses the implications for riparian zone management of recognizing the importance of hot phenomena in annual solute budgets at the watershed scale. Examples are presented to show that biogeochemical process-driven hot spots and moments occur along the stream/riparian zone/upland interface for a wide variety of constituents. A basic understanding of the possible co-occurrence of hot spots and moments for a variety of contaminants in riparian systems will increase our understanding of the influence of riparian zones on water quality and guide management strategies to enhance nutrient or pollutant removal at the landscape scale.     

Evaluation of Methods for Delineating Riparian Zones in a Semi‐Arid Montane Watershed

Riparian zones in semi‐arid, mountainous regions provide a disproportionate amount of the available wildlife habitat and ecosystem services. Despite their importance, there is little guidance on the best way to map riparian zones for broad spatial extents (e.g., large watersheds) when detailed maps from field data or high‐resolution imagery and terrain data are not available. Using well‐established accuracy metrics (e.g., kappa, precision, computational complexity), we evaluated eight methods commonly used to map riparian zones. Focusing on a semi‐arid, mountainous watershed, we found that the most accurate and robust method for mapping riparian zones combines data on upstream drainage area and valley topography. That method performed best regardless of stream order, and was most effective when implemented with fine resolution topographic and stream line data. Other commonly used methods to model riparian zones, such as those based on fixed‐width buffers, yielded inaccurate results. We recommend that until very‐high resolution (<1 m) elevation data are available at broad extents, models of riparian zones for semi‐arid mountainous regions should incorporate drainage area, valley topography, and quantify uncertainty.     

GROUND WATER NITRATE REDUCTION IN GIANT CANE AND FOREST RIPARIAN BUFFER ZONES1

ABSTRACT: Ground water contamination by excess nitrate leaching in row‐crop fields is an important issue in intensive agricultural areas of the United States and abroad. Giant cane and forest riparian buffer zones were monitored to determine each cover type's ability to reduce ground water nitrate concentrations. Ground water was sampled at varying distances from the field edge to determine an effective width for maximum nitrate attenuation. Ground water samples were analyzed for nitrate concentrations as well as chloride concentrations, which were used as a conservative ion to assess dilution or concentration effects within the riparian zone. Significant nitrate reductions occurred in both the cane and the forest riparian buffer zones within the first 3.3 m, a relatively narrow width. In this first 3.3 m, the cane and forest buffer reduced ground water nitrate levels by 90 percent and 61 percent, respectively. Approximately 40 percent of the observed 99 percent nitrate reduction over the 10 m cane buffer could be attributed to dilution by upwelling ground water. Neither ground water dilution nor concentration was observed in the forest buffer. The ground water nitrate attenuation capabilities of the cane and forest riparian zones were not statistically different. During the spring, both plant assimilation and denitrification were probably important nitrate loss mechanisms, while in the summer nitrate was more likely lost via denitrification since the water table dropped below the rooting zone.     

Persistent elevated nitrate in a riparian zone aquifer

Streamside vegetated buffer strips (riparian zones) are often assumed to be zones of ground water nitrate (NO3(-)) attenuation. At a site in southwestern Ontario (Zorra site), detailed monitoring revealed that elevated NO3(-) -N (4-93 mg L(-1)) persisted throughout a 100-m-wide riparian floodplain. Typical of riparian zones, the site has a soil zone of recent river alluvium that is organic carbon (OC) rich (36 +/- 16 g kg(-1)). This material is underlain by an older glacial outwash aquifer with a much lower OC content (2.3 +/- 2.5 g kg(-1). Examination of NO3(-), Cl(-), SO4(2-), and dissolved organic carbon (DOC) concentrations; N/Cl ratios; and NO3(-) isotopic composition (delta15N and delta18O) provides evidence of four distinct NO3(-) source zones within the riparian environment. Denitrification occurs but is incomplete and is restricted to a narrow interval located within ~0.5 m of the alluvium-aquifer contact and to one zone (poultry manure compost zone) where elevated DOC persists from the source. In older ground water close to the river discharge point, denitrification remains insufficient to substantially deplete NO3(-). Overall, denitrification related specifically to the riparian environment is limited at this site. The persistence of NO3(-) in the aquifer at this site is a consequence of its Pleistocene age and resulting low OC content, in contrast to recent fluvial sediments in modern agricultural terrain, which, even if permeable, usually have zones enriched in labile OC. Thus, sediment age and origin are additional factors that should be considered when assessing the potential for riparian zone denitrification.     

Minimizing disruptions caused by damming and ensuring the health of downstream ecosystems

In view of the Brazilian Ten‐Year Energy Expansion Plan 2021, this article presents a discussion on environmental flow (e‐flow). The authors analyze the literature to show the evolution of publications concerning e‐flow releases from the perspective of ecosystems services preservation considering results from different case studies from throughout the world. Finally, two main recommendations are drawn regarding e‐flow are: (1) performing a holistic approach to e‐flow planning, including hydrological, hydraulic, water quality, habitat, and riparian zone considerations; and (2) installing in new structures adequate bottom outlets to allow a range of adjustable e‐flow from reservoir dams to reproduce natural flow variations.     

Testing a Simple Field Method for Assessing Nitrate Removal in Riparian Zones1

Abstract: Being able to identify riparian sites that function better for nitrate removal from groundwater is critical to using efficiently the riparian zones for water quality management. For this purpose, managers need a method that is quick, inexpensive, and accurate enough to enable effective management decisions. This study assesses the precision and accuracy of a simple method using three ground water wells and one measurement date for determining nitrate removal characteristics of riparian buffer zones. The method is a scaled‐down version of a complex field research method that consists of a large network of wells and piezometers monitored monthly for over two years. Results using the simplified method were compared to those from the reference research method on a date‐by‐date basis on eight sites covering a wide range of hydrogeomorphic settings. The accuracy of the three‐well, 1 day measurement method was relatively good for assessing nitrate concentration depletion across riparian zones, but poor for assessing the distance necessary to achieve a 90% nitrate removal and for estimating water and nitrate fluxes compared to the reference method. The simplified three‐well method provides relatively better estimates of water and nitrate fluxes on sites where ground‐water flow is parallel to the water table through homogeneous aquifer material, but such conditions may not be geographically widespread. Despite limited overall accuracy, some parameters that are estimated using the simplified method may be useful to water resource managers. Nitrate depletion information may be used to assess the adequacy of existing buffers to achieve nitrate concentration goals for runoff. Estimates of field nitrate runoff and buffer removal fluxes may be adequate for prioritizing management toward sites where riparian buffers are likely to have greater impact on stream water quality.     

Denitrification in a coastal plain riparian zone contiguous to a heavily loaded swine wastewater spray field

Riparian zones are recognized as landscape features that buffer streams from pollutants, particularly nitrogen. The objectives of this experiment were to (i) assess denitrification activity within a riparian zone and (ii) determine the influence of physical, chemical, and landscape features on denitrification. This experiment was conducted from 1994 to 1997 in North Carolina on a riparian zone contiguous to a spray field that was heavily loaded with swine lagoon wastewater. Denitrification enzyme activity (DEA) was measured on soils collected from (i) the soil surface, (ii) midway between the soil surface and water table, and (iii) above the water table. The DEA ranged from 3 to 1660 microg N(2)O-N kg(-1) soil h(-1). The DEA was highest next to the stream and lowest next to the spray field. Nitrate was found to be the limiting factor for denitrification. The DEA generally decreased with soil depth; means for the surface, middle, and bottom depths were 147, 83, and 67 microg N(2)O-N kg(-1) soil h(-1), respectively. These DEA values are higher than those reported for riparian zones adjoining cropland of the southeastern United States, but are lower than those reported for a constructed wetland used for treatment of swine wastewater. Regression analysis indicated that soil total nitrogen was the highest single factor correlated to DEA (r(2) = 0.65). The inclusion of water table depth, soil depth, and distance from the spray field improved the R(2) to 0.86. This riparian zone possessed sufficient soil area with high denitrifying conditions to be a significant factor in the removal of excess nitrogen in the ground water.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

RIPARIAN ZONE CLASSIFICATION FOR MANAGEMENT OF STREAM WATER QUALITY AND ECOSYSTEM HEALTH1

ABSTRACT: Riparian zones perform a variety of biophysical functions that can be managed to reduce the effects of land use on instream habitat and water quality. However, the functions and human uses of riparian zones vary with biophysical factors such as landform, vegetation, and position along the stream continuum. These variations mean that “one size fits all” approaches to riparian management can be ineffective for reducing land use impacts. Thus riparian management planning at the watershed scale requires a framework that can consider spatial differences in riparian functions and human uses We describe a pilot riparian zone classification developed to provide such a framework for riparian management in two diverse river systems in the Waikato region of New Zealand. Ten classes of riparian zones were identified that differed sufficiently in their biophysical features to require different management. Generic “first steps” and “best practical” riparian management recommendations and associated costs were developed for each riparian class. The classification aims to not only improve our understanding of the effectiveness of riparian zone management as a watershed management tool among water managers and land owners, but to also provide a basis for deciding on management actions.     

Seasonal dynamics of denitrification along topohydrosequences in three different riparian wetlands

We investigated the seasonal patterns of denitrification rates and potentials in soil profiles along the topohydrosequence formed at the upland-wetland interface in three riparian wetlands with different vegetation cover (i.e., forest, understory vegetation, and grass). Denitrification was measured using the acetylene inhibition method on soil cores and slurries, which provided a means of comparing the relative activity of this process in different locations. We evaluated, on a seasonal basis, the respective importance of the vegetative cover and the hydromorphic gradient as factors limiting denitrification. Regardless of the season, vegetation type, or lateral position along each topohydrosequence in the riparian wetlands, strong significant gradients of both in situ and potential denitrification rates were measured within a soil profile. Results confirm that the upper organic soil horizon is the most active, when in contact with the ground water. In deeper soil horizons, denitrification activity was low (from 0.004 to 0.5 mg N kg(-1) dry soil d(-1)), but contributed significantly to the reduction of ground water NO3- load along the riparian ground water flowpath (from 9.32 to 0.98 mg NO3-N L(-1)). Along the soil topohydrosequence, the denitrifying community of the upper soil horizons did not vary significantly on a seasonal basis despite the large seasonal ground water fluctuations. Along each topohydrosequence, the denitrification-limiting factor gradually shifted from anaerobiosis to NO3- supply. In situ denitrification rates in the forested, understory vegetation and grass sites were not significantly different. This result emphasizes the importance of the topography of the valley rather than the vegetation cover in controlling denitrification activity in riparian wetlands.     

DIN Retention‐Transport Through Four Hydrologically Connected Zones in a Headwater Catchment of the Upper Mississippi River1

Abstract: Dissolved inorganic nitrogen (DIN) retention‐transport through a headwater catchment was synthesized from studies encompassing four distinct hydrologic zones of the Shingobee River Headwaters near the origin of the Mississippi River. The hydrologic zones included: (1) hillslope ground water (ridge to bankside riparian); (2) alluvial riparian ground water; (3) ground water discharged through subchannel sediments (hyporheic zone); and (4) channel surface water. During subsurface hillslope transport through Zone 1, DIN, primarily nitrate, decreased from ～3 mg‐N/l to <0.1 mg‐N/l. Ambient seasonal nitrate:chloride ratios in hillslope flow paths indicated both dilution and biotic processing caused nitrate loss. Biologically available organic carbon controlled biotic nitrate retention during hillslope transport. In the alluvial riparian zone (Zone 2) biologically available organic carbon controlled nitrate depletion although processing of both ambient and amended nitrate was faster during the summer than winter. In the hyporheic zone (Zone 3) and stream surface water (Zone 4) DIN retention was primarily controlled by temperature. Perfusion core studies using hyporheic sediment indicated sufficient organic carbon in bed sediments to retain ground water DIN via coupled nitrification‐denitrification. Numerical simulations of seasonal hyporheic sediment nitrification‐denitrification rates from perfusion cores adequately predicted surface water ammonium but not nitrate when compared to 5 years of monthly field data (1989‐93). Mass balance studies in stream surface water indicated proportionally higher summer than winter N retention. Watershed DIN retention was effective during summer under the current land use of intermittently grazed pasture. However, more intensive land use such as row crop agriculture would decrease nitrate retention efficiency and increase loads to surface water. Understanding DIN retention capacity throughout the system, including special channel features such as sloughs, wetlands and floodplains that provide surface water‐ground water connectivity, will be required to develop effective nitrate management strategies.     

Explaining landholders' decisions about riparian zone management: the role of behavioural, normative, and control beliefs

Water quality is a key concern in the current global environment, with the need to promote practices that help to protect water quality, such as riparian zone management, being paramount. The present study used the theory of planned behaviour as a framework for understanding how beliefs influence decisions about riparian zone management. Respondents completed a survey that assessed their behavioural, normative, and control beliefs in relation to intentions to manage riparian zones on their property. The results of the study showed that, overall, landholders with strong intentions to manage their riparian zones differed significantly in terms of their beliefs compared to landholders who had weak intentions to manage their riparian zones. Strong intentions to manage riparian zones were associated with a favourable cost-benefit analysis, greater perceptions of normative support for the practice and lower perceptions of the extent to which barriers would impede management of riparian zones. It was also evident that willingness to comply with the recommendations of salient referents, beliefs about the benefits of riparian zone management and perceptions of the extent to which barriers would impede riparian zone management were most important for determining intentions to manage riparian zones. Implications for policy and extension practice are discussed.     

Variations in Base‐Flow Nitrate Flux in a First‐Order Stream and Riparian Zone1

Abstract: Nonpoint source pollution, which contributes to contamination of surface waters, is difficult to control. Some pollutants, particularly nitrate (), are predominantly transmitted through ground water. Riparian buffer zones have the potential to remove contaminants from ground water and reduce the amount of that enters surface water. This is a justification for setting aside vegetated buffer strips along waterways. Many riparian zone hydrologic models assume uniform ground‐water flow through organic‐rich soil under reducing conditions, leading to effective removal of ground‐water prior to discharge into a stream. However, in a small first‐order stream in the mid‐Atlantic coastal plain, base‐flow generation was highly variable (spatially and temporally). Average base‐flow loads were greater in winter than summer, and higher during a wetter year than in dryer years. Specific sections of the stream consistently received greater amounts of high ground water than others. Areas within the riparian zone responsible for most of the exported from the watershed are termed “critical areas.” Over this 5‐year study, most of the exported during base flow originated from a critical area comprising less than 10% of the total riparian zone land area. Allocation of resources to address and improve mitigation function in critical areas should be a priority for continued riparian zone research.     

Early Vegetation Development on an Exposed Reservoir: Implications for Dam Removal

The 4-year drawdown of Horsetooth Reservoir, Colorado, for dam maintenance, provides a case study analog of vegetation response on sediment that might be exposed from removal of a tall dam. Early vegetation recovery on the exposed reservoir bottom was a combination of (1) vegetation colonization on bare, moist substrates typical of riparian zones and reservoir sediment of shallow dams and (2) a shift in moisture status from mesic to the xeric conditions associated with the pre-impoundment upland position of most of the drawdown zone. Plant communities changed rapidly during the first four years of exposure, but were still substantially different from the background upland plant community. Predictions from the recruitment box model about the locations of Populus deltoides subsp. monilifera (plains cottonwood) seedlings relative to the water surface were qualitatively confirmed with respect to optimum locations. However, the extreme vertical range of water surface elevations produced cottonwood seed regeneration well outside the predicted limits of drawdown rate and height above late summer stage. The establishment and survival of cottonwood at high elevations and the differences between the upland plant community and the community that had developed after four years of exposure suggest that vegetation recovery following tall dam removal will follow a trajectory very different from a simple reversal of the response to dam construction, involving not only long time scales of establishment and growth of upland vegetation, but also possibly decades of persistence of legacy vegetation established during the reservoir to upland transition.     

SUBSURFACE RETURN FLOW AND GROUND WATER RECHARGE OF TERRACE FIELDS IN NORTHERN TAIWAN1

ABSTRACT: This study estimates subsurface return flow and effective ground water recharge in terraced fields in northern Taiwan. Specifically, a three dimensional model, FEMWATER, was applied to simulate percolation and lateral seepage in the terraced fields under various conditions. In the terraced paddy fields, percolation mainly moves vertically downward in the central area, while lateral seepage is mainly focused around the bund. Although the simulated lateral seepage rate through the bund exceeded the percolation rate in the central area of the paddy field, annual subsurface return flow at Pei‐Chi and Shin‐Men was 0.17 × 10⁶ m³ and 0.37 × 10⁶ m³, representing only 0.17 percent and 0.21 percent of the total irrigation water required for rice growth at Pei‐Chi and Shin‐Men, respectively. For upland fields, the effective ground water recharge rate during the second crop period (July to November) exceeded that during the first crop period (January to May) because of the wet season in the second crop period. Terraced paddy fields have the most efficient ground water recharge, with 21.2 to 23.4 percent of irrigation water recharging to ground water, whereas upland fields with a plow layer have the least efficient ground water recharge, with only 4.8 to 6.6 percent of irrigation water recharging to ground water. The simulation results clearly revealed that a substantial amount of irrigation water recharges to ground water in the terraced paddy, while only a small amount of subsurface return flow seeps from the upstream to the downstream terraced paddy. The amounts of subsurface flow and ground water recharge determined in the study are useful for the irrigation water planning and management and provide a scientific basis to reevaluate water resources management in the terrace region under irrigated rice.