Modelling snowpack surface temperature in the Canadian Prairies using simplified heat flow models

Three practical schemes for computing the snow surface temperature T_s, i.e. the force–restore method (FRM), the surface conductance method (SCM), and the Kondo and Yamazaki method (KYM), were assessed with respect to T_s retrieved from cloud‐free, NOAA‐AVHRR satellite data for three land‐cover types of the Paddle River basin of central Alberta. In terms of R², the mean T_s, the t‐test and F‐test, the FRM generally simulated more accurate T_s than the SCM and KYM. The bias in simulated T_s is usually within several degrees Celsius of the NOAA‐AVHRR T_s for both the calibration and validation periods, but larger errors are encountered occasionally, especially when T_s is substantially above 0 °C. Results show that the simulated T_s of the FRM is more consistent than that of the SCM, which in turn was more consistent than that of the KYM. This is partly because the FRM considers two aspects of heat conduction into snow, a stationary‐mean diurnal (sinusoidal) temperature variation at the surface coupled to a near steady‐state ground heat flux, whereas the SCM assumes a near steady‐state, simple heat conduction, and other simplifying assumptions, and the KYM does not balance the snowpack heat fluxes by assuming the snowpack having a vertical temperature profile that is linear. Copyright © 2005 John Wiley & Sons, Ltd.     

Verification of land surface evapotranspiration estimation from remote sensing spatial contextual information

Estimation of evapotranspiration (ET) is of great significance in modeling the water and energy interactions between land and atmosphere. Negative correlation of surface temperature (T_s) versus vegetation index (VI) from remote sensing data provides diagnosis on the spatial pattern of surface soil moisture and ET. This study further examined the applicability of T_s–VI triangle method with a newly developed edges determination technique in estimating regional evaporative fraction (EF) and ET at MODIS pixel scale through comparison with large aperture scintillometer (LAS) and high‐level eddy covariance measurements collected at Changwu agro‐ecological experiment station from late June to late October, 2009. An algorithm with merely land and atmosphere products from MODIS onboard Terra satellite was used to estimate the surface net radiation (R_n) and soil heat flux. In most cases, the estimated instantaneous R_n was in good agreement with surface measurement with slight overestimation by 12 W/m². Validation results from LAS measurement showed that the root mean square error is 0.097 for instantaneous EF, 48 W/m² for instantaneous sensible heat flux, and 30 W/m² for daily latent heat flux. This paper successfully presents a miniature of the overall capability of T_s–VI triangle in estimating regional EF and ET from limited number of data. For a thorough interpretation, further comprehensive investigation needs to be done with more integration of remote sensing data and in‐situ surface measurements. Copyright © 2011 John Wiley & Sons, Ltd.     

Factors affecting longitudinal dispersion in estuaries of different scale

Traditionally, the overall diluting capacity of an estuary is characterized using a coefficient of longitudinal dispersion, K_xe, which is given by the ratio of the dispersive flux of a dissolved substance to its tidally averaged longitudinal gradient. A steady-state model, which assumes a balance between the dispersive and non-dispersive fluxes and an exponential increase in estuary cross-sectional area towards the sea, has been used to derive expressions for the axial salinity distribution and the dispersive flux of salt. The model was set up assuming either a constant dispersion coefficient along the estuary or one that increased with distance towards the sea. By comparing salinity predictions with data from five UK estuaries, estimates were made of the maximum dispersive salt flux and the corresponding maximum salinity gradient for each system. The results indicated that there was an approximately linear relation between the fluxes and gradients, and the slope of a line plotted through the origin provided an estimate of a common K_xe for all five estuaries. The magnitude of K_xe was found to be about 90 m² s^–1 with a standard deviation of approximately ±32 m² s^–1. It is concluded that a representative value of 100 m² s^–1 for K_xe is a reasonable first choice when setting up a cross-sectionally averaged estuary model. The results also showed that larger systems, such as the Thames, had lower salinity gradients and lower dispersive salt fluxes, whilst smaller estuaries displayed the opposite characteristics. The model was used to predict the variation in the non-advective flux of salt along an estuary. The distribution was found to be similar to the corresponding flux distribution estimated from observations at the seaward end of the Tees estuary, despite appreciable spatial variations in the individual flux components. Allowing for a small decrease in freshwater flow, the model indicated that there was a decrease in the maximum dispersive flux between neap and spring tides. It is argued that such a reduction in flux can result in a seaward shift in the salinity distribution to a region of greater cross-section, where the freshwater transport per unit area again balances the reduced upstream dispersive flux, as found in the neap to spring response in the Tees estuary.Responsible Editor: Hans Burchard     

