Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

Remote estimation of thermal inertia and soil heat flux for bare soil

A method is presented which allows thermal inertia (the soil heat capacity times the square root of the soil thermal diffusivity, C(h)rootD(h)), to be estimated remotely from micrometeorological observations. The method uses the drop in surface temperature, T-s, between sunset and sunrise, and the average night-time net radiation during that period, for clear, still nights. A Fourier series analysis was applied to analyse the time series of T-s . The Fourier series constants, together with the remote estimate of thermal inertia, were used in an analytical expression to calculate diurnal estimates of the soil heat flux, G. These remote estimates of C(h)rootD(h) and G compared well with values derived from in situ sensors. The remote and in situ estimates of C(h)rootD(h) both correlated well with topsoil moisture content. This method potentially allows area-average estimates of thermal inertia and soil heat flux to be derived from remote sensing, e.g. METEOSAT Second Generation, where the area is determined by the sensor's height and viewing angle. (C) 2003 Elsevier B.V. All rights reserved.

Spatio-temporal surface soil heat flux estimates from satellite data: results for the AMMA experiment at the Fakara (Niger) supersite

This paper describes a method that employs Earth Observation (EO) data to calculate spatiotemporal estimates of soil heat flux, G, using a physically-based method (the Analytical Method). The method involves a harmonic analysis of land surface temperature (LST) data. It also requires an estimate of near-surface soil thermal inertia; this property depends on soil textural composition and varies as a function of soil moisture content. The EO data needed to drive the model equations, and the ground-based data required to provide verification of the method, were obtained over the Fakara domain within the African Monsoon Multidisciplinary Analysis (AMMA) program. LST estimates (3 km × 3 km, one image 15 min−1) were derived from MSG-SEVIRI data. Soil moisture estimates were obtained from ENVISAT-ASAR data, while estimates of leaf area index, LAI, (to calculate the effect of the canopy on G, largely due to radiation extinction) were obtained from SPOT-HRV images. The variation of these variables over the Fakara domain, and implications for values of G derived from them, were discussed. Results showed that this method provides reliable large-scale spatiotemporal estimates of G. Variations in G could largely be explained by the variability in the model input variables. Furthermore, it was shown that this method is relatively insensitive to model parameters related to the vegetation or soil texture. However, the strong sensitivity of thermal inertia to soil moisture content at low values of relative saturation (<0.2) means that in arid or semi-arid climates accurate estimates of surface soil moisture content are of utmost importance, if reliable estimates of G are to be obtained. This method has the potential to improve large-scale evaporation estimates, to aid land surface model prediction and to advance research that aims to explain failure in energy balance closure of meteorological field studies.

Radiation fluxes, soil heat flux, air temperature, soil temperature, soil moisture and vegetation at dwarf shrubs and wet sedges in Kytalyk, NE Siberia

Vegetation changes, such as shrub encroachment and wetland expansion, have been observed in many Arctic tundra regions. These changes feed back to permafrost and climate. Permafrost can be protected by soil shading through vegetation as it reduces the amount of solar energy available for thawing. Regional climate can be affected by a reduction in surface albedo as more energy is available for atmospheric and soil heating. Here, we compared the shortwave radiation budget of two common Arctic tundra vegetation types dominated by dwarf shrubs (Betula nana) and wet sedges (Eriophorum angustifolium) in North-East Siberia. We measured time series of the shortwave and longwave radiation budget above the canopy and transmitted radiation below the canopy. Additionally, we quantified soil temperature and heat flux as well as active layer thickness. The mean growing season albedo of dwarf shrubs was 0.15 ± 0.01, for sedges it was higher (0.17 ± 0.02). Dwarf shrub transmittance was 0.36 ± 0.07 on average, and sedge transmittance was 0.28 ± 0.08. The standing dead leaves contributed strongly to the soil shading of wet sedges. Despite a lower albedo and less soil shading, the soil below dwarf shrubs conducted less heat resulting in a 17 cm shallower active layer as compared to sedges. This result was supported by additional, spatially distributed measurements of both vegetation types. Clouds were a major influencing factor for albedo and transmittance, particularly in sedge vegetation. Cloud cover reduced the albedo by 0.01 in dwarf shrubs and by 0.03 in sedges, while transmittance was increased by 0.08 and 0.10 in dwarf shrubs and sedges, respectively. Our results suggest that the observed deeper active layer below wet sedges is not primarily a result of the summer canopy radiation budget. Soil properties, such as soil albedo, moisture, and thermal conductivity, may be more influential, at least in our comparison between dwarf shrub vegetation on relatively dry patches and sedge vegetation with higher soil moisture.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.

Patterns of water and heat flux across a biome gradient from tropical forest to savanna in Brazil

We investigated the seasonal patterns of water vapor and sensible heat flux along a tropical biome gradient from forest to savanna. We analyzed data from a network of flux towers in Brazil that were operated within the Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA). These tower sites included tropical humid and semideciduous forest, transitional forest, floodplain (with physiognomies of cerrado), and cerrado sensu stricto. The mean annual sensible heat flux at all sites ranged from 20 to 38 Wm(-2), and was generally reduced in the wet season and increased in the late dry season, coincident with seasonal variations of net radiation and soil moisture. The sites were easily divisible into two functional groups based on the seasonality of evaporation: tropical forest and savanna. At sites with an annual precipitation above 1900 mm and a dry season length less than 4 months (Manaus, Santarem and Rondonia), evaporation rates increased in the dry season, coincident with increased radiation. Evaporation rates were as high as 4.0 mm d(-1) in these evergreen or semidecidous forests. In contrast, ecosystems with precipitation less than 1700 mm and a longer dry season (Mato Grosso, Tocantins and Sao Paulo) showed clear evidence of reduced evaporation in the dry season. Evaporation rates were as low as 2.5 mm d(-1) in the transitional forests and 1 mm d(-1) in the cerrado. The controls on evapotranspiration seasonality changed along the biome gradient, with evaporative demand (especially net radiation) playing a more important role in the wetter forests, and soil moisture playing a more important role in the drier savannah sites.

