Next Article in Journal
Analysis of Tillage Depth and Gear Selection for Mechanical Load and Fuel Efficiency of an Agricultural Tractor Using an Agricultural Field Measuring System
Next Article in Special Issue
Agricultural Robotics for Field Operations
Previous Article in Journal
Resource Allocation for Machine-Type Communication of Energy-Harvesting Devices in Wi-Fi HaLow Networks
Previous Article in Special Issue
An All-Solid-State Nitrate Ion-Selective Electrode with Nanohybrids Composite Films for In-Situ Soil Nutrient Monitoring
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 20
Issue 9
10.3390/s20092444
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal)

by
William Frodella
1,*,
Giacomo Lazzeri
1,
Sandro Moretti
1,
Jacob Keizer
2 and
Frank G. A. Verheijen
2
1
Department of Earth Sciences, University of Firenze, Via La Pira 4, 50121 Firenze, Italy
2
Department of Environment and Planning, University of Aveiro, Campus Universitário de Santiago, 3810-193 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(9), 2444; https://doi.org/10.3390/s20092444
Submission received: 27 March 2020 / Revised: 14 April 2020 / Accepted: 20 April 2020 / Published: 25 April 2020
(This article belongs to the Special Issue Advanced Sensors in Agriculture)
Browse Figures
Versions Notes

Abstract

:
The soil surface albedo decreases with an increasing biochar application rate as a power decay function, but the net impact of biochar application on soil temperature dynamics remains to be clarified. The objective of this study was to assess the potential of infrared thermography (IRT) sensing by monitoring soil surface temperature (SST) with a high spatiotemporal and thermal resolution in a scalable agricultural application. We monitored soil surface temperature (SST) variations over a 48 h period for three treatments in a vineyard: bare soil (plot S), 100% biochar cover (plot B), and biochar-amended topsoil (plot SB). The SST of all plots was monitored at 30 min intervals with a tripod-mounted IR thermal camera. The soil temperature at 10 cm depth in the S and SB plots was monitored continuously with a 5 min resolution probe. Plot B had greater daily SST variations, reached a higher daily temperature peak relative to the other plots, and showed a faster rate of T increase during the day. However, on both days, the SST of plot B dipped below that of the control treatment (plot S) and biochar-amended soil (plot SB) from about 18:00 onward and throughout the night. The diurnal patterns/variations in the IRT-measured SSTs were closely related to those in the soil temperature at a 10 cm depth, confirming that biochar-amended soils showed lower thermal inertia than the unamended soil. The experiment provided interesting insights into SST variations at a local scale. The case study may be further developed using fully automated SST monitoring protocols at a larger scale for a range of environmental and agricultural applications.

Graphical Abstract

1. Introduction

In earth science applications, land surface temperature (LST) is one of the key physical parameters of land surface processes, from local through to global scales [1,2]. Soil surface temperature (SST) in particular is an important factor that controls the energy exchange processes between the land surface and the atmosphere, representing a fundamental input parameter in meteorological and climatological studies, as well as in hydrological and agricultural analyses [3,4,5,6,7]. SST measures are usually carried out using thermocouples and radiometers [3,7,8]; nevertheless, these devices can present long-term issues regarding logistics, access, and technician costs. In the last few decades, there has been a strong interest in developing methodologies to measure LST and SST using remote sensing (RS) techniques via spaceborne [9,10,11,12,13], airborne [14,15,16,17,18], or ground-based sensors [19,20,21,22].
Since the 1980s, thermal sensors from terrestrial, airborne, and spaceborne platforms have been used in a variety of geo-environmental applications [21,22,23,24,25,26,27]. Although the advantages of satellite remote sensing in thermal applications appear to be many, the calibration and validation of the retrieved data is complex and requires complicated correction procedures [28,29,30,31]. To alleviate these issues, ground-based thermal sensors can offer a proper solution for the analysis of the local energy balance and micro-meteorological conditions that can be used as a proxy in large-scale land surface and hydrology experiments in the context of a growing interest for global change research [32,33]. Thermal imaging, or infrared thermography (IRT), offers non-contact, wide-area detection of surface temperatures [34,35,36] and can be used as an alternative or complement to conventional inspection technologies in a wide variety of geo-environmental and agricultural applications [37,38,39,40,41,42,43,44]. Thermal imaging ground-based sensors, such as thermal cameras, have undergone a significant increase in use in recent years, thanks to the technological development of focal plane array uncooled microbolometer sensors, which in turn have allowed for improvement in the accuracy, spatial resolution, and fast measurement of thermal data [45].
Soils provide crucial ecosystem services, e.g., supporting human health and habitation, storing carbon, and regulating emissions of greenhouse gases [8]. Climate change and increasing human populations are threatening soil resources and agro-ecological balances, with important implications for food security. Sustainable biochar application to soils [46] has been suggested as a means to mitigate climate change via C sequestration, as well as to improve other ecosystem services, such as crop yields [47,48,49], greenhouse gas fluxes [50,51,52], soil water storage [53,54], soil erosion [55], and soil remediation [56]. However, the effects of biochar application on soil temperature dynamics are not well understood. Yet, soil temperature drives many physicochemical and especially biological soil processes, e.g., soil organic carbon decomposition [57], with climate feedbacks. Therefore, considering the expected changes in global air temperature [58] and precipitation [59], there is an urgent need for a better understanding of both the impacts of increasing soil temperatures in general and the application of biochar to soils as a mitigation and/or adaptation measure. Some authors have shown that soil surface albedo decreases with increasing biochar application rate as a power decay function, with stronger decreases in albedo for dry soil than wet soil [60]. At 50 t/ha, a common biochar application rate, the surface albedo of dry soil decreased by around 30% and that of wet soil by around 18%. All other factors remaining equal, biochar application can be expected to increase SSTs. However, this has not been confirmed by experimental studies so far, with existing studies finding no effect [61], an effect at the surface but not at a 7.5 cm depth [62], or a decrease in soil thermal conductivity but not in soil temperature [63]. Therefore, the impact of biochar application on soil temperature dynamics remains to be clarified [64,65,66]. The lack of available data and their inconsistency is, at least in part, caused by the difficulty in monitoring soil temperature with soil probes in a representative manner given the typically pronounced spatial heterogeneity in key soil variables. The objective of this study was to assess the potential of IRT sensing to overcome this bottleneck, thanks to the scalability of its results, by monitoring SST with high spatial and thermal resolutions in an agricultural application. More specifically, a tripod-mounted IRT camera was used to measure SST with a high temporal accuracy (every 30 min over a 48 h surveying campaign) in a vineyard test site where a biochar field experiment had been set up. Three selected plots with different treatments were analyzed, two of which had sub-surface T probe instruments (5 min measuring interval). A two-spatial-scale approach was adopted: (i) the single-plot scale (given by a 42.5 × 32 cm IRT image); and (ii) the within-plot scale, given by biochar-covered and bare soil spots (from centimetric to millimetric-scale). For this second approach, a semi-automatic procedure in a geographical information system (GIS) environment was adopted for the analysis and the integration of the large IR image dataset. The presented approach can be expanded to the wider agricultural and geo-environmental domain, as well as to larger spatial scales by using IRT sensors based on drone Unmanned Aerial Vehicles platforms. Furthermore, general operative recommendations can be elaborated upon for a fully automated protocol based on an innovative approach using high-resolution spatiotemporal IRT data.

2. Materials and Methods

2.1. IRT: Theoretical Principles

