Multi-objective conditioning of a simple SVAT model

Air and biomass heat storage fluxes in a forest canopy: Calculation within a soil vegetation atmosphere transfer model

The effect of a thin and light greening system with bamboo charcoal layer for water retention on heat fluxes, in particular latent heat flux (evapotranspiration rate), under no irrigation condition, on the rooftop of a building in Higashi-Hiroshima, West Japan, was investigated. In April 2019, lawn seeds (Zoysia tenuifolia) were sown which were germinated, reached a height of 70 mm by May when 100% of the vegetation area was covered. The air temperature and humidity at two different heights (0.3 m and 1.8 m) above greening soil surface, latent, and sensible heat fluxes were estimated. Bowen ratio was employed to collect the data on surface heat balance and soil water content during the summer season (June to September) in 2019 on the rooftop of a building in Higashi-Hiroshima, West Japan. The latent heat during daytime for a week without rainfall in each month was compared with the evapotranspiration rate. Owning to the vegetation development, the ground heat flux on greening soil surface decreased from -400 W/m2 to -200 W/m2 (flux from air to soil) during sunny daytime in July, and it was less than -100 W/m2 in August, although net radiation was maintained around 800 W/m2 over the season except in September. The monthly net radiation flux for an entire day (daytime and nighttime) ranged between 55 and 125 W/m2 (average: 95 W/m2) for the summer season of which 32-66% (average: 48%) was occupied by latent heat. Evapotranspiration from greening soil ranged between 1.24 and 1.82 mm/day, averaged at 1.51 mm/day throughout the season, which corresponded to about 26% of total rainfall over the season ( r 2 = 0.88 , p < 0.01 ; S . E = 0.06 ) between the estimated and measured values. These observations suggested that the thin and bamboo coal light soil layer greening system, even without constant irrigation, could maintain the development of lawn grass and transformed more than half of net radiation to latent heat, i.e., evapotranspiration, insulating most ground heat in midsummer, which may be mostly due to bamboo charcoal sublayer.

Thermal infrared imagery provides a measure of surface temperature by quantifying the thermal emission of radiation from objects. In this work, we propose a novel approach for computing surface sensible heat fluxes from a time series of thermal imagery sampled at 20 Hz. The approach consists of using a surface energy budget, where the ground heat flux is easily computed from limited measurements using a force–restore‐type methodology, the latent heat fluxes are neglected, and the energy storage is computed using a lumped capacitance model. Preliminary validation of the method is presented using experimental data acquired at a desert playa in western Utah from an intensive observation period during the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) Program. The sensible heat fluxes computed from thermal imagery show encouraging agreement (≈4% difference for 1 h averaged fluxes computed with the force–restore method) compared to heat fluxes obtained using the eddy‐covariance technique, however additional evaluation is required to confirm the method's validity. Further, the technique illustrates the richness of surface heat flux spatio‐temporal structure even over very smooth surfaces.

New satellite-derived latent and sensible heat fluxes are performed by using WindSat wind speed, WindSat sea surface temperature, the European Centre for Medium-range Weather Forecasting (ECMWF) air humidity, and ECMWF air temperature from 2004 to 2014. The 55 moored buoys are used to validate them by using the 30 min and 25 km collocation window. Furthermore, the objectively analyzed air-sea heat fluxes (OAFlux) products and the National Centers for Environmental Prediction-National Center for Atmospheric Research reanalysis 2 (NCEP2) products are also used for global comparisons. The mean biases of sensible and latent heat fluxes between WindSat flux results and buoy flux data are–0.39 and–8.09 W/m2, respectively. In addition, the root-mean-square (RMS) errors of the sensible and latent heat fluxes between them are 5.53 and 24.69 W/m2, respectively. The RMS errors of sensible and latent heat fluxes are observed to gradually increase with an increasing buoy wind speed. The difference shows different characteristics with an increasing sea surface temperature, air humidity, and air temperature. The zonal average latent fluxes have some high regions which are mainly located in the trade wind zones where strong winds carry dry air in January, and the maximum value centers are found in the eastern waters of Japan and on the US east coast. Overall, the seasonal variability is pronounced in the Indian Ocean, the Pacific Ocean, and the Atlantic Ocean. The three sensible and latent heat fluxes have similar latitudinal dependencies; however, some differences are found in some local regions.

