应用EOS/MODIS卫星资料反演与估算西北雨养农业区地表能量和蒸散量

应用定西地区的陆面过程野外观测资料和四次EOS_MODIS卫星资料,对典型的西北半干旱雨养农业区的基本地表特征参数进行反演,在此基础上对各能量通量和蒸散量进行估算,进而分别分析其季节和空间变化.结果表明：应用卫星数据估算的各参量的相对误差基本都在20%以内,其空间和时间变化基本反映了当地的实际情况;区域各通量和日蒸散量分布极不均匀,表现为四季相差比较大的单峰型变化特征,冬季最小、其次是秋季、夏春两季最大;与面积较大的中、低植被覆盖度区域相比,小范围的高植被覆盖度区域的地表净辐射、潜热和日蒸散量相对较高,土壤热通量和显热通量相对较低,并且都在春、夏季节表现得更加显著.低、中植被覆盖度区的各能量通量的季节变化不显著的特征反映了半干旱西北雨养农业区土壤的干土层相对比较厚的特征.     

Evaporation and land surface energy budget at the Salar de Atacama, Northern Chile

Playa systems are driven by evaporation processes, yet the mechanisms by which evaporation occurs through playa salt crusts are still poorly understood. In this study we examine playa evaporation as it relates to land surface energy fluxes, salt crust characteristics, groundwater and climate at the Salar de Atacama, a 3000 km² playa in northern Chile containing a uniquely broad range of salt crust types. Land surface energy budget measurements were taken at eight representative sites on this playa during winter (August 2001) and summer (January 2002) seasons. Measured values of net all-wave radiation were highest at vegetated and rough halite crust sites and lowest over smooth, highly reflective salt crusts. Over most of the Salar de Atacama, net radiation was dissipated by means of soil and sensible heat fluxes. Dry salt crusts tended to heat and cool very quickly, whereas soil heating and cooling occurred more gradually at wetter vegetated sites. Sensible heating was strongly linked to wind patterns, with highest sensible heat fluxes occurring on summer days with strong afternoon winds. Very little energy available at the land surface was used to evaporate water. Eddy covariance measurements could only constrain evaporation rates to within 0.1 mm d⁻¹, and some measured evaporation rates were less than this margin of uncertainty. Evaporation rates ranged from 0.1 to 1.1 mm d⁻¹ in smooth salt crusts around the margin of the salar and from 0.4 to 2.8 mm d⁻¹ in vegetated areas. No evaporation was detected from the rugged halite salt crust that covers the interior of the salar, though the depth to groundwater is less than 1 m in this area. These crusts therefore represent a previously unrecorded end member condition in which the salt crusts form a practically impermeable barrier to evaporation.     

Optical depth,size distribution and flux of dust from Owens Lake,California

A joint United States/Russian/French collaborative experiment was undertaken in March 1993 and March 1996. Projects LODE I and II (Lake Owens Dust Experiments) took place on the anthropogenically desertified playa (dry lakebed) and surrounding regions of Owens Lake, in east-central California. One of the five parts of Project LODE was to determine relationships between optical depth and flux of dust emitted from the dry lake. Project LODE II included subsequent dust plume measurements and size distributions obtained through April 1996, to further refine the flux measurements for distinct mineral aerosol source regions at Owens Lake. Size distributions of dust aerosol were determined and aerosol optical depths were calculated from sunphotometer solar extinction measurements taken downwind in plumes coming from the emissive areas of Owens Lake. This source was visually observed for 10 measured dust storms. The plume mass was calculated to be 1·5 × 10⁹ g using ground-based measurements and ≥1·6 × 10⁹ g from satellite data. Project LODE II results were found to be consistent with LODE I results for the south end of the playa, but flux values were found to be reduced for the northeastern portion of the playa by comparison. Vertical flux values estimated by sunphotometry were found to be consistent with values estimated via a micrometeorological method. Copyright © 1999 John Wiley & Sons, Ltd.     

Sensitivity of model parameterizations for simulated latent heat flux at the snow surface for complex mountain sites