An Analysis of Turbulent Heat Fluxes and the Energy Balance During the REFLEX Campaign

Three eddy covariance stations were installed at the Barrax experimental farm during the Land-Atmosphere Exchanges (REFLEX) airborne training and measurement campaign to provide ground truth data of energy balance fluxes and vertical temperature and wind profiles. The energy balance closure ratio (EBR) was 105% for a homogeneous camelina site, 86% at a sparse reforestation site, and 73% for a vineyard. We hypothesize that the lower closure in the last site was related to the limited fetch. Incorporating a vertical gradient of soil thermal properties decreased the RMSE of the energy balance at the camelina site by 16 W m^?2. At the camelina site, eddy covariance estimates of sensible and latent heat fluxes could be reproduced well using mean vertical profiles of wind and temperature, provided that the Monin—Obukhov length is known. Measured surface temperature and sensible heat fluxes suggested high excess resistance for heat (kB^?1 = 17).     

CO2 absorption of sandy soil induced by rainfall pulses in a desert ecosystem

Soil CO₂ flux is strongly influenced by precipitation in many ecosystem types, yet knowledge of the effects of precipitation on soil CO₂ flux in semi‐arid desert ecosystems remains insufficient, particularly for sandy soils. To address this, we investigated the response of sandy soil CO₂ flux to rainfall pulses in a desert ecosystem in northern China during August–September 2011. Significant changes (P < 0.05) were found in diel patterns of soil CO₂ flux induced by small (2.1 mm), moderate (12.4 mm) and large (19.7 mm) precipitation events. Further analysis indicated that rainfall pulses modified the response of soil CO₂ flux to soil temperature, including hysteresis between soil CO₂ flux and soil temperature, with F_s higher when T_s was increasing than when T_s was decreasing, and the linear relationship between them. Moreover, our results showed that rainfall could result in absorption of atmospheric CO₂ by soil, possibly owing to mass flow of CO₂ induced by a gradient of gas pressure between atmosphere and soil. After each precipitation event, soil CO₂ flux recovered exponentially to pre‐rainfall levels with time, with the recovery times exhibiting a positive correlation with precipitation amount. On the basis of the amounts of precipitation that occurred at our site during the measurement period (August–September), the accumulated rain‐induced carbon absorption evaluated for rainy days was 1.068 g C m^?2; this corresponds approximately to 0.5–2.1% of the net primary production of a typical desert ecosystem. Thus, our results suggest that rainfall pulses can strongly influence carbon fluxes in desert ecosystems. Copyright © 2014 John Wiley & Sons, Ltd.     

Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes

Abstract

Reliable estimation of sensible heat flux (H) is important in energy balance models for quantifying evapotranspiration (ET). This study was conducted to evaluate the value of adding the Priestley-Taylor (PT) equation to the METRIC (Mapping Evapotranspiration at high Resolution with Internalized Calibration) model. METRIC was used to estimate energy fluxes for 10 Landsat images from the 2005, 2006 and 2007 crop growing seasons in south-central Nebraska, USA, where each image owing to recent rainfall exhibited high residual moisture content even at the hot pixel. The METRIC model performed satisfactorily for net radiation (R_n ) and soil heat flux (G) estimation with a root mean square error (RMSE) of 52 and 24 W m^-2, respectively. A RMSE of 122 W m^-2 for H indicated the limitation of the METRIC model in estimating H for high residual moisture content of the hot pixel (Alfalfa reference ET fraction, ET _r F > 0.15). The modified METRIC model (wet METRIC or wMETRIC) incorporating the PT equation was applied to calculate H at the anchor pixels (hot and cold) for high residual moisture content of the hot pixel. The α coefficient of the PT equation was locally calibrated using hourly meteorological data from an automatic weather station and R_n and G data from a Bowen ratio flux tower. The mean α coefficient value was 1.14. The wMETRIC model reduced the RMSE of H from 122 to 64 W m^-2 and that of latent heat flux, LE, from 163 to 106 W m^-2. The RMSE of daily ET decreased from 1.7 to 1.1 mm d^-1 with wMETRIC. The results indicate that treatment of anchor pixels for high residual moisture content with the PT approach gives improved estimation of H, LE and daily ET. It is recommended that the wMETRIC model be used for estimating ET if the hot pixel has high residual moisture (i.e. reference ET fraction > 0.15).

