Estimation of potential evapotranspiration in the mountainous Panama Canal watershed

Spatially distributed hydrometeorological and plant information within the mountainous tropical Panama Canal watershed is used to estimate parameters of the Penman–Monteith evapotranspiration formulation. Hydrometeorological data from a few surface climate stations located at low elevations in the watershed are complemented by (a) typical wet‐ and dry‐season fields of temperature, wind, water vapour and pressure produced by a mesoscale atmospheric model with a 3 × 3 km² spatial and hourly temporal resolution, and (b) leaf area index fields estimated over the watershed during a few years using satellite data with two different spatial and temporal resolutions. The mesoscale model estimates of spatially distributed surface hydrometeorological variables provide the basis for the extrapolation of the surface climate station data to produce input for the Penman–Monteith equation. The satellite information and existing digital spatial databases of land use and land cover form the basis for the estimation of Penman–Monteith spatially distributed parameter values. Spatially distributed 3 × 3 km² potential evapotranspiration estimates are obtained for the 3300 km² Panama Canal watershed. Estimates for Gatun Lake within the watershed are found to reproduce well the monthly and annual lake evaporation obtained from submerged pans. Sensitivity analysis results of potential evapotranspiration estimates with respect to cloud cover, dew formation, leaf area index distribution and mesoscale model estimates of surface climate are presented and discussed. The main conclusion is that even the limited spatially distributed hydrometeorological and plant information used in this study contributes significantly toward explaining the substantial spatial variability of potential evapotranspiration in the watershed. These results also allow the determination of key locations within the watershed where additional surface stations may be profitably placed. Copyright © 2006 John Wiley & Sons, Ltd.     

Distributed modelling of future changes in hydrological processes of Spencer Creek watershed

The hydrologic impact of climate change has been largely assessed using mostly conceptual hydrologic models. This study investigates the use of distributed hydrologic model for the assessment of the climate change impact for the Spencer Creek watershed in Southern Ontario (Canada). A coupled MIKE SHE/MIKE 11 hydrologic model is developed to represent the complex hydrologic conditions in the Spencer Creek watershed, and later to simulate climate change impact using Canadian global climate model (CGCM 3·1) simulations. Owing to the coarse resolution of GCM data (daily GCM outputs), statistical downscaling techniques are used to generate higher resolution data (daily precipitation and temperature series). The modelling results show that the coupled model captured the snow storage well and also provided good simulation of evapotranspiration (ET) and groundwater recharge. The simulated streamflows are consistent with the observed flows at different sites within the catchment. Using a conservative climate change scenario, the downscaled GCM scenarios predicted an approximately 14–17% increase in the annual mean precipitation and 2–3 °C increase in annual mean maximum and minimum temperatures for the 2050s (i.e., 2046–2065). When the downscaled GCM scenarios were used in the coupled model, the model predicted a 1–5% annual decrease in snow storage for 2050s, approximately 1–10% increase in annual ET, and a 0·5–6% decrease in the annual groundwater recharge. These results are consistent with the downscaled temperature results. For future streamflows, the coupled model indicated an approximately 10–25% increase in annual streamflows for all sites, which is consistent with the predicted changes in precipitation. Overall, it is shown that distributed hydrologic modelling can provide useful information not only about future changes in streamflow but also changes in other key hydrologic processes such as snow storage, ET, and groundwater recharge, which can be particularly important depending on the climatic region of concern. The study results indicate that the coupled MIKE SHE/MIKE 11 hydrologic model could be a particularly useful tool for understanding the integrated effect of climate change in complex catchment scale hydrology. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparison of pre‐ and post‐remediation water quality,Mineral Creek,Colorado

