Evaporation from boreal reservoirs: A comparison between eddy covariance observations and estimates relying on limited data

Hydrological models used for reservoir management typically lack an accurate representation of open-water evaporation and must be run in a scarce data context. This study aims to identify an accurate means to estimate reservoir evaporation with simple meteorological inputs during the open-water season, using long-term eddy covariance observations from two boreal hydropower reservoirs with contrasting morphometry as reference. Unlike the temperate water bodies on which the majority of other studies have focused, northern reservoirs are governed by three distinct periods: ice cover in the cold season, warming in the summer and energy release in the fall. The reservoirs of interest are Eastmain-1 (52°N, mean depth of 11 m) and Romaine-2 (51°N, mean depth of 42 m), both located in eastern Canada. Four approaches are analysed herein: a combination approach, a radiation-based approach, a mass-transfer approach, and empirical methods. Of all the approaches, the bulk transfer equation with a constant Dalton number of 1.2 x 10⁻³ gave the most accurate estimation of evaporation at hourly time steps, compared with the eddy covariance observations (RMSE of 0.06 mm h⁻¹ at Eastmain-1 and RMSE of 0.04 mm h⁻¹ at Romaine-2). The daily series also showed good accuracy (RMSE of 1.38 mm day⁻¹ at Eastmain-1 and RMSE of 0.62 mm day⁻¹ at Romaine-2) both in the warming and energy release phases of the open-water season. The bulk transfer equation, on the other hand, was incapable of reproducing condensation episodes that occurred soon after ice breakup. Basic and variance-based sensitivity analyses were conducted, in particular to measure the variation in performance when the bulk transfer equation was applied with meteorological observations collected at a certain distance (~10–30 km) from the reservoir. This exercise illustrated that accurate estimates of open water evaporation require representative measurements of wind speed and water surface temperature.     

River evaporation and corresponding heat fluxes in forested catchments

River water temperature is a very important variable in ecological studies, especially for the management of fisheries and aquatic resources. Temperature can impact on fish distribution, growth, mortality and community dynamics. River evaporation has been identified as an important heat loss and a key process in the thermal regime of rivers. However, its quantification remains a challenge, mainly because of the difficulty of making direct measurements. The objectives of this study were to characterize the evaporative heat flux at different scales (brook vs river) and to improve the estimation of the evaporative heat flux in a stream temperature model at the hourly timescale. Using a mass balance approach with floating minipans, we measured river evaporation at an hourly timescale in a medium‐sized river (Little Southwest Miramichi) and a small brook (Catamaran Brook) in New Brunswick, Canada. With these direct measurements of evaporation, we developed mass transfer equations to estimate hourly evaporation rates from microclimate conditions measured 2 m above the stream. During the summer 2012, river evaporation was more important for the medium‐sized river with a mean daily evaporation rate of 3.0 mm day^?1 in the Little Southwest Miramichi River compared with that of 1.0 mm day^?1 in Catamaran Brook. Evaporation was the main heat loss mechanism in the two studied streams and was responsible for 42% of heat losses in the Little Southwest Miramichi River and 34% of heat losses in Catamaran Brook during the summer. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessing the significance of wet-canopy evaporation from forests during extreme rainfall events for flood mitigation in mountainous regions of the United Kingdom

There is increased interest in the potential of tree planting to help mitigate flooding using nature-based solutions or natural flood management. However, many publications based upon catchment studies conclude that, as flood magnitude increases, benefit from forest cover declines and is insignificant for extreme flood events. These conclusions conflict with estimates of evaporation loss from forest plot observations of gross rainfall, through fall and stem flow. This study explores data from existing studies to assess the magnitudes of evaporation and attempts to identify the meteorological conditions under which they would be supported. This is achieved using rainfall event data collated from publications and data archives from studies undertaken in temperate environments around the world. The meteorological conditions required to drive the observed evaporation losses are explored theoretically using the Penman–Monteith equation. The results of this theoretical analysis are compared with the prevailing meteorological conditions during large and extreme rainfall events in mountainous regions of the United Kingdom to assess the likely significance of wet canopy evaporation loss. The collated dataset showed that event Ewc losses between approximately 2 and 38% of gross rainfall (1.5 to 39.4 mm day⁻¹) have been observed during large rainfall events (up to 118 mm day⁻¹) and that there are few data for extreme events (>150 mm day⁻¹). Event data greater than 150 mm (reported separately) included similarly high percentage evaporation losses. Theoretical estimates of wet-canopy evaporation indicated that, to reproduce the losses towards the high end of these observations, relative humidity and the aerodynamic resistance for vapour transport needed to be lower than approximately 97.5% and 0.5 to 2 s m⁻¹ respectively. Surface meteorological data during large and extreme rainfall events in the United Kingdom suggest that conditions favourable for high wet-canopy evaporation are not uncommon and indicate that significant evaporation losses during large and extreme events are possible but not for all events and not at all locations. Thus the disparity with the results from catchment studies remains.     