Citation Singh, R. K. & Irmak, A. (2011) Treatment of anchor pixels in the METRIC model for improved estimation of sensible and latent heat fluxes. Hydrol. Sci. J. 56(5), 895–906.     

It goes both ways: measurements of simultaneous evapotranspiration and fog droplet deposition at a montane cloud forest

Fluxes of latent heat, sensible heat, and water vapor, including turbulent deposition of fog droplets, were measured for two months in autumn 2005 within a subtropical montane cypress forest in Taiwan. The goal of the study was to determine whether significant evapotranspiration can occur during foggy conditions. Water vapor fluxes, Q_W, as determined with the Bowen Ratio method, were compared to those simultaneously measured with the eddy covariance method. The median Bowen Ratio was 1.06, and the median Q_W flux was 5 · 2 × 10^?5 kg m^?2 s^?1. The vertical gradients of temperature and specific humidity over the forest, ΔT and Δq, peaked around noon during days without fog, and were reduced during foggy conditions. For 66% of the data points, ΔT and Δq were negative, corresponding to positive (upward) fluxes of sensible heat Q_H and latent heat Q_E. A Monte Carlo simulation proved that statistically significant evapotranspiration rates, i.e., upward water vapor fluxes, occurred during fog. At the same time, deposition fluxes of fog droplets occurred. Our results show that even during fog events, significant evapotranspiration may occur. Copyright © 2008 John Wiley & Sons, Ltd.     

Estimation of mass and energy balance of glaciers using a distributed energy balance model over the Chandra river basin (Western Himalaya)

The ongoing glacier shrinking in the Himalayan region causes a significant threat to freshwater sustainability and associated future runoff. However, data on the spatial climatic contribution of glacier retreat is scanty in this region. To investigate the spatially distributed glacier surface energy and mass fluxes, a two-dimensional mass balance model was developed and applied to the selected glaciers of the Chandra basin, in the Upper Indus Basin, Western Himalaya. This model is driven by the remote sensing data and meteorological variables measured in the vicinity of the Chandra basin for six hydrological years (October 2013 to September 2019). The modelled variables were calibrated/validated with the in-situ observation from the Himansh station in the Chandra basin. We have derived air temperature (T_a ) spatially using the multivariate statistical approach, which indicates a relative error of 0.02–0.05°C with the observed data. Additionally, the relative error between the modelled and observed radiation fluxes was <10.0 W m⁻². Our study revealed that the Chandra basin glaciers have been losing its mass with a mean annual mass balance of −0.59 ± 0.12 m w.e. a⁻¹ for the six hydrological years. Results illustrated that the mean surface melt rate of the selected glaciers ranged from −5.1 to −2.5 m w.e. a⁻¹ that lies between 4500 and 5000 m a.s.l. The study revealed that the net radiation (R_N) contributes ~75% in total energy (F_M ) during the melt season while sensible heat (H_S) , latent heat (H_l) , and ground heat (H_G) fluxes shared 15%, 8%, and 2%, respectively. Sensitivity analysis of the energy balance components suggested that the mass balance is highly sensitive to albedo and surface radiations in the study area. Overall, the proposed model performed well for glacier-wide energy and mass balance estimation and confirms the utility of remote sensing data, which may help in reducing data scarcity in the upper reaches of the Himalayan region.     

Application of a spreadsheet hydrological model for computing the long-term water balance of Lake Awassa,Ethiopia

Abstract

The water balance of the closed freshwater Lake Awassa was estimated using a spreadsheet hydrological model based on long-term monthly hydrometeorological data. The model uses monthly evaporation, river discharge and precipitation data as input. The net groundwater flux is obtained from model simulation as a residual of other water balance components. The result revealed that evaporation, precipitation, and runoff constitute 131, 106 and 83 × 10⁶ m³ of the annual water balance of the lake, respectively. The annual net groundwater outflow from the lake to adjacent basins is 58 × 10⁶ m³. The simulated and recorded lake levels fit well for much of the simulation period (1981–1999). However, for recent years, the simulated and recorded levels do not fit well. This may be explained in terms of the combined effects of land-use change and neotectonism, which have affected the long-term average water balance. With detailed long-term hydrogeological and meteorological data, investigation of the subsurface hydrodynamics, and including the effect of land-use change and tectonism on surface water and groundwater fluxes, the water balance model can be used efficiently for water management practice. The result of this study is expected to play a positive role in future sustainable use of water resources in the catchment.     