Pre‐ and post‐remediation data sets are used herein to assess the effectiveness of remedial measures implemented in the headwaters of the Mineral Creek watershed, where contamination from hard rock mining has led to elevated metal concentrations and acidic pH. Collection of pre‐ and post‐remediation data sets generally followed the synoptic mass balance approach, in which numerous stream and inflow locations are sampled for the constituents of interest and estimates of streamflow are determined by tracer dilution. The comparison of pre‐ and post‐remediation data sets is confounded by hydrologic effects and the effects of temporal variation. Hydrologic effects arise due to the relatively wet conditions that preceded the collection of pre‐remediation data, and the relatively dry conditions associated with the post‐remediation data set. This difference leads to a dilution effect in the upper part of the study reach, where pre‐remediation concentrations were diluted by rainfall, and a source area effect in the lower part of the study reach, where a smaller portion of the watershed may have been contributing constituent mass during the drier post‐remediation period. A second confounding factor, temporal variability, violates the steady‐state assumption that underlies the synoptic mass balance approach, leading to false identification of constituent sources and sinks. Despite these complications, remedial actions completed in the Mineral Creek headwaters appear to have led to improvements in stream water quality, as post‐remediation profiles of instream load are consistently lower than the pre‐remediation profiles over the entire study reach for six of the eight constituents considered (aluminium, arsenic, cadmium, copper, iron, and zinc). Concentrations of aluminium, cadmium, copper, lead, and zinc remain above chronic aquatic‐life standards, however, and additional remedial actions may be needed. Future implementations of the synoptic mass balance approach should be preceded by an assessment of temporal variability, and modifications to the synoptic sampling protocol should be made if necessary. Published in 2009 by John Wiley & Sons, Ltd.     

Long-term ecosystem and biogeochemical research in Loch Vale watershed,Rocky Mountain National Park,Colorado

Loch Vale watershed was instrumented in 1983 with initial support from the National Acid Precipitation Assessment Program to ask whether ecosystems of Rocky Mountain National Park (RMNP) were affected by acidic atmospheric deposition. Research and monitoring activities were expanded in 1991 by the U.S. Geological Survey Water, Energy, and Biogeochemical Budgets program to understand the processes, and their interactions, controlling water, energy, and biogeochemical fluxes. With help from many collaborators we have characterized trends and patterns in atmospheric deposition, climate, and hydrology, including glaciers and other ice features. Instead of acidic deposition, we documented high atmospheric inputs of reactive nitrogen (Nr), and have studied the ecological consequences in soils, surface water, and vegetation. Using paleolimnology, we documented the onset of human-caused change to lake primary producers ca. 1950 in response to increased Nr deposition and warming. Our results provided the basis for the Colorado Nitrogen Deposition Reduction Plan, a state policy that aims to reduce Nr emissions to protect resources in RMNP by 2032. Carbon cycle research revealed mountain wetlands now release more carbon than they store, and respiration and methane flux occurs even during winter through deep snow packs. Trend analyses found export of Nr to be closely tied to atmospheric inputs, but can lag in response to drought. Current research explores consequences of the combination of warming, changes in precipitation dynamics, and atmospheric deposition of Nr and dust on stream and lake CO₂ dynamics, lake biology and trophic state, and soil carbon composition. Dramatic increases in park visitors have prompted studies on the effects of recreational use on water quality. New tools such as remote sensing and high frequency instream water quality sensors are being applied to lake and stream studies. Monitoring, combined with experiments, models, and spatial comparisons is an essential foundation for science-based resource management.     

Sources of variability in springwater chemistry in Fool Creek,a high-elevation catchment of the Rocky Mountains,Colorado, USA

Springs are the point of origin for most headwater streams and are important regulators of their chemical composition. We analysed solute concentrations of water emerging from 57 springs within the 3 km² Fool Creek catchment at the Fraser Experimental Forest and considered sources of spatial variation among them and their influence on the chemical composition of downstream water. On average, calcium and acid neutralizing capacity (bicarbonate-ANC) comprised 50 and 90% of the cation and anion charge respectively, in both spring and stream water. Variation in inorganic chemical composition among springs reflected distinct groundwater sources and catchment geology. Springs emerging through glacial deposits in the upper portion of the catchment were the most dilute and similar to snowmelt, whereas lower elevation springs were more concentrated in cations and ANC. Water emerging from a handful of springs in a geologically faulted portion of the catchment were more concentrated than all others and had a predominant effect on downstream ion concentrations. Chemical similarity indicated that these springs were linked along surface and subsurface flowpaths. This survey shows that springwater chemistry is influenced at nested spatial scales including broad geologic conditions, elevational and spatial attributes and isolated local features. Our results highlight the role of overlapping factors on solute export from headwater catchments.     