IRT (or thermal imaging) is the branch of RS dealing with measuring the thermal radiation of the Earth’s surface features from a distance [45]. Nowadays, most thermal cameras use a focal plane array (FPA) uncooled microbolometer type of detector, mainly because cooling systems make thermal cameras more expensive and less manageable [66]. Microbolometers sensors operate through their response to radiant energy, through changes in their temperature and electrical conductivity, depending on the intensity of absorption of incident radiation [67]. IRT is accomplished by using IR-calibrated cameras (thermal cameras), whose sensors are capable of mapping and monitoring the spatial distribution of thermal radiation [68,69]. An IRT survey produces a matrix of pixel values acquired by the thermal camera array sensor (called a “thermogram” or “thermographic image”), which then undergoes a correction of the sensitive parameters (object emissivity, path length, air temperature, and humidity) by the built-in processor to create a radiant temperature map of the investigated scenario [37,45]. Emissivity can be defined as the ratio of the object’s efficiency at radiating energy compared to a black body (i.e., ε, which is the radiant exitance of an object at a given temperature over the radiant exitance of a black body at the same temperature [70]). While according to Kirchoff’s laws on radiation, black bodies have an emissivity equal to 1, in most remote sensing applications, the objects in question are real bodies (also called grey bodies), e.g., in the 6–14 µm wavelength range, water bodies have very high emissivity values (0.98 < ε < 0.99), whereas soils and rocks have slightly lower values (ε ≥ 0.8) [71]. The concept of emissivity is crucial in IRT because it allows us to relate the IR radiation to an object’s temperature based on Stefan–Bolzmann’s law [45]:
W = ε σ T 4
where W is the total radiant exitance from an object’s surface, expressed in W·m2; σ is the Stefan–Boltzmann constant (5.6697 × 10−8 W·m−2·K−4); and T is the absolute temperature (expressed in K) of the emitting material. Thermal sensors measure the radiant temperatures of objects. The true kinetic temperature of an object (Tkin), generated by the vibration of molecules) can be estimated from its radiant temperature (Trad), provided the object’s emissivity is known [66]:
Trad = ε1/4 Tkin
Thermal inertia is a key thermal property that controls the diurnal and seasonal surface temperature variations. It can be thought of as a surface’s ability to conduct and store heat during the day and re-radiate it away at night in the thermal infrared portion of the spectrum [27,66]. Thermal inertia is defined as [72]:
I = (KρC)1/2
where K is the thermal conductivity, ρ is the bulk density, and C is the specific heat capacity of the surface. For rocks and soils, the bulk density (ρ) does not vary significantly, and neither does the specific heat capacity (C) for a dry surface; therefore, variations in thermal inertia are primarily due to variations in thermal conductivity [72]. In remote sensing applications, thermal inertia represents a complex combination of particle size, rock abundance, bedrock outcropping, and the degree of induration [73]; a rough approximation to thermal inertia is sometimes obtained from the amplitude of the diurnal temperature curve (i.e., maximum minus minimum surface temperature) [74]. The thermal inertia has primary effects on the shape of a surface’s diurnal temperature curve: a low thermal inertia surface will experience higher peak temperatures and lower minimum temperatures than a high thermal inertia surface because the low thermal inertia surface becomes heated and cools down faster [75,76]. This is in line with Newton’s law of cooling, which states that the rate of change in the temperature of a body is proportional to the difference between the body’s temperature and the ambient temperature (meaning that the cooling of the body is faster when the temperature difference is higher) [77].
The built-in camera processor and the IRT software (e.g., that from FLIR Systems Inc. [69]) allows for visualizing the thermal pattern image using a false-color temperature scale. The radiant temperature values acquired through an IRT survey do not correspond exactly to the temperatures emitted by the targeted object due to the attenuation and scattering by air molecules, reflection, and emissivity [78,79]. The correction of IRT data requires the calibration of object emissivity values, air temperature and humidity, and the distance between the sensor and the target. In most applications, the emissivity of targeted objects is estimated based on values from the literature; alternatively, there are empirical methods that can be used to estimate emissivity values with a reasonable precision [79,80,81,82,83,84].

2.2. Study Area and Biochar Treatment Plots

The study area was a vineyard located within the urban area of the Anadia municipality (Aveiro District, central Portugal) at 40°26′25.91″ N and 8°26′22.90″ W. The climate is warm temperate and mesic [85], with a mean annual precipitation, air temperature, and relative humidity of 700 mm, 15 °C, and ≈ 60%, respectively. At the time of this study (from 1–3 August 2013), the vineyard consisted of mature, five-year-old vines of the Sauvignon Blanc variety under conventional management and without shading. Within the vineyard, an experimental site had been set-up in March 2012 to test the impacts of biochar amendment on key topsoil properties, as well as vine hydric stress. The setup consisted of a total of 12 experimental plots, divided over three blocks, with one control and two biochar application treatments per block. A soil pit was dug in the middle of the experimental site to determine the soil type and to sample the topsoil at a 0–15 cm depth for determining the soil texture (the pipette method), organic matter content (loss-on-ignition method), dry bulk density (using a 5 cm diameter core), and pH (using an electrode in a 1:10 soil:water suspension). For more details, please see Verheijen et al. [55]. The soil was a well-draining Cambisol of 60–80 cm depth with a sandy loam topsoil texture (sand-silt-clay contents of 70%-16%-14%, respectively). The organic matter content of the topsoil was 3.9% (≈2.3% Soil Organic Carbon) [86,87]. The topsoil bulk density was 1.25 g·cm−3 and the topsoil pH was 6.0. An HCl test revealed no presence of CaCO3, and a molarity of an ethanol droplet (MED) test [88,89] showed the air-dried soil to be wettable. The biochar that was applied in the experimental site was produced from mixed wood sievings in a Pyreg® 500 III pyrolysis unit (Dörth, Germany), 620 °C maximum temperature, 20 min duration, and a H/C ratio of 0.18 (Swiss-Biochar GmbH, Basel, Switzerland). The median particle size was 29.5 µm with the following particle size distribution (w/w): 4% (<0.1 mm), 25% (0.1–0.5 mm), 34% (0.5–2 mm), 37% (>2 mm). The pH was 10.4 and the ash content was 28% (full characterization in Prodana et al. [90]).

2.3. Experimental Setup

The rationale behind the present study was to assess the effect of a biochar amendment on the vineyard SST and to relate the SST patterns to simultaneous records of air temperature and relative humidity data from a standard weather station, as well as with sub-surface soil temperature data from automatic sensor probes. For this study, the three plots of one of the vineyard blocks were selected (Figure 1). These plots were characterized by different treatments: the control plot was dominated by bare soil (namely plot “S”); one plot had biochar (a soil amendment) applied to the surface, providing a 100% soil cover (namely plot “B”); while the third plot had biochar incorporated into the topsoil, providing a soil cover of 10% (namely plot “SB”). The original biochar application rates were 40 and 4 t·ha−1 for plots B and SB, respectively. A tripod-mounted IRT camera was used to measure the SST with a high temporal accuracy (30 min) on all three plots over a 48 h surveying campaign, starting with plot B and ending with plot SB, from 08:00 on 1 August 2013 to 09:00 on 3 August 2013. The instrument employed in this work was a tripod-mounted SC620 thermal camera (FLIR Systems, Wilsonville, OR, USA) [68]. This device can measure electromagnetic radiation in the Long Wawelenght InfraRed band, with a 640 × 480 pixel detector size and a spatial resolution of 0.65 mrad (see Table 1 for further specifications). The IR camera was mounted on a tripod at a 1.3 m height above the soil surface, leading to an image geometrical resolution of 0.85 mm (Figure 1 and Figure 2). Sunrise and sunset during this period were at 06:24 and 20:41, respectively. Furthermore, plots S and SB were had two T-probe instruments at a 10 cm depth and were ≈ 15 cm apart, including (i) a Degacon 5TE probe (Pullman, WA, USA), which combines soil moisture content v/v (capacitance sensor) with electrical conductivity and temperature (0.1 °C resolution), and (ii) a Degacon MPS2 probe, which combines soil water potential (dielectric water potential) with temperature (0.1 °C resolution). The readings of the probes were recorded at 5 min intervals using Decagon EM50 loggers (Figure 2). In this study, the mean values of the two probes (ST10) per plot were used. Local weather data (air temperature, relative humidity, rainfall) were acquired using a pocket thermohygrometer (OM-EL-USB-2-LCD datalogger; Omega Engineering, Manchester, UK) located at a 1 m height above the soil and acquiring with 5 min intervals (Figure 2), and using the Instituto Português do Mar e da Atmosfera (IPMA) Anadia weather station [91] located approximately 190 m away from the test site (Table 2).
The acquired SST dataset was compared to the soil sub-surface temperature data set from the soil probes, as well as to local weather conditions. Based on these weather data, the monitoring period was divided into two sub-periods: a first one from 08:00 on 1 August to 02:00 on 2 August, which was characterized by dry atmospheric conditions and clear sky, and a second one from 02:00 on 2 August to 09:00 on 3 August, which started with a low-intensity precipitation event that lasted approximately 2 h and produced 0.9 mm of rainfall (Figure 3).
Emissivity values of the topsoil and biochar were retrieved using laboratory tests following the “reference emissivity material method,” as given in References [82,83,84]. The method briefly consists of attaching a reference emitter to the samples and heating both to 20 °C above room temperature; the unknown emissivity of the samples is then measured through the IRT software using a calibration procedure based on the known reference value. The obtained results (ε = 0.96 for the soil and ε = 0.98 for the biochar) agreed well with the values in the existing literature [79,80,92,93,94,95]. These emissivity values were then used to correct the plot thermograms using the FLIR Tools + software [69]. The latter was also used for the single-plot scale analysis to measure the maximum, minimum, and mean SST for the acquired IRT dataset by using the square tool function on the plot area of the thermograms (Figure 1f).
To analyze the effect of biochar grains on the SST of the SB plots, a “within-plot scale” analysis was performed on the whole IRT dataset using a semi-automatic procedure in a GIS environment. Regions of interest (ROIs) from B1 to B6 for biochar and from S1 to S3 for soil were manually selected on the SB thermogram. The criterion used to establish these areas was determined by the necessity of having SST only representing the investigated material. Therefore, biochar clusters were chosen to be represented by a coalescence of three or more grains. Similarly, the bare soil sectors were selected in a way that avoided the presence of rock clasts grains/pebbles or vegetation. These criteria were used to avoid the thermal influence caused by different materials (e.g., the cooling effect of vegetation due to evapotranspiration). To analyze the SST in detail, the time series were analyzed on three selected biochar-clustered grains and bare soil sectors each within the SB plot, and the thermograms were manually converted to a raster format, spatially co-registered, and analyzed using the ESRI ArcGIS software “zonal statistics tool” [96] (Figure 4). The tool was run in batch mode to automatically extract the SST values from the plot’s entire thermogram dataset, both for the thermograms and the single polygon representing the ROIs. For the pixels located in the same area, the tool computes the mean temperature, reports the minimum and maximum temperatures, and their range. As these values are extracted from all the thermograms, they were recorded in a table for each polygon separately to allow for the monitoring of SST data for the specific ROI.