Mining lakes as groundwater‐dominated hydrological systems: assessment of the water balance of Mining Lake Plessa 117 (Lusatia,Germany) using stable isotopes

In the present study, the stable isotopes δ¹⁸O and δ²H were used for assessment of the water balance in a heterogeneously structured catchment area in the Lusatian Lignite Mining District, in particular, for estimation of the annual groundwater inflow and outflow (I_GW and O_GW) of Mining Lake Plessa 117. The application of stable isotopes was possible since the water exchange in the catchment area had reached steady‐state conditions after the abandonment of mining activities in 1968 and the filling of the voids and aquifers by re‐rising groundwater in the years thereafter. Diverging slopes of the Evaporation Line and the Global Meteoric Water Line manifested as evaporation from the lake catchment area. The calculated isotope water balance was compared with the commonly used surface water balance, which is unable to differentiate between I_GW and O_GW, and with a local groundwater model. The groundwater model calculated an I_GW of about 811 000 m³ yr^?1 and an O_GW close to zero, whereas the isotope water balance showed fluxes of about 914 000 and 140 000 m³ yr^?1, respectively. Considering the contribution of the groundwater inflow to the total annual input into the lake (ΔI_T) and the mean residence time (τ), where the groundwater model and the isotope water balance calculated 42 and 47% for ΔI_T and 4·3 and 3·9 years for τ, respectively, it was shown that both water balance calculation methods led to comparable results despite the differences in I_GW and O_GW. Copyright © 2008 John Wiley & Sons, Ltd.     

大型浅水湖泊与大气之间的动量和水热交换系数——以太湖为例

湖泊水面与大气之间垂直方向的动量、水汽和热量通量与风速、湿度和温度梯度之间存在比例关系,因此在湖泊水-气相互作用研究中,这比例系数(交换系数)是关键因子.在以往的研究中,交换系数通常直接采用水面梯度观测法或海洋大气近地层的参数化方案进行计算.本文采用涡度相关系统和小气候系统仪器在太湖平台上直接观测的通量和气象要素,对上述交换系数(最小均方差原则)进行优化,结果为:动量交换系数C_D10N=1.52×10^-3、水汽交换系数C_E10N=0.82×10^-3、热量交换系数CH10N=1.02×10^-3,与其他内陆湖泊涡度相关观测数据的推导结果一致.本文的研究结果表明:与海洋参数化方案相比,在相同的风速条件下,湖面的空气动力学粗糙度比海洋高,这可能是由于受到水深的影响;如果采用海洋参数化方案,会导致湖泊年蒸发量的估算值偏大40%.太湖的动量、水汽和热量交换系数可以视为常数,可以不考虑稳定度和风速的影响.这是因为本文中83%的数据为近中性条件.敏感性分析表明:如果考虑稳定度的影响,LE模拟值的平均误差降低了0.5 W/m²,H的平均误差降低了0.4 W/m²,u*的计算值没有变化;如果考虑风速的影响,u*模拟值的平均误差降低了0.004 m/s,LE的平均误差升高了1.3 W/m²,H的模拟结果几乎不受影响.这一结果能为湖气相互作用研究提供参考.     

A simple method using climatic variables to estimate canopy temperature,sensible and latent heat fluxes in a winter wheat field on the North China Plain

Estimation of evapotranspiration from a crop field is of great importance for detecting crop water status and proper irrigation scheduling. The Penman–Monteith equation is widely viewed as the best method to estimate evapotranspiration but it requires canopy resistance, which is very difficult to determine in practice. This paper presents a simple method simplified from the Penman–Monteith equation for estimating canopy temperature (T_c). The proposed method is a biophysically‐sound extended version of that proposed by Todorovic. The estimated canopy temperature is used to calculate sensible heat flux, and then latent heat flux is calculated as the residual of the surface energy balance. An eddy covariance (EC) system and an infrared thermometer (IRT) were installed in an irrigated winter wheat field on the North China Plain in 2004 and 2005, to measure T_c, and sensible and latent heat fluxes were used to test the modified Todorovic model (MTD). The results indicate that the original Todorovic model (TD) severely underestimates T_c and sensible heat flux, and hence severely overestimates the latent heat flux. However, the MTD model has good capability for estimating T_c, and gives acceptable results for latent heat flux at both half‐hourly and daily scales. The MTD model results also agreed well with the evapotranspiration calculated from the measured T_c. Copyright © 2008 John Wiley & Sons, Ltd.     