Approaching four decades of forest watershed research at Upper Penticton Creek,British Columbia: A synthesis

Over the past 35 years, the Upper Penticton Creek (UPC) Watershed Experiment has supported forest hydrology research in south-central British Columbia (BC), Canada. This paper provides a synthesis of research results, highlights the challenges facing UPC and identifies new research directions. Clearcutting approximately 50% of two small, snow-dominated (Dfb Koppen classification) watersheds advanced the timing of snowmelt-generated high flows and decreased late-summer low flows, relative to predictions based on pre-treatment regressions. Changes in high flows did not have a significant effect on stream channels due to low stream power, coarse substrate, and limited riparian disturbance. Changes in summer low flows reduced modelled useable fish habitat by 20%–50%. Evaporation averaged 52% of the annual precipitation in the mature forest, was reduced to 30% in a clearcut, and recovered to 40% and 47% in a 10 and 25 year-old stand, respectively. Groundwater recharge to the bedrock was estimated at 19% of annual precipitation, indicating that, even with the large uncertainty associated with this estimate, deep groundwater should not be ignored in the water balance. Suspended sediment, turbidity, and colour increased post-logging; however, chemical surface water quality did not change. Aquatic community structure changed post-logging; and although this affected the processing of organic matter, the effects on habitat quality were considered minimal. The information gained at UPC has supported provincial policies, management guidelines, forest stewardship plans and watershed risk assessments. The undisturbed control watershed, re-growing treatment watersheds and ongoing long-term hydrometric monitoring continue to provide opportunities for future research addressing issues such as the effects of young forests on streamflow and hydrologic recovery, and the influence of climate change on the hydrologic regime.     

Urban irrigation in the modeling of a semi-arid urban environment: Ballona Creek watershed,Los Angeles,California

ABSTRACT

Ballona Creek watershed in Los Angeles, California provides a unique combination of heterogeneous urban land cover, a semi-arid environment, and a large outdoor water-use flux that presents a challenge for physically-based models. We ran simulations using the Noah Land Surface Model and Parflow-Community Land Model and compared to observations of evapotranspiration (ET), runoff, and land surface temperature (LST) for the entire 11-year study period. Both models were systematically adjusted to test the impact of land cover and urban irrigation on simulation results. Monthly total runoff and ET results are greatly improved when compared to an in-situ stream gauge and meteorological tower data: from 0.64 to 0.81 for the Nash–Sutcliffe efficiency (NSE) for runoff and from a negative NSE to 0.82 for ET. The inclusion of urban irrigation in semi-arid urban environments is found to be vital, but not sufficient, for the accurate simulation of variables in the studied models.     

Fifty-eight years and counting of watershed science at the Caspar Creek Experimental Watersheds in northern California

The Caspar Creek Experimental Watersheds are the site of a long-term paired watershed study in the northern Coast Ranges of California. The watersheds are predominately forested with coast redwood and Douglas-fir. Old-growth forest was logged between 1860 and 1904. Two harvesting experiments have been completed since then and a third experiment is currently underway. Caspar Creek data are split into three phases corresponding to three experiments: Phase 1 (1962–1985) reports on a selection harvest (1971–1973) and initial recovery in the South Fork watershed; Phase 2 (1985–2017) includes clearcut harvesting of ~50% of the North Fork watershed (1985–1992) and recovery; and Phase 3 (2017 onward) corresponds to a second selection harvest in the South Fork watershed with a range of subwatershed harvest intensities (2017–2019) and recovery. All three experiments included harvest-related road-building and relied primarily on measurements of streamflow and sediment delivery from both treated and reference watersheds. Major findings include modest increases in post-harvest peak flows and cumulative flow volumes, post-harvest low flows that initially increased and then decreased 12 to 15 years after harvesting, and the consequences of different yarding techniques and road design on sediment yields. Some of the data for Phase 1 and Phase 2 are available in a USDA Forest Service online archive. The archived data include precipitation, streamflow, suspended sediment concentrations, turbidity, accumulated weir pond sediment volumes, bedload transport rates, water stable isotope data, and geospatial data. Archiving activities are ongoing. Phase 3 data are currently being collected and will be archived after a post-harvest monitoring period.     

