Effects of land cover change on streamflow in the interior Columbia River Basin (USA and Canada)

An analysis of the hydrological effects of vegetation changes in the Columbia River basin over the last century was performed using two land cover scenarios. The first was a reconstruction of historical land cover vegetation, c. 1900, as estimated by the federal Interior Columbia Basin Ecosystem Management Project (ICBEMP). The second was current land cover as estimated from remote sensing data for 1990. Simulations were performed using the variable infiltration capacity (VIC) hydrological model, applied at one‐quarter degree spatial resolution (approximately 500 km² grid cell area) using hydrometeorological data for a 10 year period starting in 1979, and the 1900 and current vegetation scenarios. The model represents surface hydrological fluxes and state variables, including snow accumulation and ablation, evapotranspiration, soil moisture and runoff production. Simulated daily hydrographs of naturalized streamflow (reservoir effects removed) were aggregated to monthly totals and compared for nine selected sub‐basins. The results show that, hydrologically, the most important vegetation‐related change has been a general tendency towards decreased vegetation maturity in the forested areas of the basin. This general trend represents a balance between the effects of logging and fire suppression. In those areas where forest maturity has been reduced as a result of logging, wintertime maximum snow accumulations, and hence snow available for runoff during the spring melt season, have tended to increase, and evapotranspiration has decreased. The reverse has occurred in areas where fire suppression has tended to increase vegetation maturity, although the logging effect appears to dominate for most of the sub‐basins evaluated. Predicted streamflow changes were largest in the Mica and Corralin sub‐basins in the northern and eastern headwaters region; in the Priest Rapids sub‐basin, which drains the east slopes of the Cascade Mountains; and in the Ice Harbor sub‐basin, which receives flows primarily from the Salmon and Clearwater Rivers of Idaho and western Montana. For these sub‐basins, annual average increases in runoff ranged from 4·2 to 10·7% and decreases in evapotranspiration ranged from 3·1 to 12·1%. In comparison with previous studies of individual, smaller sized watersheds, the modelling approach used in this study provides predictions of hydrological fluxes that are spatially continuous throughout the interior Columbia River basin. It thus provides a broad‐scale framework for assessing the vulnerability of watersheds to altered streamflow regimes attributable to changes in land cover that occur over large geographical areas and long time‐frames. Copyright © 2000 John Wiley & Sons, Ltd.     

Simulated water budget of a small forested watershed in the continental/maritime hydroclimatic region of the United States

Annual streamflows have decreased across mountain watersheds in the Pacific Northwest of the United States over the last ~70 years; however, in some watersheds, observed annual flows have increased. Physically based models are useful tools to reveal the combined effects of climate and vegetation on long‐term water balances by explicitly simulating the internal watershed hydrological fluxes that affect discharge. We used the physically based Simultaneous Heat and Water (SHAW) model to simulate the inter‐annual hydrological dynamics of a 4 km² watershed in northern Idaho. The model simulates seasonal and annual water balance components including evaporation, transpiration, storage changes, deep drainage, and trends in streamflow. Independent measurements were used to parameterize the model, including forest transpiration, stomatal feedback to vapour pressure, forest properties (height, leaf area index, and biomass), soil properties, soil moisture, snow depth, and snow water equivalent. No calibrations were applied to fit the simulated streamflow to observations. The model reasonably simulated the annual runoff variations during the evaluation period from water year 2004 to 2009, which verified the ability of SHAW to simulate the water budget in this small watershed. The simulations indicated that inter‐annual variations in streamflow were driven by variations in precipitation and soil water storage. One key parameterization issue was leaf area index, which strongly influenced interception across the catchment. This approach appears promising to help elucidate the mechanisms responsible for hydrological trends and variations resulting from climate and vegetation changes on small watersheds in the region. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparative analysis of relationships between NLDAS‐2 forcings and model outputs