3. Results

3.1. SSTs from IRTs

The analyzed IRT dataset showed a double-humped daily SST distribution, characterized by four heating and cooling phases during the 48 h monitoring period (Figure 5 and Figure 6). The air temperature (air T) peaked around 14:00 on day 1 (36.5 °C), while the maximum SSTs for all the plots peaked about an hour later, i.e., at 15:04 (73.5 °C for plot B, 62.4 °C for plot S, and 63.4 °C for Plot SB). After these peaks, the SSTs dropped quickly for all the plots from around 16:00 to 19:00 (≈ 30 °C), when the differences with air T were smallest (this was because direct sunlight stopped hitting the soil surface from 17:00). The SST peaks continued to drop during the first night down to 19 °C at 06:30 on day 2, followed by a sharp increase when direct sunlight started reaching the plots again, starting from 08:30. During day 2, the effects of rainfall from 02:30 to 04:30, as well as the increased air humidity (Figure 3) and cloud cover, were clearly visible in the cooling of all the treatment plots.
The air T peaked at 12:54 (32.5 °C) and the peak SSTs of the treatments were staggered, i.e., plots S and SB at 14:34 (55.4 and 57.3 °C, respectively) and plot B an hour later at 15:34 (66.3 °C). Moderate and short T drops were visible in both daily cycles from 10:30 to around 12:30, resulting from the vine canopy shading. These drops were more pronounced in the maximum SSTs than in the minimum SSTs (T min) and were the least marked in the mean values (T mean). On day 2, this drop lasted longer than on day 1, i.e., from 11:30 to 13:00 as opposed to from 12:30 to 13:30. This was probably due to a cooling effect from the early morning rainfall.
During the nighttime of day 2, the difference between the air T and SST increased again from 19:00 to 05:00, while the rain during the first night quickly equalizes the air and SSTs (from 03:00 to 08:30). By comparing the shape of the distribution for the air T and SST, it can be seen how the former had a more regular pattern of heating and cooling, and how the variations in sun heating required a specific time to be observed on the distribution. Compared with the air T, the latter was evidently sharper and less homogeneous, presenting steep sides and peak values when related to the distribution of the air T. The minimum values of both the air T and SST occurred between 06:30 and 07:00 on both surveying days (day 1: 17.9 °C for the air T and 21.5, 21.4, and 21.1 °C for plots S, B, and SB, respectively; day 2: 13 °C for the air T and 17.1, 17.7, and 17.6 °C for plots S, B, and SB, respectively). This coincided with the time just after sunrise and before sunlight hit the plots again. The only time when the SST was lower than the air T was from 11:30 to 12:00, i.e., during the time when the vine canopy shaded the soil.

3.2. Soil Temperatures at a 10 cm Depth

At a 10 cm depth, the SB plot reached a T max that was 0.9 °C higher than the S plot around 18.00 on day 1 and 19.00 for day 2 (3 h after the air T and 1.5 h after the soil surface T; Figure 7). This indicates that plot SB warmed up and cooled down more quickly. From its peak, the T of plot SB seemed to decrease more quickly, leading to a slightly lower T from 20:00 onward. The stronger cooling seemed even more obvious for the canopy shading time in the morning, particularly on day 2. Additionally, by comparing the elapsed time between the SST and the air T, on day 2, the time interval was greater, and the peaks had wider bases compared to day 1. Figure 8 shows the differences in SST, as well as the soil temperature at a 10 cm soil depth between the control plot (S probe) and the plot with biochar mixed into the topsoil at an agronomically relevant 40 t·ha−1 application rate (SB probe). Positive values indicate a heating effect, with negative values indicating a cooling effect.
It can be observed that the temperature changes at a 10 cm soil depth occurred with a slight lag to SST and were several times smaller in magnitude. The air temperature did not appear to drive biochar effects on soil temperatures.

3.3. SST Raster Analysis for a Spot in Plot SB

The SST data from the six biochar (grains/clusters) and three bare soil spots in the SB plot are shown in Figure 9. The biochar areas (B label in Figure 9) show a similar thermal pattern, as well as the bare soil sectors (S label in Figure 8). The different thermal inertias of the materials (lower in the biochar, higher in the soil) were clearly visible in the daily peaks of the SST pattern (around 13 °C of variation in the 13:00–15:00 interval) and at night time (around 5 °C in the 01:30–05:00 interval), when the lowest temperatures were measured.
This was confirmed by analyzing the mean SST variations of biochar grains/clusters and bare soil sectors (Figure 10), where biochar showed stronger variations in the heating and cooling cycles than bare soil; the biochar’s lower thermal inertia was due to the material’s coarser texture and grain size, influencing the pore size distribution, and therefore, the air circulation. It can be seen how its temperature rose more quickly around the zenith, i.e., during the strongest solar radiation hours. Similarly, the temperature dropped after the solar irradiation decreased because of its lower thermal capacity compared to the soil. To provide a descriptive distribution of the temperatures for plot SB, box plots representing each ROI in the plot are also reported (Figure 11). Generally, both plot types had a denser distribution of the values in the lower portion, as highlighted by the median, showing the non-normal distribution of temperatures (as shown by the mean represented using the symbol “×” in the boxes). Biochar boxplots showed higher maximums and generally a wider range of minimum values compared to the soil boxplots, as well as a greater interquartile range, testifying a general tendency toward extreme temperatures, which is also evidenced by the error bars (that depict a difference between soil and biochar areas of up to 10 °C).

4. Discussion