Fuzzy logic model approaches to daily pan evaporation estimation in western Turkey / Estimation de l’évaporation journalière du bac dans l’Ouest de la Turquie par des modèles à base de logique floue

Abstract

Abstract Evaporation is one of the fundamental elements in the hydrological cycle, which affects the yield of river basins, the capacity of reservoirs, the consumptive use of water by crops and the yield of underground supplies. In general, there are two approaches in the evaporation estimation, namely, direct and indirect. The indirect methods such as the Penman and Priestley-Taylor methods are based on meteorological variables, whereas the direct methods include the class A pan evaporation measurement as well as others such as class GGI-3000 pan and class U pan. The major difficulty in using a class A pan for the direct measurements arises because of the subsequent application of coefficients based on the measurements from a small tank to large bodies of open water. Such difficulties can be accommodated by fuzzy logic reasoning and models as alternative approaches to classical evaporation estimation formulations were applied to Lake Egirdir in the western part of Turkey. This study has three objectives: to develop fuzzy models for daily pan evaporation estimation from measured meteorological data, to compare the fuzzy models with the widely-used Penman method, and finally to evaluate the potential of fuzzy models in such applications. Among the measured meteorological variables used to implement the models of daily pan evaporation prediction are the daily observations of air and water temperatures, sunshine hours, solar radiation, air pressure, relative humidity and wind speed. Comparison of the classical and fuzzy logic models shows a better agreement between the fuzzy model estimations and measurements of daily pan evaporation than the Penman method.     

用稳定同位素方法估算大型浅水湖泊蒸发量——以太湖为例

湖泊蒸发量的准确估算对于水文学、气象学和湖泊学等研究有重要的意义.基于2013-2015年太湖水量收支资料、气象观测数据和稳定同位素观测资料,采用稳定同位素质量守恒模型、水量平衡法和Priestley-Taylor模型估算太湖蒸发量,分析太湖蒸发量的季节变化和年际变化特征,并以Priestley-Taylor模型结果为参考值,评价水量平衡法和同位素质量守恒方程的计算精度.结果表明:5-9月太湖蒸发量较高,冬季最低.2013-2015年太湖年总蒸发量分别为1069、894和935 mm,蒸发量的年际变化受到天气条件的影响.2013年12月2014年11月期间,用Priestley-Taylor模型计算的湖泊蒸发量为885 mm;同位素质量守恒模型的估算结果较一致,为893 mm;而水量平衡方程的估算结果明显偏高,为1247 mm.     

Effects of fire and vegetation cover on hydrological characteristics of a Mediterranean shrubland soil