This study reports results from an analysis of the relationship between atmospheric forcing and model‐simulated water and energy fluxes for the North American Land Data Assimilation System Project Phase 2 (NLDAS‐2). The relationships between mean monthly precipitation and total runoff are stronger in the Sacramento (SAC) and variable infiltration capacity (VIC) models, which grew out of the hydrological community, than in the Noah and Mosaic models, which grew out of the soil‐vegetation‐atmosphere transfer (SVAT) community. The reverse is true for the relationship between mean monthly precipitation and evapotranspiration. In addition, surface energy fluxes in VIC are less sensitive to model forcing (except for air temperature) than those in the Noah and Mosaic model. Notwithstanding these general conclusions, the relationships between forcings and model‐simulated water and energy fluxes for all models vary for different seasons, variables, and regions. These findings will ultimately inspire a combination of SVAT‐type model energy components with hydrological model water components to develop a SVAT‐hydrology model to improve both evapotranspiration and total runoff simulations. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of soil moisture aggregation on surface evaporative fluxes

The effects of small-scale heterogeneity in land surface characteristics on the large-scale fluxes of water and energy in the land-atmosphere system has become a central focus of many climatology research experiments. The acquisition of high resolution land surface data through remote sensing and intensive land-climatology field experiments (like HAPEX, FIFE, and BOREAS) has provided data to investigate the interactions between microscale land-atmosphere interactions and macroscale models. To determine the effect of small scale heterogeneities, the spatially averaged evaporative fraction is analytically derived for spatially variable soil moisture and soil-atmospheric controls on evaporation at low soil moisture. This average evaporative fraction is compared with the evaporative fraction determined using the spatially averaged soil moisture, as if from a lumped, or aggregated, land surface model. Results show that the lumped-model based evaporation will over estimate evaporation during periods of low atmospheric demands (early morning/late afternoon, Winter periods, etc.) and under estimate evaporation during periods of high demand (midday Summer periods.) The accuracy of using ‘effective’ parameters in lumped macroscale models depends on the variability of soil moisture and the sensitivity of the soil-vegetation system to low soil moisture.     

Appropriate vertical discretization of Richards' equation for two‐dimensional watershed‐scale modelling

A number of watershed‐scale hydrological models include Richards' equation (RE) solutions, but the literature is sparse on information as to the appropriate application of RE at the watershed scale. In most published applications of RE in distributed watershed‐scale hydrological modelling, coarse vertical resolutions are used to decrease the computational burden. Compared to point‐ or field‐scale studies, application at the watershed scale is complicated by diverse runoff production mechanisms, groundwater effects on runoff production, runon phenomena and heterogeneous watershed characteristics. An essential element of the numerical solution of RE is that the solution converges as the spatial resolution increases. Spatial convergence studies can be used to identify the proper resolution that accurately describes the solution with maximum computational efficiency, when using physically realistic parameter values. In this study, spatial convergence studies are conducted using the two‐dimensional, distributed‐parameter, gridded surface subsurface hydrological analysis (GSSHA) model, which solves RE to simulate vadose zone fluxes. Tests to determine if the required discretization is strongly a function of dominant runoff production mechanism are conducted using data from two very different watersheds, the Hortonian Goodwin Creek Experimental Watershed and the non‐Hortonian Muddy Brook watershed. Total infiltration, stream flow and evapotranspiration for the entire simulation period are used to compute comparison statistics. The influences of upper and lower boundary conditions on the solution accuracy are also explored. Results indicate that to simulate hydrological fluxes accurately at both watersheds small vertical cell sizes, of the order of 1 cm, are required near the soil surface, but not throughout the soil column. The appropriate choice of approximations for calculating the near soil‐surface unsaturated hydraulic conductivity can yield modest increases in the required cell size. Results for both watersheds are quite similar, even though the soils and runoff production mechanisms differ greatly between the two catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

A remote sensing model for monitoring soil evaporation based on differential thermal inertia and its validation

To develop geosciences quantification and multi-dimensional researches will be an inevitable trend in the 21st century. The interaction between the land surface and the atmosphere not only serves as an important component in geosciences quantification, bu…     