Soil pipe collapses in a loess pasture of Goodwin Creek watershed,Mississippi: role of soil properties and past land use

Little is known about the association of soil pipe collapse features with soil properties or land use history. Three loess covered catchments in northern Mississippi, USA were characterized to investigate these relationships. Soil pipe collapses were characterized for their size, type feature and spatial location along with soil properties across the three catchments. Although mapped as the same soil, one of the catchments did not contain pipe collapse features while the other two had 29.4 and 15.4 pipe collapses per hectare. These loess soils contained fragipan layers that are suspected of perching water, thereby initiating the piping processes. Pipe collapses associated with subsurface flow paths were not always consistent with surface topography. The surface layer tended to be non‐erodible while layers below, even the upper fragipan layers, were susceptible to erosion by pipeflow. Soil properties of the lowest fragipan layer were highly variable but tended to prevent further downward erosion of soil pipes and thus formed a lower boundary for gullies. Middle to lower landscape positions in one of the piped catchments contained anthropic soils that were highly erodible. These anthropic soils were previously gullies that were filled‐in in the 1950s when forested areas, assumed to have been established when land was previously converted from crop to forest land, were converted to pasture. Three decades after this land use change from forest to pasture, pipe collapses became evident. In contrast, the adjacent catchment that does not exhibit pipe collapse features experienced severe sheet and rill erosion prior to the 1930s while in cotton production. The surface horizons above the lower fragipan layer were completely removed during this period, thus the top‐soil layer that tends to form a bridge above soil pipes in the more erodible subsoil layers was removed. This study showed that knowledge of soil characteristics or topography alone do not explain the distribution of soil pipe collapses as past land use can play a definitive role. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater monitoring at the watershed scale: An evaluation of recharge and nonpoint source pollutant loading in the Clear Creek Watershed,Iowa

Determining the groundwater contribution of nonpoint source pollution at a watershed scale is a challenging issue. In this study, we utilized a top‐down approach to characterize representative groundwater response units (GRUs) based on land use and landscape position (e.g., upland, sideslope, or floodplain) in the 275‐km² Clear Creek Watershed, Iowa. Groundwater monitoring wells were then established along downslope transects in representative GRUs. This unique combination of top‐down/bottom‐up approaches allowed us to estimate groundwater pollutant loads at the watershed scale with minimal monitoring. For the 2015 study period, results indicated that more groundwater recharge occurred in the floodplain (404 mm) compared to the uplands or sideslopes (281 and 165 mm, respectively), irrespective of land use. Recharge in the floodplains consisted of 37% of the annual precipitation, whereas upland wells averaged 26% and sideslopes averaged 15% of the annual precipitation. Less recharge was found to occur beneath perennial grass compared to row crop and urbanized areas. Baseflow discharge accounted for 69% of the total NO₃‐N exported from the Clear Creek Watershed, with row crop areas contributing approximately 95% of the annual load. Orthophosphorus (OP) yields were approximately 0.72 kg/ha beneath urban and suburban areas, three times higher than those in row crop or perennial areas. Urban and suburban areas accounted for 21.4% of groundwater orthophosphorus and chloride loads in the watershed compared to only 8.5% of the land area. Overall, the groundwater load allocation model for baseflow nutrient discharge to Clear Creek can be used to target future nonpoint source load reduction strategies at the watershed scale. The use of GRUs can pinpoint better areas of concern for controlling nutrient loads.     

Analysis of temporal variation and scaling of hydrological variables based on a numerical model of the Sagehen Creek watershed

