Canopy Morphology Research Articles

Enhancement of ecosystem carbon uptake in a dry shrubland under moderate warming: The role of nitrogen-driven changes in plant morphology.

Net ecosystem CO2 exchange is the result of net carbon uptake by plant photosynthesis and carbon loss by soil and plant respiration. Temperature increases due to climate change can modify the equilibrium between these fluxes and trigger ecosystem‐climate feedbacks that can accelerate climate warming. As these dynamics have not been well studied in dry shrublands, we subjected a Mediterranean shrubland to a 10‐year night‐time temperature manipulation experiment that analyzed ecosystem carbon fluxes associated with dominant shrub species, together with several plant parameters related to leaf photosynthesis, leaf morphology, and canopy structure. Under moderate night‐time warming (+0.9°C minimum daily temperature, no significant reduction in soil moisture), Cistus monspeliensis formed shoots with more leaves that were relatively larger and denser canopies that supported higher plant‐level photosynthesis rates. Given that ecosystem respiration was not affected, this change in canopy morphology led to a significant enhancement in net ecosystem exchange (+47% at midday). The observed changes in shoot and canopy morphology were attributed to the improved nutritional state of the warmed plants, primarily due to changes in nitrogen cycling and higher nitrogen resorption efficiency in senescent leaves. Our results show that modifications in plant morphology triggered by moderate warming affected ecosystem CO2 fluxes, providing the first evidence for enhanced daytime carbon uptake in a dry shrubland ecosystem under experimental warming.

Simulation on the impacts of tree canopy morphology on thermal comfort at the pedestrian level

To explore the role of tree canopy morphology on the urban thermal environment and human comfort, this paper considers the variation of four factors (Top Canopy Height (TCH), Leaf Area Density (LAD), Bottom Canopy Height (BCH) and Canopy Diameter (CD)) and uses an orthogonal experiment method (OEM) to build 25 simulation scenarios. In an idealized model planted with trees, the effects of changes in canopy morphological parameters on air temperature and Physical Equivalent Temperature (PET) at a pedestrian height of 1.5m are simulated and discussed using ENVI-met software. The results show that the air temperatures for 25 different groups differed significantly during the day and less at night. At the hottest time of the day (15:00), there was a difference of 7.8°C between the coolest (ID 11) and hottest (ID 16) groups. TCH, LAD and CD are all negatively correlated with PET, while the BCH is positively correlated with PET. Among the four canopy morphology parameters, CD has the most significant effect on PET at a pedestrian level of 1.5m.

Modeling Maize Canopy Morphology in Response to Increased Plant Density.

Increased plant density markedly affects canopy morphophysiological activities and crop productivity. This study aims to model maize canopy final morphology under increased interplant competition by revising a functional–structural plant model, i.e., ADEL-Maize. A 2-year field experiment was conducted at Mengcheng, Anhui Province, China, in 2016 and 2018. A randomized complete block design of five plant densities (PDs), i.e., 4.5, 6, 7.5, 9, and 15 plants m–2, with three replications was applied using a hybrid, i.e., Zhengdan 958. Canopy morphology at different PDs was measured with destructive samplings when maize canopy was fully expanded. The relationship of changes of organ morphology in relation to increased plant density was analyzed based on 2016 data. The ADEL-Maize was first calibrated for the hybrid at 4.5 plants m–2 and then revised by introducing relationships identified from 2016 data, followed by independent validation with 2018 field data. A heatmap visualization was shown to clearly illustrate the effects of increased plant density on final morphology of laminae, sheaths, and internodes. The logarithmic + linear equations were found to fit changes for the organ size versus increased plant density for phytomers excluding ear position or linear equations for the phytomer at ear position based on 2016 field data. The revision was then further tested independently by having achieved satisfactory agreements between the simulations and observations in canopy size under different PDs with 2018 field data. In conclusion, this study has characterized the relationship between canopy morphology and increased interplant competition for use in the ADEL-Maize and realized the simulations of final size of laminae, sheaths, and internodes, as affected by increased plant density, laying a foundation to test an ideotype for maize withstanding high interplant competition.

Mobile LiDAR Scanning System Combined with Canopy Morphology Extracting Methods for Tree Crown Parameters Evaluation in Orchards.