Comparative hydrology: analysis of a semiarid and a humid tropical watershed

This paper analyses measured data from two small tropical watersheds: one in a semiarid (Aiuaba, Brazil, 12·0 km², 5 years of measurements) and another in a humid environment (Jaruco, Cuba, 43·5 km², 21 years of measurements). The watersheds are similar with respect to catchment area (tens of km²), potential evaporation (2·1–2·6 m year^?1), temperature (22–30 °C) and relief (mild hillslope steepness); but show considerable hydrological discrepancies: average precipitation in the humid watershed is two times higher; average river discharge (mm year^?1) is five times higher; and surface water availability (mm year^?1) is 14 times higher than in the semiarid watershed. Long‐term operation of hypothetical surface reservoirs in both basins is simulated. The analysis shows that 73% of the average river discharge are available (with 90% annual reliability) in the humid watershed, against only 28% in the semiarid. The main cause of this difference is the excess evaporation, which consumes 55% of the stored water in the semiarid reservoir, but only 12% in the humid one. The research concludes that: (1) although precipitation indicators are higher in the humid area, they are of the same order of magnitude as in the semiarid; and (2) fluvial‐regime and water‐availability variables are more than one order of magnitude higher in the humid basin, which shows a multiplication effect of these hydrological processes. Such major hydrological differences, despite the similarities between the two tropical watersheds, show the importance of further investigations in the field of comparative hydrology. Copyright © 2009 John Wiley & Sons, Ltd.     

Calibration of paired watersheds: Utility of moving sums in presence of externalities

Historically, paired watershed studies have been used to quantify the hydrological effects of land use and management practices by concurrently monitoring 2 similar watersheds during calibration (pretreatment) and post‐treatment periods. This study characterizes seasonal water table and flow response to rainfall during the calibration period and tests a change detection technique of moving sums of recursive residuals (MOSUM) to select calibration periods for each control–treatment watershed pair when the regression coefficients for daily water table elevation were most stable to minimize regression model uncertainty. The control and treatment watersheds were 1 watershed of 3–4‐year‐old intensely managed loblolly pine (Pinus taeda L.) with natural understory, 1 watershed of 3–4‐year‐old loblolly pine intercropped with switchgrass (Panicum virgatum), 1 watershed of 14–15‐year‐old thinned loblolly pine with natural understory (control), and 1 watershed of switchgrass only. The study period spanned from 2009 to 2012. Silvicultural operational practices during this period acted as external factors, potentially shifting hydrologic calibration relationships between control and treatment watersheds. MOSUM results indicated significant changes in regression parameters due to silvicultural operations and were used to identify stable relationships for water table elevation. None of the calibration relationships developed using this method were significantly different from the classical calibration relationship based on published historical data. We attribute that to the similarity of historical and 2010–2012 leaf area index on control and treatment watersheds as moderated by the emergent vegetation. Although the MOSUM approach does not eliminate the need for true calibration data or replace the classic paired watershed approach, our results show that it may be an effective alternative approach when true data are unavailable, as it minimizes the impacts of external disturbances other than the treatment of interest.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Development of a three‐dimensional watershed modelling system for water cycle in the middle part of the Heihe rivershed,in the west of China

A three‐dimensional numerical modelling system is developed to study transformation processes of water resources in alluvial fan and river basin along the middle reaches of the Heihe River Basin, Northwest China, an arid and semi‐arid region. Integrating land utilization, remote sensing and geographic information systems, we have developed a numerical modelling system that can be used to quantify the effects of land use and anthropogenic activities on the groundwater system as well as to investigate the interaction between surface water and groundwater. Various hydraulic measurements are used to identify and calibrate the hydraulic boundary conditions and spatial distributions of hydraulic parameters. In the modelling study, various water exchanges and human effects on the watershed system are considered. These include water exchange between surface water and groundwater, groundwater pumping, lateral water recharges from mountain areas, land utilization, and infiltration and evaporation in the irrigation and non‐irrigation areas. The modelling system provides a quantitative method to describe spatial and temporal distributions and transformations between various water resources, and it has application to other watersheds in arid and semi‐arid areas. Copyright © 2011 John Wiley & Sons, Ltd.     

