Geometry of the hemispherical radiometric footprint over plant canopies

Radiometric measurements of hemispherical surface reflectance and long-wave irradiance are required to quantify the broadband albedo and the outgoing thermal radiation. These observations are typically integrated with eddy covariance measurements of sensible and latent heat fluxes to characterize the surface energy budget. While the aerodynamic footprint has been widely investigated, the geometry of the hemispherical radiometric footprint over plant canopies has been rarely tackled. In the present work, the size and shape of the hemispherical radiometric footprint are formalized for a bare surface and in presence of a vegetation cover. For this purpose, four idealized canopies are analyzed and the dependency of the radiometric footprint on leaf area index and canopy height is explored. Besides, the radiometric footprint is compared with the aerodynamic footprint in conditions of neutral stability. It was observed that almost 100% of the hemispherical radiometric signal originates within a distance of a few radiometer heights, while only about 50–80% of the cumulative aerodynamic signal is generated within a distance of about 20 sensor heights. In order to achieve comparable extensions of the footprint areas, hemispherical radiometric measurements should therefore be taken about 6–15 times higher than turbulent flux ones, depending on the vegetation type. The analysis also highlights that the size of the radiative footprint decreases at increasing leaf area index, whereas the aerodynamic footprint shows an opposite behavior. For the abovementioned reasons, this work may support the interpretation of energy flux measurements and the optimal design of eddy covariance stations located in heterogeneous sites.


Introduction
Radiometric measurements of broadband upwelling shortwave and long-wave radiation are fundamental in micrometeorology to quantify surface albedo, net radiation and the surface energy balance Eklundh et al. 2011;Cescatti et al. 2012;Stoy et al. 2013). In addition, in situ observations of spectral reflectance are increasingly applied to link surface measurements of energy and carbon fluxes with Earth observations (Gamon et al. 2006). The measurements of hyper-spectral surface reflectance at eddy covariance sites is an important step for the spatial extrapolation of in situ measurements with remote sensing retrievals from satellite platforms (Balzarolo et al. 2011) and to parameterize light use efficiency models for the up-scaling of gross primary productivity (GPP) (Sims et al. 2006;Meroni et al. 2009;Meroni et al. 2011;Peñuelas et al. 2011;Rossini et al. 2012).
To fulfill the objectives listed above, radiometric measurements have to be integrated with observations of surface fluxes performed with the eddy covariance technique. This methodology is nowadays applied for the quantification of the turbulent exchange of carbon, water and energy between vegetation and the atmosphere at hundreds of experimental sites organized in continental networks (Baldocchi 2003). Given the typical inhomogeneity of the land surface in topography, soil properties and vegetation cover, an evaluation of the spatial footprint of the observations is required for the proper integration of hemispherical radiometric measurements with turbulent fluxes. In fact, all components in a measurement system should be well matched to each other in order to make their retrievals consistent and usable in a common contest (Schmid 1997).
The source area of eddy covariance measurements has been largely investigated with both modeling and experimental approaches in order to characterize the dependence of the footprint size on measurement height, canopy structure and atmospheric conditions Rannik et al. 2003;Vesala et al. 2004;Marcolla and Cescatti 2005). On the contrary, the spatial footprint of hemispherical radiometric measurements over vegetated surfaces has not yet been formalized nor geometrically described. The quantitative description of the hemispherical radiometric footprint and its comparison with the aerodynamic footprint are ultimately required for the correct interpretation of flux measurements and for optimizing the setup of eddy covariance stations located in spatially heterogeneous landscapes. In fact, considering that the sources of both turbulent and radiative fluxes are often spatially inhomogeneous, the inability to match the source areas of turbulent fluxes with that of the upwelling radiative fluxes introduces a random error in the energy balance closure at flux sites . As far as carbon fluxes are concerned, the mismatch of radiometric and turbulent fluxes may lead to uncertainties and biases in the parameterization of models for the prediction of GPP based on spectral indexes.
The objective of the present work is to provide a formal description of the hemispherical radiometric footprint over a vegetation cover. As example, four idealized plant types are investigated to show the dependence of the radiometric footprint on major structural properties of the canopy, like leaf area index (LAI) and canopy height. Finally, the hemispherical radiometric footprint of the idealized cases is compared with the aerodynamic footprint, in order to highlight their differences in size, shape and dependence on vegetation properties.
The formalization of the radiative footprint over plant canopies is currently lacking in the literature, whereas a considerable number of articles have focused on the formalization of the footprint of turbulent fluxes. The quantitative description of the radiative footprint is a prerequisite to perform a proper comparison with the turbulent footprint and, therefore, to detect the occurrence of footprint mismatch and to elaborate alternative solutions, either to avoid the problem (e.g. with a proper choice of measurement sites) or to face it with alternative technical solutions.