An experimental study based on the effects of fire on soil hydrology was developed at the Experimental Station of ‘La Concordia’ (Valencia, Spain). It is located on a calcareous hillside facing SSE and composed of nine erosion plots (4 × 20 m). In summer 2003, after eight years of soil and vegetation recovery from previous fires in 1995 (with three fire treatments: T1 high‐intensity fire, T2 moderate intensity, and T3 not burnt), experimental fires of low intensity were again conducted on the plots already burnt, to study the effects of repeated fires on the soil water infiltration, soil water content and runoff. Infiltration rates and capacities were measured by the mini‐disk infiltrometer method (MDI), assessing the effects of vegetation cover by comparing the under‐canopy microenvironment (UC) and its absence on bare soil (BS), immediately before and after the fire experiments. Soil properties like water retention capacity (SWRC) and water content (SWC) were also determined for the different fire treatments (T1, T2 and T3) and microsites (UC and BS). Hydrological parameters, such as runoff and infiltration rate, were monitored at plot scale from July 2002 to July 2004. In the post‐fire period, data displayed a 20% runoff increase and a decrease in infiltration (18%). Differences in the steady‐state infiltration rate (SSI) and infiltration capacity (IC) were tested with the MDI on the different treatments (T1, T2 and T3), and between the UC and BS microsites of each treatment. After fire, the SSI of the UC soil declined from 16 mm h⁻¹ to 12 mm h⁻¹ on T1, and from 24 mm h⁻¹ to 19 mm h⁻¹ on T2. The IC was reduced by 2/3 in the T1 UC soil, and by half on T2 UC soil. On the BS of T1 and T2, the fire effect was minimal, and higher infiltration rates and capacities were reached. Therefore, the presence/absence of vegetation when burnt influenced the post‐burnt infiltration patterns at soil microscale. On the T3, different rates and capacities were obtained depending on the microsites (UC and BS), with higher SSI (25 mm h⁻¹) and IC (226 mm h⁻¹) on BS than on UC (SSI of 18 mm h⁻¹ and IC of 136 mm h⁻¹). The SWRC and SWC were recovered from 1995 to 2003 (prior to the fires). The 2003 fire promoted high variability on the SWC at pF 0·1, 2 and 2·5, and the SWRC on burnt soils were reduced. To summarize, the IC and SSI post‐fire decreases were related to the lower infiltration rate at plot scale, the significant differences in the SWRC between burnt and control treatments, and the increase in the runoff yield (20%). According to the results, the MDI was a useful tool to characterize the soil infiltration on the vegetation patches of the Mediterranean maquia, and contrary to other studies, on the UC soil, the infiltration rate and IC, when soil was dry, were lower than that obtained on BS. Once the soil gets wet, similar values were found on both microenvironments. Copyright © 2010 John Wiley & Sons, Ltd.     

Estimating fog deposition at a Puerto Rican elfin cloud forest site: comparison of the water budget and eddy covariance methods

The deposition of fog to a wind‐exposed 3 m tall Puerto Rican cloud forest at 1010 m elevation was studied using the water budget and eddy covariance methods. Fog deposition was calculated from the water budget as throughfall plus stemflow plus interception loss minus rainfall corrected for wind‐induced loss and effect of slope. The eddy covariance method was used to calculate the turbulent liquid cloud water flux from instantaneous turbulent deviations of the surface‐normal wind component and cloud liquid water content as measured at 4 m above the forest canopy. Fog deposition rates according to the water budget under rain‐free conditions (0·11 ± 0·05 mm h^?1) and rainy conditions (0·24 ± 0·13 mm h^?1) were about three to six times the eddy‐covariance‐based estimate (0·04 ± 0·002 mm h^?1). Under rain‐free conditions, water‐budget‐based fog deposition rates were positively correlated with horizontal fluxes of liquid cloud water (as calculated from wind speed and liquid water content data). Under rainy conditions, the correlation became very poor, presumably because of errors in the corrected rainfall amounts and very high spatial variability in throughfall. It was demonstrated that the turbulent liquid cloud water fluxes as measured at 4 m above the forest could be only ～40% of the fluxes at the canopy level itself due to condensation of moisture in air moving upslope. Other factors, which may have contributed to the discrepancy in results obtained with the two methods, were related to effects of footprint mismatch and methodological problems with rainfall measurements under the prevailing windy conditions. Best estimates of annual fog deposition amounted to ～770 mm year^?1 for the summit cloud forest just below the ridge top (according to the water budget method) and ～785 mm year^?1 for the cloud forest on the lower windward slope (using the eddy‐covariance‐based deposition rate corrected for estimated vertical flux divergence). Copyright © 2006 John Wiley & Sons, Ltd.     

Evaluating the performance of topobathymetric LiDAR to support multi-dimensional flow modelling in a gravel-bed mountain stream