Effect of watershed subdivision on simulation of water balance components

The present effect of watershed subdivision on simulated water balance components using the thoroughly tested Soil and Water Assessment Tool (SWAT) model has been evaluated for the Nagwan watershed in eastern India. Observed meteorological and hydrological data (daily rainfall, temperature, relative humidity and runoff) for the years 1995 to 1998 were collected and used. The watershed and sub‐watershed boundaries, slope and soil texture maps were generated using a geographical information system. A supervised classification method was used for land‐use/cover classification from satellite imagery of 1996. In order to study the effect of watershed subdivision, the watershed was spatially defined into three decomposition schemes, namely a single watershed, and 12 and 22 sub‐watersheds. The simulation using the SWAT model was done for a period of 4 years (1995 to 1998). Results of the study showed a perfect water balance for the Nagwan watershed under all of the decomposition schemes. Results also revealed that the number and size of sub‐watersheds do not appreciably affect surface runoff. Except for runoff, there was a marked variation in the individual components of the water balance under the three decomposition schemes. Though the runoff component of the water balance showed negligible variation among the three cases, variations were noticed in the other components: evapotranspiration (5 to 48%), percolation (2 to 26%) and soil water content (0·30 to 22%). Thus, based on this study, it is concluded that watershed subdivision has a significant effect on the water balance components. Copyright © 2005 John Wiley & Sons, Ltd.     

Implementation of an automatic calibration procedure for HYDROTEL based on prior OAT sensitivity and complementary identifiability analysis

Efficiency of hydrological models mostly depends on the quality of the calibration performed prior to use. In this paper, an automatic calibration framework for the distributed hydrological model HYDROTEL is proposed. The calibration procedure was performed for three watersheds characterized with different hydroclimatological conditions: the Sassandra located in Ivory Coast, Africa, and the Montmorency and Beaurivage watersheds located in Quebec (Canada). Results of one‐a‐time (OAT) sensitivity analysis showed that the order of the most sensitive parameters differs for each watershed. Thus, the sensitivity depends on the hydroclimatic and physiographic characteristics of the watersheds. Co‐linearity indices showed that all model parameters were identifiable, that is, none of the studied parameters could be explained by a combination of the other parameters. Following these findings, an automatic calibration was run. Results indicated there was good agreement between simulated and measured streamflows at the outlet of each watershed; Nash–Sutcliffe efficiency (NSE) ranging between 0.77 and 0.92 and R² ranging from 0.87 to 0.97. When comparing NSE and R² values obtained using a process‐oriented, multiple‐objective, manual calibration strategy, a slight increase in model efficiency was reached with the automatic calibration procedure (4.15% for NSE and 2.95% for R²) improving predictions of peak flows for the Montmorency and Beaurivage watersheds (temperate climate conditions) and flows beyond the rainfall season in the Sassandra watershed. The proposed automatic calibration procedure introduced in this paper may be applied to other distributed hydrological model. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of urbanization on stream hydrology in hillslope watersheds in central Texas

This study examined the effect of urbanization on stream hydrology in hillslope watersheds. Ten streams (seven in hillslope and three in gentle slope watersheds) around Austin, Texas were selected for analysis. For each stream, we compared parameters of transfer function (TF) models estimated from daily rainfall and streamflow data collected in two study periods (October 1988–September 1992 and October 2004–September 2008) representing different degrees of watershed urbanization. As expected, the streams became more intermittent as the watersheds were more urbanized in all the study streams. However, the effect of urbanization on peakflow differs between hillslope and gentle slope watersheds. After watershed urbanization, peakflow increased in gentle slope watersheds, but decreased in hillslope watersheds. Based on the results of the TF models, we found that urbanization made stream not flashier but drier in hillslope watersheds. Overpumpage of aquifer has been recognized as a problem that leads to the stream dryness in the study area. However, the overpumpage alone cannot explain the differences in hydrological changes between the two types of watersheds. We attributed the reduced peakflow and stream dryness in the hillslope watersheds to land grading for construction forming stair‐stepped or terraced landscape. Compared with natural hillslope, a stair‐stepped landscape could infiltrate more stormwater by slowing down surface runoff on tread portions of the stair. Our findings suggest that a watershed management scheme should take into account local hydrogeologic conditions to mitigate the stream dryness resulting from urbanization in hillslope watersheds. Copyright © 2010 John Wiley & Sons, Ltd.     