The snowcover energy balance is typically dominated by net radiation and sensible and latent heat fluxes. Validation of the two latter components is rare and often difficult to undertake at complex mountain sites. Latent heat flux, the focus of this paper, is the primary coupling mechanism between the snow surface and the atmosphere. It accounts for the critical exchange of mass (sublimation or condensation), along with the associated snowcover energy loss or gain. Measured and modelled latent heat fluxes at a wind‐exposed and wind‐sheltered site were compared to evaluate variability in model parameters. A well‐tested and well‐validated snowcover energy balance model, Snobal, was selected for this comparison because of previously successful applications of the model at these sites and because of the adjustability of the parameters specific to latent heat transfer within the model. Simulated latent heat flux and snow water equivalent (SWE) were not sensitive to different formulations of the stability profile functions associated with heat transfer calculations. The model parameters of snow surface roughness length and active snow layer thickness were used to improve latent heat flux simulations while retaining accuracy in the simulation of the SWE at an exposed and sheltered study site. Optimal parameters for simulated latent heat flux and SWE were found at the exposed site with a shorter roughness length and thicker active layer, and at the sheltered site with a longer roughness length and thinner active layer. These findings were linked to physical characteristics of the study sites and will allow for adoption into other snow models that use similar parameters. Physical characteristics of wind exposure and cover could also be used to distribute critical parameters in a spatially distributed modelling domain and aid in parameter selection for application to other watersheds where detailed information is not available. Copyright © 2012 John Wiley & Sons, Ltd.     

Modelling the water balance of a mesoscale catchment basin using remotely sensed land cover data

Hydrological modelling of mesoscale catchments is often adversely affected by a lack of adequate information about specific site conditions. In particular, digital land cover data are available from data sets which were acquired on a European or a national scale. These data sets do not only exhibit a restricted spatial resolution but also a differentiation of crops and impervious areas which is not appropriate to the needs of mesoscale hydrological models. In this paper, the impact of remote sensing data on the reliability of a water balance model is investigated and compared to model results determined on the basis of CORINE (Coordination of Information on the Environment) Land Cover as a reference. The aim is to quantify the improved model performance achieved by an enhanced land cover representation and corresponding model modifications. Making use of medium resolution satellite imagery from SPOT, LANDSAT ETM+ and ASTER, detailed information on land cover, especially agricultural crops and impervious surfaces, was extracted over a 5-year period (2000–2004). Crop-specific evapotranspiration coefficients were derived by using remote sensing data to replace grass reference evapotranspiration necessitated by the use of CORINE land cover for rural areas. For regions classified as settlement or industrial areas, degrees of imperviousness were derived. The data were incorporated into the hydrological model GROWA (large-scale water balance model), which uses an empirical approach combining distributed meteorological data with distributed site parameters to calculate the annual runoff components. Using satellite imagery in combination with runoff data from gauging stations for the years 2000–2004, the actual evapotranspiration calculation in GROWA was methodologically extended by including empirical crop coefficients for actual evapotranspiration calculations. While GROWA originally treated agricultural areas as homogeneous, now a consideration and differentiation of the main crops is possible. The accuracy was determined by runoff measurements from gauging stations. Differences in water balances resulting from the use of remote sensing data as opposed to CORINE were analysed in this study using a representative subcatchment. Resulting Nash–Sutcliff model efficiencies improved from 0.372 to 0.775 and indicate that the enhanced model can produce thematically more accurate and spatially more detailed local water balances. However, the proposed model enhancements by satellite imagery have not exhausted the full potential of water balance modelling, for which a higher temporal resolution is required.     

Modelling the spatial–temporal variability of spring snowmelt in an arctic catchment

