Received for publication: 11 April, 2012. Accepted for publication: 30 October, 2012. 1 Consultant. Bristol (UK). 2 Oil Palm Biology and Breeding Research Program, Colombian Oil Palm Research Center, Cenipalma. Bogota (Colombia). 3 Department of Biology, Universidad Nacional de Colombia. Bogota (Colombia). hmromeroa@unal.edu.co Agronomía Colombiana 30(3), 359-369, 2012 The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon offsets Balance de los gases de efecto invernadero de la industria de la palma de aceite en Colombia: análisis preliminar. I. Secuestro de carbono y créditos de carbono Ian E. Henson1, Rodrigo Ruiz R.2, and Hernán Mauricio Romero2, 3 ABSTRACT RESUMEN Colombia is currently the world’s fifth largest producer of palm oil and the largest producer in South and Central America. It has substantial areas of land that could be used for additional oil palm production and there is considerable scope for increasing yields of existing planted areas. Much of the vegetation on land suitable for conversion to oil palm has a low biomass, and so establishing oil palm plantations on such land should lead to an increase in carbon stock, thereby counteracting greenhouse gas (GHG) emissions responsible for global warming. The first part of this study examines changes in carbon stock in Colombia resulting from expansion of oil palm cultivation together with factors (offsets) that act to minimize carbon emissions. The results are subsequently used to construct a net GHG balance. Colombia es el quinto productor de aceite de palma en el mun- do y el más grande en América Central y del Sur. Tiene áreas sustanciales de tierra que podrían ser usadas para producción adicional de palma de aceite y existe una brecha interesante para incrementar los rendimientos en las áreas existentes. La mayoría de la vegetación de las áreas adecuadas para la con- versión a siembra de palma tienen una biomasa baja, de tal manera que al establecer plantaciones de palma de aceite en esas tierras debería incrementarse el carbono almacenado, de tal manera que haya una reducción en las emisiones de gases de efecto invernadero (GHG) responsables del calentamiento global. La primera parte de este estudio examina los cambios en el carbono almacenado como resultado de la expansión del cultivo de la palma de aceite junto con aquellos factores (crédi- tos) que actuan para minimizar las emisiones. Los resultados son subsecuentemente usados para construir el balance neto de GHG. Key words: carbon footprint, CO2 balance, carbon stock, land use change. Palabras clave: huella de carbono, balance de CO2, almace- namiento de carbono, cambio de uso de tierra. Introduction Colombia is currently the world’s fifth largest producer of palm oil, the second largest producer outside South- east Asia, and the largest producer in South and Central America. In a previous paper (Henson et al., 2011), we described the growth of the industry in Colombia from the early 1960’s, when detailed records began, up until 2009, when the total planted area in the country was just under 361,000 ha. This is small for a country with a total land area of over 103 million ha, even though much of this lies at altitudes inhospitable to growth of tropical lowland species such as oil palm. Studies on the area of land in Colombia suitable for oil palm (Romero et al., 1999; Mejia, 2000; Rubiano et al., 2009) have indicated a considerable potential for expanding the oil palm planted area, based on climatic, edaphic and topographical considerations. However, there are several additional factors that need to be addressed that may serve to limit such potential growth, such as impacts on food production, biodiversity, carbon stock and greenhouse gas (GHG) emission. The net impact of a process on climate change is determined by the balance between carbon uptake from the atmo- sphere, mainly by photosynthesis, and the emission to the atmosphere of GHGs that contribute to global warming. Subtracting sequestered carbon from total GHG emissions expressed in carbon or CO2 equivalents gives the net GHG 360 Agron. Colomb. 30(3) 2012 emission. The main GHGs associated with the cultivation of oil palm and the processing of its products are CO2, methane (CH4) and nitrous oxide (N2O). Assuming a resi- dence time of 100 years, the latter two gases have global warming potentials (GWPs) 25 and 298 times greater per unit mass than CO2. In this paper, we examine the impact of the oil palm industry in Colombia on the sequestration of carbon in oil palm biomass, palm products and by-products, and evaluate factors that can act to offset or minimize GHG emissions. An accompanying paper will examine emissions and produce a net GHG budget for the industry throughout its development. Materials and methods Information on different aspects of the oil palm industry from 1959 to 2009 was gathered from different informa- tion sources including the Colombian Federation of Palm Growers (Fedepalma), the Instituto Geográfico Agustín Codazzi, and the Colombian Ministry of Environment. The raw data were used to estimate the total GHG balance from the cultivation of oil palm and the milling of fresh fruit bunches (FFB) in Colombia. The construction of the carbon budget for Colombia closely followed procedures for Malaysia used in an earlier study (Henson, 2008, 2009). We constructed separate budgets for the four regions of Colombia where oil palm is grown (Eastern, Northern, Central and Western), using data for oil palm planted areas and production gathered in a previous paper (Henson et al., 2011). The main components of a carbon budget are: i) carbon sequestration in the oil palm plantation. This includes C in oil palm biomass, ground cover vegeta- tion and litter, including frond piles, detached petiole bases and shed male inflorescences. ii) carbon sequestration in mill products (palm oil, palm kernel oil and palm kernel cake) and by-products (empty fruit bunches (EFB), palm oil mill effluent (POME), fruit fibre and shell). iii) carbon emission due to land use conversion to oil palm (loss of carbon during clearance of vegetation). iv) carbon emission due to direct use of fossil fuels (diesel consumed in plantation vehicles and machinery). v) carbon emission due to indirect use of fossil fuels con- sumed in production and transport of inputs such as fertilizers and crop protection chemicals. vi) carbon equivalent emissions in the form of N2O (mainly produced from inorganic N sources). vii) carbon equivalent emissions in the form of CH4 pro- duced by POME ponds, mill boilers and plantation livestock. viii) carbon offsets that serve to reduce emissions such as the use of fibre and shell as fuel for mill boilers, the use of other mill by-products and biomass residues as nutrient sources (thus reducing inorganic fertilizer requirements), and the long-term storage of carbon in any timber harvested during land clearing. Note that other offsets may become important in the near future, notably conversion of CH4 from POME to CO2 by flaring, the generation of electricity from CH4, the use of palm bio-diesel as a fossil fuel substitute, and the production of bio-char from felled trunks, reduc- ing the rate of on-site carbon loss following oil palm replanting. For this study, carbon sequestration in plantations was calculated based on a) annual values of standing biomass derived using representative curves for oil palm, ground cover and litter (Fig. 1), b) biomass-age distributions such as those illustrated for oil palm in Fig. 2, and c) conversion of biomass to carbon assuming a carbon content of 45% (Henson, 2009). There are few, if any, direct measurements of oil palm standing biomass in Colombia but based on mean FFB yields and taking into account the generally conservative nature of oil palm vegetative growth, trends found in Malaysia were assumed to apply and calculations were performed using two contrasting data sets (Fig. 1), firstly, the Malaysian ‘national average’ (NA) oil palm based on long-term country-averaged yield performance (Henson, 2003), and secondly, a data set representing a higher yield- ing crop termed OP2, generated using the OPRODSIM model (Henson, 2005). Both sets of data include below-as well as above-ground biomass. To account for all carbon sequestered in the plantations, carbon in the ground cover and litter biomass was added to that in the oil palms. The amounts of carbon in these components were calculated in the same way as that for oil palm, based on biomass changes with age of planting shown in Fig. 1. 361Henson, Ruiz R., and Romero: The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon... The production data on which estimates were based were obtained from Fedepalma (2010 and earlier editions). It should be noted that amounts of PKO and PKC given by Fedepalma add up to less than the amounts of PK from which they are derived (ca. 89%), representing losses due to incomplete extraction and loss of residual moisture. Amounts of by-products were derived as fractions of bunch dry matter remaining after extraction of CPO and PK. The proportions of each by-product were based on median val- ues given by Wood and Corley (1993). Remaining residual amounts were assigned as unidentified losses. Carbon sequestration in mill products and by-products was calculated using carbon contents, estimated lifetimes and, in the case of fibre, shell and EFB, the proportions of materials burnt, given by Henson (2004, 2009). It was assumed that 22% of EFB was burnt to produce bunch ash and 46.8% of both fibre and shell were used in the mill as boiler fuel. The mean quantities present during any year were calculated assuming linear rates of utilization or decay. Burnt material was assumed to have a lifetime one-tenth that of un-burnt. No allowance was made for possible field application of by-products that might affect their decomposition rates and the amounts of carbon sequestered. Results and discussion Carbon sequestration in oil palm plantations Changes in biomass and carbon stock during the oil palm life cycle During development, oil palms accumulate carbon both in living biomass and in the litter they produce, which mainly comprises cut fronds, shed frond bases and male inflorescences. Carbon is also stored in the ground cover vegetation and in soil organic matter (SOM). Little is known of changes in SOM concentration during the establishment and growth of oil palm. There is evidence for both consi- derable soil carbon losses at clearance and for increases in SOM as plantations age (Henson, 2004). Because of such uncertainty, this aspect is omitted in the present analysis. The palms themselves generally contain about 85% of the total plantation non-SOM carbon store (Henson, 2009). Standing biomass of the NA oil palm (Fig. 1A) declines in later years due to stand losses resulting from age and disease (mainly basal stem rot in Malaysia and bud rot in Colombia), augmented by loss of frond bases. The stand density of OP2 oil palm (Fig. 1B) is maintained until clearance, which in some cases may be up to 30 years after planting the palms. Changes in mean age and biomass of oil palm plantations From annual data on planted areas combined with chan- ges in biomass with palm age, area-age and biomass-age distributions were derived annually for the different oil palm regions and for Colombia as a whole. Examples for the whole country are shown in Fig. 2. The distributions are affected by both expansion into newly planted palm areas and by the areas of old palms that are replanted. Two methods were devised to estimate replanted areas (Henson, 2003). These depend on the age at replanting termed the replanting cycle time (RCT), for which there are generally no records, on either a regional or a national basis. RCT can be quite variable but usually ranges from 20 to 30 years. Because of this, the calculations were performed for alternative periods of 20, 25 and 30 years. The mean palm ages resulting from these assumptions are shown in Tab. 1. In addition to RCT, the oil palm growth curves used (NA and OP2) and the method of estimating annually re- planted and newly planted areas (Henson, 2003) affect the FIgURE 1. Oil palm, ground cover, and litter standing biomass curves used to calculate oil palm plantation biomass in Colombia. (A) ‘National ave- rage’ oil palm biomass; (B) OprOdSIm Op2 modelled biomass. Sources: Henson (2003, 2005, 2009). 0 2 4 6 8 10 12 14 16 0 20 40 60 80 100 120 0 5 10 15 20 25 30 Li tt er a nd g ro un d co ve r st an di ng b io m as s (t h a- 1 ) O il pa lm s ta nd in g bi om as s (t h a- 1 ) Years after planting Oil palm Litter Ground cover 0 2 4 6 8 10 12 14 16 0 20 40 60 80 100 120 0 5 10 15 20 25 30 Li tt er a nd g ro un d co ve r st an di ng b io m as s (t h a- 1 ) O il pa lm s ta nd in g bi om as s (t h a- 1 ) Years after planting Oil palm Litter Ground cover A B 362 Agron. Colomb. 30(3) 2012 FIgURE 2. Changes in oil palm biomass-age and area-age distributions at 10-year intervals for the NA oil palm, assuming a 30-year replanting cycle using themeans of two methods to determine replanted areas. 1989 0 2 4 6 8 10 12 5 10 15 20 25 30 P er ce nt o f t ot al Biomass 14 Area 1999 0 2 4 6 8 10 12 P er ce nt o f t ot al 5 10 15 20 25 30 0 2 4 6 8 10 5 10 15 20 25 30 P er ce nt o f t ot al Palm age (years after planting) 2009 TABlE 1. Effect of replanting cycle time (rCT) on mean oil palm age, 1959 to 2009, in different regions of Colombia1. RCT (years) Eastern Northern Central Western Whole country Mean palm age (years after planting) 20 6.99 6.92 7.63 7.07 7.20 25 7.65 7.66 8.71 7.75 7.98 30 7.92 8.09 9.53 8.20 8.44 Mean 7.52 7.56 8.62 7.67 7.87 1 Data are means of two methods used to determine replanted areas (Fig.3C). were smaller compared with differences between 25 and 30-year cycles (Fig. 3B). The assessment method options for the replanted area (Fig. 3C) also had an effect on the oil palm standing biomass, with the two methods giving divergent estimates in later years. Differences in total standing biomass between regions (Fig. 3D) reflected both differences in planted areas and in mean oil palm biomass density (t ha-1). The latter changed substantially over time (Fig. 4). In the early years, biomass density was low since most palms were young. After reach- ing a peak in the mid-seventies, density in most regions declined before increasing again to a peak prior to a second decline. The Central region showed the most stability. Declines indicate periods of increased replanting while increases in density resulted from reduced replanting. In the last ten years or so, biomass density fell in all regions, consistent with the decrease in the mean palm age noted previously (Henson et al., 2011). Changes in carbon sequestration in oil palm plantations Mean sequestration rates given for different RCTs and regions in Tab. 3 reveal large regional differences, with the highest rates in the Eastern and Northern regions. estimation of biomass density (t ha-1) and the total stand- ing biomass in each region, needed to calculate carbon sequestration. The effects of such factors on total standing biomass are shown in Fig. 3. The biomass curve (Fig. 1A) had a substantial effect on total standing biomass in later years, with the OP2 oil palm hav- ing greater total biomass (Fig. 3A), and hence greater carbon sequestration (Fig. 5A) than the NA palm. Both mean palm age (Tab. 1) and palm biomass density (Tab. 2) increased in all regions as RCT was extended, with the increase in biomass being somewhat greater with the OP2 than with the NA oil palm. Differences between 25 and 30-year RCTs TABlE 2. Effect of oil palm growth curve and replanting cycle time (rCT) on mean oil palm biomass density, 1959 to 2009, in different regions of Colombia1. Biomass curve RCT (years) Eastern Northern Central Western Whole country Biomass (t ha-1) National average (NA) 20 36.93 37.75 40.33 39.99 38.65 25 40.06 41.27 45.78 43.33 42.58 30 40.65 42.18 47.41 44.55 43.56 Mean 39.21 40.40 44.51 42.62 41.60 OpRODsIM (Op2) 20 38.65 36.88 41.98 40.84 39.53 25 42.48 40.63 48.67 44.72 44.10 30 43.75 42.34 52.51 46.91 46.17 Mean 41.62 39.95 47.72 44.16 43.27 1 Data are means of two methods used to determine new plus replanted areas. 363Henson, Ruiz R., and Romero: The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon... 0 4,000 8,000 12,000 To ta l o il pa lm b io m as s (1 03 t) OP2 16,000 NA 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 0 4,000 8,000 12,000 O il pa lm b io m as s (1 03 t) 20 years 16,000 25 years 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 30 years 0 4,000 8,000 12,000 O il pa lm b io m as s (1 03 t) Method 1 16,000 Method 2 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years Mean 0 1,000 2,000 3,000 O il pa lm b io m as s (1 03 t) East 5,000 North 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years Central 4,000 West FIgURE 3. Effect of (A) oil palm growth curve, (B) replanting cycle time (rCT), (C) method of estimating replanted areas, and (d) oil palm region, on estimates of oil palm standing biomass in Colombia from 1959 to 2009. data in (A) are for all regions combined and are means of three rCTs and two methods of estimating replanted areas. data in (B) are for the NA oil palm, for all regions combined, and are means of two methods of estimating replanted areas. data in (C) are for the NA oil palm, for all regions combined, and are means of three rCTs. data in (d) are for the NA oil palm and are means of three rCTs and two methods of estimating replanted areas. FIgURE 4. regional differences in mean oil palm biomass density in Colombia from 1959 to 2009. data are for the NA oil palm and are means of 20, 25 and 30-year replanting cycle times and two methods of estimating replanted areas. 0 20 40 60 O il pa lm b io m as s (t h a- 1 ) East North 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years Central 80 West TABlE 3. mean rates of carbon sequestration in oil palm standing bio- mass in different regions of Colombia, 1959 to 2009, calculated for three replanting cycle times (rCTs)1. RCT (yr) Eastern Northern Central Western Whole country Rate of carbon sequestration (103 t year-1) 20 30.68 34.12 27.17 14.06 106.03 25 48.52 42.11 33.76 15.82 140.21 30 51.57 45.18 34.36 17.35 148.46 Mean 43.59 40.47 31.76 15.74 131.57 1 Based on means of two methods used to estimate replanted areas. TABlE 4. mean rates of carbon sequestration in plantation litter and ground cover biomass in different regions of Colombia, 1959 to 2009, calculated for three replanting cycle times (rCTs)1. RCT (yr) Eastern Northern Central Western Whole country Rate of carbon sequestration (103 t year-1) 20 10.40 9.67 7.79 2.99 30.84 25 11.89 10.53 8.52 3.23 34.17 30 12.32 11.07 8.73 3.50 35.67 Mean 11.54 10.42 8.34 3.24 33.54 1 Based on means of two methods used to estimate replanted areas. A B C D 364 Agron. Colomb. 30(3) 2012 FIgURE 5. Effect of (A) oil palm growth curve, (B) replanting cycle length, (C) method of estimating replanted areas, and (d) oil palm region, on esti- mates of annual carbon sequestration in oil palm standing biomass in Colombia from 1959 to 2009. data in (A) are for all regions combined and are means of three replanting cycle times and two methods of estimating replanted areas. data in (B) are for the NA oil palm, for all regions combined, and are means of two methods of estimating replanted areas. data in (C) are for the NA oil palm, for all regions combined, and are means of three re- planting cycle times. data in (d) are for the NA oil palm and are means of three planting cycle times and two methods of estimating replanted areas. 0 100 200 300 C s eq ue st er ed in o il pa lm b io m as s ( 10 3 t y ea r-1 ) OP2 500 NA 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 400 0 100 200 300 C s eq ue st ra tio n in o il pa lm b io m as s (1 03 t ye ar -1 ) 600 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 400 500 -100 20 years 25 years 30 years Mean 0 100 200 400 Method 1 500 Method 2 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years Mean C s eq ue st ra tio n in o il pa lm b io m as s (1 03 t ye ar -1 ) 300 600 0 40 80 160 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years C s eq ue st ra tio n in o il pa lm b io m as s (1 03 t ye ar -1 ) 120 200 East North Central West 0 20 60 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years East North Central West C s eq ue st ra tio n in li tt er a nd g ro un d co ve r bi om as s (1 03 t ye ar -1 ) 40 40 C s eq ue st ra tio n in li tt er a nd g ro un d co ve r bi om as s (1 03 t ye ar -1 ) 160 1959 1964 1969 1974 1979 1984 1989 1994 1999 2004 2009 Years 80 120 0 20 years 25 years 30 years Mean FIgURE 6. Carbon sequestration in plantation litter and ground cover in Colombia from 1959 to 2009 showing (A), the effect of replanting cycle time (rCT), and (B), differences in carbon sequestration between regions (means of 20, 25 and 30-year rCT). All data are means of two methods of assessing replantedareas. Annual changes in carbon sequestration are shown in Fig. 5. Carbon loss occurred in 1980 and 2000 with 20-year cycles but otherwise carbon was sequestered in each year, mainly due to continued increase in planted area. Results for the NA oil palm are given in Tab. 4 and Fig. 6. Carbon sequestration in ground cover and litter averaged c. 25% of that sequestered in the palms, higher than the 19% found for Malaysia (Henson, 2009). Carbon sequestration related to previous land use Existing oil palm areas: oil palm-to-oil palm replanting Additional, but albeit short-term, carbon sequestration occurs in the residual palm biomass left on-site after fe- lling the old crop. The amount of carbon sequestered in the felled material was calculated as described by Henson (2009), based on previously observed rates of biomass decomposition. The majority of residues disappear within A B C D A B 365Henson, Ruiz R., and Romero: The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon... one to two years under moist tropical conditions (Khalid et al., 2000). The calculations ignore the possible use of fire for clearing, now generally discouraged. Burning would cause a rapid release of carbon and result in emissions of other GHGs besides CO2. The amount of residual material in any one year depends strongly on RCT (Henson, 2009), with long cycles result- ing in the least carbon sequestration. Thus, there is con- siderable uncertainty in the estimated amounts of carbon sequestration in felled material, which for the NA oil palm varied from ten to <1% of the carbon sequestered in the living palms when RCT was changed from 20 to 30 years (Tab. 5, Fig. 7). One cause of the lower felled biomass with longer RCT is the decline in biomass with age that is a marked feature of the NA oil palm (Fig. 1A). However, this only partly ex- plains the reduction, as there was a similar, but albeit less pronounced,decline with the OP2 oil palm, even though there was no marked decline in its biomass with age (data not presented). New oil palm areas: clearance of previous vegetation While most of the carbon in biomass present in areas newly cleared for oil palm planting will, as for oil palm, be rapidly lost, some will be retained on-site in amounts that vary with the initial quantity of biomass, its chemical composition, and the conditions during decomposition. 100 200 300 C in p al m b io m as s at fe lli ng ( 10 -3 t) 600 1978 1983 1988 1993 1998 2003 2008 400 500 0 20 years RCT 25 years RCT 30 years RCT Mean 100 200 300 C in f el le d pa lm b io m as s re s di ue s (1 0- 3 t) 1978 1983 1988 1993 1998 2003 2008 400 0 -40 0 C s eq ue st er ed in f el le d pa lm b io m as s re si du es ( 10 -3 t ye ar -1 ) 120 1978 1983 1988 1993 1998 2003 2008 Year 40 80 -80 FIgURE 7. Changes from 1978 to 2009 (the years during which palms were replaced) in (A) carbon in oil palm biomass at felling, (B) carbon in felled biomass after loss by decomposition, and (C) carbon sequeste- red annually in the felled residues. data (3-year running means) shown for the three replanting cycle times (rCTs) are for the NA oil palm. TABlE 5. Carbon stocks in oil palm biomass and biomass residues at and after felling, and amounts sequestered in different regions in Co- lombia, 1959 to 2009, calculated for the NA oil palm for three replanting cycle times (rCTs)1. RCT (year) Eastern Northern Central Western Whole country Carbon in biomass at felling (103 t) 20 37.56 28.74 27.07 9.80 103.17 25 10.21 13.23 13.10 5.53 42.07 30 1.41 2.24 4.59 0.70 8.95 Mean 16.39 14.74 14.92 5.34 51.40 Carbon in biomass residues (103 t) 20 26.40 20.62 19.43 7.05 73.50 25 7.22 9.51 9.43 3.98 30.14 30 1.01 1.57 3.26 0.50 6.33 Mean 11.54 10.57 10.71 3.84 36.66 Rate of carbon sequestration in residues (103 t year-1) 20 3.69 1.07 0.93 0.30 5.99 25 0.98 0.69 0.53 0.23 2.43 30 0.10 0.27 0.37 0.07 0.81 Mean 1.59 0.68 0.61 0.20 3.08 1Data are averages of 3-year running means over the study period. In the case of forests, carbon may also be stored in timber collected during felling. In Malaysia, the carbon sequestered in biomass residues and timber was small in relation to total carbon sequestra- tion (5.1% in residues and 3.0% in timber; Henson, 2009), despite a high proportion of forest being present in the areas cleared. In Colombia, the proportion of total carbon sequestered this way is likely to be even lower, given that much of the land converted to oil palm contained low initial biomass. Due to the sparseness of data, a detailed analysis was not possible but rough calculations based on likely A B C 366 Agron. Colomb. 30(3) 2012 previous land use (Tabs. 10 and 12) suggest sequestration in residues to amount to no more than 1.46·103 t year-1 of car- bon or c. 0.7% of total carbon sequestration. For harvested timber (assuming as for Malaysia a 10% timber yield and exponential decay over a 30-year period), sequestration was calculated to be 0.487·103 t year-1 of carbon, 89% of which was from the Western region. Carbon sequestration in mill products and by-products Carbon is sequestered or stored temporally in both the products of FFB milling; crude palm oil (CPO), palm kernel (PK), palm kernel oil (PKO) and palm kernel cake (PKC), and the by-products; empty fruit bunches (EFB), fibre, shell and palm oil mill effluent (POME). Yearly production of FFB, mill products and by-products are given in Fig. 8. Changes in standing stocks of carbon in mill materials, after allowing for decomposition and consumption, are shown in Fig. 9. Tab. 6 gives the mean annual values over the period studied for the four regions and for Colombia as a whole. Total sequestration was roughly proportional to land areas cultivated, but there were some differences bet- ween regions on a per hectare basis (Tab. 7). Approximately two-thirds of the carbon sequestered was in the products. TABlE 6. Carbon sequestered in palm oil mill products and by-products in Colombia, 1959 to 20091. Products/by-products Eastern Northern Central Western All Rate of carbon sequestration (103 t year-1) palm oil 1.923 2.146 2.000 0.135 6.204 palm kernel oil 0.154 0.148 0.171 0.011 0.484 palm kernel cake 0.092 0.088 0.102 0.007 0.288 Total products 2.169 2.382 2.272 0.153 6.976 EFB-burnt 0.010 0.009 0.010 0.001 0.029 EFB-un-burnt 0.337 0.303 0.352 0.034 1.027 Fibre-burnt 0.024 0.021 0.031 0.002 0.078 Fibre-un-burnt 0.269 0.242 0.285 0.027 0.823 shell-burnt 0.015 0.013 0.027 0.001 0.057 shell-un-burnt 0.338 0.304 0.368 0.034 1.043 palm oil mill effluent 0.086 0.050 0.153 0.009 0.298 Losses 0.014 0.012 0.014 0.001 0.042 Total by-products 1.093 0.955 1.241 0.109 3.397 Grand total 3.262 3.336 3.513 0.261 10.373 1Note, there were no FFB harvested before 1963. Carbon emission offsets Henson (2009), considered three factors that could offset the emission of CO2 and other GHGs resulting from palm oil production: i) the use of fibre and shell as boiler fuel to generate energy for running the mill, ii) the use of mill by-products as a source of nutrients, thus reducing the need for inorganic fertilizers, and iii) the generation of energy by burning methane (biogas) emitted from POME. 1,000 2,000 FF B , P O a nd P K p ro du ct io n 4,000 1959 1969 1979 1989 1999 2009 3,000 0 FFB PO PK P K , P K O a nd P K C p ro du ct io n 160 1959 1969 1979 1989 1999 2009 PK PKC PKO 40 80 120 0 100 200 M ill b y- pr od uc t p ro du ct io n 400 1959 1969 1979 1989 1999 2009 300 0 Fibre EFB Shell Year POME Losses (1 0- 3 t y ea r-1 ) (1 0- 3 t y ea r-1 ) (1 0- 3 t y ea r-1 ) FIgURE 8. Annual production (dry mass) of FFB (A), mill products (A, B) and mill by-products (C) in Colombia from 1959 to 2009. Biogas production is yet to be actively pursued in Colombia, although plans to promote it have been initiated (FEDE- PALMA, 2009). The present analysis is therefore confined to the other two factors. Offsets arising from mill by-product use The principle offset (Tab. 8) is the use of fibre and shell as boiler fuel. This serves to provide most, if not all, of the energy needed to run the mill. Even then, not all the po- tential energy from these sources is recovered. Other uses of mill by-products, such as field application of EFB and POME to supply nutrients and organic matter, contribute only a relatively small extent towards offsetting carbon losses. If EFB is burnt, the produced bunch ash can also A B C 367Henson, Ruiz R., and Romero: The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon... serve as a fertilizer, although this practice is being phased out based on broader environmental grounds. Potential for further offsetting fossil fuel use There are several other possible means of reducing fossil fuel use that have yet to be fully exploited, and while not presently included as items in the carbon budget, they may well become important in the future. Four such items are considered briefly here. Use of nutrients from biomass residues resulting from land-use change Biomass residues of previous vegetation contain nutrients that are released during decomposition, which could serve as sources for newly planted palms, thus reducing the need for inorganic fertilizers in early years of growth. In practice, this factor has seldom been taken into account due to the uncertainty of the uptake efficiency of young palms with restricted root systems. Research is required to quantify this potential nutrient source and devise means of handling residues, so as to maximize nutrient recovery. Collection and use of methane emitted from POME ponds There is considerable potential for trapping biogas and using it to generate electricity and at the same time reduce methane emissions to the atmosphere. Several schemes to do this have recently been initiated in Colombia (Fede- palma, 2009). Production of biodiesel There are now several biodiesel plants in Colombia. The impact of biodiesel production on GHG balance needs, however, to be carefully analyzed, as it can be either positive or negative depending on the land use change (LUC) during establishment of plantations. In cases that involve clearance of sites with high carbon stocks, biofuel production can increase, rather that decrease, GHG emissions (Danielsen et al., 2008; Fargione et al., 2008). TABlE 7. Amounts of carbon sequestered annually per producing hec- tare in palm oil products and by-products in Colombia, 1959 to 20091. Products/by-products Eastern Northern Central Western Rate of carbon sequestration (kg ha-1 year-1) palm oil 121.6 116.8 110.4 72.9 palm kernel oil 10.2 9.5 9.5 5.6 palm kernel cake 5.2 4.9 5.0 2.7 Total products 137.0 131.2 124.9 81.2 EFB-burnt 0.7 0.7 0.6 0.4 EFB-un-burnt 26.2 23.2 22.0 14.5 Fibre-burnt 1.7 1.6 1.6 1.0 Fibre-unburnt 19.4 18.5 17.6 11.6 shell-burnt 1.1 1.0 1.2 0.6 shell-unburnt 24.7 23.2 22.3 14.5 palm oil mill effluent 6.3 5.6 6.7 3.7 Losses 1.0 1.0 0.9 0.6 total by-products 81.1 74.7 73.1 46.9 Grand total 218.1 206.0 198.0 128.1 1 Note, there was no FFB harvested before 1963. FIgURE 9. Annual stocks of carbon in (A), mill products and (B), mill by-products, in Colombia from 1959 to 2009. 50 100 P O , P K , P K O a nd P K C C s to ck s (1 03 t) 300 1959 1969 1979 1989 1999 2009 250 0 PO PK PKO Year PKC200 150 M ill b y- pr od uc t c ar bo n st oc ks ( 10 3 t) 100 1959 1969 1979 1989 1999 2009 Shell Fibre EFB Year POME Losses 20 40 60 0 80 TABlE 8. Carbon equivalent offsets arising from the use of mill by-pro- ducts in Colombia, 1959 to 2009. Source Eastern Northern Central Western Whole country Carbon eq. (103 t year-1) Fibre and shell1 7.80 7.43 7.88 2.68 25.79 EFB 2 0.04 0.04 0.04 0.01 0.14 Bunch ash 3 0.07 0.07 0.07 0.02 0.24 pOME4 0.12 0.11 0.12 0.04 0.39 total 8.03 7.65 8.11 2.76 26.55 1Used as fuel in mill boilers. 2Nutrient and organic matter source. 3produced by burning EFB (22% burning is assumed). Used as fertilizer. 4Nutrient, organic matter and water source. A B 368 Agron. Colomb. 30(3) 2012 Production of biochar It is possible to reduce the rate of carbon loss during land clearing by producing biochar, a long-lived form of carbon resistant to decomposition, from the residues. The process involves controlled, on-site, burning. Studies to evaluate and optimize the process are currently in progress (Gar- zón, 2009). Total carbon sequestration The amounts of carbon sequestered in various forms du- ring the expansion of oil palm cultivation in Colombia are shown in Tab. 9. These values were obtained using ‘default’ options, namely the NA oil palm growth curve, a mean of 20, 25 and 30-year replanting cycles, and a mean of the two methods of estimating areas of replanting. These options are considered either the most likely or the least extreme of the possible alternatives. Lower values of sequestration were obtained using a 20-year RCT and higher values using a 30-year RCT or the OP2 growth curve (Fig. 3). The effects of alternative scenarios on sequestration rates are further presented and discussed in part II of this paper in the context of the total carbon budget. The oil palms formed the largest ‘sink’ for carbon followed by the other plantation components (litter and ground cover) and the various offsets. Mill products and by-prod- ucts contributed much less. On a per hectare basis, there was little difference between the four growing regions in the amounts sequestered, so that the totals for each region mainly reflected the area under cultivation. On the other hand, partly due to yield differences (Henson et al., 2011), sequestration per tonne of CPO was lower in the Central region than in the other regions. Conclusions This report is part of a preliminary estimate of the total GHG balance that results from the cultivation of oil palm TABlE 9. mean carbon sequestration per hectare, per tonne crude palm oil, and per region, in the four palm oil producing areas of Colombia, 1959 to 2009, calculated using default options. Component Eastern Northern Central Western Whole country Ceq. (t ha-1 year-1) Oil palm standing biomass 1.427 1.492 1.285 1.585 1.425 plantation litter and ground cover 0.378 0.384 0.337 0.326 0.364 Felled oil palm residues 0.052 0.025 0.025 0.029 0.033 Other residues 0.013 0.010 0.005 0.068 0.016 Wood harvest 0.001 0.001 <0.001 0.044 0.005 Mill products 0.071 0.088 0.092 0.015 0.076 Mill by-products 0.036 0.036 0.050 0.011 0.037 Offsets 0.263 0.282 0.328 0.278 0.288 Total 2.241 2.318 2.122 2.356 2.244 (t Ceq./t per year CPO) Oil palm standing biomass 0.606 0.620 0.455 0.643 0.568 plantation litter and ground cover 0.160 0.160 0.119 0.132 0.145 Felled oil palm residues 0.022 0.010 0.009 0.008 0.013 Other residues 0.005 0.004 0.002 0.026 0.006 Wood harvest <0.001 <0.001 <0.001 0.017 0.002 Mill products 0.030 0.036 0.033 0.006 0.030 Mill by-products 0.015 0.015 0.018 0.004 0.015 Offsets 0.112 0.117 0.116 0.113 0.115 Total 0.950 0.962 0.752 0.949 0.894 Ceq. (103 t year-1) Oil palm standing biomass 43.59 40.47 31.76 15.74 131.57 plantation litter and ground cover 11.54 10.42 8.34 3.24 33.54 Felled oil palm residues 1.59 0.68 0.61 0.20 3.08 Other residues 0.40 0.26 0.12 0.68 1.46 Wood harvest 0.02 0.03 0.01 0.43 0.49 Mill products 2.17 2.38 2.27 0.15 6.98 Mill by-products 1.09 0.96 1.24 0.11 3.40 Offsets 8.03 7.65 8.11 2.76 26.55 Total 68.42 62.85 52.46 23.31 207.07 % of total 33.0 30.4 25.3 11.3 100 369Henson, Ruiz R., and Romero: The greenhouse gas balance of the oil palm industry in Colombia: a preliminary analysis. I. Carbon sequestration and carbon... and the milling of FFB in Colombia. It sought to identify the relative importance of the various forms of sequestra- tion and highlight the uncertainties of the analysis. Hope- fully, this should promote further investigations aimed at improving the estimation of carbon sequestration. The quantification of sequestration, however, is only part of the process of producing a full account of the GHG balance. To complete the process, we also need to quantify GHG emissions. This is the subject of the second paper in this series, leading to the presentation of the full GHG budget. literature cited Danielsen, F., H. Beukema, N.D. Burgess, F. Parish, C.A. Bruhl, P.F. Donald, D. Murdiyarso, B. Phalan, L. Reijnders, M. Struebig, and E.B. Fitzherbert. 2008. Biofuel plantations on forested lands: double jeopardy for biodiversity and climate. Conserv. Biol. 23, 348-358. Fargione, J., J. Hill, D. Tilman, S. Polasky, and P. Hawthorne. 2008. Land clearing and the biofuel carbon debt. Science 319, 1235-1238. Fedepalma, Federación Nacional de Cultivadores de Palma de Aceite. 2009. Leaflet. What oil palm has to offer. Bogota. Fedepalma, Federación Nacional de Cultivadores de Palma de Aceite. 2010. Statistical Yearbook and earlier editions. Bogota. Garzón, E.M. 2009. Potential use of bio-coal as a soil improvement measure. p. 42. In: Proceedings of XVI International Oil Palm Conference. Fedepalma; Cenipalma, Cartagena, Colombia. Henson, I.E. 2003. The Malaysian national average oil palm: concept and evaluation. Oil Palm Bull. 46, 15-27. Henson, I.E. 2004. Modelling carbon sequestration and emissions due to oil palm cultivation and associated land use change in Malaysia. MPOB Technol. 27, 51. Henson, I.E. 2005. OPRODSIM, a versatile, mechanistic simulation model of oil palm dry matter production and yield. pp. 801-832. In: Proceedings of PIPOC 2005 Agriculture, Biotechnology and Sustainability Conference. Malaysian Palm Oil Board, Kuala Lumpur. Henson, I.E. 2008. The carbon cost of palm oil production in Ma- laysia. The Planter 84, 445-464. Henson, I.E. 2009. Modelling carbon sequestration and greenhouse gas emissions associated with oil palm cultivation and land-use change in Malaysia. A re-evaluation and a computer model. MPOB Technol. 31, 116. Henson, I.E., R. Ruiz, and H.M. Romero. 2011. The growth of the oil palm industry in Colombia. J. Oil Palm Res. 23, 1121-1128. Khalid, H., Z.Z. Zin, and J.M. Anderson. 2000. Decomposition processes and nutrient release patterns from oil palm residues. J. Oil Palm Res. 12, 46-63. Mejia, J.O. 2000. Cenipalma internal report. Centro de Investigación en Palma de Aceite (Cenipalma), Bogota. Romero, C., A. Moreno, and F. Munevar. 1999. Evaluación edafo- climática de las tierras del trópico bajo colombiano para el cultivo de la palma de aceite. Cenipalma; Corpoica, Bogota. Rubiano, D.J., G. Corzo, M. Fadul, A. Palacios, C. Martin, E. Boohorquez, S. Ramirez, L. Gomez, C. Sarmiento, and K. Ar- cilaburgos. 2009. Zonas ambientalmente aptas para el cultivo de la palma de aceite en Colombia. Ministerio de Ambiente, Vivienda y Desarrollo Rural, Bogota. Wood, B.J. and R.H.V. Corley. 1993. The energy balance of oil palm cultivation. pp. 130-143. In: Proceedings of 1991 PORIM International Palm Oil Conference, Agriculture. Palm Oil Research Institute of Malaysia, Kuala Lumpur.