To meet the demand for canopy morphological parameter measurements in orchards, a mobile scanning system is designed based on the 3D Simultaneous Localization and Mapping (SLAM) algorithm. The system uses a lightweight LiDAR-Inertial Measurement Unit (LiDAR-IMU) state estimator and a rotation-constrained optimization algorithm to reconstruct a point cloud map of the orchard. Then, Statistical Outlier Removal (SOR) filtering and European clustering algorithms are used to segment the orchard point cloud from which the ground information has been separated, and the k-nearest neighbour (KNN) search algorithm is used to restore the filtered point cloud. Finally, the height of the fruit trees and the volume of the canopy are obtained by the point cloud statistical method and the 3D alpha-shape algorithm. To verify the algorithm, tracked robots equipped with LIDAR and an IMU are used in a standardized orchard. Experiments show that the system in this paper can reconstruct the orchard point cloud environment with high accuracy and can obtain the point cloud information of all fruit trees in the orchard environment. The accuracy of point cloud-based segmentation of fruit trees in the orchard is 95.4%. The R2 and Root Mean Square Error (RMSE) values of crown height are 0.93682 and 0.04337, respectively, and the corresponding values of canopy volume are 0.8406 and 1.5738, respectively. In summary, this system achieves a good evaluation result of orchard crown information and has important application value in the intelligent measurement of fruit trees.

The Effects of Canopy Morphology on Flow Over a Two‐Dimensional Isolated Ridge

AbstractMomentum and mass exchanges between the atmosphere and forests situated on complex terrain continue to draw significant research attention primarily because of their significance to a plethora of applications. In this paper, we investigated flows behavior on the leeward side of a two‐dimensional forested ridge under neutrally stratified conditions using large‐eddy simulations (LESs). The goal is to understand how variations in leaf area index (LAI), vertical canopy foliage distributions, and forest edge positions affect mean/turbulent flow statistics, momentum fluxes, and onset of recirculation patterns. Although pressure perturbations are dominated by the hill shape, it is demonstrated here that changes in canopy foliage distribution modulate intensities and patterns of the leeward adverse pressure gradients. Such changes in the adverse pressure gradients alter the mean velocity streamlines including the patterns and magnitudes of the leeward downward mean vertical velocity and the velocity variances and momentum flux in the wake region. While a downwind recirculation zone develops in all cases, the details regarding the incipient location and recirculation zone size vary including positions of the separation and reattachment points. Furthermore, changes in the strength and depth of the zone occur due to canopy‐induced changes in adverse pressure gradients, advection, and canopy drag. Because the recirculation zone impacts the local mean advective terms in momentum and scalar exchanges, the simulations here indicate that canopy morphology‐induced changes in the leeward flows have significant implications to both measurements and models of biosphere‐atmosphere exchange over complex terrain.

Dynamics of canopy development of Cunninghamia lanceolata mid-age plantation in relation to foliar nitrogen and soil quality influenced by stand density

It has been generally accepted that different silvicultural practices affect the forest canopy morphology and structure. During forest establishment, many natural sites were converted to coniferous plantations in southern China. Retention of the canopy during stand conversion may be desirable to promote ecological function and meet conservation objectives. We tested the impact of planting density, foliar nitrogen and soil chemical properties on the canopy development of Chinese fir (Cunninghamia lanceolata) mid-age monoculture stands. Low density (1450 trees hm−2 with planting spacing of 2.36 × 2.36 m), intermediate-density (2460 trees hm−2 with planting spacing of 1.83 × 1.83 m) and high density (3950 trees hm−2 with planting spacing of 1.44 × 1.44 m) stands were selected in Xinkou forest plantations in Sanming City, China. Canopy characteristics such as leaf area index (LAI), mean tilt angle of the leaf (MTA) and average canopy openness index (DIFN) were measured. Measurements were taken using LAI-2200 PCA. The results illustrated that stand density was the primal factor responsible in canopy structuring while soil chemical properties seem to play a secondary role for canopy dynamics. LAI increased from 3.974 m2 m-2 to 5.072 m2 m-2 and MTA increases from 34.8° to 48.7° as the stand density increased while the DIFN decreased from 0.1542 to 0.0902 with the increasing stand density but it was no significantly different in intermediate and high-density stands. Additionally, LAI and MTA were positively correlated to foliar nitrogen while the DIFN was negatively correlated. In general, soil available nitrogen, available phosphorus and soil pH were not significant to canopy parameters. The results presented provide guiding principles about the canopy dynamics distribution in varying stand densities from LICOR measurements in mid-age Chinese fir monoculture. Furthermore, this provides a base to study canopy dynamics at mature stage forests because of more senescence activities.