城市景观湖泊不同时间尺度CH4通量特征及影响因素分析——以南京市莫愁湖为例

城市景观湖泊对温室气体的收支发挥着重要作用。本文以南京市莫愁湖为研究对象,采用静态箱—温室气体分析仪法实时监测湖泊水—气界面CH₄通量,分析湖泊主要温室气体CH₄在日尺度和季节尺度上因冒泡和扩散排放方式不同对其通量的影响,探究影响湖泊CH₄通量的因素。结果表明:(1)在日尺度上,四季24 h内CH₄均呈排放状态,受白天冒泡影响,四季CH₄总通量均存在白天高于夜间的日变化特征。(2)在季节尺度上,莫愁湖CH₄排放通量呈现显著的时空异质性,受冒泡通量的影响夏季CH₄通量明显高于春、秋、冬三季;B区的CH₄总通量(6.04 nmol/(m²·s))显著高于A区(3.82 nmol/(m²·s)),水体的营养化程度和离岸距离是空间变化的主要影响因素。A、B两区CH₄排放夏季以冒泡排放为主,春、秋、冬以扩散排放为主。(3)在日尺度上,CH₄     

Seasonal climate patterns and their influence on calibration of the Hargreaves-Samani equation

Abstract

Annual patterns in climate parameters were studied to evaluate how these influence the quality of reference evapotranspiration (ETo) estimates obtained from the Hargreaves-Samani (HS) equation, since the method only uses the measured temperature directly. The work evaluates how these patterns can be used to improve the HS ETo estimates. Ten-year moving averages from a set of California Irrigation Management Information System (CIMIS) stations were used to evaluate the relationships between solar radiation (R_s), temperature (T) and ETo. The results indicate that T treads behind solar radiation and its value peaks some 25 days later. Thus, the main irrigation season in the Mediterranean climate (1 May–30 September) can be divided into three phases: increasing R_s and T; decreasing R_s with increasing T; and decreasing R_s and T. Non-univocal annual cycles were observed between R_s and T, ETo and R_s, and ETo and T. These annual patterns result in important seasonal changes in the ratio between the HS and Penman-Monteith (FAO PM) ETo estimates. The changes are particularly important during the irrigation season, where the FAO PM initially calculates greater ETo values than the HS methodology, and from the end of May to early September, where the HS equation overestimates the ETo values (by 17 mm, or 3%). These patterns obtained from 2000–2009 data were used to calibrate and improve HS ETo estimates at new sites for the 2010–2011 period. Calibration based on the proposed seasonal region-wide FAO PM/HS ETo ratios improved both the bias (decreased from 0.40 to 0.36 mm d^-1) and r² (increased from 0.67 to 0.87) of the ETo estimates for the irrigation season. The proposed methodology can be easily applied to other regions, even when the existing weather stations are sparse.

Editor Z.W. Kundzewicz     

Pulsating of the generalized ion and neutral polar winds

A three-dimensional, time-dependent model of the ion and neutral polar winds was used to study their dynamic evolution during the May 4, 1998 magnetic storm. The simulation tracked the dynamics of five species (O⁺, H⁺, H_s, O_s, and electrons) and covered a 9-h period. During the storm, D_st decreased to −210 nT, Ap reached 300, and K_p was elevated. The IMF B_z component was southward at the start of the storm and for several hours thereafter and then turned northward. However, the magnetospheric energy input to the ionosphere exhibited a 50-min oscillation, with the plasma convection and particle precipitation patterns expanding and contracting in a periodic manner. As a consequence, the ion and neutral polar winds pulsated with an approximate 50-min period. The H⁺ and O⁺ ions displayed cyclic upflows and downflows in the topside ionosphere as well as a highly structured spatial distribution that varied with time. The vertical flux of the neutral H_s atoms was upward at the top of the ionosphere, but the magnitude varied in a cyclic manner in response to the oscillating stormtime energy input. The vertical flux of neutral O_s atoms was downward at the top of the ionosphere and varied significantly with the stormtime energy input. For H⁺, O⁺, and H_s, the maximum total (integrated) vertical flux during the storm was upward at the top of the ionosphere, with values of 8–9×10²⁵ particles/s for H⁺, 2–4×10²⁶ particles/s for O⁺, and 2–3×10²⁷ particles/s for H_s. The corresponding total vertical O_s flux was predominately downward, with only localized areas with positive fluxes.     

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.     