Analyzing the SST daily variations can have important implications for agronomy and environmental studies since the soil temperature drives soil physicochemical and biological processes. Traditionally, soil temperature data is gathered using soil probes at varying soil depths; however, the spatial heterogeneity of key soil properties at various scales hinders representative soil monitoring using point data. IRT can potentially overcome this limitation thanks to its scalability, and thereby improve our understanding of land and soil management via understanding the soil temperature.
IRT was used similarly by previous authors, but rather than to assess the SST, it was used in laboratory experiments to estimate the soil surface microrelief and rill morphology [97], to map soil surface permeability and identify preferential flow [98], to map macroporosity at the soil surface [99], to assess cooling the soil surface with cold water for assessing soil water repellency [100]. This paper proposed a novel, scalable approach for applying IRT monitoring of SST in environmental applications. In particular, the SST patterns of selected plots (biochar plot B), bare soil (plot S), and biochar-amended soil (plot SB) of a vineyard test site in a warm temperate environmental zone (Anadia, Portugal) were analyzed over a 48 h period. This was characterized by dry (day 1) and wet conditions (day 2) (Figure 5), which were due to a minor rainfall event that occurred during 2:00–4:00 on day 2 (Figure 3). The IRT data were acquired at high spatial and thermal resolutions (sub-millimeter and hundredth of a degree Celsius, respectively) with 30 min intervals using a tripod-mounted infrared thermal camera equipped with a focal plane array uncooled microbolometer. The purpose was to test the effect of biochar on the soil thermal properties since biochar can be applied to soils to retain more water and to improve the quality and fertility of the soil [47,48,49,50,51,52,53,54,55,56,101]. Due to its cost-effectiveness, portability, and easy battery power supply, this type of sensor was more suitable for this specific application as a cooled system.
The results of the experiments showed that with the 100% biochar cover (plot B), the SST reached a higher daily temperature peak relative to the other plots, and showed a faster rate of T increase during the day (especially in the last part of the rise toward the peak). This effect may be explained by the reduction in soil surface albedo, which according to the dose–response study of Verheijen et al. [61], was expected to be in the range of 80%. However, on both days, the SST of plot B dipped below the T of the control treatment (plot S) and the biochar-amended soil (plot SB) from about 1800 onward and throughout the night. This thermal behavior was likely connected to the biochar’s higher emissivity and a lower thermal inertia due to coarser and loose grains that allowed for a lower resistance to air circulation relative to soil [47,97,98]. This thermal behavior was further supported by analyzing the range in the SST graph for all three plots (Figure 12). Plot B had a larger range compared to plots S and SB, which was caused by the higher thermal inertia of the soil in both plots. This is also supported by the dose–response study of Verheijen et al. [61], which showed that the biochar concentration in plot SB was expected to result in a soil surface albedo reduction of 30%–40%. This effect on SST variation was also evident when biochar was mixed with soil (plot SB), even at a 10 cm depth (Figure 7 and Figure 8), as noted by the fast response of the T probe. The interpretation of the subsurface soil temperature data was limited by the probe’s resolution (i.e., 0.1 °C) and potential minor differences in the installation depth. The latter was compensated for to some extent by using the mean of two side-by-side probes for each plot. The effect size of the higher T peaks was several times the size of the probe resolution, which strengthened our confidence in the results. The temporal resolution of this case study is also relevant to experiments in soil science aimed at determining whether, and by how much, biochar application may increase the soil T, which is an important driver of soil biological processes, including relevant implications for soil respiration [102] and evaporation [103] at the effect size on soil T that we found here. At the same time, regarding the interpretation of the soil SSTs, the adopted thermal camera thermal sensitivity was 0.04 °C, while the observed effects on the measured values were many times larger. The confidence in this would also increase through upscaling in future studies. A drawback of this exploratory study is that we did not have physical replicates; therefore, new experiments in a similar climatic context will be planned to enhance confidence in the interpretation of the SST effects, also through upscaling to a larger scale of observation. Further investigations with more and higher-resolution T probes are recommended to further our understanding of how the SST regimes relate to temperature regimes within the soil profile, and in particular, how biochar impacts this relationship. The effect of rainfall on the SST of all the graphs was particularly evident in day 2 when the soil water content increased, and therefore, all of the temperature peaks decreased, the diurnal temperature fluctuation range was reduced, and the SST patterns tended to homogenize following the second peak.
The “within-plot scale” analysis in plot SB allowed for assessing the thermal behavior of single picture elements over the entire monitoring period, and consequently, getting a better understanding of the thermal differences between bare and biochar-covered spots of the soil surface.
The approach used to extract the spatial information from the IRTs showed sharp peaks in the cooling and heating cycles for biochar, confirming its different thermal behavior relative to bare soil, and suggesting that biochar grains can modify the thermal properties of the topsoil by lowering its thermal inertia relative to bare soil. Indeed, the thermal behavior of the single bare soil sectors in the Plot SB showed the same thermal behavior (SST heating/cooling phases evolution and maximum and minimum values) as plot S. The same thermal behavior was not found for the biochar grains because individual biochar grains had the same thermal behavior but a different SST range compared to the whole B plot.

5. Conclusions

The adopted methodology has shown how its high temporal and thermal resolutions can provide useful insights for the understanding of daily microclimate variations in a warm temperate environmental zone; it also allowed for observing differences related to one specific type of soil management practice. The abilities regarding data collection, as well as management and soil properties determination, teamed with agronomists’ knowledge of how the influence of such parameters on vineyards, allows for developing special microculture methods to increase the productivity of plantations. IRT has shown the potential to upscale soil T data to the crop scale and as input for precision farming operations. For biochar field experiments, it may be especially useful to improve the spatial resolution of soil T data and support furthering our understanding of biochar’s effects on soil T dynamics. Although the proposed methodology was applied to soil plots, in future research, the raster approach allows for an upscaling of the procedure through the application of the sensors on a drone platform. This, combined with the high resolution of the infrared sensor. will allow for extending the analysis of larger areas, such as the whole vineyard/plantation. To achieve such a goal, a proper number of thermohygrometer sensors should be placed in the investigated area to collect the necessary parameters for the thermogram image corrections. This procedure will also require the placing of georeferenced optical markers that are necessary for the composition of the acquired images mosaic. The SST analysis will require the creation of mask features to analyze only the soil areas by excluding the vegetated portion; this could be achieved through the application of a supervised classification of the area in a GIS environment. Through the acquisition of multiple images with a designated time step and thanks to the batch analysis, it can be possible to obtain a time-dependent variation of the SST (theoretically, by computing a raster difference, we could obtain a map showing the highest soil T for the analyzed time). Future developments of the procedure can lead to a fully automated protocol that can monitor the interested areas on a daily basis.

Author Contributions

W.F. and F.G.A.V. planned and performed both the experiment and the methodological approach; W.F. carried out the IR surveys, with the help of F.G.A.V., and performed the software elaboration with G.L. Writing—review and editing was done by W.F., F.G.A.V., G.L., and J.K. Supervision, funding acquisition, and overall review was done by S.M. and J.K. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the post-graduate internship of William Frodella at the University of Aveiro, funded by the University of Florence, and the Fundação para a Ciencia e Tecnologia FCT-funded research contract of Frank G. A. Verheijen (CEECIND/02509/2018) and Jan J. Keizer (IF/01465/2015). Thanks are due to FCT/MCTES (Ministério da Ciência, Tecnologia e Ensino Superior) for the funding of the MASSC project (“Mediterranean Agricultural Soils Conservation under Global Change”; ref. ARIMNET2/0006/2015) in the framework of ARIMNET2, as well as for the financial support to Centre for Environmental And Marine Studies (CESAM) (UIDP/50017/2020+UIDB/50017/2020), in both cases through national funds.

Acknowledgments