Effect of Plant Morphological Traits on Throughfall, Soil Moisture, and Runoff

Plant canopy morphology plays an important role in water balance, peculiarly in semiarid environments. Through a field experiment, the impact of plant morphology of two native plant species, Artemisia sacrorum Ledeb (ASL) and Spiraea pubescens Turcz (SPT), on partitioning rainwater was revealed. The results indicated that a fragmented leaf shape and apparently high stem density of ASL reduced the throughfall and its intensity effectively but facilitated soil moisture replenishment. Although SPT has a greater canopy cover and canopy depth than ASL does, larger throughfall, a smaller throughfall threshold, and stronger throughfall intensity were observed in the SPT site. Moreover, the SPT site has a higher soil steady infiltration rate, but produced more surface runoff and caused lower soil moisture content. This study highlights the importance of plant morphological traits, peculiarly leaf morphology, in partitioning rainwater in this semiarid region. We argue that plant morphological traits should be considered when selecting plant species for revegetation and assessing water balance.

Uncertainty in Parameterizing Floodplain Forest Friction for Natural Flood Management, Using Remote Sensing

One potential Natural Flood Management (NFM) option is floodplain reforestation or manage existing riparian forests, with a view to increasing flow resistance and attenuate flood hydrographs. However, the effectiveness of floodplain forests as resistance agents, during different magnitude overbank floods, has yet to be appropriately parameterized for hydraulic models. Remote sensing offers high-resolution datasets capable of characterizing vegetation structure from a variety of platforms, but they contain uncertainty. For the first time, we demonstrate uncertainty propagation in remote sensing derivations of complex vegetation structure through roughness prediction and floodplain flow for extreme flows and different forest types (young and old Poplar plantations, young and old Pine plantations, and an unmanaged riparian forest). The lowest uncertainties resulted from terrestrial and airborne lidar, where airborne lidar is currently best at defining canopy leaf area, but more research is needed to determine wood area. Mean literature uncertainties in stem density, trunk diameter, wood, and leaf area indices (20, 10, 30, 20%, respectively) resulted in a combined Manning’s n uncertainty from 11–13% to 11–17% at 2 m to 8 m flow depths. This equates to 7–8% roughness uncertainty per 10% combined forest structure uncertainty. Individually, stem density and trunk diameter uncertainties resulted in the largest Manning’s n uncertainty at all flow depths, especially for flow though Pine plantations. For deeper flows, leaf and woody areas become much more important, especially for unmanaged riparian forests with low canopy morphology. Forest structure errors propagated to flow depth demonstrate that even small flows can change by a decimeter, while deeper flows can change by 40 cm or more. For flow depth, errors in canopy structure are deemed more severe in flows depths beyond 4–6 m. This study highlights the need for lower uncertainty in all forest structure components using remote sensing, to improve roughness parameterization and flood modeling for NFM.

Individual Tree Crown Segmentation of a Larch Plantation Using Airborne Laser Scanning Data Based on Region Growing and Canopy Morphology Features

The detection of individual trees in a larch plantation could improve the management efficiency and production prediction. This study introduced a two-stage individual tree crown (ITC) segmentation method for airborne light detection and ranging (LiDAR) point clouds, focusing on larch plantation forests with different stem densities. The two-stage segmentation method consists of the region growing and morphology segmentation, which combines advantages of the region growing characteristics and the detailed morphology structures of tree crowns. The framework comprises five steps: (1) determination of the initial dominant segments using a region growing algorithm, (2) identification of segments to be redefined based on the 2D hull convex area of each segment, (3) establishment and selection of profiles based on the tree structures, (4) determination of the number of trees using the correlation coefficient of residuals between Gaussian fitting and the tree canopy shape described in each profile, and (5) k-means segmentation to obtain the point cloud of a single tree. The accuracy was evaluated in terms of correct matching, recall, precision, and F-score in eight plots with different stem densities. Results showed that the proposed method significantly increased ITC detections compared with that of using only the region growing algorithm, where the correct matching rate increased from 73.5% to 86.1%, and the recall value increased from 0.78 to 0.89.

Location Distribution Characteristics in Leaf Lateral Asymmetry of Hybrid Indica Rice

Rice leaves show lateral asymmetry. Differences in width, thickness, chlorophyll meter readings, nitrogen content, density, and yield on each asymmetrical have been reported. However, limited information is available on the distribution characteristics of the smooth side (SS) and rough side (RS) of rice leaves. The SS of rice leaves can be divided into two types according to its location in leaf: one is the smooth side appearing on the left (SSAL) of the leaf, while the smooth side of the other type appears on the right (SSAR). Through the investigation of four rice varieties, we identified that the proportion of SSAL or SSAR of flag leaves on main stems or tillers was close to 0.50 in the rice population. On the same tiller stem, SSAL and SSAR present alternately, and the SS always stretches ahead of the RS in the curled state. Among the tiller flag leaves of a single rice plant, the proportion of tiller flag leaves with SSAL or SSAR to the total number of tillers is a discrete random variable. Frequency distribution analysis of these proportions showed that the percent of frequency of some class limits was close to each other, and the ratios of percent of frequency of SSAL to SSAR were close to 1:1. To our knowledge, this is a novel finding in rice canopy morphology. According to this study, we can divide rice population into two groups in four hybrid rice varieties.

Notice of release of seven little bluestem populations selected for canopy morphology: selected classes of germplasm

Notice of release of seven little bluestem populations selected for canopy morphology: selected classes of germplasm

Estimating Canopy Parameters Based on the Stem Position in Apple Trees Using a 2D LiDAR

Data of canopy morphology are crucial for cultivation tasks within orchards. In this study, a 2D light detection and range (LiDAR) laser scanner system was mounted on a tractor, tested on a box with known dimensions (1.81 m × 0.6 m × 0.6 m), and applied in an apple orchard to obtain the 3D structural parameters of the trees (n = 224). The analysis of a metal box which considered the height of four sides resulted in a mean absolute error (MAE) of 8.18 mm with a bias (MBE) of 2.75 mm, representing a root mean square error (RMSE) of 1.63% due to gaps in the point cloud and increased incident angle with enhanced distance between laser aperture and the object. A methodology based on a bivariate point density histogram is proposed to estimate the stem position of each tree. The cylindrical boundary was projected around the estimated stem positions to segment each individual tree. Subsequently, height, stem diameter, and volume of the segmented tree point clouds were estimated and compared with manual measurements. The estimated stem position of each tree was defined using a real time kinematic global navigation satellite system, (RTK-GNSS) resulting in an MAE and MBE of 33.7 mm and 36.5 mm, respectively. The coefficient of determination (R2) considering manual measurements and estimated data from the segmented point clouds appeared high with, respectively, R2 and RMSE of 0.87 and 5.71% for height, 0.88 and 2.23% for stem diameter, as well as 0.77 and 4.64% for canopy volume. Since a certain error for the height and volume measured manually can be assumed, the LiDAR approach provides an alternative to manual readings with the advantage of getting tree individual data of the entire orchard.

The effects of tree characteristics on rainfall interception in urban areas

Trees in urban areas have significant effects on the urban ecosystem. They can be used to improve the water cycle in urban areas by increasing evaporation and reducing runoff through rainfall interception. Street trees placed in planters on impervious areas reduce runoff by intercepting rainfall and by temporarily storing raindrops on leaves. Therefore, understanding tree canopy geometry and the effect of rainfall interception is important in urban hydrology. In this study, we assessed the effect of tree canopy morphology on rainfall interception using four major street tree species, Sophora japonica L., Ginkgo biloba L., Zelkova serrata (Thunb.) Makino, and Aesculus turbinata Blume, in Seoul, South Korea. We measured throughfall for each tree and also derived three-dimensional data of tree canopy morphology with a terrestrial laser scanner. Tree height, canopy crown width, leaf area index (LAI), leaf area density, mean leaf area, and mean leaf angle were used to determine canopy morphology. The interception rate was mostly affected by the LAI; a higher LAI tended to result in a higher interception rate. Leaf area affected the rainfall interception rate when trees had similar LAIs. These findings on individual tree canopy rainfall interception can help us to understand the importance of rainfall interception in hydrology and for ecological restoration when planning urban green spaces.

Identification, Distribution and Abundance of Scale Insect Associate with Several Clones of Durian (Durio zibethinus Merr.)

Durian contributes a high economic value in agribusiness of horticultural product. However, recently the attack of scale insect has devastated most of pre-production trees in Indonesia. A research with objective to know the basic information about the correct species, distribution and population abundance was done at Aripan Research Station of Indonesian Tropical Fruits Research Institute from July to August 2018. A proper identification and observation method were used to know the species identity, distribution and abundance of the insect. Biotic and abiotic factors in each clone such as nutrients content, tree morphological characteristic, light intensity, temperature and relative humidity under canopy were analysed. The result showed that Aulacaspis vitis Green was a scale insect identified associating with durian leaves. Clone Bintana and Sikapal were the most and the less preferred with the intensity of damage was 88.60 % and 38.30 %, respectively. Nutrient content especially nitrogen and water as well as canopy morphology were simultaneously the most contributed in the preference of a certain clone. The pest was preferred to clump in the East cardinal direction and no significant difference between upper and lower part of the tree including within and outside of the canopy.

Friction factor for turbulent open channel flow covered by vegetation

The need for operational models describing the friction factor f in streams remains undisputed given its utility across a plethora of hydrological and hydraulic applications concerned with shallow inertial flows. For small-scale roughness elements uniformly covering the wetted parameter of a wide channel, the Darcy-Weisbach f = 8(u*/Ub)2 is widely used at very high Reynolds numbers, where u* is friction velocity related to the surface kinematic stress, Ub = Q/A is bulk velocity, Q is flow rate, and A is cross-sectional area orthogonal to the flow direction. In natural streams, the presence of vegetation introduces additional complications to quantifying f, the subject of the present work. Turbulent flow through vegetation are characterized by a number of coherent vortical structures: (i) von Karman vortex streets in the lower layers of vegetated canopies, (ii) Kelvin-Helmholtz as well as attached eddies near the vegetation top, and (iii) attached eddies well above the vegetated layer. These vortical structures govern the canonical mixing lengths for momentum transfer and their influence on f is to be derived. The main novelty is that the friction factor of vegetated flow can be expressed as fv = 4Cd(Uv/Ub)2 where Uv is the spatially averaged velocity within the canopy volume, and Cd is a local drag coefficient per unit frontal area derived to include the aforemontioned layer-wise effects of vortical structures within and above the canopy along with key vegetation properties. The proposed expression is compared with a number of empirical relations derived for vegetation under emergent and submerged conditions as well as numerous data sets covering a wide range of canopy morphology, densities, and rigidity. It is envisaged that the proposed formulation be imminently employed in eco-hydraulics where the interaction between flow and vegetation is being sought.

Characterization of the Effect of Increased Plant Density on Canopy Morphology and Stalk Lodging Risk.

Plants react to the environment and to management interventions by undergoing architectural and structural modifications. A field trial was conducted in China in 2016 to study the effects of the plant population on morphological development of the maize canopy. The main objectives of the current study were (i) to characterize the effects of increased plant density on canopy morphology and stalk lodging and (ii) to explore the relationships between organ morphology and stalk lodging. The field experiment was composed of five plant densities (4.5, 6, 7.5, 9, and 15 plants m−2) of three cultivars: Zhengdan 958 (lodging-resistant cultivar), Longping 206 and Jinqiu 119 (lodging-susceptible cultivars). In response to plant densities of all the three cultivars, the lamina and sheath lengths increased in lower phytomers but decreased in upper phytomers. The lamina width and internode diameter decreased for all phytomers in response to plant densities for all the cultivars. The correlation between organ morphology, plant density and stalk lodging was linear. Data obtained from characterization used in this study (that is, canopy morphology, correlation of organ morphology with stalk lodging traits in response to various plant densities for different cultivars, etc.) will be useful in future modeling studies to predict the morphology characteristics of the canopy affected by interplant competition and stalk lodging.

Effects of leaf area index and density on ultrafine particle deposition onto forest canopies: A LES study

A framework to describe transport and deposition of ultrafine particles (UFP) within forests using large eddy simulation (LES) is presented. Comparison with measurements collected within and above a Scots pine stand in Southern Finland are used to explore the plausibility of the simulations. The numerical model is then employed to quantify the effects of canopy morphology (leaf area index and leaf area density), turbulence intensity, and particle size on the partitioning between upper canopy and subcanopy deposition and the overall deposition velocities. Results show a complex interplay between canopy morphology and turbulence, which is reflected on the particle flux profiles within the canopy. However, mean particle concentration profiles, total deposition, and deposition velocities at the canopy top are insensitive to the leaf area density profile but show dependence mostly on leaf area index, turbulence levels, and particle size. Finally, with the goal of understanding the sensitivity of the deposition velocity to all these parameters, an analytical model is developed that shows good agreement with LES results (within ±20%) for all conditions simulated here.

Partial Least Square Discriminant Analysis Based on Normalized Two-Stage Vegetation Indices for Mapping Damage from Rice Diseases Using PlanetScope Datasets.

In recent decades, rice disease co-epidemics have caused tremendous damage to crop production in both China and Southeast Asia. A variety of remote sensing based approaches have been developed and applied to map diseases distribution using coarse- to moderate-resolution imagery. However, the detection and discrimination of various disease species infecting rice were seldom assessed using high spatial resolution data. The aims of this study were (1) to develop a set of normalized two-stage vegetation indices (VIs) for characterizing the progressive development of different diseases with rice; (2) to explore the performance of combined normalized two-stage VIs in partial least square discriminant analysis (PLS-DA); and (3) to map and evaluate the damage caused by rice diseases at fine spatial scales, for the first time using bi-temporal, high spatial resolution imagery from PlanetScope datasets at a 3 m spatial resolution. Our findings suggest that the primary biophysical parameters caused by different disease (e.g., changes in leaf area, pigment contents, or canopy morphology) can be captured using combined normalized two-stage VIs. PLS-DA was able to classify rice diseases at a sub-field scale, with an overall accuracy of 75.62% and a Kappa value of 0.47. The approach was successfully applied during a typical co-epidemic outbreak of rice dwarf (Rice dwarf virus, RDV), rice blast (Magnaporthe oryzae), and glume blight (Phyllosticta glumarum) in Guangxi Province, China. Furthermore, our approach highlighted the feasibility of the method in capturing heterogeneous disease patterns at fine spatial scales over the large spatial extents.

Universal scaling-law for flow resistance over canopies with complex morphology

Flow resistance caused by vegetation is a key parameter to properly assess flood management and river restoration. However, quantifying the friction factor or any of its alternative metrics, e.g. the drag coefficient, in canopies with complex geometry has proven elusive. We explore the effect of canopy morphology on vegetated channels flow structure and resistance by treating the canopy as a porous medium characterized by an effective permeability, a property that describes the ease with which water can flow through the canopy layer. We employ a two-domain model for flow over and within the canopy, which couples the log-law in the free layer to the Darcy-Brinkman equation in the vegetated layer. We validate the model analytical solutions for the average velocity profile within and above the canopy, the volumetric discharge and the friction factor against data collected across a wide range of canopy morphologies encountered in riverine systems. Results indicate agreement between model predictions and data for both simple and complex plant morphologies. For low submergence canopies, we find a universal scaling law that relates friction factor with canopy permeability and a rescaled bulk Reynolds number. This provides a valuable tool to assess habitats sustainability associated with hydro-dynamical conditions.

EQUATIONS TO ESTIMATE TREE GAPS IN A PRECISION FOREST MANAGEMENT AREA THE AMAZON BASED ON CROWN MORPHOMETRY

ABSTRACT The precision forest management technique still has much to be improved with the incorporation of forest biometric techniques and forest profiling with airborne LIDAR. When planning the cutting of a tree in forest management, the volume to be produced for industry is estimated but not the area impacted by removal of the tree. The objective of the present study was to develop equations for the Amazon rainforest that are able to estimate the impact area of gaps from harvesting individual dominant and co-dominant trees based on the canopy morphology obtained through forest profiling. On two separate occasions profiles were made in an annual forest-production unit in the Antimary State Forest (FEA) in the state of Acre, Brazil. The first was done a few days before the start of logging in 2010 and the second was done after completion of harvest activities in 2011. With field measurements and processing of the cloud of LIDAR points, dendrometric and morphometric variables were obtained for the canopy in order to develop equations for estimating gap areas. After evaluation of the explanatory variables with the highest correlation with gap area, the method used considered all possible models and included 2-4 parameters. The explanatory variables that best represent the impact of clearings are volume of the crown (VCop) and crown-projection area (APC). Ten equations were selected, of which two were chosen for use; these had R2 aj > 75% and Syx <23%. The good fit of the equations demonstrates the potential use of LIDAR to obtain information for estimating in advance the gaps in the forest cover that will be created from harvesting trees of different sizes.