Radiometric footprint on bare soil
According to Lambert's cosine law of radiation, the radiation flux (F R ) on an inverted flat plate radiometer received from a solid angle dω can be expressed as (Schwerdtfeger 1976;Schmid 1997): where I o is the normal component of the specific radiation intensity originated from a differential source area dA and dω is the differential solid angle from which dA is seen from the center of the radiometer plate.
Considering that the radiation flux originated from a solid angle dω on a sensor with hemispherical field of view is independent of the azimuthal angle φ between the sensor and the area for the specific design of the sensor itself, it is convenient to integrate the radiance over concentric annulus centered at the sensor location ( Fig. 1) and to express the signal as a mere function of the nadir angle θ.
For this purpose, the solid angle dω subtended by a differential annulus of radial thickness dθ at a nadir angle θ can be calculated as: Assuming that the observed surface is homogeneous and Lambertian (i.e. I o is independent of (θ, φ)), it is possible to model the contribution to the total signal of the radiance coming from a differential annulus by combining Eqs. (1) and (2) (Schwerdtfeger 1976;Meroni et al. 2011): where cosθ represents the cosine correction of the radiation intensity due to the vertical angle between the radiometer plate and the source, and sinθ accounts for the area of the differential annulus at varying nadir angles (Fig. 2a). It follows that the total flux from the entire hemisphere can be obtained integrating Eq. (3) from 0 to π/2 obtaining F R,tot = π I o ; hence, the total integrated footprint function for radiation is: while the normalized cumulative footprint function F for a given θ can be expressed as ( Fig. 2c): Given that θ ¼ tan −1 x h s ; where x is the radial distance and h s is the sensor height ( Fig. 1), the cumulative footprint function can be expressed in terms of radial distance x as follows ( Fig. 2d): whereas its derivative with respect to x is equal to (Fig. 2d): Equation (7) is therefore the analytic formulation of the hemispherical radiometric contribution to a radiative flux measurement from a differential annulus over a flat surface as a function of the radial distance and can be expressed in terms of the nadir angle θ as ( Fig. 2c): The proposed formalization focuses on hemispherical radiometric sensors for their importance in the measurement of surface energy fluxes (e.g. short-wave and long-wave upwelling radiation, net radiation, albedo). However, this analytical framework can be adapted to any narrower fore-optics by changing the integration domain.