Seasonal and diurnal variation of energy budget components in coniferous forests

基于三江平原稻田2005~2007年5~10月涡度相关通量观测数据, 分析了该区稻田能量通量的日变化、季节变化和能量分配特征以及能量平衡状况。结果表明: 三江平原稻田净辐射和潜热通量日变化均表现为明显的单峰特征, 感热通量日变化在水稻发育进入成熟期后才较明显, 而土壤热通量在水稻整个发育期内日变化特征都不明显。稻田净辐射季节变化特征显著, 6月下旬至7月上旬达到最大值18~20 MJ·m^-2·d^-1。潜热通量季节变化与净辐射同步, 最大值为13~19 MJ·m^-2·d^-1。相比之下感热通量较小, 观测期间变化于-3.90~ 3.94 MJ·m^-2·d^-1, 且没有明显的季节变化。5~10月土壤热通量呈下降趋势, 变化于-2.67~3.62 MJ·m^-2·d^-1。三江平原地区稻田能量分配特征明显, 潜热通量占净辐射的比例(<i>LE</i>/<i>R</i>n) 5~10月平均值为0.67, 表明净辐射大部分以潜热通量形式所消耗, 但生长旺季<i>LE</i>/<i>R</i>n略大于生长季初期和末期。感热通量占净辐射的比例(<i>H</i>s/<i>R</i>n)的季节变化特征与<i>LE</i>/<i>R</i>n比值相反, 观测期间平均值为0.10。这导致波文比在水稻生长旺季较小而在初期和末期较大。5~10月土壤热通量占净辐射的比例(<i>G</i>/<i>R</i>n)呈逐渐下降趋势, 其月平均值由5月的0.14下降到10月的-0.08。线性回归法和能量平衡比率均表明三江平原稻田能量明显不闭合, 2005、2006年5~10月能量不闭合度分别为22%和16%, 而2007年能量"过闭合", 能量平衡比率平均值为1.07。

Soil Vegetation Atmosphere Transfer (SVAT) models have been implemented to simulate energy and mass fluxes between soil, vegetation and atmosphere of various ecosystems. Usually, these models are simple, but they use realistic descriptions of radiative, turbulent and water transfers. These include description of stomatal control of transpiration and CO/sub 2/ fluxes. They can be used for assimilating remote sensing data and derive vegetation canopy evapotranspiration or photosynthesis. Various remote sensing data may provide useful information to drive SVAT models. Surface temperature may be used through inversion procedures to retrieve parameters related to stomatal conductance or root zone soil moisture. Parameters related to vegetation structure (LAI, vegetation height) may be retrieved from reflectance measurements in the solar domain, either through direct relationships with some vegetation index or by inverting radiative transfer formulation against spectral reflectance measurements. The microwave data contribution has not been studied very often in the case of vegetation canopies, but they were proposed for estimating surface soil moisture. In this paper, inversions of the ALiBi model were performed to retrieve canopy evapotranspiration from thermal infrared, spectral reflectances and microwave data on two water stressed soybean crops. In a previous study, thermal infrared data alone were used to invert the model on plant water status parameters, while other parameters, related to canopy structure and soil surface water status, were prescribed from in situ measurements. In the present study, spectral reflectance and radar measurements were used to retrieve canopy structure parameters (LAI and vegetation height) and surface soil moisture by inverting radiative transfer models.

This research parameterized sensible and latent heat fluxes using Bowen Ratio (BR) also known as (Bowen Ratio Energy Balance), Aerodynamic Gradient (AG) and Aerodynamic Resistance (AR) techniques. Energy Residuum as well as Closure fractions were obtained and were used to assess closures for both direct and indirect fluxes. The data used were sourced from the Nigeria Micrometeorological Experiment (NIMEX-1) conducted within the transition period, 2004. The research site is located at the agricultural farm of Obafemi Awolowo University, Ile-Ife, (lat. 7.55o N, Long. 4.56o E), Nigeria. The results from each of the techniques were correlated with the measured heat fluxes in order to decide the best theoretical approach for the estimation of fluxes over the region of study. BR and AG compared moderately with the measured fluxes. BR scheme overestimated latent heat flux but underestimated sensible heat flux; AG technique slightly underestimated sensible heat flux but strongly underestimated latent heat fluxes with only few days of overestimation of sensible heat flux. AR, on the other hand, underestimated both latent and sensible heat fluxes. The obtained closure for AG proved much closer to the closure obtained from the direct measurement, BR has a better theoretical closure while AR failed as a whole. The trends as well as statistical analysis proved that AG performed best for sensible heat flux while BR presented good results for both latent heat flux and energy closure.

