Comparison of Sensible Heat Fluxes Measured by a Large Aperture Scintillometer and Eddy Covariance System over a Heterogeneous Farmland in East China

The sensible heat is an important component in surface energy partitioning over the land surface. This paper compared the sensible heat fluxes measured by a large aperture scintillometer system (LAS) and an eddy covariance system (EC) over a rice paddy with a patch of mulberry seedlings in the east China coastal region during the period from 13 September–11 October 2015. During the observation period, easterlies and northerlies prevailed, and 96% easterlies and northerlies had a speed of 0–6 m s−1. The sensible heat fluxes measured by the two systems reflected that the value of HLAS generally was inclined to be larger than HEC with the average difference of 20.30 W m−2, and the uncertainty for two instruments was less than 17 W m−2. Analysis of the average footprint resulted that the mulberry seedling field always had a higher contribution to LAS than that to EC, which could be the reason that HLAS was always larger than HEC. During the days when the contributions of the mulberry seedling field to the two systems were close to each other, the sensible heat flux measurements of the two instruments were similar. The case analysis on typical sunny days showed that there would be larger sensible heat fluxes over the mulberry seedling field than in the rice paddy field especially under larger net radiation conditions.


Introduction
Farmland is a typical underlying surface and widely distributed around the world.In China, due to the family contract responsibility system, farmlands are divided into small pieces allocated to different families, and different crops may be planted within different farmland pieces upon the decision of the owners (e.g., Figure 1).Measurement of the sensible heat fluxes over these mixed-crop farmland fields can help to understand the energy distribution, which plays an important role in the local weather, water cycle and climate [1][2][3][4].
Traditionally, eddy covariance (EC), which is based on direct measurements of the product of vertical velocity fluctuations (w') and scalar concentration fluctuations (c'), is considered as the standard method for sensible heat fluxes measurement [5][6][7][8][9].However, recently, a new instrument, the large aperture scintillometer (LAS), was used in a number of field experiments and gradually widely used [10][11][12][13].For LAS, sensible heat fluxes are derived from the line-integral of the structure parameter of the refraction index; the average of sensible heat fluxes is obtained over an area formed by the path length of the light beam of the scintillometer and a line in the upwind direction.In contrast to the EC system, the measurement range of LAS is extended to a few kilometers, which coincides with the common resolution of atmospheric numerical models and satellite remote sensing retrievals (i.e., several kilometers).The application of LAS over different surfaces has been evaluated in a number of experiments [11,[14][15][16][17][18][19][20][21][22][23], and in some of the experiments, the LAS-derived sensible heat fluxes were compared with EC-derived sensible heat fluxes [14][15][16]20,[24][25][26][27][28].Over either a homogenous or a heterogeneous surface, these studies demonstrated that the fluxes measured by these two methods had a good consistency, and LAS was a reliable method for deriving the sensible heat fluxes over heterogeneous complex topography.
Atmosphere 2017, 8, 101 2 of 21 over these mixed-crop farmland fields can help to understand the energy distribution, which plays an important role in the local weather, water cycle and climate [1][2][3][4].
Traditionally, eddy covariance (EC), which is based on direct measurements of the product of vertical velocity fluctuations (w') and scalar concentration fluctuations (c'), is considered as the standard method for sensible heat fluxes measurement [5][6][7][8][9].However, recently, a new instrument, the large aperture scintillometer (LAS), was used in a number of field experiments and gradually widely used [10][11][12][13].For LAS, sensible heat fluxes are derived from the line-integral of the structure parameter of the refraction index; the average of sensible heat fluxes is obtained over an area formed by the path length of the light beam of the scintillometer and a line in the upwind direction.In contrast to the EC system, the measurement range of LAS is extended to a few kilometers, which coincides with the common resolution of atmospheric numerical models and satellite remote sensing retrievals (i.e., several kilometers).The application of LAS over different surfaces has been evaluated in a number of experiments [11,[14][15][16][17][18][19][20][21][22][23], and in some of the experiments, the LAS-derived sensible heat fluxes were compared with EC-derived sensible heat fluxes [14][15][16]20,[24][25][26][27][28].Over either a homogenous or a heterogeneous surface, these studies demonstrated that the fluxes measured by these two methods had a good consistency, and LAS was a reliable method for deriving the sensible heat fluxes over heterogeneous complex topography.Field experiments were conducted over different farmlands, such as arid grassland [8] and rice paddy [6,[29][30][31], with EC systems over the world.In this study, combined LAS and EC measurements were conducted over the selected heterogeneous farmland in Jiangsu province of China during the period from 13 September-11 October 2015.First, LAS sensible heat fluxes were compared with the EC system here to investigate their reliability and characteristics.Second, average footprints were calculated to analyze the statistical characteristic.Third and the last, the footprint model was used to assess LAS and EC measurements to understand the influence of the heterogeneity in typical sunny days.It is expected that the results found here will provide a reference for sensible heat fluxes' analysis and/or LAS experiments over heterogeneous farmlands.
The study is organized as follows: Section 2 describes the observational site and micrometeorological measurement systems; Section 3 presents the data processing methods; Section 4 analyzes the results; and Section 5 summarizes and concludes the paper.

Observation Site
The experimental site is at a farmland about 45 km from the east China coast (32.76 • N, 120.47 • E; 4 m above sea level).The site was generally flat, and around the site, there are some farmhouses (Figure 1).An LAS transmitter unit (marked A in Figure 1), a receiver unit (marked B in Figure 1), EC instruments (marked C in Figure 1) and an automatic meteorological station (AWS marked D in Figure 1) were installed at the site.At three meters east of the LAS optical path (the A-B line in Figure 1), there is a patch of mulberry seedlings with a north, west and east side length of 175 m, 175 m and 125 m, respectively.Additionally, there is a 70 m × 62 m poplar grove at the west of the LAS optical path (and at the north of the EC system).The rest of the areas are rice paddy fields.The experiment was carried out from 13 September-11 October 2015, and continuous sensible heat fluxes were derived by both the LAS and EC systems.During the observation period, the canopy height increased from 0.5 m-0.6 m for the rice paddy, from 0.15 m-0.4 m for the mulberry seedling field, and it remained about 12 m for the poplars in the grove, respectively.

Large Aperture Scintillometer
The LAS instrument is manufactured by Rainroot Scientific Limited China (Beijing, China), and it has a similar observation principle to the LAS from Kipp&Zonen, but independently designed [32], and has an aperture diameter (D) of 0.15 m, operating at a near-infrared wavelength of 880 nm.The LAS receiver (marked B in Figure 1) and transmitter unit (marked A in Figure 1) were deployed at a height of 6 m on the tower over the paddy.The electromagnetic radiation was transmitted from north to south with an optical path of 470 m.Such a short path is caused by the limitation of the position of local villagers' houses.

EC System
The EC instruments (marked C in Figure 1) were about 100 m to the west of the LAS optical path, including a H 2 O analyzer (LI-7500, LI-COR Biotechnology, Lincoln, NE, USA) and a three-dimensional sonic anemometer (CSAT3, Campbell Scientific Inc., Edmonton, AB, Canada).The instruments were installed on a tower at the height of 10 m, and they measure three-dimensional wind speed components u, v, w, theta, the density of air, etc., at a sampling rate of 10 Hz.

Automatic Meteorological Station
The automatic meteorological station (AWS, marked D in Figure 1) was about 185 m to the east of the LAS optical path.The air temperature, humidity and wind speed sensors were placed at each height (3 m, 5 m, 8 m and 10 m) on the tower.A wind direction sensor was deployed at the height of 10 m.The radiometer, measuring upward and downward short and long wave radiations, was installed at a height of 3 m.Each of these sensors sampled at a rate of 1 Hz, and the measurements were averaged at 30-min intervals.

Parameter Set
The aerodynamic roughness length (z 0 ) and the zero plane displacement height (d) were calculated through the gradient method with EC and AWS data.Based on the Monin-Obukhov (M-O) similarity hypothesis, the vertical wind profile is expressed as follows: where z is the height with respect to wind speed u, L is the Obukhov length derived by EC systems, κ(= 0.4) is the von Karman constant and u * is the friction velocity.The integrated stability correction function Ψ m was given by Dyer [33]: (3) where Equations ( 2) and ( 3) are valid under unstable conditions and Equation ( 4) is valid under stable conditions.Then, the fitting process was conducted on AWS-measured wind speed data, and the results showed that the mean of z 0 was 0.025 m; and d varied from 0.43 m on 13 September to 0.55 m on 11 October 2015.Besides, owing to the flat underlying surface, the LAS effective height (z eff ) was estimated to be its installation height of 6 m.

Data Process Scheme of EC
In this paper, we used EddyPro 5.2.1 (software developed by LI-COR Biotechnology) to process data measured by the EC system.The EddyPro software applied the following corrections: despiking algorithm [34], spectral corrections [35,36], compensation for density fluctuations [37], time lag compensation, double rotation for tilt correction, block averaging, statistical tests [34], etc.Then, the half-hour averages of heat fluxes were derived.As the sensible and latent heat fluxes were calculated, the Bowen ratio β could be derived as: where H EC is EC-derived sensible heat fluxes, and LE EC is EC-derived latent heat fluxes.For the missing data or the results that β > 3 (unstable phase) during the observation period, the Bowen ratios were set up as the average β (=0.22).Then β can be used in the calculation of sensible heat fluxes from LAS.

Scintillometry Method
The intensity of fluctuations in the refractive index of air, which can be converted to the structure parameter of the refractive index (C 2 n ) turbulence [38], could indicate the variation of atmosphere turbulence.LAS measures the intensity of optical fluctuations as the laser goes through a turbulent atmosphere.C 2 n can be expressed as the function of the structure parameter of temperature (C 2 T ) and humidity (C 2 q ) [39].Additionally, it is sensitive to the fluctuations of temperature more than humidity at a near-infrared wavelength [40].Further, Odhiambo and Savage [41] verified that the sensible heat fluxes derived from the surface layer scintillometry (one kind of scintillometer that is similar to LAS in sensible heat fluxes' measurements between 50 and 350 m) with the β-corrected C 2 T has greater agreement with the EC-derived sensible heat fluxes than that without β correction.Following Odhiambo and Savage [41], C 2 n is expressed as: where T is the air temperature (K), P is the atmospheric pressure (Pa) and β is the Bowen ratio derived from EC. Normally, the Bowen ratios were smaller than 1 (mean of 0.22 for this study) under rice paddy.The calculated C 2 T could overestimate about 30% without the Bowen ratio correction.Therefore, this correction could not be neglect in this paper.The half-hourly average C 2 n was calculated to match the averaging period of EC, and the noise was eliminated by using a criterion of X(t) < X − 4σ or X(t) > X + 4σ, where X(t) is the measurement, X is the mean over the interval and σ is the standard deviation.Besides, the mathematical relationship of C 2 n and the variance of the natural logarithm of light intensity was not established, and the measurement should be eliminated when the LAS signal is saturated (with higher turbulence intensity).Based on Ochs and Hill [42], the threshold value of C 2 n for LAS signal saturation is set to 5.87 × 10 −12 m −2/3 in our experiment.
Then, according to the M-O similarity hypothesis, the following equation can be derived: where T * is the surface layer temperature scale parameter, z eff is the LAS effective height, d is the zero plane displacement, L is the Obukhov length and f T is a universal dimensionless function of the temperature structure parameter.Several forms of f T have been proposed from different experiments [43][44][45][46][47][48][49].This paper adopts the function of Andreas (1988) [45]: for the unstable condition, (8) ) for the stable condition.
The friction velocity u * and the Obukhov length L were calculated as: where κ(= 0.4) is the von Karman constant, and g(= 9.8 m s −2 ) is the acceleration of gravity.u is the wind speed at the height of z eff .Owing to the different heights of LAS and EC, u was obtained through neutral-condition log-law correction of the EC-measured wind speed.The surface layer is considered as unstable when these following inequalities were satisfied: where R n is net radiation.The surface layer parameters T * , u * and L are solved using an iterative calculation of Equations ( 2)-( 4) and ( 6)-( 11), then area-averaged H LAS can be calculated through the T * (based on C 2 n ) and u * (sensitive to the local roughness): where ρ is the air density and C P is the specific heat capacity of air at constant pressure.Similar to Gruber and Fochesatto [50], the convergence of this iterative calculation was identified through the fixed-point method [50].Combining Equations ( 2)-( 4) and ( 6)-( 11), the recursive function was: where under unstable conditions: Tκu 2 ( 15) Equation ( 14) using fixed-point recursion under unstable conditions is seen in Figure 2. As the relationship of M and z e f f − d /L is monotonic, the iterative calculation was applicable to derive sensible heat fluxes from LAS measurement under unstable conditions.
where ρ is the air density and P C is the specific heat capacity of air at constant pressure.
Similar to Gruber and Fochesatto [50], the convergence of this iterative calculation was identified through the fixed-point method [50].Combining Equations ( 2)-( 4) and ( 6)-( 11), the recursive function was: where under unstable conditions: Equation ( 14) using fixed-point recursion under unstable conditions is seen in Figure 2. As the relationship of M and ( ) is monotonic, the iterative calculation was applicable to derive sensible heat fluxes from LAS measurement under unstable conditions.

Calculation Uncertainty
The uncertainty of EC-derived sensible heat fluxes was also calculated through Finkelstein and Sims' [51] method; the results under unstable conditions are shown in Figure 3a with the increase of H EC ; the uncertainty of H EC could increase from almost zero to about 17 W m −2 ; and the mean uncertainty was 5.46 W m −2 .
Even if turbulence is being sampled above an extremely flat field, for LAS, the uncertainty in z will still be present.The uncertainty of LAS should be assessed before comparison.The sensitivity function for the sensible heat flux H LAS under unstable conditions was given by Gruber et al. [52]: where z(u) is the height of the beam along the relative path position u.The weight function from Hartogensis et al. [53] is as follow: where L is the optical path length, K = 2π/λ is the optical wavenumber, k the turbulent spatial wavenumber, φ n (k) = 0.033k 11/3 the three-dimensional spectrum of the refractive index in the inertial range and J 1 (x 1 ) and J 2 (x 2 ) are Bessel functions of the first kind, first order with x 1 = kDu/2 and , where D is the aperture diameter.Combining Equations ( 16)-( 19), the uncertainty of LAS-derived sensible heat fluxes could be calculated through:

Footprint Analysis
The flux footprint provides a means of estimating both the source area and the relative contribution of each surface element to the measured fluxes.In this paper, the flux footprint model was calculated through [54]:

Footprint Analysis
The flux footprint provides a means of estimating both the source area and the relative contribution of each surface element to the measured fluxes.In this paper, the flux footprint model was calculated through [54]: where G(u) is the weight function and x 1 and x 2 are the position of transmitter and receiver of LAS, respectively.(y) and x (y ) are the point in the optical path and in the up-wind areas, respectively.f is the point-flux footprint model from Kormann and Meixner [55].

Meteorological Conditions
Temporal variations of wind speed, air temperature, relative humidity, surface atmospheric pressure and precipitation from 12 September to 11 October 2015 was showed in Figures 4 and 5.During the observation period, the wind speed (WS) showed obvious diurnal variation, and the maximum WS occurred at noon every day with a value varying from 1.95-10 m s −1 on different days (Figure 4a).Between the height of 10 m and 3 m, the WS difference was always less than 1 m s −1 .The overall average WS was 2.43 m s −1 at 10 m, 2.21 m s −1 at 8 m, 1.94 m s −1 at 5 m and 1.72 m s −1 at 3 m, respectively.The highest value of WS was 10 m s −1 at the height of 10 m at 10:30 on 29 September.Combined with Figure 5, it is easy to find that the prevailing winds were easterlies at the observation height of 10 m in this site.During the daytime, more than 62 percent of winds were easterlies (including southeasterlies and northeasterlies), and only 21 percent of winds were westerlies (including southwesterlies and northwesterlies).At night, about 51 percent of winds were from the east, and 26 percent were from the west.The air temperature (Figure 4b) also significantly varied with a diurnal cycle.Most of the values were between 283 K and 298 K, and the overall average air temperature was 291.62 K at 10 m, 291.75 K at 8 m, 291.95K at 5 m and 292.03K at 3 m, respectively.As the height increased, the air temperature slightly decreased, but the temperature difference between the highest and lowest levels was less than 0.5 K at daylight.Driven by the variation of the surface net radiation, the highest and lowest temperature often occurred at 14:00 and 2:00, respectively.The maximum and minimum temperatures were 298.17K and 283.88 K (caused by a strong precipitation process) at 14:30 on 18 September and at 5:30 on 2 October, respectively.Relative humidity (RH) also varied largely between mostly 50% and 100% during a day.Generally, the maximum RH (Figure 4c) followed the minimum air temperature at about 2:00, while the minimum RH corresponded to the maximum air temperature at about 14:00.It is noteworthy that RH was also influenced by other factors such as precipitation, advective flow, transpiration and turbulent intensity.The variation of atmospheric pressure is presented in Figure 4d.During the observation period, the value of atmospheric pressure varied from 1020.92 hPa (at 9:30 on 13 September) to 1006.44 hPa (at 12:30 on 23 September) and then increased to 1019.86 hPa (at 10:30 a.m. on 29 September).Later, because the air mass transited from northwest, the pressure went sharply down to 1005.18 hPa (at 4:00 a.m. on 1 October) and finally picked up to 1025.19 hPa (at 9:30 a.m. on 4 October).Combined with the wind speed and wind direction from 29 September-1 October, it could be found that a low pressure moved from west to the east of the site.Meanwhile, the values of pressure varied significantly and had two maximums and minimums during the diurnal cycle.The first low value and high value were usually at about 3:00 a.m. and 10:00 a.m.; the second low value and high value were usually at about 14:00 and 20:00.The temporal variation of precipitation is shown in Figure 4e; there were four obvious precipitation processes: 22-23 September (32.9 mm), 29 September (24.3 mm), 30 September-1 October (34.3 mm) and 7-8 October (4.5 mm).
from west to the east of the site.Meanwhile, the values of pressure varied significantly and had two maximums and minimums during the diurnal cycle.The first low value and high value were usually at about 3:00 a.m. and 10:00 a.m.; the second low value and high value were usually at about 14:00 and 20:00.The temporal variation of precipitation is shown in Figure 4e; there were four obvious precipitation processes: 22-23 September (32.9 mm), 29 September (24.3 mm), 30 September-1 October (34.3 mm) and 7-8 October (4.5 mm).from west to the east of the site.Meanwhile, the values of pressure varied significantly and had two maximums and minimums during the diurnal cycle.The first low value and high value were usually at about 3:00 a.m. and 10:00 a.m.; the second low value and high value were usually at about 14:00 and 20:00.The temporal variation of precipitation is shown in Figure 4e; there were four obvious precipitation processes: 22-23 September (32.9 mm), 29 September (24.3 mm), 30 September-1 October (34.3 mm) and 7-8 October (4.5 mm).3, the uncertainty (<17 W m −2 for EC) was always smaller than the difference between LAS and EC shown in Figure 6; this may indicate that the differences are not caused by the uncertainty of two instruments.The maximums of H LAS and H EC were 125.36 W m −2 (at 10:30 on 6 October) and 93.36 W m −2 (at 14:00 on 10 October), respectively; and the averages of the daily maximum of H LAS and H EC were 80.14 W m −2 and 59.84 W m −2 , respectively.The same as other experiments [22,56], the sensible heat fluxes had a similar diurnal variation pattern with net radiation.The LE also showed an evident diurnal variation pattern with the average daily maximum of 383.71 W m −2 .
The relationship of H LAS and H EC is shown in Figure 7 with negative values omitted.The red points and the corresponding bars represent the median values binned at 5 W m −2 of the H EC interval and the interquartile ranges, respectively.The median of H LAS was 13.09 W m −2 greater than the median of H EC on average when H EC was at the range of 0-80 W m −2 .For the data where H EC exceeded 80 W m −2 , the sample numbers were too small to exactly estimate the relationship between H LAS and H EC .Overall, H EC shows that the sensible heat fluxes varied from 0 to about 80 W m −2 on the rice paddy fields in eastern China during the maturity period of rice, and it is similar to the measurement during the rice growth period [31].
always larger than the corresponding HEC, and the difference is mostly experienced in the central part of the day, while towards the beginning and end of the diurnal cycle, the fluxes converged.Considering the uncertainty of the two systems' derived sensible heat fluxes shown in Figure 3, the uncertainty (<17 W m −2 for EC) was always smaller than the difference between LAS and EC shown in Figure 6; this may indicate that the differences are not caused by the uncertainty of two instruments.The maximums of HLAS and HEC were 125.36 W m −2 (at 10:30 on 6 October) and 93.36 W m −2 (at 14:00 on 10 October), respectively; and the averages of the daily maximum of HLAS and HEC were 80.14 W m −2 and 59.84 W m −2 , respectively.The same as other experiments [22,56], the sensible heat fluxes had a similar diurnal variation pattern with net radiation.The LE also showed an evident diurnal variation pattern with the average daily maximum of 383.71 W m −2 .
The relationship of HLAS and HEC is shown in Figure 7 with negative values omitted.The red points and the corresponding bars represent the median values binned at 5 W m −2 of the HEC interval and the interquartile ranges, respectively.The median of HLAS was 13.09 W m −2 greater than the median of HEC on average when HEC was at the range of 0-80 W m −2 .For the data where HEC exceeded 80 W m −2 , the sample numbers were too small to exactly estimate the relationship between HLAS and HEC.Overall, HEC shows that the sensible heat fluxes varied from 0 to about 80 W m −2 on the rice paddy fields in eastern China during the maturity period of rice, and it is similar to the measurement during the rice growth period [31].

The Effects of Winds and Associated Footprint on HLAS and HEC
To clarify the reason that HLAS tended to be larger than HEC during the observation period, the effects of wind direction, speed and the associated footprint (the model of Kormann and Meixner [55]) are investigated in this subsection.

The Effects of Winds and Associated Footprint on H LAS and H EC
To clarify the reason that H LAS tended to be larger than H EC during the observation period, the effects of wind direction, speed and the associated footprint (the model of Kormann and Meixner [55]) are investigated in this subsection.
Figure 8a-d  It should be noted that the H LAS − H EC regression slope (1.46) was specifically large when easterlies prevailed, while when winds were from other directions, the H LAS − H EC regression slopes were closer to one (1.17,1.17 and 1.1 for westerlies, southerlies and northerlies, respectively).Therefore, the heterogeneity caused by the mulberry field, which was at the east side of the LAS and EC observation instruments, could play a significant role.
Since the easterlies and northerlies conditions dominated the observation period, the observed sensible heat fluxes and associated mulberry seedling field contributions under the two conditions are analyzed in detail here.First, it is found that when the wind direction switched from northerlies to easterlies, the average H LAS increased from 33.72-55.09W m −2 associated with the average mulberry seedling field contribution from 27.23-40.14%.Second, the average H EC increased from 25.66-31.96W m −2 associated with the average mulberry seedling field contribution from 1.57%-26.66%.Third, similar average mulberry seedling field contributions to LAS under northerlies (27.23%) and to EC under easterlies (26.66%) yielded similar average sensible heat flux observations (33.72 W m −2 for LAS and 31.96W m −2 for EC).Based on this evidence, it can be concluded that, under the same meteorological condition and for the same area, the mulberry seedling field generally produces more sensible heat flux than the rice paddy; and the different mulberry seedling field contribution to LAS and EC is the substantial cause of the difference between the H LAS and the H EC .
To clarify the physical mechanism giving the origin to the difference of two instrument under unstable conditions, the factors derived from sensible heat fluxes (shown in Equation ( 13)) were compared here under north wind and east wind conditions.Under unstable conditions, the LAS-EC regression slope (k), determination coefficient (Rs), mean square error (RMSE) and mean air density (ρ), mean friction velocity (u * ), mean temperature scale parameter (T * ) were calculated (Table 1).First, it is clear that the air density measured from two instruments has a small bias here both under east (1.20 kg m −3 for LAS and 1.18 kg m −3 for EC) and north wind conditions (1.20 kg m −3 for LAS and 1.17 kg m −3 for EC).Additionally, the differences were small such that they had a slight influence (less than 3%) on sensible heat fluxes.Second, it is found that when the wind direction switched from northerlies to easterlies, the average u * LAS (u * EC ) increased from 0.20 (0.37)-0.26 (0.38) m s −1 , and the average T * LAS (T * EC ) decreased from −0.16 (−0.067)-−0.21(−0.079)K associated with the average mulberry seedling field contribution from 27.23% (1.57%)-40.14%(26.66%).Third, under the similar average mulberry seedling field contributions (27.23% for LAS under northerlies and 26.66% for EC under easterlies), the absolute value of average T * LAS was significantly larger than the absolute value of T * EC while u * LAS was smaller than u * EC .This results in a similar value on the average sensible heat flux for two instruments.Therefore, it can be concluded that the mulberry seedling field had significant influence not only on u * , but also on T * .The difference of sensible heat fluxes derived from two instruments could be attributed to the different u * and T * from the two instruments.
evidence, it can be concluded that, under the same meteorological condition and for the same area, the mulberry seedling field generally produces more sensible heat flux than the rice paddy; and the different mulberry seedling field contribution to LAS and EC is the substantial cause of the difference between the HLAS and the HEC.To clarify the physical mechanism giving the origin to the difference of two instrument under unstable conditions, the factors derived from sensible heat fluxes (shown in Equation ( 13)) were compared here under north wind and east wind conditions.Under unstable conditions, the LAS-EC regression slope (k), determination coefficient (Rs), mean square error (RMSE) and mean air density ( ρ ), mean friction velocity ( * u ), mean temperature scale parameter ( * T ) were calculated (Table 1).First, it is clear that the air density measured from two instruments has a small bias here both under east (1.20 kg m −3 for LAS and 1.18 kg m −3 for EC) and north wind conditions (1.20 kg m −3 for LAS and 1.17 kg m −3 for EC).Additionally, the differences were small such that they had a slight influence (less than 3%) on sensible heat  Further, measurements under east and north wind conditions were selected to analyze the influence of wind speeds and the associated footprint.The wind speeds are categorized by 0-2, 2-4 and 4-6 m s −1 , while the condition is which wind speeds are larger than 6 m s −1 is not considered here due to their infrequent occurrence.For different wind speed conditions under east (north) wind directions, Figure 9a,b shows the observed H LAS and H EC , and Figure 9c-h presents the corresponding averaged footprints.Under easterlies of 0-2, 2-4 and 4-6 m s −1 , the footprint extended more eastward when the wind speed increased, and the corresponding contributions of the mulberry seedling field were 53.51% (19.21%), 41.07%(28.25%) and 35.23% (26.71%) for LAS (EC), respectively.Generally, with the increase of the speed of easterlies, the mulberry seedling field contributed less for LAS, but more for EC.When the east winds were at speed of 0-2, 2-4 and 4-6 m s −1 ; the average H LAS (H EC ) were 33.87 (21.50), 47.76 (30.17) and 69.91 (37.69)W m −2 , respectively.Following the increasing of wind speed category, the contributions from the mulberry seedling field to LAS and EC became closer to each other, but the difference between the average H LAS and H EC became larger, which was probably caused by the larger sensible heat fluxes' difference over the mulberry seedling field and rice paddies under higher wind speeds and larger average net radiation (the average net radiation was 320.05, 363.24 and 395.75 W m −2 for 0-2, 2-4 and 4-6 m s −1 ).Under northerlies, the footprint extended northward with the wind speed increase.The contributions of the mulberry seedling field were 36.41%(1.03%), 27.66% (1.24%) and 22.09% (2.28%) for LAS (EC) with wind speed increasing from 0-2 to 4-6 m s −1 .When the north winds were at speeds of 0-2, 2-4 and 4-6 m s −1 , the average net radiation was 281.92, 291.74 and 362.09W m −2 , respectively.Additionally, the average HLAS (HEC) were 32.39 (23.08), 27.78 (23.81) and 44.54 (30.91)W m −2 , respectively.The results show that the larger mean net radiation also corresponds to the higher wind speed.words, on 2 October, H LAS and H EC were generally from the homogenous rice paddies; while on 5 October, the mulberry seedling field affected H LAS and H EC greatly, and the influence was different for the two instruments.Generally, under the same meteorological condition (i.e., wind, downward short and long wave radiation), the sensible heat fluxes were homogeneous over a uniform underlying surface such as a rice paddy field [18], which is supported by the observation on 2 October when the footprints of LAS and EC were both over the rice paddy field.However, when the footprints of LAS covered more the mulberry seedling field than that of EC on 5 October, H LAS and H EC diverged drastically.On 5 October, based on the fact that the mulberry seedling field accounted for more weight in footprints for LAS than that for EC, it can be concluded that the mulberry seedling field contributed more sensible heat fluxes to LAS than that to EC.On the other hand, on 5 October, the H LAS was larger than H EC , while it can be further deduced that the mulberry seedling field generated more sensible heat fluxes than that rice paddy field.Therefore, due to the influence of the mulberry seedling field at the east side of the instruments, the sensible heat fluxes obtained by LAS were generally larger on the days with easterlies than on the days with westerlies.Alfieri and Blanken [57] found that a single point measurement (EC systems) is not sufficient to accurately describe this sand-sagebrush environment where the components of the surface energy budget can differ substantially across short distances.Meanwhile, EC could not capture the large-scale motions, which were beyond its footprint [58].The difference between LAS and EC might have been caused by the large-scale turbulence at the experimental site.with easterlies than on the days with westerlies.Alfieri and Blanken [57] found that a single point measurement (EC systems) is not sufficient to accurately describe this sand-sagebrush environment where the components of the surface energy budget can differ substantially across short distances.Meanwhile, EC could not capture the large-scale motions, which were beyond its footprint [58].The difference between LAS and EC might have been caused by the largescale turbulence at the experimental site.2).Table 2 shows that HLASavg was larger than HECavg not only on 28 September but also on 5 October, and the sensible heat fluxes obtained by the two instruments were larger on 28 September than on 5 October.Figure 11a,b further confirms that HLAS was always larger than HEC on 28 September and 2 October, but their differences decreased with larger wind speeds.During the two days (28 September and 5  2).Table 2 shows that H LASavg was larger than H ECavg not only on 28 September but also on 5 October, and the sensible heat fluxes obtained by the two instruments were larger on 28 September than on 5 October.Figure 11a,b further confirms that H LAS was always larger than H EC on 28 September and 2 October, but their differences decreased with larger wind speeds.During the two days (28 September and 5 October), the calculated footprint areas are both at the east of the instruments (Figures 10d and 11b).Still, there was the mulberry seedling field at the east of the optical path (green trapezoid in Figures 10d and 11b).During 10:00-15:00 on 28 September and 5 October, the mulberry seedling field contributed 33.78% and 45.15% to H LAS , 36.05% and 21.92% to H EC , respectively.Generally, the sensible heat fluxes were homogeneous over the rice paddy field (discussed above).Meanwhile, for EC, the mulberry seedling field contributed more sensible heat fluxes on 28 September than on 5 October, and EC obtained larger values on 28 September.Therefore, it can be reconfirmed that the mulberry seedling field generates more sensible heat fluxes than that rice paddy field.Due to the influence of the mulberry seedling field at the east side of the instruments, the sensible heat fluxes obtained by EC were increased while LAS-derived sensible heat fluxes were almost invariable under greater wind speeds conditions when the weather was dominated by prevailing easterly wind.Due to the influence of the mulberry seedling field at the east side of the instruments, the sensible heat fluxes obtained by EC were increased while LAS-derived sensible heat fluxes were almost invariable under larger wind speeds conditions when the weather was dominated by prevailing easterly wind.The different phenomenon of LAS and EC might be caused by the difference of the measurement theory between two instruments.The LAS-derived sensible heat fluxes are based on M-O similarity equations; while EC measured turbulence directly and required a certain level of mechanical turbulence present in the flow.The error by using M-O similarity equations under a heterogeneous underlying surface may affect the accuracy of the LAS-derived sensible heat fluxes, and the EC measurement under low wind conditions may also be erroneous.

Conclusions
This study compared sensible heat fluxes measured by a large aperture scintillometer and an eddy covariance system over a rice paddy in the coastal region of east China during the period from 13 September-11 October 2015.Through the comparison of the two sets of data, it

Conclusions
This study compared sensible heat fluxes measured by a large aperture scintillometer and an eddy covariance system over a rice paddy in the coastal region of east China during the period from 13 September-11 October 2015.Through the comparison of the two sets of data, it can be found that the majority of the values of H LAS were larger than H EC , and the daily maximum of H LAS (H EC ) during the period was about 100 (80) W m −2 .
First, the daytime sensible heat fluxes and corresponding average footprints under different wind direction conditions were compared.During the observation period, east wind and north wind prevailed, and the footprint analysis showed that under larger contribution conditions of the mulberry seedling field, the average H LAS was always larger than H EC .The consistent measurements of sensible heat fluxes by the two instruments often correspond to the consistent contributions from the mulberry seedling field for the two instruments' source areas.That is, the mulberry seedling field has significantly affected the sensible heat fluxes' observation of the two systems, and the influence is not only on u * , but also on T * .The difference of the sensible heat fluxes derived from the two instruments could be attributed to the different u * and T * from the two instruments.
Under easterlies and northerlies, with increasing wind speed, which always corresponded to larger net radiation, the footprint extended more eastward and northward, respectively; and the mulberry seedling field's contributions to the two instruments became more similar to each other; but the difference between H LAS and H EC became larger.It is also concluded that the difference of the sensible heat fluxes over the mulberry seedling field and over rice paddies increased with the average net radiation conditions.
Finally, the H LAS was generally larger under easterlies than under westerlies conditions.Footprint analysis showed that, under easterlies, the source areas of instrument were rice paddies together with a mulberry seedling field.Additionally, under westerlies, the source areas were basically rice field.This confirmed that the sensible heat fluxes over the mulberry seedling field were larger than the value over the rice paddies field.For LAS, the contribution of the mulberry seedling field was larger, and that made the LAS-derived sensible heat fluxes larger than those that were EC derived on east wind conditions.Through the comparison of sensible heat fluxes measured by the two instruments under different east wind speed conditions, it is found that under east wind conditions, the enlargement of H EC is consistent with the increasing contribution of the mulberry seedling field for EC.This confirms that there were larger sensible heat fluxes over the mulberry seedling field than over the rice paddy again.
This paper found that during the observation period, the sensible heat fluxes over the mulberry seedling field were larger than over the rice paddy field especially under larger net radiation conditions.Additionally, the heterogeneity caused by the different crops in nearby fields played an important role in the sensible heat fluxes' observation by the LAS and EC systems.This paper shows an example of using the two observation systems (LAS and EC) to expose the effect of surface heterogeneity.More detailed analysis of the variation of sensible heat fluxes over each kind of crop in east China will be carried out with further observation experiments.

Figure 1 .
Figure 1.(a) Map of eastern China; the location of the experimental site is shown as a red star.(b) The experimental site.Points A and B are the transmitter and receiver of the large aperture scintillometer (LAS) with an optical path of 470 m, respectively; C is the eddy covariance (EC) system; and D is the automatic meteorological station (AWS).The mulberry seedling field and grove field are also shown with green and blue frames, respectively.Concrete roads are shown as white bands.

Figure 1 .
Figure 1.(a) Map of eastern China; the location of the experimental site is shown as a red star; (b) The experimental site.Points A and B are the transmitter and receiver of the large aperture scintillometer (LAS) with an optical path of 470 m, respectively; C is the eddy covariance (EC) system; and D is the automatic meteorological station (AWS).The mulberry seedling field and grove field are also shown with green and blue frames, respectively.Concrete roads are shown as white bands.

Figure 2 .
Figure 2. Solution of Equation (14) using fixed-point recursion, and the relationship of M and ( ) eff z d L − is monotonic.

Figure 2 .
Figure 2. Solution of Equation (14) using fixed-point recursion, and the relationship of M and z e f f − d /L is monotonic.

21 Figure 3 .
Figure 3. Uncertainty of (a) EC-and (b) LAS-derived sensible heat fluxes from measurements during this observation period.

2 xFigure 3 .
Figure 3. Uncertainty of (a) EC-and (b) LAS-derived sensible heat fluxes from measurements during this observation period.

Figure 5 .
Figure 5. Distributions of (a) wind directions at 10-m height during the daytime and (b) wind directions at 10-m height during nighttime.

Figure 5 .
Figure 5. Distributions of (a) wind directions at 10-m height during the daytime and (b) wind directions at 10-m height during nighttime.

Figure 5 .
Figure 5. Distributions of (a) wind directions at 10-m height during the daytime and (b) wind directions at 10-m height during nighttime.

4. 2 .
Turbulent Heat Fluxes The variations of H LAS , H EC , LE EC and net radiation (Rn) are shown in Figure 6.The maximum of Rn varied from 200 W m −2 (on cloudy days) to near 700 W m −2 (on sunny days).The maximum of Rn in this period was 682.3 W m −2 at 11:30 on 16 September, and the average of the daily maximum of Rn was 489.3 W m −2 .Because of the occurrence of the clouds, there were always some fluctuations near the peak of Rn.The H LAS and H EC varied significantly from very low values to about 130 W m −2 and about 90 W m −2 , respectively.The values of H LAS were always larger than the corresponding H EC , and the difference is mostly experienced in the central part of the day, while towards the beginning and end of the diurnal cycle, the fluxes converged.Considering the uncertainty of the two systems' derived sensible heat fluxes shown in Figure

Figure 6 .
Figure 6.Variations of sensible heat fluxes (H), latent heat fluxes (LE) and net radiation (Rn) from 13 September-11 October 2015.The three days that are discussed in Section 5.3 are marked in yellow.

Figure 7 .
Figure 7. Variations of H LAS and H EC when the sensible heat fluxes were greater than zero from 12 September-11 October 2015.The red points represent the median values binned at 5 W m −2 of the interval.The corresponding bars represent the interquartile ranges.The sample numbers are also shown in the figure.The regression lines are plotted in blue.The slopes and determination coefficients are also shown in the picture.

Figure 8 .
Figure 8. Relationship between LAS-derived and EC-derived sensible heat fluxes under (a) east wind, (b) west wind, (c) south wind and (d) north wind conditions is shown in the left panels.The slopes and determination coefficients are also shown in the picture.The regression lines are plotted in red.The corresponding average footprints for LAS and EC under (e) east wind, (f) west wind, (g) south wind and (h) north wind conditions are shown in the right panels.The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.The numbers (n) of data points are also given in the figure.

Figure 8 .
Figure 8. Relationship between LAS-derived and EC-derived sensible heat fluxes under (a) east wind; (b) west wind; (c) south wind and (d) north wind conditions is shown in the left panels.The slopes and determination coefficients are also shown in the picture.The regression lines are plotted in red.The corresponding average footprints for LAS and EC under (e) east wind; (f) west wind; (g) south wind and (h) north wind conditions are shown in the right panels.The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.The numbers (n) of data points are also given in the figure.

Figure 9 .
Figure 9. Relationship between LAS-derived and EC-derived sensible heat fluxes under (a) east wind and (b) north wind conditions.The different wind speed was plotted in black point (0-2 m s −1 ), green plus sign (2-4 m s −1 ) and blue asterisk (4-6 m s −1 ).The slopes and determination coefficients were also shown in the picture.The regression lines were plotted in red.The corresponding average footprints when east wind was 0-2 m s −1 (c), 2-4 m s −1 (d) and 4-6 m s −1 (e).The corresponding average footprints when north wind was 0-2 m s −1 (f), 2-4 m s −1 (g) and 4-6 m s −1 (h).The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.The numbers (n) of data points are also given in the figure.

Figure 9 .
Figure 9. Relationship between LAS-derived and EC-derived sensible heat fluxes under (a) east wind and (b) north wind conditions.The different wind speed was plotted in black point (0-2 m s −1 ), green plus sign (2-4 m s −1 ) and blue asterisk (4-6 m s −1 ).The slopes and determination coefficients were also shown in the picture.The regression lines were plotted in red.The corresponding average footprints when east wind was 0-2 m s −1 (c), 2-4 m s −1 (d) and 4-6 m s −1 (e).The corresponding average footprints when north wind was 0-2 m s −1 (f), 2-4 m s −1 (g) and 4-6 m s −1 (h).The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.The numbers (n) of data points are also given in the figure.

Figure 10 .
Figure 10.Variations of HLAS and HEC on 2 (a) and 5 October (b).Footprints of sensible heat fluxes measured by EC and LAS systems on 2 October (c), and on 5 October (d).The blue rectangle represents the grove field, and the green trapezoid represents the mulberry seedling field.Additionally, the corresponding time is marked in the lower right corner of picture.

Figure 10 .
Figure 10.Variations of H LAS and H EC on 2 (a) and 5 October (b).Footprints of sensible heat fluxes measured by EC and LAS systems on 2 October (c), and on 5 October (d).The blue rectangle represents the grove field, and the green trapezoid represents the mulberry seedling field.Additionally, the corresponding time is marked in the lower right corner of picture.
field.Due to the influence of the mulberry seedling field at the east side of the instruments, the sensible heat fluxes obtained by EC were increased while LAS-derived sensible heat fluxes were almost invariable under greater wind speeds conditions when the weather was dominated by prevailing easterly wind.Due to the influence of the mulberry seedling field at the east side of the instruments, the sensible heat fluxes obtained by EC were increased while LAS-derived sensible heat fluxes were almost invariable under larger wind speeds conditions when the weather was dominated by prevailing easterly wind.The different phenomenon of LAS and EC might be caused by the difference of the measurement theory between two instruments.The LAS-derived sensible heat fluxes are based on M-O similarity equations; while EC measured turbulence directly and required a certain level of mechanical turbulence present in the flow.The error by using M-O similarity equations under a heterogeneous underlying surface may affect the accuracy of the LAS-derived sensible heat fluxes, and the EC measurement under low wind conditions may also be erroneous.

Figure 11 .
Figure 11.Variations of HLAS and HEC on 28 September were shown in (a,b) Footprints of sensible heat fluxes measured by the EC and LAS systems on 28 September.The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.Additionally, the corresponding time is marked in the lower right corner of picture.

Figure 11 .
Figure 11.Variations of H LAS and H EC on 28 September were shown in (a,b) Footprints of sensible heat fluxes measured by the EC and LAS systems on 28 September.The blue rectangle represents the grove, and the green trapezoid represents the mulberry seedling field.Additionally, the corresponding time is marked in the lower right corner of picture.

Table 1 .
Comparison of the LAS and EC measured factors derived for sensible heat fluxes under unstable conditions.The LAS-EC regression slope (k), determination coefficient (Rs), mean square error (RMSE), mean air density (ρ), mean friction velocity (u * ), mean temperature scale parameter (T * ) are shown in this table.
RMSE: square error of factors' regression; LAS: averaged measurement from LAS; EC: averaged measurement from EC;