Temporal variations of the main hydrological variables over 16 years were systematically investigated based on the results from an integrated hydrological modeling at the Sagehen Creek watershed in northern Sierra Nevada. Temporal scaling of these variables and damping effects of the hydrological system as well as its subsystems, i.e., the land surface, unsaturated zone, and saturated zone, were analyzed with spectral analyses. It was found that the hydrological system may act as a cascade of hierarchical fractal filters which sequentially transfer a non-fractal or less correlated fractal hydrological signal to a more correlated fractal signal. The temporal scaling of simulated infiltration (SI), simulated actual evapotranspiration (SET), simulated recharge (SR), measured baseflow (MBF), measured streamflow (MSF) exist and the temporal autocorrelation of these variables increase as water moves through the system. The degree of the damping effect of the subsystems is different and is strongest in the unsaturated zone compared with that of the land surface and saturated zone. The temporal scaling of the simulated groundwater levels (Sh) also exists and is strongly affected by the river: the temporal autocorrelation of Sh near the river is similar to that of the river stage fluctuations and increases away from the river. There is a break in the temporal scaling of Sh near the river at low frequencies due to the effect of the river. Temporal variations of the simulated soil moisture (Sθ) is more complicated: the value of the scaling exponent (β) for Sθ increases with depth as water moves downwards and its high-frequency fluctuations are damped by the unsaturated zone. The temporal fluctuations of measured precipitation and SI are fractional Gaussian noise, those of SET, SR, MBF, and MSF are fractional Brownian motion (fBm), and those of Sh away from the river are 2nd-order fBm based on the values of β obtained in this study.     

Dynamic evolution of nutrient discharge under stormflow and baseflow conditions in a coastal agricultural watershed in Ishigaki Island,Okinawa, Japan

Excessive terrestrial nutrient loadings adversely impact coral reefs by primarily enhancing growth of macroalgae, potentially leading to a phase‐shift phenomenon. Hydrological processes and other spatial and temporal factors affecting nutrient discharge must be examined to be able to formulate effective measures for reducing nutrient export to adjacent reefs. During storm events and baseflow periods, water samples were obtained from the tropical Todoroki River, which drains an intensively agricultural watershed into Shiraho coral reef. In situ nutrient analyzers were deployed for 6 months to hourly measure dissolved nutrient (NO₃⁻‐N and PO₄³⁻‐P) concentrations. Total phosphorus (TP) and suspended solid concentration (TSS) were increased by higher rainfall intensity (r = 0·94, p < 0·01) and river discharge Q (r = 0·88, p < 0·01). In contrast, NO₃⁻‐N concentration tends to decrease drastically (e.g. from 3 to 1 mg l⁻¹) during flood events. When base flow starts to dominate afterwards, NO₃⁻‐N manifested an increasing trend, but decreases when baseflow discharge becomes low. This counter‐clockwise hysteresis for NO₃⁻‐N highlights the significant influence of groundwater discharge. N delivery can therefore be considered a persistent process compared to sediment and P discharge, which are highly episodic in nature. Based on GIS analysis, nutrient concentration along the Todoroki River was largely affected by the percentage of sugarcane/bare areas and bedrock type. The spatial distribution of N concentration in the river reflects the considerable influence of subsurface geology—higher N levels in limestone‐dominated areas. P concentrations were directly related to the total length of artificial drainage, which enhances sediment transport. The use of high‐resolution monitoring data coupled with GIS‐based spatial analysis therefore enabled the clarification of control factors and the difference in the spatio‐temporal discharge characteristics between N and P. Thus, although erosion‐reduction schemes would reduce P discharge, other approaches (e.g. minimize fertilizer) are needed to reduce N discharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrologic‐response simulations for the North Fork of Caspar Creek: second‐growth,clear‐cut,new‐growth,and cumulative watershed effect scenarios

This study demonstrates that comprehensive hydrologic‐response simulation can be a useful tool for studying cumulative watershed effects. The simulations reported here were conducted with the Integrated Hydrology Model (InHM). The location of the 473 ha study site is the North Fork of the Caspar Creek Experimental Watershed, near Fort Bragg, California. Existing information from a long‐term monitoring programme and new soil‐hydraulic property measurements made for this study were used to parameterize InHM. Long‐term continuous wet‐season simulations were conducted for the North Fork catchments and main stem for second‐growth, clear‐cut and new‐growth scenarios. The simulation results show that the increases and decreases, respectively, for throughfall and potential evapotranspiration related to clear‐cutting had quantifiable impacts on the simulated hydrologic response at both the catchment and watershed scales. Model performance was best for the new‐growth simulation scenarios. To improve upon the simulations reported here would require additional soil‐hydraulic property information from across the study area. Although principally focused on the integrated hydrologic response, the effort reported here demonstrates the potential for characterizing distributed responses with physics‐based simulation. The search for a comprehensive understanding of hydrologic response will require both data‐intensive discovery and concept‐development simulation, from both integrated and distributed perspectives. Copyright © 2013 John Wiley & Sons, Ltd.     