Hourly lysimetric and micrometeorological data taken over a grass surface at the Meteorological Research Unit, Cardington U.K. have been analysed. Measured quantities include sensible heat flux and friction velocity, obtained directly by the eddy correlation method, wind speeds and temperatures at various heights, screen dry and wet-bulb temperatures, net radiation, soil heat flux at several soil depths and lysimetric data taken hourly by a weighing type lysimeter. A temperature and wind difference allowed sensible heat and momentum fluxes determination by the profile method on an hourly basis. Estimation of latent and sensible heat flux by the resistance method is also obtained. Use is made of the energy balance at the surface and of the Monteith-Penman formula for estimation of the latent heat flux. Equations based on the Monin-Obukhov similarity theory are used to account for stability effects. Aerodynamic and surface resistances, necessary for the Penman-Monteith equation are also calculated from routinely measured meteorological data. The comparison of the various methods leads to the establishment of empirical relationships which correlate various quantities such as soil heat flux, resistances, evapotranspiration etc. to routinely measured meteorological data. The experiment was performed under near neutral or slightly unstable conditions. It appears that, under these environmental conditions the profile method for estimation of sensible heat flux and friction velocity is found to work excellently. Surface latent and sensible heat fluxes can also be described very well by the resistance method. The results also substantiate the accuracy of the Penman-Monteith equation for estimating latent heat flux and its insensitivity to adjustments for atmospheric stability.

The urban heat island (UHI) is formed due to higher local temperatures in urban areas than surrounding rural areas. Impervious surfaces and vegetation are the two major land covers that have the capability to modify energy fluxes. The Latent Heat flux, Sensible Heat Flux and the Ground Heat flux was obtained for Bangalore city and ratio of fluxes to net radiation was calculated. The mean of ratio of each flux to net radiation was calculated for each land cover and found that the urban and bare soil region shows higher values of <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\mathrm{G}/\mathrm{R}_{\mathrm{n}}$</tex> as compared to the vegetated land covers. The urban land cover showed higher values of <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\mathrm{H}/\mathrm{R}_{\mathrm{n}}$</tex> as compared to the natural land covers and the vegetation and waterbody showed higher values of <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\text{LE}/\mathrm{R}_{\mathrm{n}}$</tex> than other land covers. The relation of sensible heat flux with imperviousness and latent heat flux with vegetation was determined. The plot of <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\mathrm{H}/\mathrm{R}_{\mathrm{n}}$</tex> with ISF and <tex xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">$\text{LE}/\mathrm{R}_{\mathrm{n}}$</tex> with NDVI, both showed good positive relation.

The Tibetan Plateau (TP) has become a focus of strong scientific interest due to its role in the global water cycle and its reaction to climate change. Regional flux estimates of sensible and latent heat are important variables for linking the energy and hydrological cycles at the TP’s surface. Within this framework, a 3-year dataset (2008–2010) of eddy covariance measured turbulent fluxes was compiled from four stations on the TP into a standardised workflow: corrections and quality tests were applied using an internationally comparable software package. Second, the energy balance closure (C EB) was determined and two different closure corrections applied. The four stations (Qomolangma, Linzhi, NamCo and Nagqu) represent different locations and typical land surface types on the TP (high altitude alpine steppe with sparse vegetation, a densely vegetated alpine meadow, and bare soil/gravel, respectively). We show that the C EB differs between each surface and undergoes seasonal changes. Typical differences in the turbulent energy fluxes occur between the stations at Qomolangma, Linzhi and NamCo, while Nagqu is quite similar to NamCo. Specific investigation of the pre-monsoon, the Tibetan Plateau summer monsoon, post-monsoon and winter periods within the annual cycle reinforces these findings. The energy flux of the four sites is clearly influenced by the Tibetan Plateau monsoon. In the pre-monsoon period, sensible heat flux is the major energy source delivering heat to the atmosphere, whereas latent heat flux is greater than sensible heat flux during the monsoon season. Other factors affecting surface energy flux are topography and location. Land cover type also affects surface energy flux. The energy balance residuum indicates a typically observed overall non-closure in winter, while closure (or ‘turbulent over-closure’) is achieved during the Tibetan Plateau summer monsoon at the Nagqu site. The latter seems to depend on ground heat flux, which is higher in the wet season, related not only to a larger radiation input but also to a thermal decoupling of dry soils. Heterogeneous landscape modelling using a MODIS product is introduced to explain energy non-closure.