Soil CO2 flux: a method comparison of closed static chambers in a sugarcane field

A large variety of techniques have been used to measure soil CO2 released from the soil surface, and much of the variability observed between locations must be attributed to the different methods used by the investigators. Therefore, a minimum protocol of measurement procedures should be established. The objectives of this study were (a) to compare different absorption areas, concentrations and volumes of the alkali trapping solution used in closed static chambers (CSC), and (b) to compare both, the optimized alkali trapping solution and the soda-lime trapping using CSC to measure soil respiration in sugarcane areas. Three CO2 absorption areas were evaluated (7; 15 and 20 % of the soil emission area or chamber); two volumes of NaOH (40 and 80 mL) at three concentrations (0.1, 0.25 and 0.5 mol L-1). Three different types of alkaline traps were tested: (a), 80 mL of 0.5 mol L-1 NaOH in glass containers, absorption area 15 % (V0.5); (b) 40 mL of 2 mol L-1 NaOH retained in a sponge, absorption area 80 % (S2) and (c) 40 g soda lime, absorption area 15 % (SL). NaOH concentrations of 0.5 mol L-1 or lower underestimated the soil CO2-C flux or CO2 flux. The lower limit of the alkali trap absorption area should be a minimum of 20 % of the area covered by the chamber. The 2 mol L-1 NaOH solution trap (S2) was the most efficient (highest accuracy and highest CO2 fluxes) in measuring soil respiration.

Design of an amplifier and an instrumentation setup for a harsh environment heat flux measurement application – case axial flux electric machine

Measurement is a tool for researching. Therefore, it is important that the measuring process is carried out correctly, without distorting the signal or the measured event. Researches of thermoelectric phenomena have been focused more on transverse thermoelectric phenomena during recent decades. Transverse Seebeck effect enables to produce thinner and faster heat flux sensor than before. Studies about transverse Seebeck effect have so far focused on materials, so in this Master’s Thesis instrumentation of transverse Seebeck effect based heat flux sensor is studied, This Master’s Thesis examines an equivalent circuit of transverse Seebeck effect heat flux sensors, their connectivity to electronics and choosing and design a right type amplifier. The research is carried out with a case study which is Gradient Heat Flux Sensors and an electrical motor. In this work, a general equivalent circuit was presented for the transverse Seebeck effect-based heat flux sensor. An amplifier was designed for the sensor of the case study, and the solution was produced for the measurement of the local heat flux of the electric motor to improve the electromagnetic compatibility.

Heat flux measurement inside internal combustion engine with gradient heat flux sensor.

This master’s thesis is devoted to study different heat flux measurement techniques such as differential temperature sensors, semi-infinite surface temperature methods, calorimetric sensors and gradient heat flux sensors. The possibility to use Gradient Heat Flux Sensors (GHFS) to measure heat flux in the combustion chamber of compression ignited reciprocating internal combustion engines was considered in more detail. A. Mityakov conducted an experiment, where Gradient Heat Flux Sensor was placed in four stroke diesel engine Indenor XL4D to measure heat flux in the combustion chamber. The results which were obtained from the experiment were compared with model’s numerical output. This model (a one – dimensional single zone model) was implemented with help of MathCAD and the result of this implementation is graph of heat flux in combustion chamber in relation to the crank angle. The values of heat flux throughout the cycle obtained with aid of heat flux sensor and theoretically were sufficiently similar, but not identical. Such deviation is rather common for this type of experiment.

Thermalhydraulic behavior of electrically heated rod during a critical heat flux transient

In nuclear reactors, the occurrence of critical heat flux leads to fuel rod overheating with clad fusion and radioactive products leakage. To predict the effects of such phenomenon, experiments are performed using electrically heated rods to simulate operational and accidental conditions of nuclear fuel rods. In the present work, it is performed a theoretical analysis of the drying and rewetting front propagation during a critical heat flux experiment, starting with the application of an electrical power step from steady state condition. After the occurrence of critical heat flux, the drying front propagation is predicted. After a few seconds, a power cut is considered and the rewetting front behavior is analytically observed. Studies performed with various values of coolant mass flow rate show that this variable has more influence on the drying front velocity than on the rewetting one.

Understanding electron heat flux signatures in the solar wind

Suprathermal electrons (E > 80 eV) carry heat flux away from the Sun. Processes controlling the heat flux are not well understood. To gain insight into these processes, we model heat flux as a linear dependence on two independent parameters: electron number flux and electron pitch angle anisotropy. Pitch angle anisotropy is further modeled as a linear dependence on two solar wind components: magnetic field strength and plasma density. These components show no correlation with number flux, reinforcing its independence from pitch angle anisotropy. Multiple linear regression applied to 2 years of Wind data shows good correspondence between modeled and observed heat flux and anisotropy. The results suggest that the interplay of solar wind parameters and electron number flux results in distinctive heat flux dropouts at heliospheric features like plasma sheets but that these parameters continuously modify heat flux. This is inconsistent with magnetic disconnection as the primary cause of heat flux dropouts. Analysis of fast and slow solar wind regimes separately shows that electron number flux and pitch angle anisotropy are equally correlated with heat flux in slow wind but that number flux is the dominant correlative in fast wind. Also, magnetic field strength correlates better with pitch angle anisotropy in slow wind than in fast wind. The energy dependence of the model fits suggests different scattering processes in fast and slow wind.