Arctic spring landscapes are usually characterized by a mosaic of coexisting snow‐covered and bare ground patches. This phenomenon has major implications for hydrological processes, including meltwater production and runoff. Furthermore, as indicated by aircraft observations, it affects land‐surface–atmosphere exchanges, leading to a high degree of variability in surface energy terms during melt. The heterogeneity and related differences when certain parts of the landscape become snow free also affects the length of the growing season and the carbon cycle. Small‐scale variability in arctic snowmelt is addressed here by combining a spatially distributed end‐of‐winter snow cover with simulations of variable snowmelt energy balance factors for the small arctic catchment of Trail Valley Creek (63 km²). Throughout the winter, snow in arctic tundra basins is redistributed by frequent blowing snow events. Areas of above‐ or below‐average end‐of‐winter snow water equivalents were determined from land‐cover classifications, topography, land‐cover‐based snow surveys, and distributed surface wind‐field simulations. Topographic influences on major snowmelt energy balance factors (solar radiation and turbulent fluxes of sensible and latent heat) were modelled on a small‐scale (40 m) basis. A spatially variable complete snowmelt energy balance was subsequently computed and applied to the distributed snow cover, allowing the simulation of the progress of melt throughout the basin. The emerging patterns compared very well visually to snow cover observations from satellite images and aerial photographs. Results show the relative importance of variable end‐of‐winter snow cover, spatially distributed melt energy fluxes, and local advection processes for the development of a patchy snow cover. This illustrates that the consideration of these processes is crucial for an accurate determination of snow‐covered areas, as well as the location, timing, and amount of meltwater release from arctic catchments, and should, therefore, be included in hydrological models. Furthermore, the study shows the need for a subgrid parameterization of these factors in the land surface schemes of larger scale climate models. Copyright © 2005 John Wiley & Sons, Ltd.     

Vegetation physiological parameter setting in the Simple Biosphere model 2(SiB2) for alpine meadows in the upper reaches of Heihe river

Land surface process modeling of high and cold area with vegetation cover has not yielded satisfactory results in previous applications. In this study, land surface energy budget is simulated using a land surface model for the A’rou meadow in the upper-reach area of the Heihe River Basin in the eastern Tibetan Plateau. The model performance is evaluated using the in-situ observations and remotely sensed data. Sensible and soil heat fluxes are overestimated while latent heat flux is underestimated when the default parameter setting is used. By analyzing physical and physiological processes and the sensitivities of key parameters, the inappropriate default setting of optimum growth and inhibition temperatures is identified as an important reason for the bias. The average daytime temperature during the period of fastest vegetation growth(June and July) is adopted as the optimum growth temperature, and the inhibition temperatures were adjusted using the same increment as the optimum temperature based on the temperature acclimation. These adjustments significantly reduced the biases in sensible, latent, and soil heat fluxes.     

SESTET: A spatially explicit stream temperature model based on equilibrium temperature

Stream-water temperature is a key variable controlling chemical, biological, and ecological processes in freshwater environments. Most models focus on a single river cross-section; however, temperature gradients along stretches and tributaries of a river network are crucial to assess ecohydrological features such as aquatic species suitability, growth and feeding rates, or disease transmission. We propose SESTET, a deterministic, spatially explicit stream temperature model for a whole river network, based on water and energy budgets at a reach scale and requiring only commonly available spatially distributed datasets, such as morphology and air temperature, as input. Heat exchange processes at the air–water interface are modelled via the widely used equilibrium temperature concept, whereas the effects of network structure are accounted for through advective heat fluxes. A case study was conducted on the prealpine Wigger river (Switzerland), where water temperatures have been measured in the period 2014–2018 at 11 spatially distributed locations. The results show the advantages of accounting for water and energy budgets at the reach scale for the entire river network, compared with simpler, lumped formulations. Because our approach fundamentally relies on spatially distributed air temperature fields, adequate spatial interpolation techniques that account for the effects of both elevation and thermal inversion in air temperature are key to a successful application of the model. SESTET allows the assessment of the magnitude of the various components of the heat budget at the reach scale and the derivation of reliable estimates of spatial gradients of mean daily stream temperatures for the whole catchment based on a limited number of conveniently located (viz., spanning the largest possible elevation range) measuring stations. Moreover, accounting for mixing processes and advective fluxes through the river network allows one to trust regionalized values of the parameters controlling the relationship between equilibrium and air temperature, a key feature to generalize the model to data-scarce catchments.     

Spatio‐temporal variability of water and energy fluxes – a case study for a mesoscale catchment in pre‐alpine environment