This paper describes and interprets the 1987 data from Cabauw, the Netherlands, which can be used to test land surface schemes in stand-alone mode. The data are available from the authors for model development and research. It consists of half-hour averages of forcing data (wind, temperature, specific humidity at 20-m height, downward solar and thermal radiation, and precipitation) and validation data (net radiation, sensible heat flux, latent heat flux, ground heat flux, and soil temperature). To obtain a continuous time series of the forcing parameters and the surface energy fluxes, it was necessary to use a model to fill in the missing observations. The quality of the observations and the reliability of model data are assessed by exploiting the redundancy in the observations and by comparing the model output with the data when both are available. The monthly averages of sensible heat flux are believed to be accurate to within ±5 W m−2 and the monthly means of net radiation and latent heat flux...

The surface energy balance is a key factor influencing the ground thermal regime. With ongoing climate change, it is crucial to understand the interactions of the individual heat fluxes at the surface and within the subsurface layers as well as their relative impacts on permafrost thermal regime. A unique set of high-altitude meteorological measurements has been analysed to determine the energy balance at three mountain permafrost sites in the Swiss Alps (Murtèl-Corvatsch, Schilthorn and Stockhorn), where data is being collected since the late 1990s in collaboration with the Swiss Permafrost Monitoring Network (PERMOS). All stations are equipped with sensors for four-component radiation, air temperature, humidity, wind speed and direction as well as ground temperatures and snow height. The three sites differ considerably in their surface and ground material composition as well as their ground ice contents. The energy fluxes are calculated based on two decades of field measurements. While the determination of the radiation budget and the ground heat flux is comparatively straightforward (by the four-component radiation sensor and thermistor measurements within the boreholes), larger uncertainties exist for the determination of turbulent sensible and latent heat fluxes. Our results show that mean air temperature at Murtèl-Corvatsch (1997–2018, 2600 m asl.) is −1.66 °C and has increased by about 0.7 °C during the measurement period. At the Schilthorn site (1999–2018, 2900 m asl.) a mean air temperature of −2.60 °C with a mean increase of 1.0 °C was measured. The Stockhorn site (2003–2018, 3400 m asl.) recorded lower air temperatures with a mean of −6.18 °C and an increase of 0.7 °C. Measured net radiation, as the most important energy input at the surface, shows substantial differences with mean values of 30.59 W m−2 for Murtèl-Corvatsch, 32.40 W m−2 for Schilthorn and 6.91 W m−2 for Stockhorn. The calculated turbulent fluxes show values of around 7 to 13 W m−2 using the Bowen ratio method and 3 to 15 W m−2 using the bulk method at all sites. Large differences are observed regarding the energy used for melting of the snow cover: at Schilthorn a value of 8.46 W m−2, at Murtèl-Corvatsch of 4.17 W m−2 and at Stockhorn of 2.26 W m−2 is calculated reflecting the differences in snow height at the three sites. In general, we found considerable differences in the energy fluxes at the different sites. These differences may help to explain and interpret the causes of the varying reactions of the permafrost thermal regime at the three sites to a warming atmosphere. We recognize a strong relation between the net radiation and the ground heat flux. Our results further demonstrate the importance of long-term monitoring in order to better understand the impacts of changes in the surface energy balance components on the permafrost thermal regime. The dataset presented can be used to improve permafrost modelling studies aiming at e.g. advancing knowledge about permafrost thaw processes. The data presented and described in this study is available for download at the following site http://dx.doi.org/10.13093/permos-meteo-2021-01 (Hoelzle et al., 2021).