Morphology of Red Creek,Wyoming, an arid-region anastomosing channel system

Red Creek, in the Red Desert area of the Great Divide Basin, Wyoming, is an arid-region anastomosing stream. The narrow, deep, and sinuous main channel is flanked by anastomosing flood channels, or anabranches. Most anabranches are initiated at meander bends. The primary mechanism of anabranch initiation is avulsion during overbank floods. Anabranch enlargement occurs by headward erosion. Anabranches act as distributary channels during floods, when water and sediment from overbank flows are transported to and deposited on the floodplain via the anabranches. During periods of low discharges, the anabranches act as tributaries to the main channel, transporting runoff from the floodplain and surrounding hillslopes to the main channel of Red Creek. Aggradation is occurring in the main channel and on the floodplain throughout the study reach. Infilling of the main channel occurs primarily by lateral accretion, while the floodplain accretes vertically through deposition of overbank sediment from the main channel and anabranches. Infilling of the main channel may cause avulsion of the main channel into an anabranch. The abandoned main channel segment may then fill completely or act as an anabranch. Because lateral migration of channels is inhibited by the high cohesion of the silt and clay channel sediment, periodic avulsion is the primary form of lateral mobility in the system.     

Spatial distribution of pipe collapses in Goodwin Creek Watershed,Mississippi

The internal erosion of soil pipes can induce pipe collapses that affect soil erosion processes and landform evolution. The objective of this study was to determine the spatial distribution of pipe collapses in agricultural fields of Goodwin Creek watershed. Ground survey was carried out to detect pipe collapses, and the location, size and surface elevation was measured with differential GPS. A total of 143 of the 145 pipe collapses were found in cropland, and the density was approximately 0.58 collapses per hectare. The spatial distribution of pipe collapses was not uniform as pipe collapses were concentrated in the flat alluvial plains where the land use was dominated by cropland. One of the four parcels had 90% of the pipe collapses with a density of 7.7 collapses per hectare. The mean depth, area and volume of these pipe collapses were 0.12 m, 0.34 m² and 0.02 m³, respectively, and all these properties exhibited a skewed distribution. The drainage area–slope gradient equation, which has been widely used for erosion phenomenon prediction, did not represent pipe collapses in this study as the coefficient of determination was <0.01. This is clear evidence that subsurface flow is not represented by surface topographic characteristics. The pipe collapses were found to intercept runoff, thereby reducing the slope length factor by 6% and the drainage area by 7%. Both of these factors can reduce the sheet and rill erosion; however, the increased subsurface flow could enhance ephemeral gully erosion. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.     

The permeability of the Elkhorn fault zone,South Park,Colorado

The purposes of this study are to use both field and modeling approaches to characterize the permeability of a fault and to assess the role of the fault on regional ground water flow. The study subject is the Elkhorn fault, a low-angle reverse fault that brings Precambrian crystalline rocks over the sediments of Colorado's South Park Basin. The fault is hypothesized to act as a low-permeability barrier to flow, restricting interaction between the crystalline aquifer and the basin sediments. To test this hypothesis and to better predict the permeability structure of the fault, we synthesized geologic data to create a geologic model of the fault, conducted aquifer tests to estimate the hydrogeologic properties of the fault zone, and used ground water modeling to test the influence of a range of hydraulic properties for the fault zone on ground water flow in the region. Our study suggests that the fault is a low-permeability feature. Estimated heads are best matched to observations by modeling the fault as a 10-foot-thick interval of low-permeability fault gouge. Steady-state flow models show that much of the flow in the study area is topographically driven near land surface. Flow rates decrease with depth in the aquifers. In the footwall, ground water moves updip in the Michigan-San Isabel syncline to discharge in the South Park Basin. In the hanging wall, ground water moves east to a regional ground water divide. Sensitivity analyses indicate that hydraulic heads are most sensitive to changes in hydraulic conductivity and recharge.