The Estação Vitivinícola da Bairrada—Regional Ministry of Agriculture (DRAPC) is thanked for making available the field site and logistical support. We wish to thank the Earth Surface Processes (ESP) team at the University of Aveiro for useful feedback and support. The authors would like to Thank Federico Di Traglia for his accurate revisions

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van De Griend, A.A.; Owe, M. On the relationship between thermal emissivity and the normalized difference vegetation index for natural surfaces. Int. J. Remote Sens. 1993, 14, 1119–1131. [Google Scholar] [CrossRef]
  2. Bastiaanssen, W.; Pelgrum, H.; Wang, J.; Ma, Y.; Moreno, J.; Roerink, G.; Van Der Wal, T. A remote sensing surface energy balance algorithm for land (SEBAL). J. Hydrol. 1998, 212, 213–229. [Google Scholar] [CrossRef]
  3. Heusinkveld, B.G.; Jacobs, A.; Holtslag, A.; Berkowicz, S. Surface energy balance closure in an arid region: Role of soil heat flux. Agric. For. Meteorol. 2004, 122, 21–37. [Google Scholar] [CrossRef]
  4. Niu, G.-Y.; Yang, Z.; Mitchell, K.E.; Chen, F.; Ek, M.B.; Barlage, M.; Kumar, A.; Manning, K.; Niyogi, D.; Rosero, E.; et al. The community Noah land surface model with multiparameterization options (Noah-MP): 1. Model description and evaluation with local-scale measurements. J. Geophys. Res. Space Phys. 2011, 116, 116. [Google Scholar] [CrossRef] [Green Version]
  5. Sellers, P.J.; Dickinson, R.E.; Randall, D.A.; Betts, A.K.; Hall, F.G.; Berry, J.A.; Collatz, G.J.; Denning, A.S.; Mooney, H.A.; Nobre, C.A.; et al. Henderson-Sellers A. Modeling the exchanges of energy, water, and carbon between continents and the atmosphere. Science 1997, 275, 502–509. [Google Scholar] [CrossRef] [Green Version]
  6. Mira, M.; Valor, E.; Boluda, R.; Caselles, V.; Coll, C. Influence of soil water content on the thermal infrared emissivity of bare soils: Implication for land surface temperature determination. J. Geophys. Res. Space Phys. 2007, 112, 112. [Google Scholar] [CrossRef] [Green Version]
  7. Steenpass, C.; VanderBorght, J.; Herbst, M.; Šimůnek, J.; Vereecken, H. Estimating Soil Hydraulic Properties from Infrared Measurements of Soil Surface Temperatures and TDR Data. Vadose Zone J. 2010, 9, 910–924. [Google Scholar] [CrossRef] [Green Version]
  8. Rossel, R.V.; Webster, R. Predicting soil properties from the Australian soil visible–near infrared spectroscopic database. Eur. J. Soil Sci. 2012, 63, 848–860. [Google Scholar] [CrossRef]
  9. McMillin, L.M. Estimation of sea surface temperatures from two infrared window measurements with different absorption. J. Geophys. Res. 1975, 80, 5113–5117. [Google Scholar] [CrossRef]
  10. Hulley, G.C.; Hook, S.J.; Baldridge, A.M. Investigating the effects of soil moisture on thermal infrared land surface temperature and emissivity using satellite retrievals and laboratory measurements. Remote Sens. Environ. 2010, 114, 1480–1493. [Google Scholar] [CrossRef]
  11. Tomlinson, C.J.; Chapman, L.; Thornes, J.E.; Baker, C. Remote sensing land surface temperature for meteorology and climatology: A review. Meteorol. Appl. 2011, 18, 296–306. [Google Scholar] [CrossRef] [Green Version]
  12. Li, Z.L.; Tang, B.H.; Wu, H.; Ren, H.; Yan, G.; Wan, Z.; Trigo, I.F.; Sobrino, J.A. Satellite-derived land surface temperature: Current status and perspectives. Remote Sens. Environ. 2013, 131, 14–37. [Google Scholar] [CrossRef] [Green Version]
  13. Abrams, M.; Tsu, H.; Hulley, G.; Iwao, K.; Pieri, D.; Cudahy, T.; Kargel, J. The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) after fifteen years: Review of global products. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 292–301. [Google Scholar] [CrossRef]
  14. Anderson, J.M.; Wilson, S.B. The physical basis of current infrared remote-sensing techniques and the interpretation of data from aerial surveys. Int. J. Remote Sens. 1984, 5, 1–18. [Google Scholar]
  15. Kirkland, L.; Herr, K.; Keim, E.; Adams, P.; Salisbury, J.; Hackwell, J.; Treiman, A. First use of an airborne thermal infrared hyperspectral scanner for compositional mapping. Remote Sens. Environ. 2002, 80, 447–459. [Google Scholar] [CrossRef] [Green Version]
  16. Sobrino, J.A.; Jiménez-Muñoz, J.-C.; Zarco-Tejada, P.J.; Sepulcre-Canto, G.; De Miguel, E. Land surface temperature derived from airborne hyperspectral scanner thermal infrared data. Remote. Sens. Environ. 2006, 102, 99–115. [Google Scholar] [CrossRef]
  17. Eccel, E.; Arman, G.; Zottele, F.; Gioli, B. Thermal infrared remote sensing for high resolution minimum temperature mapping. Ital. J. Agrometeorol. 2008, 13, 52–61. [Google Scholar]
  18. Pasquier, C.; Bourennane, H.; Cousin, I.; Séger, M.; Dabas, M.; Thiesson, J. Comparison between thermal airborne remote sensing, multi-depth electrical resistivity profiling, and soil mapping: An example from Beauce (Loiret, France). Near Surf. Geophys. 2016, 14, 345–356. [Google Scholar] [CrossRef] [Green Version]
  19. Kant, Y.; Badarinath, K.V.S. Ground-based method for measuring thermal infrared effective emissivities: Implications and perspectives on the measurement of land surface temperature from satellite data. Int. J. Remote Sens. 2002, 23, 2179–2191. [Google Scholar] [CrossRef]
  20. Qin, Z.; Berliner, P.; Karnieli, A. Ground temperature measurement and emissivity determination to understand the thermal anomaly and its significance on the development of an arid environmental ecosystem in the sand dunes across the Israel–Egypt border. J. Arid. Environ. 2005, 60, 27–52. [Google Scholar]
  21. Pfister, L.; McDonnell, J.J.; Hissler, C.; Hoffmann, L. Ground-based thermal imagery as a simple, practical tool for mapping saturated area connectivity and dynamics. Hydrol. Process. 2010, 24, 3123–3132. [Google Scholar] [CrossRef]
  22. Price, J.C. The potential of remotely sensed thermal infrared data to infer surface soil moisture and evaporation. Water Resour. Res. 1980, 16, 787–795. [Google Scholar] [CrossRef]
  23. Casagli, N.; Frodella, W.; Morelli, S.; Tofani, V.; Ciampalini, A.; Intrieri, E.; Raspini, F.; Rossi, G.; Tanteri, L.; Lu, P. Spaceborne, UAV and ground-based remote sensing techniques for landslide mapping, monitoring and early warning. Geoenviron. Disasters 2017, 4, 9. [Google Scholar] [CrossRef]
  24. Di Traglia, F.; Nolesini, T.; Ciampalini, A.; Solari, L.; Frodella, W.; Bellotti, F.; Fumagalli, A.; De Rosa, G.; Casagli, N. Tracking morphological changes and slope instability using space-borne and ground-based SAR data. Geomorphology 2018, 300, 95–112. [Google Scholar] [CrossRef]
  25. Moran, M.S.; Peters-Lidard, C.D.; Watts, J.M.; McElroy, S. Estimating soil moisture at the watershed scale with satellite-based radar and land surface models. Can. J. Remote Sens. 2004, 30, 805–826. [Google Scholar] [CrossRef] [Green Version]
  26. Coolbaugh, M.F.; Kratt, C.; Fallacaro, A.; Calvin, W.M.; Taranik, J.V. Detection of geothermal anomalies using advanced spaceborne thermal emission and reflection radiometer (ASTER) thermal infrared images at Bradys Hot Springs, Nevada, USA. Remote Sens. Environ. 2007, 106, 350–359. [Google Scholar] [CrossRef]
  27. Gigli, G.; Frodella, W.; Garfagnoli, F.; Morelli, S.; Mugnai, F.; Menna, F.; Casagli, N. 3-D geomechanical rock mass characterization for the evaluation of rockslide susceptibility scenarios. Landslides 2014, 11, 131–140. [Google Scholar] [CrossRef]
  28. Kerr, Y.H.; Lagouarde, J.P.; Imbernon, J. Accurate land surface temperature retrieval from AVHRR data with use of an improved split window algorithm. Remote Sens. Environ. 1992, 41, 197–209. [Google Scholar] [CrossRef]
  29. Prata, A.J.; Caselles, V.; Coll, C.; Sobrino, J.A.; Ottle, C. Thermal remote sensing of land surface temperature from satellites: Current status and future prospects. Remote Sens. Rev. 1995, 12, 175–224. [Google Scholar] [CrossRef]
  30. Wan, Z.; Zhang, Y.; Zhang, Q.; Li, Z.L. Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data. Remote Sens. Environ. 2002, 83, 163–180. [Google Scholar] [CrossRef]
  31. Zhou, D.K.; Larar, A.M.; Liu, X. On the relationship between land surface infrared emissivity and soil moisture. J. Appl. Remote Sens. 2018, 12, 016030. [Google Scholar] [CrossRef] [Green Version]
  32. Shaver, G.R.; Canadell, J.; Chapin, F.S.; Gurevitch, J.; Harte, J.; Henry, G.; Ineson, P.; Jonasson, S.; Melillo, J.; Pitelka, L.; et al. Global Warming and Terrestrial Ecosystems: A Conceptual Framework for Analysis: Ecosystem responses to global warming will be complex and varied. Ecosystem warming experiments hold great potential for providing insights on ways terrestrial ecosystems will respond to upcoming decades of climate change. Documentation of initial conditions provides the context for understanding and predicting ecosystem responses. Aibs Bull. 2000, 50, 871–882. [Google Scholar]
  33. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Seidl, R.; Delzon, S.; Corona, P.; Kolström, M.; et al. Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
  34. Meola, C.; Boccardi, S.; Carlomagno, G.M. Infrared Thermography in the Evaluation of Aerospace Composite Materials: Infrared Thermography to Composites; Woodhead Publishing: Cambridge, UK, 2016; pp. 57–83. [Google Scholar] [CrossRef]
  35. Jones, H.G.; Leinonen, I. Thermal imaging for the study of plant water relations. J. Agric. Meteorol. 2003, 59, 205–217. [Google Scholar] [CrossRef] [Green Version]
  36. Grant, O.M.; Tronina, Ł.; Jones, H.G.; Chaves, M.M. Exploring thermal imaging variables for the detection of stress responses in grapevine under different irrigation regimes. J. Exp. Bot. 2006, 58, 815–825. [Google Scholar] [CrossRef]
  37. Spampinato, L.; Calvari, S.; Oppenheimer, C.; Boschi, E. Volcano surveillance using infrared cameras. Earth Sci. Rev. 2011, 106, 63–91. [Google Scholar] [CrossRef]
  38. Reinert, S.; Bögelein, R.; Thomas, F.M. Use of thermal imaging to determine leaf conductance along a canopy gradient in European beech (Fagus sylvatica). Tree Physiol. 2012, 32, 294–302. [Google Scholar] [CrossRef] [Green Version]
  39. Clark, M.R.; McCann, D.M.; Forde, M.C. Application of infrared thermography to the non-destructive testing of concrete and masonry bridges. NDT Int. 2003, 36, 265–275. [Google Scholar] [CrossRef]
  40. De Lima, J.L.M.P.; Abrantes, J.R.C.B. Can infrared thermography be used to estimate soil surface microrelief and rill morphology? Catena 2014, 113, 314–322. [Google Scholar] [CrossRef]
  41. Frodella, W.; Gigli, G.; Morelli, S.; Lombardi, L.; Casagli, N. Landslide Mapping and Characterization through Infrared Thermography (IRT): Suggestions for a Methodological Approach from Some Case Studies. Remote Sens. 2017, 9, 1281. [Google Scholar] [CrossRef] [Green Version]
  42. Nolesini, T.; Frodella, W.; Bianchini, S.; Casagli, N. Detecting Slope and Urban Potential Unstable Areas by Means of Multi-Platform Remote Sensing Techniques: The Volterra (Italy) Case Study. Remote Sens. 2016, 8, 746. [Google Scholar] [CrossRef] [Green Version]
  43. Di Traglia, F.; Nolesini, T.; Solari, L.; Ciampalini, A.; Frodella, W.; Steri, D.; Allotta, B.; Rindi, A.; Marini, L.; Monni, N.; et al. Lava delta deformation as a proxy for submarine slope instability. Earth Planet. Sci. Lett. 2018, 488, 46–58. [Google Scholar] [CrossRef]
  44. Frodella, W.; Morelli, S.; Pazzi, V. Infrared Thermographic surveys for landslide mapping and characterization: The Rotolon DSGSD (Norther Italy) case study. Ital. J. Eng. Geol. Environ. 2017, 2017. [Google Scholar] [CrossRef]
  45. Maldague, X. Theory and Practice of Infrared Technology for Non-Destructive Testing; John-Wiley & Sons: Hoboken, NJ, USA, 2001; 684p. [Google Scholar]
  46. Baldock, J.A.; Smernik, R.J. Chemical composition and bioavailability of thermally altered Pinus resinosa (Red pine) wood. Org. Geochem. 2002, 33, 1093–1109. [Google Scholar] [CrossRef]
  47. Verheijen, F.G.; Montanarella, L.; Bastos, A.C. Sustainability, certification, and regulation of biochar. Pesqui. Agropecuária Bras. 2012, 47, 649–653. [Google Scholar] [CrossRef] [Green Version]
  48. Jeffery, S.; Verheijen, F.G.; van der Velde, M.; Bastos, A.C. A quantitative review of the effects of biochar application to soils on crop productivity using meta-analysis. Agric. Ecosyst. Environ. 2011, 144, 175–187. [Google Scholar] [CrossRef]
  49. Jeffery, S.; Verheijen, F.G.; Bastos, A.C.; Van Der Velde, M. A comment on ‘Biochar and its effects on plant productivity and nutrient cycling: A meta-analysis’: On the importance of accurate reporting in supporting a fast-moving research field with policy implications. Gcb Bioenergy 2014, 6, 176–179. [Google Scholar] [CrossRef] [Green Version]
  50. Jeffery, S.; Abalos, D.; Prodana, M.; Bastos, A.C.; Van Groenigen, J.W.; Hungate, B.A.; Verheijen, F. Biochar boosts tropical but not temperate crop yields. Environ. Res. Lett. 2017, 12, 053001. [Google Scholar] [CrossRef]
  51. Jeffery, S.; Verheijen, F.G.; Kammann, C.; Abalos, D. Biochar effects on methane emissions from soils: A meta-analysis. Soil Biol. Biochem. 2016, 101, 251–258. [Google Scholar] [CrossRef]
  52. Cayuela, M.L.; Van Zwieten, L.; Singh, B.P.; Jeffery, S.; Roig, A.; Sánchez-Monedero, M.A. Biochar’s role in mitigating soil nitrous oxide emissions: A review and meta-analysis. Agric. Ecosyst. Environ. 2014, 191, 5–16. [Google Scholar] [CrossRef]
  53. Wang, J.; Xiong, Z.; Kuzyakov, Y. Biochar stability in soil: Meta-analysis of decomposition and priming effects. GCB Bioenergy 2016, 8, 512–523. [Google Scholar] [CrossRef] [Green Version]
  54. Omondi, M.O.; Xia, X.; Nahayo, A.; Liu, X.; Korai, P.K.; Pan, G. Quantification of biochar effects on soil hydrological properties using meta-analysis of literature data. Geoderma 2016, 274, 28–34. [Google Scholar] [CrossRef]
  55. Verheijen, F.G.A.; Zhuravel, A.; Silva, F.; Amaro, A.; Ben-Hur, M.; Keizer, J.J. The influence of biochar particle size and concentration on bulk density and maximum water holding capacity of sandy vs. sandy loam soil in a column experiment. Geoderma 2019, 347, 194–202. [Google Scholar] [CrossRef]
  56. Abrol, V.; Ben-Hur, M.; Verheijen, F.G.; Keizer, J.J.; Martins, M.A.; Tenaw, H.; Tchehansky, L.; Graber, E.R. Biochar effects on soil water infiltration and erosion under seal formation conditions: Rainfall simulation experiment. J. Soils Sediments 2016, 16, 2709–2719. [Google Scholar] [CrossRef]
  57. Hilber, I.; Bastos, A.C.; Loureiro, S.; Soja, G.; Marsz, A.; Cornelissen, G.; Bucheli, T.D. The different faces of biochar: Contamination risk versus remediation tool. J. Environ. Eng. Landsc. Manag. 2017, 25, 86–104. [Google Scholar] [CrossRef] [Green Version]
  58. Davidson, E.A.; Janssens, I.A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 2006, 440, 165–173. [Google Scholar] [CrossRef]
  59. Intergovernmental Panel on Climate Change. Climate Change 2014—Impacts, Adaptation and Vulnerability: Regional Aspects; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  60. Routson, C.C.; McKay, N.P.; Kaufman, D.S.; Erb, M.P.; Goosse, H.; Shuman, B.N.; Rodysill, J.R.; Ault, T. Mid-latitude net precipitation decreased with Arctic warming during the Holocene. Nature 2019, 568, 83–87. [Google Scholar] [CrossRef]
  61. Verheijen, F.G.; Jeffery, S.; Van Der Velde, M.; Penížek, V.; Béland, M.; Bastos, A.C.; Keizer, J.J. Reductions in soil surface albedo as a function of biochar application rate: Implications for global radiative forcing. Environ. Res. Lett. 2013, 8, 044008. [Google Scholar] [CrossRef] [Green Version]
  62. Nelissen, V.; Ruysschaert, G.; Manka’Abusi, D.; D’Hose, T.; De Beuf, K.; Al-Barri, B.; Boeckx, P. Impact of a woody biochar on properties of a sandy loam soil and spring barley during a two-year field experiment. Eur. J. Agron. 2015, 62, 65–78. [Google Scholar] [CrossRef]
  63. Ventura, M.; Sorrenti, G.; Panzacchi, P.; George, E.; Tonon, G. Biochar Reduces Short-Term Nitrate Leaching from A Horizon in an Apple Orchard. J. Environ. Qual. 2013, 42, 76–82. [Google Scholar] [CrossRef]
  64. Usowicz, B.; Lipiec, J.; Łukowski, M.; Marczewski, W.; Usowicz, J. The effect of biochar application on thermal properties and albedo of loess soil under grassland and fallow. Soil Tillage Res. 2016, 164, 45–51. [Google Scholar] [CrossRef]
  65. Tammeorg, P.; Bastos, A.C.; Jeffery, S.; Rees, F.; Kern, J.; Graber, E.R.; Cordovil, C.M.D.S. Biochars in soils: Towards the required level of scientific understanding. J. Environ. Eng. Landsc. Manag. 2017, 25, 192–207. [Google Scholar] [CrossRef] [Green Version]
  66. Lillesand, T.M.; Kiefer, R.W.; Chipman, J.W. Remote Sensing and Image Interpretation, 6th ed.; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2008. [Google Scholar]
  67. Dhar, N.K.; Sood, A.K.; Dat, R. Advances in Infrared Detector Array Technology. In Optoelectronics—Advanced Materials and Devices; INTECH Open Access Publisher: London, UK, 2013; Volume 1. [Google Scholar]
  68. FLIR Systems Inc. ThermaCAM SC620 Technical Specifications. 2009. Available online: http://www.flir.com/cs/emea/en/view/?id=41965 (accessed on 21 January 2020).
  69. FLIR Systems Inc. FLIR Tools+ Datasheet. 2015. Available online: https://www.infraredcamerawarehouse.com/content/FLIR%20Datasheets/FLIR%20ToolsPlus%20Datasheet.pdf (accessed on 22 January 2020).
  70. Incropera, P.F.; De Witt, D.P. Introduction to Heat Transfer; John Wiley: Hobekone, NJ, USA, 1985; pp. 654–746. [Google Scholar]
  71. Rees, W.G. Physical Principles of Remote Sensing, 2nd ed.; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  72. Gupta, R.P. Interpretation of data in the thermal-infrared region. In Remote Sensing Geology; Springer: New York, NY, USA, 2003; pp. 190–203. [Google Scholar]
  73. Verstraeten, W.W.; Veroustraete, F.; van der Sande, C.J.; Grootaers, I.; Feyen, J. Soil moisture retrieval using thermal inertia, determined with visible and thermal spaceborne data, validated for European forests. Remote Sens. Environ. 2006, 101, 299–314. [Google Scholar] [CrossRef]
  74. Hardgrove, C.; Moersch, J.; Whisner, S. Thermal imaging of sedimentary features on alluvial fans. Planet. Space Sci. 2010, 58, 482–508. [Google Scholar] [CrossRef]
  75. Van de Griend, A.A.; Camillo, P.J.; Gurney, R.J. Discrimination of soil physical parameters, thermal inertia, and soil moisture from diurnal surface temperature fluctuations. Water Resour. Res. 1985, 21, 997–1009. [Google Scholar] [CrossRef]
  76. Jakosky, B.M. On the thermal properties of Martian fines. Icarus 1986, 66, 117–124. [Google Scholar] [CrossRef]
  77. Mutlutürk, M.; Altindag, R.; Türk, G. A decay function model for the integrity loss of rock when subjected to recurrent cycles of freezing–thawing and heating–cooling. Int. J. Rock Mech. Min. Sci. 2004, 41, 237–244. [Google Scholar] [CrossRef]
  78. Ball, M.; Pinkerton, H. Factors affecting the accuracy of thermal imaging cameras in volcanology. J. Geophys. Res. 2006, 111, B11203. [Google Scholar] [CrossRef]
  79. Salisbury, J.W.; D’Aria, D.M. Emissivity of terrestrial materials in the 8–14 µm atmospheric window. Remote Sens. Environ. 1992, 42, 83–106. [Google Scholar] [CrossRef]
  80. Rubio, E.; Caselles, V.; Badenas, C. Emissivity measurements of several soils and vegetation types in the 8–14 μm wave band: Analysis of two field methods. Remote Sens. Environ. 1997, 59, 490–521. [Google Scholar] [CrossRef]
  81. Garrity, D.; Bindraban, P. A Globally Distributed Soil Spectral Library Visible near Infrared Diffuse Reflectance Spectra; ICRAF (World Agroforestry Centre)/ISRIC (World Soil Information) Spectral Library: Nairobi, Kenya, 2004. [Google Scholar]
  82. ISO 18434-1:2008 Condition Monitoring and Diagnostics of Machines—Thermography—Part 1: General Procedures. Available online: https://www.iso.org/standard/41648.html (accessed on 18 January 2019).
  83. Avdelidis, N.P.; Moropoulou, A. Emissivity considerations in building thermography. Energy Build. 2003, 35, 663–667. [Google Scholar] [CrossRef]
  84. Chicco, J.M.; Vacha, D.; Mandrone, G. Thermo-Physical and Geo-Mechanical Characterization of Faulted Carbonate Rock Masses (Valdieri, Italy). Remote Sens. 2019, 11, 179. [Google Scholar] [CrossRef] [Green Version]
  85. Metzger, M.J.; Bunce, R.G.H.; Jongman, R.H.G.; Sayre, R.; Trabucco, A.; Zomer, R. A high-resolution bioclimate map of the world: A unifying framework for global biodiversity research and monitoring. Glob. Ecol. Biogeogr. 2013, 22, 630–638. [Google Scholar] [CrossRef] [Green Version]
  86. Pereira, V.; FitzPatrick, E.A. Cambisols and related soils in north-central Portugal: Their genesis and classification. Geoderma 1995, 66, 185–212. [Google Scholar] [CrossRef]
  87. Inácio, M.; Pereira, V.; Pinto, M. The soil geochemical atlas of Portugal: Overview and applications. J. Geochem. Explor. 2008, 98, 22–33. [Google Scholar] [CrossRef]
  88. King, P.M. Comparison of methods for measuring severity of water repellence of sandy soils and assesment of some factors that affect its measurement (Australia). Aust. J. Soil Res. 1981, 19, 275–285. [Google Scholar] [CrossRef]
  89. Keizer, J.J.; Doerr, S.H.; Malvar, M.C.; Ferreira, A.J.D.; Pereira, V.M.F.G. Temporal and spatial variations in topsoil water repellency throughout a crop-rotation cycle on sandy soil in north-central Portugal. Hydrol. Process. 2007, 21, 2317–2324. [Google Scholar] [CrossRef]
  90. Prodana, M.; Bastos, A.C.; Amaro, A.; Cardoso, D.; Morgado, R.; Machado, A.L.; Loureiro, S. Biomonitoring tools for biochar and biochar-compost amended soil under viticulture: Looking at exposure and effects. Appl. Soil Ecol. 2019, 137, 120–128. [Google Scholar] [CrossRef]
  91. IPMA (Instituto Português do Mar e da Atmosfera). Available online: https://www.ipma.pt/pt/index.html (accessed on 21 January 2020).
  92. FLIR Systems Inc. ThermaCAM SC620 User’s Manual; FLIR Systems Inc.: Wilsonville, OR, USA, 2009. [Google Scholar]
  93. Rossel, R.A.V.; Behrens, T.; Ben-Dor, E.; Brown, D.; Demattê, J.; Shepherd, K.; Shi, Z.; Stenberg, B.; Stevens, A.; Adamchuk, V.; et al. A global spectral library to characterize the world’s soil. Earth-Science Rev. 2016, 155, 198–230. [Google Scholar] [CrossRef] [Green Version]
  94. Salazar, R.K.; Laurente, C.L.B. Emissivity of “Grey” Bodies, Surface Roughness and other Measures. Recoletos Multidiscip. Res. J. 2014, 2, 2. [Google Scholar]
  95. Lou, C.W.; Lin, J.H. Evaluation of bamboo charcoal/stainless steel/TPU composite woven fabrics. Fibers Polym. 2011, 12, 514. [Google Scholar] [CrossRef]
  96. ESRI ArcGis Raster Calculator ArcMap 10.3.2. Available online: http://desktop.arcgis.com/en/arcmap/10.3/tools/spatial-analyst-toolbox/raster-calculator.htm (accessed on 23 February 2020).
  97. De Lima, J.L.M.P.; Abrantes, J.R.C.B.; Silva, V.P., Jr.; Montenegro, A.A.A. Prediction of skin surface soil permeability by infrared thermography: A soil flume experiment. Quant. Infrared Thermogr. J. 2014, 11, 161–169. [Google Scholar] [CrossRef]
  98. De Lima, J.L.M.P.; Abrantes, J.R.C.B.; Silva, V.P., Jr.; De Lima, M.I.P.; Montenegro, A.A.A. Mapping soil surface macropores using infrared thermography: An exploratory laboratory study. Sci. World J. 2014, 2014, 84560. [Google Scholar] [CrossRef] [Green Version]
  99. Abrantes, J.R.C.B.; De Lima, J.L.M.P.; Prats, S.A.; Keizer, J.J. Assessing soil water repellency spatial variability using a thermographic technique: Small-scale laboratory study. Geoderma 2017, 287, 98–104. [Google Scholar] [CrossRef]
  100. Abrantes, J.R.C.B.; De Lima, J.L.M.P.; Prats, S.A.; Keizer, J.J. Field assessment of soil water repellency using infrared thermography. Forum Geogr. 2016, 15, 12–18. [Google Scholar] [CrossRef]
  101. Kirschbaum, M. The temperature dependence of soil organic matter decomposition, and the effect of global warming on soil organic C storage. Soil Biol. Biochem. 1995, 27, 753–760. [Google Scholar] [CrossRef]
  102. Usowicz, B.; Lipiec, J.; Lukowski, M.; Marczewski, W.; Usowicz, J. The effect of biochar application on thermal inertia of soil. EGU General Assembly Conference Abstracts. Available online: https://ui.adsabs.harvard.edu/abs/2015EGUGA..17.9373U/abstract (accessed on 25 April 2020).
  103. Qiu, G.Y.; Yano, T.; Momii, K. An improved methodology to measure evaporation from bare soil based on comparison of surface temperature with a dry soil surface. J. Hydrol. 1998, 210, 93–105. [Google Scholar] [CrossRef]