Water and energy fluxes at and between the land surface, the subsurface and the atmosphere are inextricably linked over all spatio‐temporal scales. Our research focuses on the joint analysis of both water and energy fluxes in a pre‐alpine catchment (55 km²) in southern Germany, which is part of the Terrestrial Environmental Observatories (TERENO). We use a novel three‐dimensional, physically based and distributed modelling approach to reproduce both observed streamflow as an integral measure for water fluxes and heat flux and soil temperature measurements at an observation location over a period of 2 years. While heat fluxes are often used for comparison of the simulations of one‐dimensional land surface models, they are rarely used for additional validation of physically based and distributed hydrological modelling approaches. The spatio‐temporal variability of the water and energy balance components and their partitioning for dominant land use types of the study region are investigated. The model shows good performance for simulating daily streamflow (Nash–Sutcliffe efficiency > 0.75). Albeit only streamflow measurements are used for calibration, the simulations of hourly heat fluxes and soil temperatures at the observation site also show a good performance, particularly during summer. A limitation of the model is the simulation of temperature‐driven heat fluxes during winter, when the soil is covered by snow. An analysis of the simulated spatial fields reveals heat flux patterns that reflect the distribution of the land use and soil types of the catchment. The water and energy partitioning is characterized by a strong seasonal cycle and shows clear differences between the selected land use types. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Estimation of spatially distributed surface energy fluxes using remotely‐sensed data for agricultural fields

Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m⁻² and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R² = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m³ ha⁻¹). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Numerical simulation and data assimilation of the water-energy cycle over semiarid northeastern China

The default fractional vegetation cover and terrain height were replaced by the estimated fractional vegetation cover, which was calculated by the Normalized Difference Vegetation Index(NDVI) of Earth Observing System Moderate-Resolution Imaging Spectroradiometer(EOS-MODIS) and the Digital Elevation Model of the Shuttle Radar Topography Mission(SRTM) system. The near-surface meteorological elements over northeastern China were assimilated into the three-dimensional variational data assimilation system(3DVar) module in the Weather Research and Forecasting(WRF) model. The structure and daily variations of air temperature, humidity, wind and energy fields over northeastern China were simulated using the WRF model. Four groups of numerical experiments were performed, and the simulation results were analyzed of latent heat flux, sensible heat flux, and their relationships with changes in the surface energy flux due to soil moisture and precipitation over different surfaces. The simulations were compared with observations of the stations Tongyu, Naiman, Jinzhou, and Miyun from June to August, 2009. The results showed that the WRF model achieves high-quality simulations of the diurnal characteristics of the surface layer temperature, wind direction, net radiation, sensible heat flux, and latent heat flux over semiarid northeastern China in the summer. The simulated near-surface temperature, relative humidity, and wind speed were improved in the data assimilation case(Case 2) compared with control case(Case 1). The simulated sensible heat fluxes and surface heat fluxes were improved by the land surface parameterization case(Case 3) and the combined case(Case 4). The simulated temporal variations in soil moisture over the northeastern arid areas agree well with observations in Case 4, but the simulated precipitation should be improved in the WRF model. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations. The assimilation datasets generated by this work can be applied to research on climate change and environmental monitoring of arid lands, as well as research on the formation and stability of climate over semiarid areas.     

Evapotranspiration from the margin and moist playa of a closed desert valley

Two long term microclimate measurement stations with Bowen ratio capability have been used to study water cycling in a closed desert basin. Microclimate variables including the temperature and vapor pressure gradients were monitored continuously and were used to estimate the Bowen ratio, sensible and latent heat fluxes during 1986 and 1987. Despite having a water table that varied between the surface and 30 cm below the surface, the playa had little evaporation except after rainfall events. The very high osmotic pressure of the soil and salt crust caused most of the absorbed radiation to be partitioned to sensible heat. In contrast, along the margin the thin grass and brush cover transpired water freely, with the latent heat flux exceeding 60% of available energy for much of the season. The higher air temperatures above the playa raised potential evapotranspiration (ET) significantly higher than along the margin throughout the summer. The annual average actual ET of the playa was only 36% of the margin. During the drier summer period (May–October), this ratio decreased to < 28%. Immediately after a rainfall event, evaporation rates of the two sites were similar, but the playa rate was quickly reduced as the osmotic potential increased. During this study, the playa lost < 229 mm of subsurface water to evaporation annually, while > 638 mm were lost from the margin groundwater supply.

The 24-h solar and net radiation correlations were 0.80 and 0.94 for the playa and margin, respectively. The lower correlation for the playa resulted from the wide variation of albedo with surface moisture changes. The annual average albedo values for the playa and margin were 0.64 and 0.46, respectively.     

Modelling runoff using optical satellite remote sensing data in a high mountainous alpine catchment of Italy