Stream biophysical processes are commonly studied using multi-dimensional numerical modelling that quantifies flow hydraulics from which parameters such as habitat suitability, stream carrying capacity, and bed mobility are derived. These analyses would benefit from accurate high-resolution stream bathymetries spanning tens of kilometres of channel, especially in small streams or where navigation is difficult. Traditional ground-based survey methods are limited by survey time, dense vegetation and stream access, and are usually only feasible for short reaches. Conversely, airborne topobathymetric LiDAR surveys may overcome these limitations, although limited research is available on how errors in LiDAR-derived digital elevation models (DEMs) might propagate through flow models. This study investigated the performance of LiDAR-derived topobathymetry in support of multi-dimensional flow modelling and ecohydraulics calculations in two gravel-bedded reaches (approximately 200 m long), one morphologically complex and one morphologically simple, and at the segment scale (32 km-long stream segment) along a 15 m-wide river in central Idaho, USA. We compared metre and sub-metre-resolution DEMs generated from RTK-GPS ground and Experimental Advanced Airborne Research LiDAR-B (EAARL-B) surveys and water depths, velocities, shear stresses, habitat suitability, and bed mobility modelled with two-dimensional (2D) hydraulic models supported by LiDAR and ground-surveyed DEMs. Residual statistics, bias (B), and standard deviation (SD) of the residuals between depth and velocity predicted from the model supported by LiDAR and ground-survey topobathymetries were up to −0.04 (B) and 0.09 m (SD) for depth and −0.09 (B) and 0.20 m s⁻¹ (SD) for velocity. The accuracy (B = 0.05 m), precision (SD = 0.09 m), and point density (1 point m⁻²) of the LiDAR topobathymetric survey (regardless of reach complexity) were sufficient to support 2D hydrodynamic modelling and derivative stream habitat and process analyses, because these statistics were comparable to those of model calibration with B = 0 m and SD = 0.04 m for water surface elevation and B = 0.05 m s⁻¹ and SD = 0.22 m s⁻¹ for velocity in our investigation. © 2020 John Wiley & Sons, Ltd.     

Quantitative effects of changes in agricultural irrigation on potential evaporation

Evaporation is a key element to the basin's water cycle. Agricultural irrigation has resulted in a significant variation of regional potential evaporation (E_pen). The spatiotemporal variation of E_pen and influencing factors in natural, agricultural, and desert areas in different developmental stages of irrigation in Heihe River Basin (HRB) from 1970 to 2017 were comparatively analysed in this study. This work focused on the correction effect of irrigation on the variation of E_pen. Agricultural water consumption in HRB significantly varied around 1998 due to agricultural development and water policy. Under the influence of irrigation, annual variations of E_pen in agricultural, natural, and desert areas were significantly different. From 1970 to 1998, the annual trend slope of E_pen in natural area only reduced by 1 mm decade⁻¹, while that in agricultural area significantly decreased by 39 mm decade⁻¹. After the implementation of water-saving irrigation, E_pen in natural and agricultural areas increased by 11 and 54 mm decade⁻¹, respectively, from 1998 to 2017. In contrast with natural and agricultural areas, E_pen in desert area decreased by 80 mm decade⁻¹ from 1970 to 1998 and continuously decreased by 41 mm decade⁻¹ from 1998 to 2017. However, the regulatory effect of irrigation on E_pen in desert area started to manifest due to the expansion of cultivated land area from 2010 to 2017. Irrigation had a significant regulatory effect on the variation of E_pen in HRB. The regulatory effect was mainly reflected on the aerodynamic term (E_aero). Results indicated that the main meteorological factors influencing E_pen in each region were wind speed, which is 2 m above the surface (U₂), and water vapour deficit (VPD).     

Evaporation estimation models for Lake Eğirdir,Turkey