Abstract. The surface energy balance is a key factor influencing the ground thermal regime. With ongoing climate change, it is crucial to understand the interactions of the individual heat fluxes at the surface and within the subsurface layers, as well as their relative impacts on the permafrost thermal regime. A unique set of high-altitude meteorological measurements was analysed to determine the energy balance at three mountain permafrost sites in the Swiss Alps (Murtèl–Corvatsch, Schilthorn and Stockhorn), where data have been collected since the late 1990s in the framework of the Swiss Permafrost Monitoring Network (PERMOS). All stations are equipped with sensors for four-component radiation, air temperature, humidity, and wind speed and direction, as well as ground temperatures and snow height. The three sites differ considerably in their surface and ground material composition, as well as their ground ice contents. The energy fluxes were calculated based on two decades of field measurements. While the determination of the radiation budget and the ground heat flux is comparatively straightforward (by the four-component radiation sensor and thermistor measurements within the boreholes), larger uncertainties exist for the determination of turbulent sensible and latent heat fluxes. Our results show that mean air temperature at Murtèl–Corvatsch (1997–2018, 2600 m a.s.l.) is −1.66 ∘C and has increased by about 0.8 ∘C during the measurement period. At the Schilthorn site (1999–2018, 2900 m a.s.l.) a mean air temperature of −2.60 ∘C with a mean increase of 1.0 ∘C was measured. The Stockhorn site (2003–2018, 3400 m a.s.l.) recorded lower air temperatures with a mean of −6.18 ∘C and an increase of 0.5 ∘C. Measured net radiation, as the most important energy input at the surface, shows substantial differences with mean values of 30.59 W m−2 for Murtèl–Corvatsch, 32.40 W m−2 for Schilthorn and 6.91 W m−2 for Stockhorn. The calculated turbulent fluxes show values of around 7 to 13 W m−2 using the Bowen ratio method and 3 to 15 W m−2 using the bulk method at all sites. Large differences are observed regarding the energy used for the melting of the snow cover: at Schilthorn a value of 8.46 W m−2, at Murtèl–Corvatsch 4.17 W m−2 and at Stockhorn 2.26 W m−2 are calculated, reflecting the differences in snow height at the three sites. In general, we found considerable differences in the energy fluxes at the different sites. These differences help to explain and interpret the causes of a warming atmosphere. We recognise a strong relation between the net radiation and the ground heat flux. Our results further demonstrate the importance of long-term monitoring to better understand the impacts of changes in the surface energy balance components on the permafrost thermal regime. The dataset presented can be used to improve permafrost modelling studies aiming at, for example, advancing knowledge about permafrost thaw processes. The data presented and described here are available for download at the following site: https://doi.org/10.13093/permos-meteo-2021-01 (Hoelzle et al., 2021).

Land surface models involve a large number of interdependent parameters that affect the physics of how surface energy fluxes are partitioned between latent heat, sensible heat, net radiative, and ground heat fluxes. The goal of an optimal parameter and uncertainty analysis of a land surface model is to identify a range of parameter sets that enable model predictions to be bounded within observational uncertainties. Here we apply Bayesian stochastic inversion (BSI) using very fast simulated annealing (VFSA) to identify parameter sets of the Chameleon surface model (CHASM) land surface model that are consistent with the uncertainty limits ascribed to a high‐quality data set collected from Cabauw, Netherlands. These results are compared to the parameter sets obtained through the multicriteria (MC) approach. All analyses evaluate model performance against daily and monthly mean observations of sensible, latent, and ground heat fluxes. BSI and MC identify similar “best fit” model parameter sets that improve CHASM performance over default parameter settings. The three most important CHASM parameters at Cabauw are minimum stomatal resistance, vegetation roughness length, and vegetation fraction cover. BSI is based on a Bayesian inference model such that that it expresses uncertainty in terms of a posterior probability density function, different moments of which provide information about parameter means and covariances. Although MC gives a range of possible optimal parameters through the concept of a Pareto set, we found that these ranges did not provide a consistent or representative view of the uncertainty within the observational data. The BSI algorithm in the current study is particularly efficient in that it only requires about double the number of model evaluations than the MC algorithm. This is a substantial saving over other more accurate methods to evaluate uncertainty such as the Metropolis/Gibbs' sampler that requires at least 40 times more computations than the BSI algorithm to obtain similar results.