Radiometric footprint over plant canopies
If a porous media, like a vegetation cover, is interposed between the sensor and the ground, the radiometric signal includes the contribution of both the canopy and the underlying soil in proportions determined by the structural characteristics and density of the vegetation layer.
Assuming a homogeneous distribution of phytoelements in the canopy space, the contributions to the radiometric signal (f l_i ) of a horizontal layer of leaves (l_i) at height z and of the soil (f s ) can be quantified by combining Eq. (8) with the probability of light interception by plants and soil as derived from a Markov light interception model (Nilson 1971). For this purpose, the canopy is discretized in n horizontal layers. The soil contribution accounts for the probability of the light reflected (or thermal radiation emitted) by the soil to reach the inverted flat plate radiometer, which is the complement to one of the canopy interception.
The interception probability at height z i can be formulated as: where LAI is the leaf area index, cumulated from the canopy top down to the height z i within the canopy, G is the leaf angular distribution and Ω is the clumping coefficient (Nilson 1971). In the present analysis, we assumed a spherical distribution of the leaf angles, and therefore, G was set constant and equal to 0.5, whereas the leaf area was uniformly distributed along the vertical profile, and therefore, the cumulated LAI increases linearly between the canopy top and the canopy bottom. The values used for Ω are reported in Table 1. The interception probability in the i-th layer is equal to: and the contribution of an horizontal annulus in the i-th canopy layer to the total source area expressed as a function of θ is: being h b and h c the canopy bottom and top heights, respectively.
Since the same θ corresponds to different horizontal distances x at different heights, in order to compute the annular footprint functions at the same radial distances in every canopy layer, Eq. (11) was interpolated using the R function approx (https://www.r-project.org/) with the linear method and recalculated at fixed x distances. The total canopy contribution was eventually calculated summing up all layer contributions at each distance: Since the total canopy interception is: where LAI TOT is the total canopy LAI, the contribution of a the concentric annulus of the soil surface to the radiometric source area can be calculated as follows: f s (x) was rescaled in order that its numerical integral over x equals the integral of f s (θ) over θ in [0, π/2], while f c (x) was rescaled in order that its integral over x equals the integral of 1 − f s (θ) over θ in [0, π/2].  Finally, the total annular footprint function is calculated summing Eqs. (12) and (14) after rescaling: For the sake of simplicity, the formulation presented in Eq. (9) is based on a 1D turbid medium analogy (Nilson 1971) and assumes that leaves are optically black, therefore ignoring the radiation scattered by vegetation. These assumptions allow the analytical formalization of the problem. In addition, the idealized 1D canopy structure is consistent with the assumptions made in the Lagrangian turbulent footprint model used in the present work. However, the numerical solution of the equation could make use of more complex radiative transfer models that account for horizontal heterogeneity in the canopy and leaf scattering (e.g. 3D turbid medium or ray tracing models; Cescatti and Niinemets 2005) to calculate the probability of light interception in canopy layers (P_l i , Eqs. (10) and (11)).
In the present analysis, four idealized plant canopies were used to explore the geometry of the radiometric footprint for different architectures, as described by the structural parameters reported in Table 1. In addition, for one short canopy (grassland, GRA) and for one tall canopy (deciduous broadleaf forest, DBF), the radiometric footprint was calculated for variable LAI in the range 1-6 and for variable canopy height in the range 0.05-0.2 h s for GRA and 0.5-0.85 h s for DBF. The positions of the 50th and 80th percentiles were calculated for every pair of parameters (LAI, canopy height).

Footprint of turbulent fluxes
For the idealized canopies described in Table 1, the footprint of turbulent fluxes was estimated for neutral condition with the Lagrangian stochastic trajectory model proposed by Thomson (1987). For this purpose, the emitted tracer was simulated as the release of a large number of passive particles, which were assumed to follow perfectly the flow. The diffusion of the scalar was described by a stochastic differential equation (a generalized Langevin equation) which determines the motion of a Lagrangian particle in space and time (Thomson 1987;Kurbanmuradov et al. 2001;Markkanen et al. 2003;Rödenbeck et al. 2003;Kljun et al. 2004;Hsieh and Katul 2009).
The analytical second-order closure model of Massman and Weil (1999) was used to predict the vertical profile of mean wind speed, second moments and the dissipation rate of turbulent kinetic energy, accounting for the structural properties of the canopy as described by Marcolla et al. (2003). The model gives an analytical expression for the second-order moments as functions of the cumulative leaf drag area, coupling the exponential wind profile after Albini (1981) with the equations for the variances of the wind components reported in Wilson and Shaw (1977), under the assumptions of horizontal homogeneity of the canopy, steady-state condition for the wind field and neutrality. Particles were released at the h = 0 to obtain the soil footprint function or uniformly released in the canopy layers to obtain the canopy footprint function.
Radiometric and flux footprint functions, calculated separately for canopy and soil fluxes, were finally compared for two of the idealized canopies (DBF and GRA) and in particular the positions of the 50th and of the 80th percentiles were calculated for LAI values ranging from 1 to 6 m 2 m −2 .
3 Results and discussion

Radiometric footprint
The theoretical case of the hemispherical radiometric footprint over a flat bare soil is presented in Fig. 2 as a function of the zenith angle (a,c), and of the radial distance (b,d). The figure shows that 50% of the total signal comes from a cone limited by a zenith angle of 45°which corresponds to a horizontal distance x = h s , and 80% of the radiometric signal originates within a distance of about 2 h s (θ~63.4°). The maximum contribution to the measured signal comes from the annulus located at a zenith angle of 45°.
When a vegetation layer is located between the sensor and the soil, the shape and size of the radiometric footprint change considerably. In Fig. 3, the contribution to the total source area by canopy layer annuli at increasing distance from the sensor is presented for the idealized ENF canopy, both as a function of the zenith angle and of the radial distance. The image clearly shows that the higher canopy layers have a larger impact on the radiometric footprint, while the position of the peak and 80% percentiles are rather similar between layers.
The radiometric annular footprint was then compared with the 1D aerodynamic footprint in neutral conditions in order to investigate the mismatch between the two (Fig. 4). The analysis was performed for the four model canopies described in Table 1. For all the canopies, the maximum contribution to the measured signal comes from a distance that is less than one radiometer height, and this distance is lower for high canopies if compared to short ones (0.25 and 0.5 h s , respectively), while in the case of the turbulent flux footprint, the order of magnitude is of few sensor heights (2 h s for the two tall canopies, and 5 and 10 h s for the idealized cropland and grassland, respectively). The mismatch between the two footprints is more evident for short canopies, whereas both the aerodynamic and radiometric footprints are more leptokurtic for tall canopies.
As far as the cumulative footprint is concerned, while 100% of the signal recorded by a radiometer comes from a distance within few sensor heights, at a distance of 20 h s , the aerodynamic signal reaches a percentage of 70% for CRO, 50% for GRA, and 88% for ENF and DBF. In addition to the size, also the shape is rather different between the aerodynamic and the radiative footprints. In fact, the footprint of turbulent fluxes is elongated in the wind direction, and therefore the cumulative footprint depends on the wind rose (Kljun et al. 2015). On the contrary, the radiative footprint in case of flat terrain is azimuthally symmetric relative to the vertical axis passing by the radiometer, whereas it becomes elliptical on slopes.
The mismatch of source areas of radiometric and aerodynamic flux measurements becomes important if the surface underlying the instruments is spatially heterogeneous, and when measurements coming from instruments characterized by different footprint are combined (e.g. in the calculation of the energy balance). This heterogeneity has been recognized as one of the potential causes of the lack of energy balance closure observed at most eddy covariance sites Stoy et al. 2013). Schmid (1997) reports a detailed discussion of the issues involved in matching scales of observations, in particular the different geometry and temporal dynamics of the source areas for radiometers and turbulent flux footprints. From the present analysis, in conditions of neutral stability, it turns out that the radiometer should be installed about 15 and 9 times higher than the eddy covariance system for GRA and CRO and about 6 times higher for ENF and DBF in order to match the distance of the radiometric and aerodynamic footprint at 80% of the total signal. Of course the problem of the temporal dynamics and directionality of turbulent flux footprint would still persist. If the wind has a dominant direction and therefore the turbulent footprint is not symmetric with respect to the tower location, an alternative method to maximize the overlap would be to center the radiometer with respect to the turbulent footprint instead of locating it on the flux tower. Alternatively, the use of multiple radiometer sensors distributed in the turbulent footprint could be envisaged.
For the radiometric footprint, the dependence of the 50th and 80th percentiles on canopy height and LAI has been investigated, and the results are reported in Fig. 5 for a short canopy (GRA) and a tall canopy (DBF). For both vegetation types, the footprint size decreases at increasing LAI and at increasing canopy height. As expected, the distances are slightly higher for the short canopy if compared to the tall one. Figure 6 reports the trend of the horizontal distance of the 50th and 80th percentiles of the radiometric (Eqs. (12) and (14)) and aerodynamic footprints as a function of LAI for GRA and DBF vegetation canopies; soil and canopy footprints are reported separately. For DBF, both the aerodynamic and the radiometric soil footprint are larger than those of the canopy. The size of the radiometric footprint shows a monotonic decreasing trend with LAI, whereas for the aerodynamic footprint, the distance of both the 50th and 80th percentiles peaks at intermediate values of LAI and then decreases. For the short canopy, given the proximity of vegetation layers to the ground, soil and canopy aerodynamic footprints are almost undistinguishable and increase with LAI. On the contrary, the size of the radiometric footprint for soil is smaller than that for the canopy and shows a monotonic decreasing trend. Hence, at increasing LAI, one should expect an increasing source area for the turbulent fluxes and a decreasing source area for radiometric measurements, which means an increasing spatial mismatch between the two set of observations. The complex interplay between LAI and aerodynamic footprint depends on several factors such as the impact of LAI on roughness length, and on the distance between the sources and the sensors. The present analysis shows that the radiometric footprint is systematically smaller and may have a different dependence on LAI changes than the turbulent footprint. For this reason, we may foresee that in the presence of a LAI seasonality, the temporal evolution of the footprint during the growing season can be different for turbulent and radiometric footprints.