Lake E?irdir is located in the Lakes District in southwestern Turkey and it is the second largest freshwater resource lake. Evaporation is an important parameter in hydrological and meteorological practical studies. This study has three objectives: (1) to develop models for the estimation of daily evaporation using measured data from the automated GroWeather meteorological station located near Lake E?irdir; (2) to compare the evaporation models with the classical Penman approach; (3) to evaluate the potential of each model. The comparisons are based on daily and monthly available data from 2001 and 2002. The evaporation estimation models (EEMs) developed in this paper have lower mean absolute errors and higher coefficient of determination R² values than the Penman method. In order to evaluate the potential of the EEMs, daily evaporation values are calculated by the Priestley–Taylor, Brutsaert–Stricker, de Bruin, Makkink and Hamon methods. The EEMs are statistically indistinguishable from the classical methods on the basis of the parameters of mean, standard deviation, etc. In the evaluation of daily and monthly values, the relative error percentage for daily evaporation has lower values than for monthly evaporation. It can be seen that the EEMs help in calculating daily evaporation rather than monthly. Final evaluation and comparison indicate that there is a good agreement between the results of EEMs and the Penman approach than with the classical methods. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of willow short‐rotation forest evaporation using a modified Shuttleworth–Wallace approach

Evaporation from a willow short‐rotation forest was analysed using a modified version of the Shuttleworth–Wallace model. The main modification consisted of a two‐layer soil module, which enabled soil surface resistance to be calculated as a function of the wetness of the top soil. Introduction of the threshold value of the leaf area index when scaling up from the leaf to the canopy resistance resulted in improvement to the simulated evaporation. The analysis was concentrated mainly on the 1988 season (May–October) when total evaporation was measured by the energy balance/Bowen ratio method throughout the growing season, covering all stages of canopy development. At the beginning of the 1994 season, soil evaporation were also measured with a ventilated chamber system. The general seasonal dynamics of the evaporation were fairly well simulated with the model. The largest deviation between measured and simulated evaporation occurred in June, when the model underestimated evaporation by about 1 mm day^?1. The model underestimated also in May but not as much as in June. In September and October the performance of the model was very good. For 130 days of the period May–October the cumulated measured evaporation was 364 mm and the simulated evaporation for the same days was 362 mm. It should be pointed out that this result was obtained without calibrating the model against the measured evaporation. The total simulated evaporation for the season was 450 mm with transpiration constituting 298 mm (66%), soil evaporation 102 mm (23%) and interception evaporation 50 mm (11%). The sensitivity analysis showed, in general, that simulated evaporation was most sensitive to changes in resistances when the leaf area index was smallest, i.e. under non‐closed canopy conditions. Changes in stomatal resistance, which is one of the most sensitive parameters, with associated changes in canopy transpiration, resulted in a negative feedback effect on soil evaporation. This reduced the total evaporation's sensitivity to stomatal resistance. This type of interaction between canopy and soil or undergrowth fluxes has been observed in other studies as well. Copyright © 2001 John Wiley & Sons, Ltd.     

The transport of aggregates associated with soil organic carbon under the rain-induced overland flow on the Chinese Loess Plateau

Organic carbon (OC) is easily enriched in sediment particles of different sizes due to aggregate breakdown and selective transport for sheet erosion. However, the transport of aggregate-associated OC has not been thoroughly investigated. To address this issue, 27 simulated rainfall experiments were conducted in a 1 m × 0.35 m box on slope gradients of 15°, 10°, and 15°and under three rainfall intensities of 45 mm h⁻¹, 90 mm h⁻¹ and 120 mm h⁻¹. The results showed that OC was obviously enriched in sediment particles of different sizes under sheet erosion. The soil organic carbon (SOC) concentrations of each aggregate size class in sediments were different from those in the original soil, especially when the rainfall intensity or slope was sufficiently low, such as 45 mm h^–1 or 5°, respectively. Under a slope of 5°, the SOC enrichment ratios (ERocs) of small macroaggregates and microaggregates were high but decreased over time. As rainfall intensity increased, OC became enriched in increasingly fine sediment particles. Under a rainfall intensity of 45 mm h^–1, the ERocs of the different aggregate size classes were always high throughout the entire erosion process. Under a rainfall intensity of > 45 mm h^–1 and slope of > 5°, the ERocs of the different aggregate size classes were close to 1.0, especially those of clay and silt. Therefore, the high ERocs in sediments resulted from the first transport of effective clay. Among total SOC loss, the proportion of OC loss caused by the transport of microaggregates and silt plus clay-sized particles was greater than 50%. We also found that low stream power and low water depth were two requirements for the high ERocs in aggregates. Stream power was closely related to sediment particle distribution. Flow velocity was significantly and positively related to the percentage of OC-enriched macroaggregates in the sediments (P > 0.01). Our study will provide important information for understanding the fate of SOC and building physical-based SOC transport models. © 2019 John Wiley & Sons, Ltd.     