Development of hydrological models for seasonal and real-time runoff forecast in rivers of high alpine catchments is useful for management of water resources. The conceptual models for this purpose are based on a temperature index and/or energy budget and can be either lumped or distributed over the catchment area. Remote sensing satellite data are most useful to acquire near real-time geophysical parameters in order to input to the distributed forecasting models. In the present study, integration of optical satellite remote sensing-derived information was made with ground meteorological and hydrological data, and predetermined catchment morphological parameters, to study the feasibility of application of a distributed temperature index snowmelt runoff model to one of the high mountainous catchments in the Italian Alps, known as Cordevole River Basin. Five sets of Landsat Multispectral Scanning System (MSS) and Thematic Mapper (TM) computer-compatible tapes (CCTs) were processed using digital image processing techniques in order to evaluate the snow cover variation quantitatively. Digital elevation model, slope and aspect parameters were developed and used during satellite data processing. The satellite scenes were classified as snow, snow under transition and snow free areas. A second-order polynomial fit has been attempted to approximate the snow depletion and to estimate daily snow cover areal extent for three elevation zones of the catchment separately. Model performance evaluation based on correlation coefficient, Nash–Sutcliffe coefficient and percentage volume deviation indicated very good simulation between measured and computed discharges for the entire snowmelt period. The use of average temperature values computed from the maximum and minimum temperatures into the model was studied and a suitable algorithm was proposed. © 1997 John Wiley & Sons, Ltd.     

Geoneutrinos and the energy budget of the Earth

The total energy loss of the Earth is well constrained by heat flux measurements on land, the plate cooling model for the oceans, and the buoyancy flux of hotspots. It amounts to 46 ± 2 TW. The main sources that balance the total energy loss are the radioactivity of the Earth's crust and mantle, the secular cooling of the Earth's mantle, and the energy loss from the core. Only the crustal radioactivity is well constrained. The uncertainty on each of the other components is larger than the uncertainty of the total heat loss. The mantle energy budget cannot be balanced by adding the best estimates of mantle radioactivity, secular cooling of the mantle, and heat flux from the core. Neutrino observatories in deep underground mines can detect antineutrinos emitted by the radioactivity of U and Th. Provided that the crustal contribution to the geoneutrino flux can be very precisely calculated, it will be possible to put robust constraints on mantle radioactivity and its contribution to the Earth's energy budget. Equally strong constraints could be obtained from a deep ocean observatory without the need of crustal correction. In the future, it may become possible to obtain directional information on the geoneutrino flux and to resolve radial variations in concentration of heat producing elements in the mantle.     

Observational study of land-atmosphere turbulent flux exchange over complex underlying surfaces in urban and suburban areas

Based on observation data from urban observation stations in Nanjing and Suzhou at two heights in the roughness sublayer above the canopy and observation data at three heights in the SORPES station at the Xianlin Campus of Nanjing University in a suburban area, the of land-atmosphere turbulent flux exchange and the energy balance over complex underlying surfaces were analyzed. The results indicated that in the roughness sublayer above the canopy, the nearsurface momentum flux, sensible heat flux, and latent heat flux increase with height, and the observation value of the surface albedo increases with height. However, the observation value of the net radiation decreases with height, thus resulting in a change in the urban surface energy budget with height. At the SORPES station in the Xianlin Campus of Nanjing University located in a hilly area, the momentum flux, sensible heat flux, and latent heat flux of the ground observation field significantly differed from those of the two heights on the tower, while the two heights on the tower were extremely close. These results indicate that the flux observation over the complex underlying surface exhibits adequate local only when it is conducted at a higher altitude above the ground. The turbulent flux observation results at a lower altitude in urban areas are underestimated, while the turbulent flux observation results near the surface produce a large deviation over the underlying hilly complex.     

Process-based snowmelt modeling: does it require more input data than temperature-index modeling?

Modeling snow hydrology for cold regions remains a problematic aspect of many hydro-environmental models. Temperature-index methods are commonly used and are routinely justified under the auspices that process-based models require too many input data. To test this claim, we used a physical, process-based model to simulate snowmelt at four locations across the conterminous US using energy components estimated from measured daily maximum and minimum temperature, i.e. using only the same data required for temperature-index models. The results showed good agreement between observed and predicted snow water equivalents, average R²>0.9. We duplicated the simulations using a simple temperature-index model best fitted to the data and results were poorer, R²<0.8. At one site we applied the process-based model without substantial parameter estimation, and there were no significant (=0.05) differences between these results and those obtained using temperature-estimated parameters, despite relatively poorly predicted specific energy budget components (R²<0.8). These results encourage the use of mechanistic snowmelt modeling approaches in hydrological models, especially in distributed hydrological models for which landscape snow distribution may be controlled by spatially distributed components of the environmental energy budget.     