Water transport through Lake Kinneret sediments traced by tritium

The transport of water through the sediments of Lake Kinneret was estimated using tritium of fallout origin as a tracer. The calculation model takes into account the diffusion of the tracer and its decay, and allows for the advective flux through the sediment pores as a free parameter. The dependence of the interstitial diffusion constant on the porosity was assumed identical to that for chloride ions. Tritium profiles in nine cores were used for estimating the advective fluxes and velocities, by finding the best fit between data and model. Water advective fluxes of 3.4 g cm^?2 yr^?1 into the lake were found. This water flux is negligible in the hydrological balance of Lake Kinneret; but it carries along 7% of the unaccounted-for chloride input to the lake.     

Short‐term stream water temperature observations permit rapid assessment of potential climate change impacts

Assessment of potential climate change impacts on stream water temperature (T_s) across large scales remains challenging for resource managers because energy exchange processes between the atmosphere and the stream environment are complex and uncertain, and few long‐term datasets are available to evaluate changes over time. In this study, we demonstrate how simple monthly linear regression models based on short‐term historical T_s observations and readily available interpolated air temperature (T_a) estimates can be used for rapid assessment of historical and future changes in T_s. Models were developed for 61 sites in the southeastern USA using ≥18 months of observations and were validated at sites with longer periods of record. The T_s models were then used to estimate temporal changes in T_s at each site using both historical estimates and future T_a projections. Results suggested that the linear regression models adequately explained the variability in T_s across sites, and the relationships between T_s and T_a remained consistent over 37 years. We estimated that most sites had increases in historical annual mean T_s between 1961 and 2010 (mean of +0.11 °C decade^?1). All 61 sites were projected to experience increases in T_s from 2011 to 2060 under the three climate projections evaluated (mean of +0.41 °C decade^?1). Several of the sites with the largest historical and future T_s changes were located in ecoregions home to temperature‐sensitive fish species. This methodology can be used by resource managers for rapid assessment of potential climate change impacts on stream water temperature. Copyright © 2014 John Wiley & Sons, Ltd.     

Surface energy balance at a tropical station

Summary Studies of various fluxes, namely net radiation, soil heat, sensible heat and latent heat observed at a tropical station are presented in this paper. The time variation of these fluxes are examined in relation to various meteorological parameters and atmospheric conditions. The turbulent transfer coefficients have been evaluated to examine the applicability of the classical theory or the non-equivalence theory for eddy transport in the lower layers of the atmosphere. The energy balance at a tropical station is evaluated. It is found over year there is a net surplus of 94,000 ly. A detailed discussion of the disposal of this energy by various consuming processes is given.Symbols and notation All the quantities represented by symbols in the text of the paper are defined below C _p specific heat at constant pressure in cal. g^–1 °K - E evaporation in g cm^–2 hr^–1 - E ^* evaporative heat flux cal. cm^–2 sec^–1 (in Eq. 10) - e vapour pressure in millibars - e _z ,e _2z vapour pressure at heightsz and 2z - g acceleration due to gravity - H sensible heat flux cal cm^–2 sec^–1 (in Eq. 12) - K _M ,K _H ,K _W coefficients of eddy diffusivities of momentum, heat, and water vapour respectively in cm² sec^–1 - k von Karman' constant=0.4 - L Monin-Obukhov length (according to Monin and Obukhov [53] the structure of the turbulent boundary layer is determined by the non-dimensional variableZ/L whereL is defined byL=–(u _* ³ C _p T)/(kgH) - ly langleys - Q _c Q—sensible heat flux in langleys (in Eqs. 3 and 4) - Q _e E—latent heat flux in langleys (in Eqs. 3 and 4) - Q _s S—soil heat flux in langleys (in Eqs. 3 and 4) - Q _i Q _c +Q _e +Q _s whenK _M K _H K _W , (in Eq. 6) - Q' _i Q' _c +Q' _e +Q _s whenK _M =K _H =K _W (in Eq. 7) - qq mean specific humidity g kg^–1