Artificial neural networks and empirical equations to estimate daily evaporation: application to Lake Vegoritis,Greece

ABSTRACT

Evaporation is one of the most important components in the energy and water budgets of lakes and is a primary process of water loss from their surfaces. An artificial neural network (ANN) technique is used in this study to estimate daily evaporation from Lake Vegoritis in northern Greece and is compared with the classical empirical methods of Penman, Priestley-Taylor and the mass transfer method. Estimation of the evaporation over the lake is based on the energy budget method in combination with a mathematical model of water temperature distribution in the lake. Daily datasets of air temperature, relative humidity, wind velocity, sunshine hours and evaporation are used for training and testing of ANN models. Several input combinations and different ANN architectures are tested to detect the most suitable model for predicting lake evaporation. The best structure obtained for the ANN evaporation model is 4-4-1, with root mean square error (RMSE) from 0.69 to 1.35 mm d^?1 and correlation coefficient from 0.79 to 0.92.

EDITOR M.C. Acreman

ASSOCIATE EDITOR not assigned     

Simulation of Atmospheric Boundary Layer Characteristics during Indian Summer Monsoon using Observations from Monsoon Trough Boundary Layer Experiment at Jodhpur, India

The diurnal structure of the boundary layer during Indian summer monsoon period is studied using a one-dimensional meteorological boundary layer model and the observations collected from the Monsoon Trough Boundary Layer Experiment conducted in 1990 at Jodhpur, India. The model was initialized with the observed temperature profiles at 0530 LST on 17 July, 1990 at Jodhpur and was run for 26 hours. The study is carried out with a geostrophic wind speed of 9.5 m s⁻¹ corresponding to the strong wind simulation. The mean thermodynamic and wind structure simulated by the model are in good agreement with those observed from 30 m tower. The computed surface layer characteristics such as the surface fluxes, TKE and standard deviations of velocity components are found to be reasonably in good agreement with those based on turbulence measurements. The shear and buoyancy budget computed from the model are also compared with the turbulence measurements. The integrated cooling budget in the nocturnal boundary layer is examined.     

Seasonal soil water storage changes beneath central Amazonian rainforest and pasture

Evaporation and infiltration were compared for tropical rainforest and pasture, near to Manaus, Brazil from October 1990 to February 1992 using measurements of soil water storage over a depth of 2 m. The soil is a clayey oxisol of low water available capacity. In both of the dry seasons studied, the maximum change in soil water storage in the forest was 154 mm and in the pasture it was 131 and 112 mm. Similar behaviour of the soil water reservoir below forest and pasture in the wet season implied that differences in evaporation and drainage were small. In the dry season, soil water storage behaviour in the upper metre of the soil was similar but there were marked differences in the second metre. The pasture took up little water from below 1.5 m but the forest appeared to utilise all of the available water in the 2 m profile in both seasons.

The water balance of the 2 m profile showed that the pasture evaporation rate was equal to that of the forest until storage had decreased 80 mm from the maximum. There was then a decline in pasture evaporation rate to 1.2 mm day⁻¹ as the storage decreased by a further 50 mm. In contrast, the forest uptake rate remained above 3.5 mm day⁻¹ until storage had decreased 140 mm from the maximum (within 15 mm of the extraction limit), before declining abruptly to less than 1.5 mm day⁻¹. There was strong evidence that the forest was able to abstract water from depths greater than 3.6 m.

Spatial variability of soil water storage was significantly greater beneath the pasture than beneath the forest, particularly following rainfall events in the dry season. This was largely the result of redistribution of rainfall as local surface runoff. There was no evidence of redistribution or runoff in the forest.     