Evaporation and energy balance of a sub‐Arctic hillslope in northern Finland

This paper presents measurements of the energy balance (radiation, sensible heat flux, evaporation) from a sub‐arctic hillside in northern Finland for a summer season. Comparisons are also made with a nearby wetland site. The hillslope measurements show an equal partition of the radiant energy into sensible and latent heat flux. The evaporative ratio of just over one half was remarkably constant throughout the season, despite very large day‐to‐day and diurnal variations of temperature, humidity deficit and radiation input. This conservative behaviour of the evaporation was caused by a strong rise in effective surface resistance to evaporation with increasing vapour pressure deficit. This suggests a strong physiological control on the evaporation, with stomata closing at times of high evaporative demand. There was no obvious impact of soil‐water stress on the evaporation. However, a comparison with the evaporation measured at a nearby mire site in 1997 suggests that the mire has a significantly lower surface resistance, even when the impact of a significantly lower humidity deficit in the earlier year is taken into account. The measurements are used to test, off‐line, the performance of MOSES (Meteorological Office Surface Exchange Scheme), a simple, but comprehensive, land surface model. The sensitivity of the energy exchanges to the thermal properties of the top soil layer (a surrogate for the upper soil/vegetation layer) is investigated with the use of the model. It is found that the evaporation is insensitive to these properties; they do, however, influence the partition of energy between the sensible heat flux and the ground heat flux (and hence the soil temperatures). It is suggested that the model needs to represent the thermal properties of the canopy more realistically. Copyright © 2002 John Wiley & Sons, Ltd.     

Impacts of spatially and temporally varying snowmelt on subsurface flow in a mountainous watershed: 1. Snowmelt simulation

Abstract

The dominant source of streamflow in many mountainous watersheds is snowmelt recharge through shallow groundwater systems. The hydrological response of these watersheds is controlled by basin structure and spatially distributed snowmelt. The purpose of this series of two papers is to simulate spatially varying snowmelt and groundwater response in a small mountainous watershed. This paper examines the spatially and temporally variable snowmelt to be used as input to the groundwater flow modelling described in the second paper. Snowmelt simulation by the Simultaneous Heat and Water (SHAW) model (a detailed process model of the interrelated heat, water and solute movement through vegetative cover, snow, residue and soil) was validated by applying the model to two years of data at three sites ranging from shallow transient snow cover on a west-facing slope to a deep snow drift on a north-facing slope. The simulated energy balances for several melt periods are presented. Snow depth, density, and the magnitude and timing of snow cover outflow were simulated well for all sites.     

Influence of a rock glacier spring on the stream energy budget and cold‐water refuge in an alpine stream

The thermal regimes of alpine streams remain understudied and have important implications for cold‐water fish habitat, which is expected to decline due to climatic warming. Previous research has focused on the effects of distributed energy fluxes and meltwater from snowpacks and glaciers on the temperature of mountain streams. This study presents the effects of the groundwater spring discharge from an inactive rock glacier containing little ground ice on the temperature of an alpine stream. Rock glaciers are coarse blocky landforms that are ubiquitous in alpine environments and typically exhibit low groundwater discharge temperatures and resilience to climatic warming. Water temperature data indicate that the rock glacier spring cools the stream by an average of 3 °C during July and August and reduces maximum daily temperatures by an average of 5 °C during the peak temperature period of the first two weeks in August, producing a cold‐water refuge downstream of the spring. The distributed stream surface and streambed energy fluxes are calculated for the reach along the toe of the rock glacier, and solar radiation dominates the distributed stream energy budget. The lateral advective heat flux generated by the rock glacier spring is compared to the distributed energy fluxes over the study reach, and the spring advective heat flux is the dominant control on stream temperature at the reach scale. This study highlights the potential for coarse blocky landforms to generate climatically resilient cold‐water refuges in alpine streams.