Conclusions
In the last two decades, the network of eddy covariance sites has rapidly expanded to sample the natural variability of landscapes. The spatial heterogeneity of the source intensity at most flux sites raises the issue of matching the source areas of the different measuring systems, given that in several scientific domains, multiple measurements need to be combined. In this context, typical examples are the surface energy balance or the use of optical proximal sensing of canopy reflectance to predict carbon fluxes with different vegetation indexes and metrics (e.g. NDVI, PRI, sun induced fluorescence) (Gamon 2015). Considerable effort has been put so far in formalizing the footprint of turbulent flux measurements. On the contrary, the issue of the radiometric footprint over vegetated canopies has been rarely tackled. Quantifying and characterizing the extent of both footprints is required for the formalized analysis of the two and for a proper comparison that could: (a) help to avoid the problem of footprint mismatch by driving the selection of the experimental site and (b) help to find proper technical solutions to maximize the overlap of the two footprints.
In this study, an analytical formulation of the hemispherical radiometric footprint in presence of vegetation has been proposed. The manuscript focuses on the hemispherical sensor because this geometry is required for several micrometeorological measurements (e.g. short-wave and long-wave upwelling radiation, net radiation, albedo). The same analytical Fig. 4 Radiometric and aerodynamic footprint as functions of the radial distance normalized by canopy height for the model canopies described in Table 1. Probability density functions (top row) and cumulative functions (bottom row) are reported Fig. 5 Dependency of the normalized radial distance of the 50th and the 80th percentiles on LAI and canopy to sensor height ratio, for one short canopy (GRA) and one tall canopy (DBF) Fig. 6 Horizontal distance of the 50th and 80th percentiles at increasing LAI for the aerodynamic and radiometric footprints of the canopy and the soil, for a short canopy (GRA) and a tall canopy (DBF) framework could be adopted for narrower fore-optics by changing the integration domain of the equations.
Four model canopies (two short and two tall canopies) were used in the analysis of the dependence of the radiometric footprint source area on the major geometrical features of the canopy. It was observed that if turbulent flux and radiation measurements are performed at the same height, radiometric footprint is always considerably smaller than the aerodynamic one. In fact, close to 100% of the signal recorded by a radiometer comes from a distance equal to few times the sensor heights, whereas at a distance of 20 sensor heights, the aerodynamic footprint reaches a percentage of the total signal of about 70% for CRO, 50% for GRA, and 88% for ENF and DBF. In order to have a comparable extension of the radiative and aerodynamic footprint, the radiometer should be installed about 15 and 9 times higher than the eddy covariance system for the idealized GRA and CRO and about 6 times higher for ENF and DBF. The size of the radiometric footprint always decreases at increasing LAI, while the aerodynamic footprint increases for GRA and shows a peak at intermediate LAI for DBF. These opposite trends are particularly relevant for deciduous canopies, for which the strong seasonality in LAI ultimately affects the mismatch between the footprint area of aerodynamic and radiometric measurements.
An important finding emerging from this study concerns with the footprint size of soil fluxes versus canopy fluxes. Our analysis shows that for short canopies the radiometric footprint of the soil is systematically smaller than that of the canopy, while for tall canopies the aerodynamic footprint of soil fluxes is typically larger than that of the canopy . The implication of these trends is that the footprint mismatch maximizes for the soil fluxes and is therefore particularly relevant for the processes characterized by large radiative or turbulent fluxes stemming from the ground (e.g. soil respiration, broadband albedo in case of snow-covered forest floors). In conclusion, the formal description of the hemispherical radiometric footprint and its comparison with the aerodynamic footprint provided in this study may help in the design of eddy covariance stations and the integration of flux measurements from multiple sensors and techniques.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.