The role of urban surfaces (permeable pavements) in regulating drainage and evaporation: development of a laboratory simulation experiment

Permeable pavements and similar stormwater control devices have not been exploited in the UK, in part because their adoption has been hindered by a lack of detailed knowledge of their hydrological performance. This paper describes a research programme that produced detailed information on the hydrological behaviour of a car park surface. The study involved the construction of full‐scale permeable pavement model car park structures and a rainfall simulator for use in the laboratory. A monitoring procedure was developed in order to measure inputs and changes in drainage, storage and evaporation over short and long time‐scales (2 hours to 3 months). A range of simulated rainfalls, which varied in intensity and duration, was applied to the model car park surfaces. Hydrological processes were monitored over an 18‐month period. Results demonstrated that evaporation, drainage and retention in the structures were strongly influenced by the particle size distribution of the bedding material and by water retention in the surface blocks. In general, an average of 55% of a one‐hour duration, 15 mm h⁻¹ rainfall event could be retained by an initially air‐dry structure. Subsequent simulations demonstrated that 30% of a one‐hour duration, 15 mm h⁻¹ rainfall event could be stored by an initially wet structure (with a minimum time interval between rainfall applications of 72 hours). Copyright © 1999 John Wiley & Sons, Ltd.     

A spreadsheet‐based (Microsoft Excel) point surface energy balance model for glacier and snow melt studies

This paper describes a point surface energy balance model which runs within the Microsoft Excel spreadsheet package. The study incorporates a large amount of previous energy balance work and presents it in a useable form. The core model calculates the net shortwave and longwave radiation fluxes, the turbulent sensible and latent heat fluxes and the surface melt rate at a point on a melting ice or snow surface, from hourly inputs of incoming shortwave radiation, vapour pressure, air temperature and wind speed data. The latitude, longitude, slope angle, aspect, elevation, local temperature lapse rate, albedo and aerodynamic roughness of the study site, and the elevation of the meteorological station, can all be specified in the model. An output file containing the hourly and daily rates, and the totals of the energy fluxes is generated. The main advantages of the model are: first, that it requires only a PC or laptop computer running standard Microsoft Windows software, enabling it to be used at a desktop or in the field; and second, that it can be adapted quickly to different sites, meteorological data formats and other application requirements. Model calculations are compared with measured surface melt rates at five points on Haut Glacier d'Arolla, Switzerland, over a 115 day ablation period. Allowing for differences in shading between the meteorological station and the glacier, the root mean square error of the calculated melt rates is 2·0 mm day⁻¹ water equivalent melt (mean error +1·2 mm day⁻¹), for measured melt rates in the range 23 to 42 mm day⁻¹ water equivalent melt. Copyright © 2000 John Wiley & Sons, Ltd.     

A stream tracer technique employing ionic tracers and specific conductance data applied to the Maimai catchment,New Zealand

The stream tracer technique and transient storage models (TSMs) have become common tools in stream solute and hyporheic exchange studies. The expense and logistics associated with water sample collection and analysis often results in limited temporal resolution of stream tracer breakthrough curves (BTCs). Samples are often collected without a priori or real‐time knowledge of BTC information, which can result in poor sample coverage of the critical shoulder (initial rise) and tail (post‐steady state fall) of the BTC. We illustrate the use of specific conductance (SC) measurements as a surrogate for conservative dissolved tracer (Br⁻) samples. The advantages of collecting SC data for use in the TSM are (1) cost, (2) ease of data collection, and (3) well‐defined breakthrough curves, which strengthen TSM parameter optimization. This method is based on developing an ion concentration (IC)–SC relationship from limited discrete tracer solute samples. SC data can be collected on a more frequent basis at no additional analysis cost. TSM simulations can then be run for the conservative tracer data derived from SC breakthrough curves and the IC–SC relationship. This technique was tested in a 120 m reach of stream (2–60 m subreaches) in the Maimai M15 catchment, New Zealand during baseflow recession. Dissolved LiBr was injected for 12·92 h, with Br⁻ as the conservative ion of interest. Four TSM simulations using the OTIS model are optimized using UCODE to fit (1) Br⁻ data derived from the Br⁻–SC relationship (n = 1307 observations at each of two stream sampling sites), (2) all stream Br⁻ data collected (n = 58 in upper reach, n = 60 in lower reach), (3) half of the stream Br⁻ data collected, and (4) 20 stream Br⁻ samples from each site. No two simulations resulted in the same optimal parameter values. Results suggest that the greater the frequency of observations, the greater the confidence in estimated parameter values. Br⁻–SC simulations resulted in the best overall model fits to the data, with the lowest calculated error variance of 6·37, narrowest 95% parameter estimate confidence intervals, and the highest correlation coefficient of 0·99 942, among the four simulations. This is largely due to the improved representation of the shoulder and tail of the BTC with this method. The IC–SC correlation method is robust in situations in which (1) changes in background SC data can be accounted for, and (2) the data used to define the IC–SC relationship are representative of the range of data collected. This method provides more efficient sample analysis, improved data resolution, and improved model results compared to the alternative stream tracer data gathering methods presented. Additionally, we describe a new parameterization of the cross‐sectional area of the stream during flow recession, as a function of discharge, based on a stream hydraulic geometry relationship. This variant of the OTIS model provides a more realistic representation of stream dynamics during unsteady discharge. Copyright © 2005 John Wiley & Sons, Ltd.     