How much does dry‐season fog matter? Quantifying fog contributions to water balance in a coastal California watershed

The seasonally‐dry climate of Northern California imposes significant water stress on ecosystems and water resources during the dry summer months. Frequently during summer, the only water inputs occur as non‐rainfall water, in the form of fog and dew. However, due to spatially heterogeneous fog interaction within a watershed, estimating fog water fluxes to understand watershed‐scale hydrologic effects remains challenging. In this study, we characterized the role of coastal fog, a dominant feature of Northern Californian coastal ecosystems, in a San Francisco Peninsula watershed. To monitor fog occurrence, intensity, and spatial extent, we focused on the mechanisms through which fog can affect the water balance: throughfall following canopy interception of fog, soil moisture, streamflow, and meteorological variables. A stratified sampling design was used to capture the watershed's spatial heterogeneities in relation to fog events. We developed a novel spatial averaging scheme to upscale local observations of throughfall inputs and evapotranspiration suppression and make watershed‐scale estimates of fog water fluxes. Inputs from fog water throughfall (10–30 mm/year) and fog suppression of evapotranspiration (125 mm/year) reduced dry‐season water deficits by 25% at watershed scales. Evapotranspiration suppression was much more important for this reduction in water deficit than were direct inputs of fog water. The new upscaling scheme was analyzed to explore the sensitivity of its results to the methodology (data type and interpolation method) employed. This evaluation suggests that our combination of sensors and remote sensing allows an improved incorporation of spatially‐averaged fog fluxes into the water balance than traditional interpolation approaches.     

Charred forests accelerate snow albedo decay: parameterizing the post‐fire radiative forcing on snow for three years following fire

As large, high‐severity forest fires increase and snowpacks become more vulnerable to climate change across the western USA, it is important to understand post‐fire disturbance impacts on snow hydrology. Here, we examine, quantify, parameterize, model, and assess the post‐fire radiative forcing effects on snow to improve hydrologic modelling of snow‐dominated watersheds having experienced severe forest fires. Following a 2011 high‐severity forest fire in the Oregon Cascades, we measured snow albedo, monitored snow, and micrometeorological conditions, sampled snow surface debris, and modelled snowpack energy and mass balance in adjacent burned forest (BF) and unburned forest sites. For three winters following the fire, charred debris in the BF reduced snow albedo, accelerated snow albedo decay, and increased snowmelt rates thereby advancing the date of snow disappearance compared with the unburned forest. We demonstrate a new parameterization of post‐fire snow albedo as a function of days‐since‐snowfall and net snowpack energy balance using an empirically based exponential decay function. Incorporating our new post‐fire snow albedo decay parameterization in a spatially distributed energy and mass balance snow model, we show significantly improved predictions of snow cover duration and spatial variability of snow water equivalent across the BF, particularly during the late snowmelt period. Field measurements, snow model results, and remote sensing data demonstrate that charred forests increase the radiative forcing to snow and advance the timing of snow disappearance for several years following fire. Copyright © 2016 John Wiley & Sons, Ltd.     

Optical remote mapping of rivers at sub‐meter resolutions and watershed extents

At watershed extents, our understanding of river form, process and function is largely based on locally intensive mapping of river reaches, or on spatially extensive but low density data scattered throughout a watershed (e.g. cross sections). The net effect has been to characterize streams as discontinuous systems. Recent advances in optical remote sensing of rivers indicate that it should now be possible to generate accurate and continuous maps of in‐stream habitats, depths, algae, wood, stream power and other features at sub‐meter resolutions across entire watersheds so long as the water is clear and the aerial view is unobstructed. Such maps would transform river science and management by providing improved data, better models and explanation, and enhanced applications. Obstacles to achieving this vision include variations in optics associated with shadows, water clarity, variable substrates and target–sun angle geometry. Logistical obstacles are primarily due to the reliance of existing ground validation procedures on time‐of‐flight field measurements, which are impossible to accomplish at watershed extents, particularly in large and difficult to access river basins. Philosophical issues must also be addressed that relate to the expectations around accuracy assessment, the need for and utility of physically based models to evaluate remote sensing results and the ethics of revealing information about river resources at fine spatial resolutions. Despite these obstacles and issues, catchment extent remote river mapping is now feasible, as is demonstrated by a proof‐of‐concept example for the Nueces River, Texas, and examples of how different image types (radar, lidar, thermal) could be merged with optical imagery. The greatest obstacle to development and implementation of more remote sensing, catchment scale ‘river observatories’ is the absence of broadly based funding initiatives to support collaborative research by multiple investigators in different river settings. Copyright © 2007 John Wiley & Sons, Ltd.     

LAND-USE IMPACT ON WATERSHED RESPONSE: THE INTEGRATION OF TWO-DIMENSIONAL HYDROLOGICAL MODELLING AND GEOGRAPHICAL INFORMATION SYSTEMS

The integration of a two-dimensional, raster-based rainfall–runoff model, CASC2D, with a raster geographical information system (GIS), GRASS, offers enhanced capabilities for analysing the hydrological impact under a variety of land management scenarios. The spatially varied components of the watershed, such as slope, soil texture, surface roughness and land-use disturbance, were characterized in GRASS at a user-specified grid cell resolution for input into the CASC2D model. CASC2D is a raster-based, single-event rainfall–runoff model that divides the watershed into grid cell elements and simulates the hydrological processes of infiltration, overland flow and channel flow in response to distributed rainfall precipitation. The five-step integration of CASC2D and GRASS demonstrates the potential for analysing spatially and temporally varied hydrological processes within a 50 square mile semi-arid watershed. By defining possible land-use disturbance scenarios for the watershed, a variety of rainfall–runoff events were simulated to determine the changes in watershed response under varying disturbance and rainfall conditions. Additionally, spatially distributed infiltration outputs derived from the simulations were analysed in GRASS to determine the variability of hydrological change within the watershed. Grid cell computational capabilities in GRASS allow the user to combine the scenario simulation outputs with other distributed watershed parameters to develop complex maps depicting potential areas of hydrological sensitivity. This GIS–hydrological model integration provides valuable spatial information to researchers and managers concerned with the study and effects of land-use on hydrological response.     

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.     

利用卫星热红外遥感技术监测长白山天池火山活动性

介绍了卫星热红外遥感技术用于火山监测的国内外研究现状,阐述了热红外遥感技术的原理,分析了卫星热红外遥感技术用于火山活动性监测的可行性。以长白山天池火山为例,基于Landsat TM/ETM影像和ASTER影像反演获得了1999—2008年的温度场,并选取了其中的3种地面覆盖类型(森林植被、土壤和植被(矮草)以及裸露岩石),从而去除了地表环境因素的影响;从每种地面覆盖类型中扣除了当日天池气象站的平均气温,去除了气象因素的影响,得到了由火山热活动可能导致的温度热异常。结果显示,从1999—2005年,由火山活动导致的温度热异常伴随着扰动发生了明显的上升,自2005年以后逐渐下降,2006—2008年趋于稳定。这些结果与测震、GPS形变以及He同位素比值变化趋势保持了较好的一致性,表明卫星热红外遥感技术用于火山活动性监测的巨大潜力和优势,可以作为一种常规的监测手段尝试性地纳入日常的火山监测工作中