Accurate estimation of vegetation biophysical variables such as the vegetation canopy height (H) is of great importance to the applications of the land surface models. It is difficult to obtain the data of H at the regional scale or larger scale, but the remote sensing provides the most useful and most effective method. The leaf area index (LAI) is closely related to the H, and we analyzed its relationship with the correlation analysis based on the dataset at 86 site-years of field measurements from sites worldwide in this study. The result indicates that there is significant positive exponent correlation between these two parameters and the change of LAI would exert great impacts on H. The higher the LAI is, the higher the H is, and vice versa. Besides, the coefficients of different land cover types are very heterogeneous, and LAI of the needleleaf forest shows strong correlation with H, while that of the cropland shows weak correlation with H. The results may provide certain reference information for the extraction of the data of H at the regional scale with the remote sensing data.
1. Introduction
The vegetation plays an important role in the transfer of heat, momentum, and substance in the Earth system [1], and the accurate estimation of vegetation biophysical variables is of great importance to the agricultural, ecological, and meteorological applications [2]. The structural factors of vegetation, such as the leaf area index (LAI) and vegetation canopy height (H), have direct influence on the surface albedo, surface roughness, surface temperature, surface moisture, and so forth, all of which are important input parameters of the models such as land surface models (LSM), regional climate model (RCM), and global climate model (GCM) [3]. Therefore, the data accuracy of the vegetation density and height may have great impacts on the uncertainties in the simulation results with these models.
Both LAI and H are important vegetation physiological parameters that have close relationship with the ecological, hydrological, and climatic models. H is an important ecological metric that can provide essential information to scientists interested in understanding or modeling a wide range of atmospheric, hydrological, biophysical, and ecological processes in the forest and shrubland [4], while LAI is also one of the most important vegetation parameters and land property indices that serve as a primary controlling factor of the exchange of energy, water, and carbon fluxes between the terrestrial ecosystems and the atmosphere [5–8]. As the primary attribute of the vertical structure, H affects the boundary layer meteorology and microclimate [9]. Besides, H plays a very important role in the interface between atmosphere and land surface in many ways. For example, since the physical processes generally change at the H, the quantities of physical parameters (e.g., wind speed, temperature and concentration of CO2) above the ground surface are often normalized at the H, and sometimes these physical parameters are also normalized with their values at the H [10]. In addition, some studies demonstrated that H also plays an important role in influencing the biogeochemical properties of the land system. For example, Desai et al. showed that H had a strong correlation with CO2 fluxes [11]. H is a fundamental variable in allometric equations that estimate the forest biomass and productivity [12, 13]. Therefore, H is a very important parameter in many research fields; however, the current land cover and vegetation classifications usually provide little and even no information about H beyond the most basic structural vegetation distinctions [4].
Land surface models such as CESM, SEBAL, and METRIC [14, 15] can calculate the surface fluxes as a function of the aerodynamic resistance for the heat transfer, which is a function of aerodynamic roughness length z0 (m). In these models zero-plane displacement d (m) and z0 may be derived from a variety of canopy structure inputs. At the most basic resolution, d and z0 are simple linear functions of site-level H, typically d≈0.66H and z0≈0.10H [16, 17]. Some studies took d and z0 as a proportion of h [18]. This increases the uncertainties of the simulation results. Pitman (1994) showed that a 30% error in the estimation of z0 within LSM could lead to a 15% error in estimation of the sensible heat flux over forested environments [19]. Besides, parameters such as d and z0 [20, 21] are often normalized by H for intersite comparison studies. Both d and z0 are primarily dependent on H, though the normalized values d/H and z0/H are known to vary with the canopy density and/or LAI [20–22].
Some researchers have found that there is a relatively high correlation between the H estimates from NASA’s active laser vegetation imaging sensor (LVIS) and maximum H values of the Howland Forest measured in the field [23, 24]. However, there is some limitation in the spatial extent of the data obtained from the airborne LIDAR sensors such as LVIS [25], and it is necessary to develop some efficient methods to estimate H using optical remote sensing [26]. Unfortunately, there have been no mature approaches yet. Currently, the light detection and ranging (LIDAR) data can provide detailed information about the forest canopy height in the vertical plane [25, 27]. For example, the current remote sensing data, for example, the Landsat data, can provide useful structural information in the horizontal plane, but these data are relatively insensitive to the H. It has become a popular and effective method to extract the H based on the remote sensing data at the regional scale [6], and it has also become an important scientific issue to analyze the relation between H and other relevant parameters of vegetation which can be easily accessed by remote sensing.
In terms of the relationship between the vegetation characteristics, other researches have evaluated the relationships between vegetation indices, canopy structure, physiology, and biomass [28, 29]. Besides, some previous researches have reported that there is close relationship between LAI and the land cover type [30]. However, there are few studies about the direct physical relationship between them. In this study, we have assumed that there is certain linear relationship between them and analyzed this relationship with the regression between H and LAI on the basis of observation data collected from 86 sites covering various vegetated surfaces. The result of this study may provide certain reference value for indirect extraction of H at the regional and global scale.
2. Data
The dataset used in this study mainly includes the data of H and LAI, both of which are considered as important parameters of the land surface characteristics in most sites [31]. In this study, we have collected and sorted out data of the vegetation and structure parameters at 86 site-years of field measurements from sites worldwide. The original vegetation cover data were sorted into 17 types [32] according to vegetation/land cover classification system standard by International Geosphere-Biosphere Program (IGBP) [33]. Figure 1 shows the spatial distribution of these sites. In Figure 1, the five site-years in the yellow circle, which do not have the coordinate information, are from the second version (Version 2) database of the Global Land Cover Characteristics.
Spatial distribution of sites.
In this study, the dataset of H and LAI comes from two data sources. The first dataset is the dataset of net carbon dioxide and water fluxes of global terrestrial ecosystems, 1969–1998. The Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC) for Biogeochemistry Dynamics organized and formatted these data for long-term archive. The data file of 133 study sites includes 110 micrometeorological studies and 23 enclosure studies. The type of study (micrometeorological or enclosure) and classification of the ecosystem given by Buchmann were added to each entry. Data used in this paper is measurement site locations and several physical and biophysical characteristics (e.g., H, LAI, and stand age) and finally sorted out valid data record 53. The second dataset comes from the research results of Cho et al. (2012) [34]. The data is measured at the peak of the leaf biomass in the specified year for each site and categorized by the IGBP vegetation classification system. In this dataset, there are 33 records.
The dataset which contains information of latitude, longitude, LAI (maximum), and H (maximum) in the sites was finally established. Descriptive statistics for the variables of this dataset are included in Table 1. For the analysis, we reclassified the original vegetation types into the types of the IGBP classification system. Then statistics of different types of record number were incorporated into NF (needleleaf forest), BF (broadleaf forest), C (croplands), G (grasslands), MF (shrublands), and savannas (Table 2). The main surface vegetation type of the site is grasslands, croplands, forests (deciduous broadleaf Forests, deciduous needleleaf forests, evergreen needleleaf forests, evergreen broadleaf forests, and open/closed shrublands), and savannas. Then the land cover classification system of vegetated surface was got after excluding nonvegetation surface types, including urban and built-up (13), snow and ice (15), water bodies (17), land cover classification types with nonsite distribution mixed forest (5), woody savannas (8), permanent wetlands (11), cropland/natural vegetation mosaic (14), and barren or sparsely vegetated (16).
Descriptive statistics of the dataset at 86 site-years of field measurements from sites worldwide.
Variable
Obs.
Mean
Std. dev.
Min
Max
H
86
12.71
11.49
0.05
60
LAI
86
3.48
2.36
0.2
10
Latitude
81
33.12
32.08
−42.2
69.47
Longitude
81
−9.18
97.61
−164.73
176
Sites statistical distribution in IGBP land cover classification system.
Code
IGBP classification
Obs.
Merge sorting
1
Evergreen needleleaf forest
16
NF
2
Evergreen broadleaf forest
2
BF
3
Deciduous needleleaf forest
19
NF
4
Deciduous broadleaf forest
21
BF
5
Mixed forest
—
—
6
Closed shrublands
1
MF
7
Open shrublands
1
MF
8
Woody savannas
—
—
9
Savannas
4
S
10
Grasslands
14
G
11
Permanent wetlands
—
—
12
Croplands
8
C
14
Cropland/natural vegetation mosaic
—
—
16
Barren or sparsely vegetated
—
—
Note. The table excluded nonvegetation surface types which are urban and built-up (13), snow and ice (15), and water bodies (17). So “—” means that such land cover types do not exist.
3. Results
The correlation analysis between H, ln(H), LAI, and ln(LAI) was implemented one by one, and finally the correlation matrix was obtained. The result indicated that the Pearson correlation value of LAI and vegetation height (H) and ln(H) reached 0.64 and 0.616, respectively, while there is a more sensitive correlation between ln(H) and LAI and ln(LAI) since their Person value was 0.609 and 0.681, respectively.
The result shows that there is a positive correlation between all of H, ln(H), LAI, and ln(LAI) (Table 3). There is the most significant correlation between ln(H) and ln(LAI), which is significant at the 1% level of significance (bilateral), and their correlation coefficient is 0.681. It is followed by that between H and LAI, the correlation coefficient of which is 0.641. The significance level of ln(H) and LAI is the weakest, and their correlation coefficient is 0.609. By inquiring the critical value table of the correlation coefficient, it is found that the correlations between H, ln(H), LAI, and ln(LAI) are all significant at the 99% level of confidence.
Correlation matrix of H, ln(H), LAI, and ln(LAI).
H
ln(H)
Pearson value
Significance (bilateral)
Pearson value
Significance (bilateral)
LAI
0.641
0.000
0.609
0.000
ln(LAI)
0.616
0.000
0.681
0.000
The regression analysis of H and LAI was implemented with the SPSS software, with H being the dependent variable and LAI being the independent variable. The results are as follows: the summarized conditions of the regression model, the estimated values of parameters, the results of curve fitting of H and LAI, and the final regression equations (Table 4, Figure 2).
Summarization of regression models and estimated values of parameters.
Models
Summarization of the model
Estimated values of parameters
Regression models
R2
F
df1
df2
Sig.
Constant
b1
Liner
0.41
58.64
1
84
0.00
1.846
3.120
y=3.120x+1.846
Logarithm
0.38
51.35
1
84
0.00
5.378
7.792
y=7.792ln(x)+5.378
Power
0.46
72.60
1
84
0.00
1.492
1.340
y=1.492x1.340
Exponential
0.37
49.51
1
84
0.00
1.057
0.461
y=1.057e0.461x
Result of fitting curve of H and LAI BF (broadleaf forest), NF (needleleaf forest), C (croplands), G (grasslands), S (savannas), and (MF) shrublands.
The result indicated that these functions all reflected the relationship between H and LAI to some degree. The result suggests that the result of curve fitting of H and LAI with the power function is the most proper, which can effectively reflect the correlation between H and LAI (y=1.4915x1.3396, R2=0.464). It is followed by the linear fitting (R2=0.411), while the result with the exponential function is the worst (R2=0.371).
Although the power function can reflect the correlation between H and LAI best, it is found that the logarithmic function can better reflect the correlation between H and LAI to some degree when the single or several land cover types are taken as the research object. So the logarithmic function has been used to analyze the correlation between H and LAI with the combined land cover types as the research object in this study, including the land cover BF, NF, C, and G.
The results of the analysis with the logarithmic function are as follows. Figures 3(a), 3(b), 3(c), and 3(d) represents the scatter plots between ln(H) and LAI of the broadleaf forest, needleleaf forest, croplands, and grasslands, respectively. The solid lines represent the fitted line of the changing trend, the fitted equations, their fitting variances (R2), and significance level; see Table 5.
Changing trends of broadleaf forest, needleleaf forest, croplands, and grasslands (with the square of correlation coefficient (R2), significance level (Sig.), and regression models).
Type
Regression model
R2
Sig.
Broadleaf forest
y=0.2348x+1.9174
0.24
0.019**
Needleleaf forest
y=0.1305x+2.0052
0.25
0.002***
Croplands
y=-0.3083x+1.192
0.49
0.064*
Grasslands
y=0.4479x-2.0488
0.26
0.055*
Notes. ***Significant at 1% level, **significant at 5% level, and *significant at 10% level.
Scatter plots between ln(H) and LAI, (a) broadleaf forest, (b) needleleaf forest, (c) croplands, and (d) grasslands.
The result indicated that the changing trends of H and LAI of broadleaf forest, needleleaf forest, and grasslands tend to converge; that is, they show a positive correlation on the whole, suggesting that the canopy heights of broadleaf forest, needleleaf forest, and grasslands all will increase with LAI. And the R2 of broadleaf forest, needleleaf forest, and grasslands are 0.24, 0.25, and 0.26. They are significant at the 10% significance level. There is a negative correlation between the changing trend of H and LAI of cropland, which is contrary to that of their land cover types and inconsistent with the common sense; it may be due to the following reasons. First, the data of croplands is greatly influenced by the types of the crop. Second, the analyzed result is not significant correlation between H and LAI since there is very limited data of this land cover type in the dataset used in this study.
By contrast, the results well reflected the correlation between H and LAI in the other three land cover types. The result indicated that there is a positive correlation between H and LAI on the whole, and the change of LAI has some influence on the change of the vegetation height; that is, the larger the LAI is, the higher the vegetation height will be, and vice versa. Besides, there is some variation in the degree of correlation between H and LAI in different land cover types.
4. Conclusions
The results show that there is significant positive correlation between H and LAI, and the power function relationship between H and LAI is especially significant. Besides, this study also analyzed the relationship between H and LAI with the vegetation type of broadleaf forest, needleleaf forest, croplands, and grasslands. The result indicates that there is significant difference in relationship between H and LAI among different land cover types. Specifically, there is a significantly positive relationship between H and LAI of broadleaf forest, needleleaf forest, and grasslands, while there is negative correlation between H and LAI of croplands. In addition, the instrument used to measure LAI in different sites is not exactly the same, which increases the uncertainties in the relationship between LAI and H.
More types of vegetation dataset should be included to improve the robustness of this study. At the same time, for the data shortage, this study only analyzed the relationship between H and LAI of broadleaf forest, needleleaf forest, grasslands, and some other vegetation cover type, so the future studies should analyze a comprehensive relationship between H and LAI and include a wider range of vegetation cover types to quantify the role of vegetation on heat transport.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was supported by the Project of the National Natural Science Foundation of China (Grant no. 41071343) and the National Key Programme for Developing Basic Science of China (Grant no. 2010CB950900).
NemaniR. R.RunningS. W.PielkeR. A.ChaseT. N.Global vegetation cover changes from coarse resolution satellite data19961013715771622-s2.0-0029751617DarvishzadehR.SkidmoreA.SchlerfM.AtzbergerC.Inversion of a radiative transfer model for estimating vegetation LAI and chlorophyll in a heterogeneous grassland20081125259226042-s2.0-4144910806710.1016/j.rse.2007.12.003TakahashiH. G.YoshikaneT.HaraM.TakataK.YasunariT.High-resolution modelling of the potential impact of land surface conditions on regional climate over Indochina associated with the diurnal precipitation cycle20103013200420202-s2.0-7794940394610.1002/joc.2119SelkowitzD. J.GreenG.PetersonB.WylieB.A multi-sensor lidar, multi-spectral and multi-angular approach for mapping canopy height in boreal forest regions20121214584712-s2.0-8485872993910.1016/j.rse.2012.02.020ChenJ. M.BlackT. A.Measuring leaf area index of plant canopies with branch architecture1991571–31122-s2.0-0026380713WuB.ZengY.HuangJ.Overview of LAI/FPAR retrieval from remotely sensed data2004104585590BonanG. B.Importance of leaf area index and forest type when estimating photosynthesis in boreal forests19934333033142-s2.0-0027335302JensenJ. L. R.HumesK. S.VierlingL. A.HudakA. T.Discrete return lidar-based prediction of leaf area index in two conifer forests200811210394739572-s2.0-5084911306110.1016/j.rse.2008.07.001RaupachM. R.Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index1994711-22112162-s2.0-0028173787KaimalJ. C.FinniganJ. J.1994New York, NY, USAOxford University PressDesaiA. R.NoormetsA.BolstadP. V.ChenJ.CookB. D.DavisK. J.EuskirchenE. S.GoughC.MartinJ. G.RicciutoD. M.SchmidH. P.TangJ.WangW.Influence of vegetation and seasonal forcing on carbon dioxide fluxes across the Upper Midwest, USA: implications for regional scaling200814822883082-s2.0-3894910117810.1016/j.agrformet.2007.08.001PatenaudeG.MilneR.DawsonT. P.Synthesis of remote sensing approaches for forest carbon estimation: reporting to the Kyoto Protocol2005821611782-s2.0-1814442316210.1016/j.envsci.2004.12.010AndersenH.ReutebuchS. E.McGaugheyR. J.A rigorous assessment of tree height measurements obtained using airborne lidar and conventional field methods20063253553662-s2.0-33846784000BastiaanssenW. G. M.MenentiM.FeddesR. A.HoltslagA. A. M.A remote sensing surface energy balance algorithm for land (SEBAL): 1. Formulation1998212-2131982122-s2.0-003241898710.1016/S0022-1694(98)00253-4AllenR. G.TasumiM.TrezzaR.Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—model200713343803942-s2.0-3454712763710.1061/(ASCE)0733-9437(2007)133:4(380)CowanI. R.Mass, heat and momentum exchange between stands of plants and their atmospheric environment19689440252354410.1002/qj.49709440208TannerC. B.PeltonW. L.Potential evapotranspiration estimates by the approximate energy balance method of Penman196065103391341310.1029/JZ065i010p03391RannikÜ.MarkkanenT.RaittilaJ.HariP.VesalaT.Turbulence statistics inside and over forest: influence on footprint prediction200310921631892-s2.0-014194022010.1023/A:1025404923169PitmanA. J.Assessing the sensitivity of a land-surface scheme to the parameter values using a single column model1994712185618692-s2.0-0028590518VerhoefA.McNaughtonK. G.JacobsA. F. G.A parameterization of momentum roughness length and displacement height for a wide range of canopy densities19971181912-s2.0-0031394322NakaiT.SumidaA.DaikokuK.MatsumotoK.van der MolenM. K.KodamaY.KononovA. V.MaximovT. C.DolmanA. J.YabukiH.HaraT.OhtaT.Parameterisation of aerodynamic roughness over boreal, cool- and warm-temperate forests200814812191619252-s2.0-5024914690910.1016/j.agrformet.2008.03.009HollingerD. Y.OllingerwS. V.RichardsonwA. D.MeyerszT. P.DailD. B.MartinwM. E.ScottN. A.ArkebauerkT. J.BaldocchiD. D.ClarkK. L.CurtisP. S.DavisK. J.DesaiA. R.DragonikkD.GouldenM. L.GuL.KatulzzzG. G.PallardyS. G.PawuK. T.SchmidH. P.StoyP. C.SuykerA. E.VermaS. B.Albedo estimates for land surface models and support for a new paradigm based on foliage nitrogen concentration20101626967102-s2.0-7814926477110.1111/j.1365-2486.2009.02028.xAndersonJ.MartinM. E.SmithM.DubayahR. O.HoftonM. A.HydeP.PetersonB. E.BlairJ. B.KnoxR. G.The use of waveform lidar to measure northern temperate mixed conifer and deciduous forest structure in New Hampshire200610532482612-s2.0-3375071936810.1016/j.rse.2006.07.001KimesD. S.RansonK. J.SunG.BlairJ. B.Predicting lidar measured forest vertical structure from multi-angle spectral data200610045035112-s2.0-3234444053510.1016/j.rse.2005.11.004LefskyM. A.CohenW. B.ParkerG. G.HardingD. J.Lidar remote sensing for ecosystem studies: lidar, an emerging remote sensing technology that directly measures the three-dimensional distribution of plant canopies, can accurately estimate vegetation structural attributes and should be of particular interest to forest, landscape, and global ecologists2002521193010.1641/0006-3568(2002)052[0019:LRSFES]2.0.CO;2ChoppingM.MoisenG. G.SuL.LaliberteA.RangoA.MartonchikJ. V.PetersD. P. C.Large area mapping of southwestern forest crown cover, canopy height, and biomass using the NASA Multiangle Imaging Spectro-Radiometer20081125205120632-s2.0-4124909981910.1016/j.rse.2007.07.024HayashiM.SaigusaN.OgumaH.YamagataY.Forest canopy height estimation using ICESat/GLAS data and error factor analysis in Hokkaido, Japan201381121810.1016/j.isprsjprs.2013.04.004GamonJ. A.FieldC. B.GouldenM. L.GriffinK. L.HartleyA. E.JoelG.PenuelasJ.ValentiniR.Relationships between NDVI, canopy structure, and photosynthesis in three Californian vegetation types19955128412-s2.0-0029132359GaoS.NiuZ.HuangN.HouX.Estimating the leaf area index, height and biomass of maize using HJ-1 and RADARSAT-22013241810.1016/j.jag.2013.02.002AsrarG.MyneniR. B.ChoudhuryB. J.Spatial heterogeneity in vegetation canopies and remote sensing of absorbed photosynthetically active radiation: a modeling study1992412-3851032-s2.0-0027010136BaldocchiD.FalgeE.GuL.OlsonR.HollingerD.RunningS.AnthoniP.BernhoferC.DavisK.EvansR.FuentesJ.GoldsteinA.KatulG.LawB.LeeX.MalhiY.MeyersT.MungerW.OechelW.PawU. K. T.PilegaardK.SchmidH. P.ValentiniR.VermaS.VesalaT.WilsonK.WofsyS.FLUXNET: a new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide, water vapor, and energy flux densities20018211241524342-s2.0-0001633245TurnerB. L.IISkoleD.SandersonS.FischerG.FrescoL. O.LeemansR.Land-use and land-cover change: science/research plan1995352-s2.0-0029529473CohenW. B.MaierspergerT. K.YangZ.GowerS. T.TurnerD. P.RittsW. D.BerterretcheM.RunningS. W.Comparisons of land cover and LAI estimates derived from ETM+ and MODIS for four sites in North America: a quality assessment of 2000/2001 provisional MODIS products20038832332552-s2.0-084230634810.1016/j.rse.2003.06.006ChoJ.MiyazakiS.YehP. J. F.KimW.KanaeS.OkiT.Testing the hypothesis on the relationship between aerodynamic roughness length and albedo using vegetation structure parameters20125624114182-s2.0-8485706254910.1007/s00484-011-0445-2