Sensors 20 02444 g001 550
Figure 1. Examples of the analyzed plots: plot B (a), plot S (b), and plot SB (c). Corresponding infrared thermography (IRT) images underneath (d–f). The black square in panel f shows the surveyed IRT sector (region of interest (ROI)).
Figure 1. Examples of the analyzed plots: plot B (a), plot S (b), and plot SB (c). Corresponding infrared thermography (IRT) images underneath (d–f). The black square in panel f shows the surveyed IRT sector (region of interest (ROI)).
Sensors 20 02444 g001
Sensors 20 02444 g002 550
Figure 2. Experimental setup in the vineyard study area: (a) tripod-mounted thermal camera, (b) Decagon Em50 data logger, and (c) vineyard control treatment soil plot.
Figure 2. Experimental setup in the vineyard study area: (a) tripod-mounted thermal camera, (b) Decagon Em50 data logger, and (c) vineyard control treatment soil plot.
Sensors 20 02444 g002
Sensors 20 02444 g003 550
Figure 3. Pocket thermohygrometer weather data and hourly precipitation from the IPMA station.
Figure 3. Pocket thermohygrometer weather data and hourly precipitation from the IPMA station.
Sensors 20 02444 g003
Sensors 20 02444 g004 550
Figure 4. Methodological workflow (a) and analyzed ROIs on the plot SB thermogram (b): bare soil sectors (red polygons) and biochar clusters (green polygons).
Figure 4. Methodological workflow (a) and analyzed ROIs on the plot SB thermogram (b): bare soil sectors (red polygons) and biochar clusters (green polygons).
Sensors 20 02444 g004
Sensors 20 02444 g005 550
Figure 5. IRT dataset of the analyzed plots showing the SST variations on the single thermograms.
Figure 5. IRT dataset of the analyzed plots showing the SST variations on the single thermograms.
Sensors 20 02444 g005
Sensors 20 02444 g006 550
Figure 6. Maximum (dashed lines), minimum (dotted lines), and mean (continuous lines) soil surface temperatures monitored during the survey. Plot B (100% biochar cover)—black lines, plot S (bare soil)—red lines, plot SB (mixed soil-biochar)—dark blue lines, and air—light blue line.
Figure 6. Maximum (dashed lines), minimum (dotted lines), and mean (continuous lines) soil surface temperatures monitored during the survey. Plot B (100% biochar cover)—black lines, plot S (bare soil)—red lines, plot SB (mixed soil-biochar)—dark blue lines, and air—light blue line.
Sensors 20 02444 g006
Sensors 20 02444 g007 550
Figure 7. Comparison between the air temperature, kinetic temperature at a 10 cm depth (S and SB probes), and mean kinetic soil surface temperature of biochar clusters and soil sectors within plot SB, obtained by converting the Trad measured by the IRT camera into Tkin.
Figure 7. Comparison between the air temperature, kinetic temperature at a 10 cm depth (S and SB probes), and mean kinetic soil surface temperature of biochar clusters and soil sectors within plot SB, obtained by converting the Trad measured by the IRT camera into Tkin.
Sensors 20 02444 g007
Sensors 20 02444 g008 550
Figure 8. Biochar effect on soil temperature, obtained by subtracting the bare soil (S) temperature from the biochar-amended soil (SB) temperature (TB − TS, read on primary y-axis), for both soil surface (orange) and a 10 cm depth (red), with air temperature (air T) in blue (read on secondary y-axis).
Figure 8. Biochar effect on soil temperature, obtained by subtracting the bare soil (S) temperature from the biochar-amended soil (SB) temperature (TB − TS, read on primary y-axis), for both soil surface (orange) and a 10 cm depth (red), with air temperature (air T) in blue (read on secondary y-axis).
Sensors 20 02444 g008
Sensors 20 02444 g009 550
Figure 9. SST comparison between the ROIs in plot SB: green lines represent biochar grains (B), while red lines represent bare soil sectors (S).
Figure 9. SST comparison between the ROIs in plot SB: green lines represent biochar grains (B), while red lines represent bare soil sectors (S).
Sensors 20 02444 g009
Sensors 20 02444 g010 550
Figure 10. Mean SST variations between the bare soil in plot S (red), biochar cluster/grains in plot SB (green), and bare soil sectors in plot SB (dashed red).
Figure 10. Mean SST variations between the bare soil in plot S (red), biochar cluster/grains in plot SB (green), and bare soil sectors in plot SB (dashed red).
Sensors 20 02444 g010
Sensors 20 02444 g011 550
Figure 11. Boxplots for the biochar (in green tones) and bare soil areas (in red tones) in plot SB showing the mean temperature distribution.
Figure 11. Boxplots for the biochar (in green tones) and bare soil areas (in red tones) in plot SB showing the mean temperature distribution.
Sensors 20 02444 g011
Sensors 20 02444 g012 550
Figure 12. Temperature range within the IRT images for plots B (black), S (red), and SB (blue).
Figure 12. Temperature range within the IRT images for plots B (black), S (red), and SB (blue).
Sensors 20 02444 g012
Table
Table 1. FLIR SC620 thermal camera main technical specifications.
Table 1. FLIR SC620 thermal camera main technical specifications.
FeatureUnitValue
Detector sizepixel640 × 480
Spectral rangeµm(7.5, 13)
Temperature range°C(−40, +500)
Thermal accuracy°C±2
Thermal sensitivitymK40
Field of view (FOV)°24 × 18
Lens°24
Spatial resolutionmrad0.65
Minimum focus distancem0.3
Image frequencyHz30
Table
Table 2. Weather data from the Instituto Português do Mar e da Atmosfera (IPMA) Anadia station. T: Temperature, RH: Relative humidity.
Table 2. Weather data from the Instituto Português do Mar e da Atmosfera (IPMA) Anadia station. T: Temperature, RH: Relative humidity.
Heading1 August2 August3 August
T med (°C)22.021.419.1
T max (°C)31.928.526.6
T min (°C)14.914.612.1
RH med (%)757575
RH max (%)9810099
RH min (%)363945
Daily Rainfall (mm)0.00.90.0

Share and Cite

MDPI and ACS Style

Frodella, W.; Lazzeri, G.; Moretti, S.; Keizer, J.; Verheijen, F.G.A. Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal). Sensors 2020, 20, 2444. https://doi.org/10.3390/s20092444

AMA Style

Frodella W, Lazzeri G, Moretti S, Keizer J, Verheijen FGA. Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal). Sensors. 2020; 20(9):2444. https://doi.org/10.3390/s20092444

Chicago/Turabian Style

Frodella, William, Giacomo Lazzeri, Sandro Moretti, Jacob Keizer, and Frank G. A. Verheijen. 2020. "Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal)" Sensors 20, no. 9: 2444. https://doi.org/10.3390/s20092444

APA Style

Frodella, W., Lazzeri, G., Moretti, S., Keizer, J., & Verheijen, F. G. A. (2020). Applying Infrared Thermography to Soil Surface Temperature Monitoring: Case Study of a High-Resolution 48 h Survey in a Vineyard (Anadia, Portugal). Sensors, 20(9), 2444. https://doi.org/10.3390/s20092444

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop