Estimation of Above Ground Biomass in Forests Using Alos Palsar Data in Kericho and Aberdare Ranges

Above Ground Biomass is one of the six pools identified in the inventory of forest resources and estimation of greenhouse gas emissions and sinks from the forestry sector. The pool varies by management practices in different agroecological or agro-climatic zones in forests. The quantification of above ground biomass (AGB) hence carbon sequestration in forests has been very difficult due to the immense costs required. This research was done to estimate AGB using ALOS PALSAR L band data (HH, HV polarisation) acquired in 2009 in relation with ground measurements data in Kericho and Aberdares ranges in Kenya. Tree data information was obtained from ground measurement of DBH and tree heights in 100 circular plots of 15 m radius, by use of random sampling technique. ALOS PALSAR image is advantageous for its active microwave sensor using L-band frequency to achieve cloud free imageries, and the ability of long wavelength cross-polarization to estimate AGB accurately for tropical forests. The variations result between Natural and plantation forest for measured and estimated biomass in Kericho HV band regression value was 0.880 and HH band was 0.520. In Aberdare ranges HV regression value of 0.708 and HH band regression value of 0.511 for measured and estimated biomass respectively. The variations can be explained by the influence of different management regimes induced human disturbances, forest stand age, density, species composition, and trees diameter distribution. However, further research is required to investigate how strong these factors affect relationship between AGB and Alos Palsar backscatters.


Introduction
Biomass refers to the material content in a forest (FAO, 2004).Its assessment allows valuation of the forest for purposes of carbon and available resource for energy.Forest Biomass in general is inclusive of the Above Ground and Below Ground living mass (Lu, 2006).These include standing trees, shrubs, roots, creepers/climbers, and the dead mass of fine and coarse litters associated with the soil.However, from most previous research (Steininger, 2000;Zheng et al., 2004;Lu et al., 2005), there is difficulty in collecting field data of Below Ground Biomass.Moreover, unlike above ground biomass, below ground biomass is not immediately affected by common deforestation and degradation processes.This therefore makes most of research to focus on biomass estimation on Above Ground Biomass (AGB).
Remote sensing approaches offer considerable potential in support of forest monitoring as it provides long-term and repetitive observations over large areas.
Satellite imagery based techniques provide an alternative to traditional methods by providing spatially explicit information and enable repeated monitoring, even in remote locations, in a cost-effective way (Patenaude et al., 2005;Rosillo-Calle et al., 2006).With the advantage of providing spatial, temporal, and spectral information (Brown, 2002), remote sensing can be used as a tool to estimate carbon.Difficulties have been experienced to acquire cloud free images in tropical regions hence a key challenge in using the widely available optical remote sensing data.
Space-borne synthetic aperture radar (SAR) sensors such as the L-band ALOS PALSAR, is an active system, transmitting microwave energy at wavelengths from 3.0 (X-band) to 23.6 cm (L-band).The major advantage of SAR systems is their weather and daylight independence.In addition, the ability to penetrate into the volume of the object (canopies) which depend on the wavelength is another important character (Lu et al., 2005).
The degree of penetration depends on the wavelength so the sensor's ability to estimate biomass.The ability to measure biomass is additionally affected by the polarization and the incidence angle of the sensor, and land cover and terrain properties (Lu et al., 2005).ALOS PALSAR is one of the most promising remote sensors to map above-ground forest carbon stock.Many studies have shown that long-wavelength radar data are sensitive to above-ground biomass (AGB).
This work aimed at estimating AGB in plantation and Natural forest through application of remote sensing technique using ALOS PALSAR L-band backscatter and forest biomass for a large range of AGB values.
The selection of the study area was based on availability of various characteristics of a typical tropical forest in Kenya, data availability and accessibility.The chosen study area for this research was Kericho that lies within longitudes 35˚02'E and 35˚40'E and between the Equator and latitude 0˚23'S in the highland west of the Great rift-valley in Kenya.The Forest Reserve occupies a total area of 32,700 ha.The area is traversed by three agro-climatic zones-humid, Semi humid and sub humid.The soils are imperfectly drained, very deep, dark grey to black vertisols (Bationo et al., 2006).
The Aberdare forest is within longitudes 36˚30'E and 36˚55'E and latitudes 0˚45'S and 0˚05'S the area is traversed by four counties-Kiambu, Murang'a, Nyeri and Nyandarua with four agroclimatic zones-humid, Semi humid, sub humid and semi humid to arid.Aberdare Ranges have inherent high fertility, being of basaltic origin.They are well drained, normally very deep, dark reddish brown, clays with a humid top-soil layer.Figure 1 shows the location of sampling plots of the study area in Kericho and Aberdare Ranges.( )

ALOS
where σ is back-scatter coefficient, DN is Digital number value and CF is Ca- libration Factor (taken as −83 in this study).
The Diameter measurement was done at Breast Height (DBH) for all trees stand inside the plots and also the height.The allometric Equation (2) was used for ground measured biomass estimation for the trees.The AGB was calculated based on mass per hectare, which gives the units of Mg/ha.Several variables were derived from ALOS PALSAR images that were correlated with the measured AGB on the ground.
( ) ( ) where B is the Biomass, D 1.3 is the Diameter at Breast height (Diameter is measured at the height of 1.3 metres from the ground level) and h is the Height.
The microwave energy transmitted that penetrates the forest canopies is largely dependent on the size and orientation of canopy structural elements, such as leaves, branches and stems.ALOS PALSAR has two polarizations, horizontal (HH) and vertical (HV) polarimetry.In this study HV responded strongly than HH in estimating AGB.However both polarizations were used in this study to investigate their correlations.By use of the strongest HV correlations coefficient factor, the empirical prediction model to estimate AGB for the study area was established.

Field Work
The purpose of the fieldwork phase was to measure the biophysical parameters such as DBH (cm), tree height (m) from the study area.To execute a forest inventory, we need to establish a relationship between directly measurable tree or stand characteristics (e.g.DBH, height) and other forest stand parameters such as volume or biomass which is impossible to be measured directly (Husch et al., 2003).
Therefore, the stand parameters of tropical forest was actually measured from the field and then used for biomass calculation by use of allometric equations.

Sampling Design
Sampling was done on two unique forest management systems.The forests lie in different agro-ecological/agro-climatic zones.Kenya is divided into 7 agro-climatic zones using a moisture index (Sombroek et al., 1982) based on annual rainfall expressed as a percentage of potential evaporation.Areas with an index greater than 50% have high potential for cropping, and are designated zones I, II, and III.These zones account for 12% of Kenya's land area.The semi-humid to arid regions (zones IV, V, VI, and VII) have indexes of less than 50% and a mean annual rainfall of less than 1100 mm.
The study areas sample plots in Kericho and Aberdares should be lie in agro-climatic zone II and III respectively.The sampling targeted all the tree species within the sampling units.The data collection was done for all the tree species for the two study areas.The sampling units were for plantation and Natural forest where an average total number of 30 species was for Natural forest and 20 tree species for plantation forest.

Sample Units
The established sampling unit was circular plot with a radius of 15 meters.A total number of fifty plots were measured in Kericho and Aberdare respectively in areas of plantation and natural forests.There were 80 plots in natural forests and 20 plots in plantation forests.The aim of the field measurements was to capture the total amount of biomass on the sample plots.All trees within 15 m radius with diameter at breast height (DBH) of at least 5 cm were measured and also the standing height as shown in Figure 2.

Sample Plots Establishment
All the plots were generated with a cell size of 25 × 25 m.The cell size corresponds to the area of a reference sample plot and the radius was chosen to be 15 m.The random forest sampling method was applied for location of the study plot in both areas.
Sample plot information for every tree species was collected.A separate copy of Sample Plot Form was filled for each sample plot.The form contained general data describing the sample plot, including the information on its identification code and whether the plot can be accessed or not.

Field Data Analysis
There have been many previous studies on the tree and stand level modelling in Figure 2. A circular sample plot with a radius (R) of 15 meters.
Kenya and other Eastern African countries.However, the focus has mainly been on plantation tree species such as Pinus patula, P. radiata, Cupressus lusitanica and Eucalyptus sp. and the main focus has been on stand growth and yield rather than individual trees.Furthermore, the models are often based on data from restricted areas and may not be relevant at the national scale.Equations to predict tree biomass has been considered only recently because of the increasing importance to estimate the forest carbon stocks.
Allometric equations were used to calculate the ABG biomass from volumetric or structural dimensions and most of them are based on tree diameter at breast height and tree height.Various allometric equations have been developed for tropical rain forest (Brown, 1997;Araújo et al., 1999;Chambers et al., 2001;Chave et al., 2005;Basuki et al., 2009).However, according to literature review, only one allometric equation to estimate AGB biomass from bio physical parameters Heights and DBH was available locally for the study area.Kinyanjui et al.
(2014) used Equation ( 2) for montane forests of Kenya.The AGB per tree was then summed over all the trees in the plot and expansion factor used to obtain a stand-level AGB estimate per plot for the two sites.

ALOS PALSAR Data
ALOS (Advanced Land Observing Satellite) is designed as a follow-on to the JERS-1 and ADEOS (Midori).Besides PRISM (for stereoscopic mapping) and an AVINIR, ALOS carries a phased array L band SAR, known as PALSAR.The SAR has a swath width of 70 km and a 2 look spatial resolution of 10 m in its observation mode, and a swath width of 250 -360 km with a spatial resolution of 100 m in a Scan SAR (wide swath width) mode (Shimada et al., 2009).
The ALOS PALSAR data used for this research was acquired in January 2009 with 2 polarisation HH and HV.It was a fine mode product with level 1.5.The preprocessing parameters used are as shown in Table 1.

Estimation of Above Ground Biomass
In this study, correlation analysis was used to assess the degree of association of biophysical parameters, AGB with radar backscatter and DN value from fusion data.Regression analysis is a common way to develop AGB estimation models (Lu, 2006).Long wavelength in SAR L-band is more reliable for AGB estimation.
Correlation analysis measures the degree of association between two or more variables.Correlation expresses the joint property or relationship between two or more variables to see how closely they are associated.When correlation exists, the size of the measurements of one variable is related to the sizes of the measurements of another variable (Husch et al., 2003).
The measure of the degree of association between two variables is called the correlation coefficient.It takes value between −1 and 1, where −1 indicate a perfect negative relation, +1 a perfect positive relation and the value of 0 indicates absence of relation (Stein et al., 1999).
The objective of regression analysis is to quantify the relationship between dependent variable and one or more independent variables.Regression implies a cause and effect relationship in which a change in the value of an independent variable will result in an expected average change in the dependent variable.The quantitative relationship is expressed by an equation and its graphic representation (Husch et al., 2003).The square value of the correlation coefficient (R 2 ) is called the coefficient of determination.It can be interpreted as indicating the percentage of variation in one variable that is associated with other variable (Husch et al., 2003).

Results
The total biomass estimation from ground data was presented by use of boxplots shown in Figure 4 for Natural and plantation forest of the two study areas respectively.The plots analysed for each study area was 40 plots for natural forest and 10 plots for plantation forest.
The estimated ground measurement biomass for Aberdare ranges forest was slightly higher than the Kericho area as represented in     In Figure 8 the results represented HV estimated biomass in Kericho, very strong with a regression of 0.880 against the HH backscattering which had a weak regression value of 0.520.
In Figure 9 the Aberdare ranges Natural forest shows a better correlation of HV in Natural forest of 0.753 while the HH backscattering coefficient was weak with a value of 0.514.
Figure 10 shows Aberdare Forest, the HV backscattering coefficient value as relatively strong with a regression value of 0.613 against the HH band with a value of 0.450.
As shown in Figure 11 the total biomass estimated correlation with the biomass measured in Aberdare was strong in HV band 0.708 unlike in the HH band which recorded a regression value of 0.511.
As represented in the box graph below in Figure 12, the above ground measured biomass in correlation with biomass estimation from L band HV backscatter is high in comparison with biomass estimation from L band HH.
The correlation analysis was employed again using Pearson correlation coeffi-     cient to see the effect of this variation in HH and HV L-band in AGB estimation.
Table 2 shows confidence level of 95% the HV band in Plantation and Natural forest that was very strong with a coefficient value of 0.880 and 0.708 respectively.
The HH band was relatively weak in both scenarios at the same confidence level with a value of 0.511 in natural forest and 0.520 in plantation.Therefore these results affirm the strength of HV L-band in estimating AGB.
The land cover classification result of the total study area of Kericho site was calculated and presented in Table 3.The highest land cover class was Forestland with an area of 17,487.Aberdare ranges land cover area had 3 classes as shown in Figure 15.The highest land cover area was forestland with an area of 14,513.8hectare and lowest was cropland.Table 4 shows the total coverage of the land cover area of the classification.
Figure 15 shows spatial representation of the above ground biomass in the Aberdare Ranges.The biomass spatial distribution (the highest to the lowest biomass value) in the Aberdares Ranges is shown in Figure 16.

Discussion
Standing biomass of a given tree species is mainly a function of DBH, tree height and density which depends on tree age, forestry practices and also environmental and genetic factors (Kasischke and Christensen, 1990).Therefore, there must be a good relationship between radar backscatter and biophysical parameters to exploit the use of SAR backscatter to determine AGB since radar energy respond  et al., 2004).
According to the results from the two study areas the HV L-band had strong correlations in all the forest and hence adopted in this research to compute biomass estimation value for the forest.More strength is alluded by the Pearson product correlation where at confidence level of 95% the HV band in Plantation and Natural forest was very strong with a coefficient value of 0.880 and 0.709 respectively as shown in Table 2.
In plantation forest, the correlation variations were due to possibly the gap in management status of sampling plots obtained from the field.This could be improved if the sampling plots could be stratified base on the different managements.However, this could not be done due to the limitation of number of plots obtained from plantation.
In another research, (Beaudoin et al., 1994) concluded that linear regression analysis between backscatter coefficients of forest stands and each forest parameter resulted in a higher correlation at HV, followed by HH with general forest parameter such as stand height, DBH, tree and stand basal area.

Conclusion
The objective of this research was to estimate and map Above Ground Biomass Figure 1.Study area-Kericho (top red circle) and Aberdare Ranges (bottom red circle).

Figure 5 .
Figure 5. Biomass estimate and biomass measured in both L-Band HH and HV.

Figure 8 .
Figure 8. Kericho Forest Biomass correlation against the Measured and Estimated.

Figure 11 .
Figure 11.Aberdare Ranges Forest Biomass correlation against the Measured and Estimated.

Figure 12 .
Figure 12.Aberdare Ranges Forest Biomass correlation against the Measured and Estimated HH and HV L-band.
Figure 14 which represents the above ground biomass distribution in the area from the highest to the lowest figure in Mg/ha.The lowest biomass was below 700 Mg/ha while the highest was 3500 Mg/ha.

(
AGB) in plantation and Natural forest using Alos Palsar data.To achieve this, correlation analysis was used to assess the relation of AGB and other stand parameter measured from the field with Palsar backscatter extracted from L-band HH and HV polarization.After that, a linear regression model was established with the chosen inputs from previous correlation analysis to estimate AGB.The AGB in the natural forest could be estimated and mapped accurately using the linear regression of L-band HH and HV.The plantation forest had a weaker correlation in HV but strong in HH backscatter.The strong correlation of HV of 0.753 developed the biomass estimation equation y = 0.0116x − 12.904.y denotes the biomass value while x is the HV value.The equation was used for biomass maps generation for the Kericho and Aberdares study areas.ALOS PALSAR images are advantageous for its active microwave sensor using L-band frequency to achieve cloud free imageries.The data has never been used in Kenya for biomass estimation.Owing to the challenges that are faced on the estimation of biomass in tropical rain forest in Kenya, the images offer a solution.

Table 4 .
Aberdare ranges land cover classification.