Roundhouse (RND) Mountain Top Research Site: Measurements and Uncertainties for Winter Alpine Weather Conditions

The objective of this work is to better understand and summarize the mountain meteorological observations collected during the Science of Nowcasting Winter Weather for the Vancouver 2010 Olympics and Paralympics (SNOW-V10) project that was supported by the Fog Remote Sensing and Modeling (FRAM) project. The Roundhouse (RND) meteorological station was located 1,856 m above sea level that is subject to the winter extreme weather conditions. Below this site, there were three additional observation sites at 1,640, 1,320, and 774 m. These four stations provided some or all the following measurements at 1 min resolution: precipitation rate (PR) and amount, cloud/fog microphysics, 3D wind speed (horizontal wind speed, U _h; vertical air velocity, w _a), visibility (Vis), infrared (IR) and shortwave (SW) radiative fluxes, temperature (T) and relative humidity with respect to water (RH_w), and aerosol observations. In this work, comparisons are made to assess the uncertainties and variability for the measurements of Vis, RH_w, T, PR, and wind for various winter weather conditions. The ground-based cloud imaging probe (GCIP) measurements of snow particles using a profiling microwave radiometer (PMWR) data have also been shown to assess the icing conditions. Overall, the conclusions suggest that uncertainties in the measurements of Vis, PR, T, and RH can be as large as 50, >60, 50, and >20 %, respectively, and these numbers may increase depending on U _h, T, Vis, and PR magnitude. Variability of observations along the Whistler Mountain slope (~500 m) suggested that to verify the models, model space resolution should be better than 100 m and time scales better than 1 min. It is also concluded that differences between observed and model based parameters are strongly related to a model’s capability of accurate prediction of liquid water content (LWC), PR, and RH_w over complex topography.     

Thermische Energieflüsse an der Wasseroberfläche: Beispiel Sempachersee

Data collected every 20 minutes for 18 months by a meteorological buoy moored on Lake Sempach in Switzerland (maximum depth 86 m, surface area 14.1 km²) are used to calculate different processes contributing to the net heat flux between water and atmosphere. The processes considered are shortwave and longwave radiation, evaporation/condensation and sensible heat transfer. The temporal resolution of the measurements allows the evaluation of the processes occurring on three different time scales: diurnal variations, weather events of a few days and yearly cycles. The heat content of the lake is calculated from quasi-continuous measurements of water temperature at different depths. The yearly amplitude of the heat content is 1100·10⁶ J/m². Short-time variations of the heat flux determined from water temperature analysis agree well with the flux variations modeled using meteorological data. However, the latter generally underestimates the measured net heat flux in the long term. Wind measurements, together with the net heat flux, are used to calculate the Flux Richardson Number and the Monin-Obukov Length. Examples are given to show the predominant influence of the wind on the stratification of the upper water column and thus on the surface water temperature.