Skip to content


  • Research
  • Open Access

Carbon budgets of potential tropical perennial grass cropping scenarios for bioenergy feedstock production

Carbon Balance and Management201813:17

  • Received: 9 November 2017
  • Accepted: 29 August 2018
  • Published:



The environmental costs of fossil fuel consumption are globally recognized, opening many pathways for the development of regional portfolio solutions for sustainable replacement fuel and energy options. The purpose of this study was to create a baseline carbon (C) budget of a conventionally managed sugarcane (Saccharum officinarum) production system on Maui, Hawaii, and compare it to three different future energy cropping scenarios: (1) conventional sugarcane with a 50% deficit irrigation (sugarcane 50%), (2) ratoon harvested napiergrass (Pennisetum purpureum Schumach.) with 100% irrigation (napier 100%), and (3) ratoon harvested napiergrass with a 50% deficit irrigation (napier 50%).


The differences among cropping scenarios for the fossil fuel-based emissions associated with agricultural inputs and field operations were small compared to the differences associated with pre-harvest burn emissions and soil C stock under ratoon harvest and zero-tillage management. Burn emissions were nearly 2000 kg Ceq ha−1 year−1 in the conventional sugarcane; whereas soil C gains were approximately 4500 kg Ceq ha−1 year−1 in the surface layer of the soil profile for napiergrass. Further, gains in deep soil profile C were nearly three times greater than in the surface layer. Therefore, net global warming potential was greatest for conventional sugarcane and least for napier 50% when deep profile soil C was included. Per unit of biomass yield, the most greenhouse gas (GHG) intensive scenario was sugarcane 50% with a GHG Index (GHGI, positive values imply a climate impact, so a more negative value is preferable for climate change mitigation) of 0.11 and the least intensive was napiergrass 50% when a deep soil profile was included (GHGI = − 0.77).


Future scenarios for energy or fuel production on former sugarcane land across the Pacific Basin or other volcanic islands should concentrate on ratoon-harvested crops that maintain yields under zero-tillage management for long intervals between kill harvest and reduce costs of field operations and agricultural input requirements. For napiergrass on Maui and elsewhere, deficit irrigation maximized climate change mitigation of the system and reduced water use should be part of planning a sustainable, diversified agricultural landscape.


  • Global warming potential
  • Greenhouse gas index
  • Napiergrass
  • Ratoon crop
  • Sugarcane


Interest in the production of renewable biofuel from lignocellulosic crops is gaining global recognition as a leading alternative energy scenario in future fuel markets. The negative environmental costs associated with fossil fuel consumption are being recognized and accounted on an international level due primarily to the adoption of the Kyoto Protocol in 2005. Since then, the development of alternative fuel sources has been a major concern for nations like the United States that rely heavily on imported fossil fuel. These pressures can result in a shift from the conventional food-crop agriculture to bioenergy systems in the United States and production of renewable biomass sources has been accelerated in recent years due to government regulations [1]. More specifically, regulations such as the Energy Independence and Security Act (EISA) of 2007 have mandated that fossil fuel sources must be mixed with at least 36 billion gallons of renewable fuel by the year 2022, meeting approximately 25% of liquid-based fuel needs by 2050 [2]. In support of national policy, Hawaii is under similar pressure resulting from the 2008 Hawaii Clean Energy Initiative (HCEI), which mandates a local, efficient, and renewable source of fuel be developed for Hawaii by the year 2030. The HCEI specifically requires that 40% of the states fossil energy be replaced with “locally generated renewable sources” over the next 17 years. These increasing fuel demands may place additional stress on already exploited agricultural lands resulting in land intensification and conversion if not managed conservatively.

Agricultural intensification can negatively impact soil carbon (C) storage, increase greenhouse gas (GHG) emissions, and offset the overall ecosystem C balance if managed incorrectly [24]. The energy sector en masse converted annual croplands into perennial biomass systems over the past decade, which was thought to provide GHG mitigation potential through an improvement in soil quality and a reduction in nutrient amendments [5, 6]. In recent years, however, conflicting results about the mitigation potential of these systems raised questions about the true system-level GHG offset [3, 7, 8]. Shortages in the world’s productive croplands and the impending scarcity of water made worse by climate change lead to additional uncertainty in the sustainability of increasing the world’s biofuel acreage [9, 10]. Due to the wide range of conflicting viewpoints on the bioenergy issue, documentation of local and regional datasets that quantify these uncertainties are of the utmost importance.

The total GHG balance of biofuel production is difficult to capture because of the large variation in cropping practices, land management, and equipment usage across these agricultural systems. Additionally, many recent studies focused on single-species scenarios that are specific to Brazilian ethanol production and are not entirely applicable to addressing regional issues outside of their study areas [3, 8, 11, 12]. Although these studies can be a powerful tool to advance the global implications of bioenergy production under a large agricultural infrastructure, additional small-scale studies are required to assess the tradeoffs between these systems and provide alternative management options on a local scale. To do this, a GHG balance needs to be created based on site-specific data that incorporates the energy inputs (fossil and non-fossil fuel) required to establish and maintain a bioenergy cropping system. For the majority of these GHG assessments, fossil energy inputs are converted using emission factors to carbon dioxide equivalents (CO2eq). Referred to as a global warming potential (GWP), these CO2eq are comparable on an international scale. This kind of GHG accounting allows for a direct comparison between the GHGs emitted by agricultural operations and the GHGs saved by the production of a renewable fuel sources [3, 8]. For the purposes of GHG accounting, all emissions were converted into a similar C equivalent (Ceq) in order to make a direct comparison with the amount of soil C stored under each cropping practice. This C budget related the GWP of each of these systems to their crop production, which allows for estimation of the overall GHG intensity (GHGI) of a production system [13].

To date, there has been no known C budget specific to the agricultural operations in Hawaii, which makes comparing future cropping scenarios against current practices impossible. With over 11,000 ha−1 of Maui being state designated and protected agricultural land, the future sustainability of this region may lie with the adoption of an energy crop scenario that is able to displace fossil fuel based GHG emissions. The purpose of this study was to create a baseline C budget for a conventionally managed sugarcane production system on (100% irrigation) Maui and compare it to three different future energy cropping scenarios: [1] conventional sugarcane with a 50% deficit irrigation (sugarcane 50%), [2] ratoon harvested napiergrass with 100% irrigation (napier 100%), and [3] ratoon harvested napiergrass with a 50% deficit irrigation (napier 50%). This comparison will help identify, in terms of their agricultural C budgets, the best-case scenario for future biofuel production in Hawaii and other Pacific Basin island nations.


Site description and experimental design

The field experiment was conducted in central Maui, Hawaii (20.89°N, 156.41°W) on Hawaiian Commercial and Sugar Company (HC&S) land. At the time of the study, HC&S was the only remaining sugarcane plantation in Hawaii. The experimental plots were within a highly weathered, very-fine, kaolinitic, isohyperthermic Typic Eutrotorrox of the Molokai series in field #609, which is approximately 100 meters above sea level and has a total commercial area of 72 ha. The soil is well drained, rocky, and has deep, well-defined horizons below the plow layer [14]. Soil pH was 7.97, and C concentration was 1.37% on average in the top 40 cm with a mean bulk density of 1.51 g cm3 as assessed by the baseline soil collection in 2011. During the trial period, average annual air temperature was 23.4 °C and annual precipitation was 241 mm, which are consistent with long-term averages for the area [15].

The experiment was a strip-plot, group-balanced design with two factors, irrigation and species with three replicates (blocks) (please see [16] for additional details). Irrigation was applied at the standard plantation rate (100%), and two deficit irrigation rates (75% and 50% of plantation standard). The original trial included four species, sugarcane (Saccharum officinarum), energycane (Saccharum officinarum x Saccharum spontaneum), napiergrass (Pennisetum purpureum Schumach.), and sweet sorghum (Sorghum bicolor (L.) Moench). For this study, two crops (sugarcane and napiergrass) were evaluated at two irrigation levels (50% and 100%). On June 26, 2011 the field plots were established in a recently harvested sugarcane field that had been in a cane-on-cane rotation for over 100 years. The sugarcane plots were planted with seed cane from an adjacent field within the HC&S plantation and napiergrass seed crop was supplied from a harvested population at the University of Hawaii’s research station in Waimanalo, Oahu.

Deficit irrigation treatments were applied to the field from November 13, 2011 depending on water availability across the plantation and were controlled with automated timers. From November 2011–October 2012, 1245 mm water ha−1 were applied to the 100% plots and 633 mm water ha−1 were applied to the 50% plots, for an actual deficit treatment of 50.8%. During the study period, the napiergrass plots were harvested every 6 months (a time interval that maximizes yield) on March 13, 2012 and September 25, 2012. The sugarcane crop was under a 2-year growth rotation and did not reach maturity during the scope of this study period and was not harvested.

Developing a baseline system for sugarcane

Sugarcane is a high-yielding, tropical C4 perennial grass of South Pacific origin. Tropical sugarcane biomass yields are up to 40 Mg dry wt ha−1 year−1 in Hawaii and 26 Mg dry wt ha−1 year−1 in Brazil [8, 17]. The species supports a drought resistant, robust root system that can improve soil structure and accumulate C on marginal lands [18, 19]. However, Hawaiian sugarcane has been grown on a 2-year crop cycle, reaching maturity after 24 months and then harvested after a low intensity burn followed by deep tillage and mechanized planting. Commercial sugarcane production existed on Maui for over 125 years. However, in January 2016 HC&S, the last remaining large-scale sugar producer, announced a wholesale transition on their 14,000 ha plantation to diversified agriculture including perennial grasses for forage, pasture, and bioenergy feedstock.

The development of a C budget or GHG analysis includes a detailed accounting of the agricultural inputs required to produce a specific crop. This includes quantifying the fossil fuel and non-fossil fuel based emissions. Fossil emissions are considered emissions resulting from fuel use during field preparation, planting, application of agrochemicals, harvesting, and maintenance [3, 5, 11, 12, 20]. Non-fossil emissions are considered biogenic GHG emissions that consider the production of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) as a result of the production system and are primarily a result of pre-harvest burn operations, soil GHG exchange, and residue management [3]. This baseline system will be referred to as the sugarcane 100% scenario for future comparisons.

Fossil emissions from fuel consumption and agricultural inputs

Field operations

Field operations (i.e., field preparation, harvest, fabrication and maintenance, seed propagation, and irrigation) often are considered to be hidden sources of emissions because of their indirect contribution to GHG flux. These emissions are caused by the burning of fossil fuel during equipment operation [12] and accounting can be challenging partly because of large variation in the descriptive energy units [20]. A standard unit of kg Ceq was used to assess the contribution of field operations to the total C budget. Emissions factors (EF) were used to convert the fuel use requirements of each operation to Ceq. The EFs used in this analysis represent a synthesis of the best-available and current values found in the literature. Specific information pertaining to type of equipment and usage has been obtained from personal communication with plantation staff.

When EFs found in the literature were inadequate to describe the operation on Maui, they were generated independently to best reflect current plantation practices. For example, irrigation emissions were calculated based on the average energy required to pump water across the plantation at a rate specific to Maui, 0.0057 L/kWh (personal communication with L. Jakeway, Chemical Engineer at HC&S, 2013). Additionally, this value was adjusted to account for the total amount of water applied to a field in the 2011–2012 year to correctly reflect the amount of renewable energy used on the plantation for these operations (approximately 3:1 renewable to fossil energy ratio). An additional input to the baseline scenario was a calculation of fossil emissions related to seed cane production, which were adjusted for the weight of seed cane used in planting operations [21].

Agricultural inputs

Emissions from agricultural inputs (i.e., fertilizer, herbicide, and lime) are a result of the energy required to produce, transport, and distribute these items [11, 20, 22]. A pre-emergence herbicide mix containing atrazine (1-chloro-3-ethylamino-5-isopropylamino-2, 4, 6-triazine), 2, 4-D (2, 4-dichlorophenoxyacetic acid), Prowl ((N-1- ethylpropyl)-3, 4-dimethyl-2, 6 dinitrobenzenamine), Rifle (3, 6-dichloro-2-methoxybenzoic acid), and Velpar (3-cyclohexyl- 6-dimethylamino-1-methyl-1, 3, 5-triazine-2, 4(1H,3H)-dione) was applied once three weeks after planting. Each plot received 345 kg N ha−1 (as liquid urea: 46-0-0) applied through the drip irrigation system. The fertilizer was applied monthly once the crops were established and concluded after 10 months. The timing and rate of urea application were optimized for the 2-year sugarcane crop and were based on current HC&S plantation practices. Deficit irrigation treatments were postponed during all fertilizer application events.

Quantification of the C emissions resulting from agricultural inputs in Maui sugarcane production was based on application rate and converted to Ceq with reported EF [11, 20]. An EF of 0.97 was used to convert the fertilizer application rate to Ceq. Lime (CaCO3) was applied at a rate of 2569 kg ha−1 prior to field planting and converted to Ceq using an EF of 0.12. Individual emission factors were identified for each chemical used in the herbicide mix reported by HC&S. These factors were averaged and a new EF (5.64) was developed for herbicide application specific to Maui.

Non-fossil emissions

Litter decomposition

Emissions from litter decomposition is a function of residue management. For example, N2O emissions can increase due to decomposition of leaf material following harvest and are greater in intact compared to burned fields [3]. The goal on Maui sugarcane fields has been to maintain 15% of total field biomass for crop residue. Using this percentage, litter Ceq were calculated based on an average biomass production of 80.4 Mg ha−1 year−1 reported by HC&S. Emission factors for litter decomposition are based on the amount of N in crop residues following harvest [3, 11, 12].

Pre-harvest burn emissions

Conventional cultivation of sugarcane in Hawaii included a pre-harvest, low intensity burn to remove unwanted leafy material prior to harvesting. Pre-harvest burning significantly increases GHG emissions through the production of CH4 and N2O and the release of black carbon (BC) to the environment [3, 12]. The values used for the cropping scenarios in Maui were based on a 15% residue retention rate of total field biomass. Current IPCC values outlined in the Guidelines for National Greenhouse Gas Inventories assessment in 2006 suggest an EF of 0.07 kg N2O per ton of dry matter burnt and an EF of 2.7 kg CH4 per ton of dry matter burnt. Black C has a GWP that is 500 times greater than CO2 on a 100-year time horizon and it is estimated that 1 kg of BC is created for every kg of trash burnt [12, 23].

Soil N2O and CH4 emissions

Field measurement of soil fluxes in CH4, and N2O began in October of 2011 and were sampled at least monthly until October of 2012 following the GRACEnet sampling protocols [24], as previously reported in detail (Pawlowski et al. in review). Mid-morning measurements were collected from sealed static PVC chambers affixed to permanent collars installed in the sugarcane and napiergrass rows and inter-rows. Samples were collected by sealing each chamber and using a 10 mL polypropylene syringe and extracting 8 mL of headspace air through a septum on the styrene lid at 0, 15, 30, 45, and 60 min after chamber closure. Each gas sample was immediately injected into an evacuated Exetainer® (Labco Limited, UK) fitted with a Doubled Wadded Teflon/Silicon septa (Labco Limited, UK) for short-term storage. Samples were analyzed using a Shimadzu GC-2014 Gas Chromatograph (Shimadzu Scientific Instruments, Inc.), which used a flame ionization detector to measure the concentration of CH4 and CO2 after methanization, and an electrical conductivity detector for N2O analysis. Flux rates were calculated by assuming a linear change in gas concentration over time [25, 26]. Cumulative annual emissions of N2O and CH4 were interpolated from daily fluxes and summer over the first year and reported in terms of kg CH4 ha−1 year−1 and kg N2O ha−1 year−1 [2, 13]. To discuss N2O and CH4 emissions in terms of a C balance, annual rates were converted into CO2eq using the IPCC 100-year horizon factors for calculating GWP. Therefore, when CO2 = 1 on a 100 year−1 time scale, then the GWP for N2O and CH4 are 298 and 25 respectively [2628]. For the purposes of making direct comparisons with soil C storage on these plots, the GWP values were converted to Ceq relative to C in CO2.

Soil C quantification

Ten baseline soil cores were collected in June 2011 using 20-cm depth increments up to a vertical depth of 2.4 m. Cores were extracted using a standard wet core diamond tipped drill bit with an internal diameter of 7 cm (Diamond Products Core Borer, Elyria, Ohio, USA). Each core barrel was inserted into the soil by a rotating hydraulic drill to minimize compaction within the barrel and to ensure accurate depth measurements. Soil samples were frozen at field moisture conditions until laboratory analysis. The cores were sieved at < 2 mm and dried for 48 h at 105 °C. Subsamples were ground to pass through a 250 micron sieve for heterogeneity, weighed, and analyzed for C concentration by combustion using a Costech ECS 4010 CNH Analyzer (Costech Analytical Technologies, Inc., Valencia, CA, USA). Soil C stock was determined with the equivalent soil mass method in increments of 3600 Mg ha−1 and a mean value for the baseline cores was determined at the 7200 and 18,000 Mg ha−1 reference masses, which represent the surface layer (25–40 cm) and deep profile (1–1.4 m) of soil [29].

Three soil cores (65-mm inner diameter bucket auger) were collected annually to a depth of 120 cm in 20 cm increments from each of the experimental plots. Samples were processed as noted above for the baseline cores for C concentration and soil C stock determination at the 18,000 Mg ha−1 reference mass. The change (∆) in soil C stock was made for each plot by subtracting the mean baseline value from the mean of the three cores for each plot and was previously reported [16]. The mean ∆ values for experimental years 1 and 2 were reported here as an annualized ∆.

Alternative-cropping scenarios

With local and national mandates for cellulosic energy sources increasing, there is opportunity to consider alternative-cropping scenarios to sugar production that may increase financial security and ameliorate environmental impacts caused by conventional management practices. One potential candidate species is napiergrass, a warm-season perennial C4 grass of African origin that produces high biomass yields under tropical conditions. In recent studies, napiergrass produced more than 45 Mg dry wt ha−1 year−1 in Florida [30] and between 40 and 53 Mg dry wt ha−1 year−1 in Hawaii [31, 32]. In contrast to sugarcane production in Hawaii, napiergrass is a ratoon-harvested crop that maintains high yields with zero-tillage management and no pre-harvest burning. Regional water shortages from drought and limitations to water access from litigation impose uncertainty on Maui’s agricultural communities. As a result, experimental scenarios include deficit irrigation to better anticipate changes in crop yields and GWP brought on by limited water. For this experiment, the napiergrass plots received the same amount of fertilizer as the sugarcane plots.

For the purpose of creating a C budget for these alternative-cropping scenarios the following assumptions were made: (1) napiergrass fossil-based emissions during the field preparation, harvest, and maintenance operations were considered equal to baseline sugarcane operations due to the similarities in production requirements between the species, (2) seed emissions for napiergrass were calculated based on actual planting rate (seed pieces per plot) relative to the sugarcane seed rate (6500 kg ha−1), because a total weight of napiergrass seed was not known, (3) both the lime and herbicide application rate were considered equal to the baseline sugarcane system, regardless of the species or irrigation treatment, (4) deficit irrigation values were adjusted to the known amount of water applied to each of the 100% and 50% treatments, (5) total biomass for each scenario was adjusted based on a 44% reduction in yield found under the deficit irrigation treatments measured in the field at the time, (6) litter decomposition rates and emissions were reduced for the deficit treatments based on residue management and napiergrass decomposition was calculated based on unburned residue N content after harvest, (7) emissions resulting from pre-harvest burn operations were adjusted for the sugarcane 50% scenario by accounting for the reduction in yield under the deficit irrigation, and (8) no SOC accumulates under the baseline scenario due to the incorporation of the pre-burn harvest followed by the associated soil disturbance of conventional tillage operations after harvest.

Calculating GHG intensity

Greenhouse gas index (GHGI) relates net GWP to crop production by dividing net GWP values by total crop biomass [13]. Net GWP was calculated based on the net C balance for each scenario (described above) and reported in kg Ceq ha−1 year−1. In the first 2 years of the trial, yields at 100% irrigation were 73.9 Mg ha−1 year−1 for sugarcane and 47.7 Mg ha−1 year−1 for napiergrass [16]. These values were similar to 2011 sugarcane crop yields data reported from HC&S as 80.4 Mg ha−1 year−1. Similarly, best-case Hawaiian napiergrass biomass averages reported from Kinoshita et al. [31] and Kinoshita and Zhou [17] were approximately 45 Mg ha−1 year−1. The measured reductions in yields under the deficit irrigation treatment during the first 2 years of the field trial were 60% for sugarcane (29.4 Mg ha−1 year−1) and 31% (32.9 Mg ha−1 year−1) for napiergrass [16]. This metric is useful for comparing scenarios with different yield potentials relative to their impact on GWP and C balance: more positive GHGI values are a stronger source of atmospheric GHGs per unit crop yield, whereas, more negative values are a stronger sink.

Results and discussion

Baseline system C balance

The majority of C emissions in the baseline scenario (sugarcane 100%) were a result of fossil fuel-based field operation inputs and BC emissions (Fig. 1a). Agricultural inputs together (778 kg Ceq ha−1 year−1) were a smaller contributor to total emissions than field operations (2181 kg Ceq ha−1 year−1) and factors associated with burning (1644 and 320 kg Ceq ha−1 year−1 for BC and CH4 + N2O emission respectively) (Table 1). Emissions associated with seed propagation comprised 89% of the field operations total input. Individually, the other field operations inputs and soil GHG flux and litter decomposition were minor, contributing less than 100 kg Ceq ha−1 year−1 each). In contrast, 1964 kg Ceq ha−1 year−1 were released through a single burn event. This system was a significant net source of C to the atmosphere, with a net GWP of 4949 kg Ceq ha−1 year−1 (Fig. 1b).
Fig. 1
Fig. 1

Carbon emission balance by agricultural input and output (a), net C balance for each energy crop scenario (b)

Table 1

Agricultural carbon (C) balance of bioenergy cropping scenarios in Maui, Hawaii

C inputs and outputs


kg Ceq ha−1 year−1

Sugarcane 100%

Sugarcane 50%

Napiergrass 100%

Napiergrass 50%

Fossil emission

 Field operations

  Field preparationa

   Lime application

John Deere 5000–7000 (140–360 HP)






John Deere 9000 (540 HP)






John Deere 9000 (540 HP)





    Strip tillage

John Deere 9000 (540 HP)






Game (250–300 HP)





    Herbicide application

John Deere 5000–7000 (140–360 HP)






   Cane harvester

John Deere 3522 (337 HP)






Cat 950 and John Deere 624 (175–200 HP)






Cat 773B-E (650–675 Hp)






   Farm machinery

Tractors, implements and trucks





  Seed propagationd












Field operations subtotal





Agricultural inputs

  Fertilizer applicationf

344.68 kg ha−1 year−1 for SC, 243.63 kg ha−1 year−1 for N






Applied at a rate 20.307 kg ha−1 year−1 for first year





  Lime (CaO)h

Applied at a rate of 2569 kg ha−1 every 2 years





 Agricultural inputs subtotal





Fossil emission subtotal





Non-fossil emissions

  Litter decompositioni






  Pre-harvest burning emissionsj

   Burn emissions of CH4 and N2O





   Black carbon






   Soil emisisonsk

    N2O emission






    CH4 emission


− 29.4

− 19.7

− 40.9

− 35.4

Non-fossil emissions subtotal





Total emissions






  Δ Soil carbon (surface layer)l



− 6828

− 6820

  Δ Soil carbon (deep profile)l



− 14,571

− 17,359

Net flux (surface layer)


− 4574.5

− 4621.1

Net flux (deep profile)


− 12,317.5

− 15,160.1

aData from Lal [20] and Macedo et al. [11] for emissions and EF related to fuel consumption from farming operations during field preparation. Equipment description and HPs gathered from personal communication with HC&S (2013). EF used 0.853

bData from Lal [20] and Macedo et al. [11] for emissions and EF related to fuel consumption from harvest operations. Equipment description and HPs gathered from personal communication with HC&S (2013). EF used 0.853

cValues for fabrication and maintenance of equipment from Samson et al. [5]

dEF 0.3 used from Six and Jastrow [21], based on total seed weight of 6500 kg ha−1 for sugarcane and 4334 kg ha−1 for napiergrass

eUsed Lal [20] conversion factor of kWh = 7.25 × 10−2 CE for conversion of total energy used to pump 1 MG of water. Total energy was reduced by 75% based on personal communication with HC&S to account for renewable energy

fEF for N fertilization used from Lal [20] and Macedo et al. [11]. EF used 0.97; SC sugarcane, N napiergrass

gIndividual herbicide EF were calculated from Lal [20], EF used 5.64

h EF used 0.12 from Macedo et al. [11]

iAssuming trash input is 15% of total biomass yield. Based on N content in residual litter of 45 kgN ha−1 in unburnt systems and 7 kgN ha−1 in burnt systems (Cerri [42]; Galdos et al. [12]; Lisboa et al. [3])

jAssuming 15% residue. EF 0.07 kg N2O DM burnt and 2.70 kgCH4 DM burnt (IPCC 2006; Macedo et al. [11]). BC based on 1.0 kg BC Mg trash burnt and a GWP of 500 relative to CO2 (Sanhueza [23])

kUsed IPCC (2007) GWP of 25 for CH4 and 298 for N2O

lChange (Δ) in SOC was measured by the equivalent soil mass method in year one and two of cultivation, reported value is the mean annual Δ

As of January 2016, HC&S, the last remaining large-scale sugar producer in the Hawaiian islands announced a transition from intensive sugarcane production to diversified agriculture with a commitment to conservation management practices and sustainability. Moving forward, this baseline system may serve as an indicator for improvements in plantation-level sustainability by tracking individual components of fossil and non-fossil fuel based emissions as well as the net GWP. Black C emissions associated with the pre-harvest burn practice already ceased and C costs associated with seed propagation should decline with the phasing-out of 2-year crop rotations. As these components contributed the majority of C emissions to the conventional sugarcane system even under deficit irrigation, the net GWP of the plantation will decrease immediately.

Comparison of crop scenarios

The differences among crop scenarios for the fossil fuel-based emissions associated with agricultural inputs and field operations were small compared to the differences associated with pre-harvest burn emissions and soil C stock under ratoon harvest and zero-tillage management (Fig. 1a). As a result, the net GWP was greatest for the sugarcane 100% followed by sugarcane 50% (Fig. 1b). The deficit irrigation treatment in sugarcane resulted in a 21.5% reduction in the net GWP to 3885 kg Ceq ha−1 year−1 (Table 1). Largely because the napiergrass scenarios do not include burning and accumulate soil C, the net GWP was negative (and therefore a sink rather than a source of C) for both the 100% and 50% scenarios.

The critical component of the net C sink under alternative scenarios was the cultivation of a ratoon-harvested energy crop that maintains high yields with zero-tillage management and requires no pre-harvest burning. In this study, we focused on the cultivation of napiergrass; however, other grasses such as sugarcane and energycane have the same potential to improve GWP from the conventional sugarcane practices. Recent studies suggest a range of soil C responses from perennial grasses depending on genotype [33], species or accession [32], and soil properties [34]. The soil C change from the initial shift to ratoon harvested napiergrass from over 100 years of intensive, monocropping cultivation was large and rapid. The annual soil C change is likely to decrease over time as the initial phase of accumulation ends and C saturation approaches [35]. How high soil C stock will go and how long the system will take to reach equilibrium remains unknown.

The difference between burned and un-burned GHG balances has been well documented in recent literature. For example, Lisboa et al. [3] found that over 17% of total emissions from sugarcane production in Brazil was due to pre-harvest trash burning and Galdos et al. [12] found that the single largest source of emissions in a traditional sugarcane system in Brazil was also due to biomass burning, resulting in over 1732 kg Ceq ha−1 year−1. Our values ranged from 938 to 1964 kg Ceq ha−1 year−1 depending on the irrigation treatment.

The sugarcane scenarios generated the largest fossil-based emissions: 2959 and 2930 kg Ceq ha−1 year−1 for the 100% and 50% respectively, whereas the 100% and 50% napiergrass scenarios were less: 2191 and 2162 kg Ceq ha−1 year−1 respectively (Table 1). The single largest contributor to fossil emissions for all scenarios was the production of seed: 1950 and 1300 kg Ceq ha−1 year−1 for sugarcane and napiergrass respectively, which accounted for more than 60% of fossil-based emissions from all of the scenarios. The seed emissions were the driving factor for the differences in fossil based emissions between the sugarcane and napiergrass scenarios. Napiergrass may be ratooned for up to 10 years without yield declines [31], so seed propagation should decline depending on management decisions concerning kill-harvest cycles.

Fertilizer contribution to GWP was over 12% of total fossil emissions for sugarcane and 11% for napiergrass, regardless of the irrigation scenario. The application of mineral fertilizers contributes a substantial amount of emissions in other sugarcane systems [3, 11]. The potential exists for the fertilizer contribution to GWP to increase over time. Many of the large inputs accounted for in this study (seed propagation, liming and field preparation) are associated with the first-year of cultivation. Their effects on overall emissions are expected to decline, especially in subsequent years that allow for ratoon harvests. Therefore, the effects of N fertilizer may become more important because it will most likely remain a constant source of inputs as nutrients are added to the system each year in order to maintain yields [6]. However, in the napiergrass scenarios, as soil quality improves over time with the accumulation of organic matter under conservation practices, the fertilizer requirement may be reduced [36].

Although the reduction in irrigation did not have a large direct effect on the fossil fuel based emissions, the impact on biomass yield and soil C in deep soil horizons in napiergrass was high. Many of the inputs for GHG accounting require biomass yields in order to calculate the associated emissions. Our estimates of the yield reduction caused by the 50% deficit resulted in a substantial reduction in litter decomposition and burn-related emissions. However, these reductions only influenced the net C balance of the sugarcane scenarios. The difference in net C balance between the napiergrass irrigation scenarios were minimal if only the surface soil horizon was considered, but substantial if the deep soil profile was considered.

Driven by a lack of pre-harvest burn, non-fossil based emissions were markedly different between the sugarcane and napiergrass systems (Fig. 1a). Even within the sugarcane scenarios, the 100% irrigation scenario produced nearly double the emissions as the 50% scenario (Table 1) due to reduction in biomass under the deficit irrigation. Although a minor contributor, litter decomposition was reduced under deficit irrigation due to lower yield, and therefore lower residual plant litter. Soil N2O and CH4 fluxes were comparable between the scenarios [16].

Soil C increased rapidly from the baseline during the 2-year crop cycle [16] resulting in substantial C credit towards the net GWP for the napiergrass scenarios (Fig. 1a). In the surface layer, the annual increase in soil C was nearly identical between the napiergrass scenarios (Table 1). However, soil C increased more in the deficit irrigation scenario if the deeper soil horizons also were considered. Because of the large contribution of soil C output to GWP in the napiergrass scenarios, the net GWPs were − 4574 and − 4621 Ceq ha−1 year−1 for the napiergrass 100% and 50% scenarios respectively if only accounting for the surface soil layer. If accounting for the deep soil profile, then the net GWP decreases further to − 12,317 Ceq ha−1 year−1 (napiergrass 100%) and − 15,160 Ceq ha−1 year−1 (napiergrass 50%).

The vast majority of studies take only the surface soil layer (typically 0–15 or 30 cm) into account for C stock assessment and inventory [37]; however, this practice ignores deeper soil processes associated with roots, the rhizosphere and dissolved organic matter transport. As pointed out above, these processes could result in soil priming (and C loss) [38], soil C accrual etc. As a tropical, perennial C4 grass, napiergrass is known to generate deep, extensive root systems in search of nutrients and water, particularly during drought [39]. Accordingly, soil C accrual as a result of root-derived inputs to the soil was even more prominent in deeper horizons of the deficit irrigation treatment compared to the 100% of plantation practice treatment.

Greenhouse gas intensity

Per unit of biomass yield, the most GHG intensive scenario was sugarcane 50% (GHGI = 0.11) and the least intensive was napiergrass 50% when the deep soil profile was considered in the net GWP calculation (GHGI = − 0.77) (Table 2). The reduction in sugarcane yield with 50% less water was not matched proportionately by reductions in C emissions, so the deficit irrigation did not provide a net benefit to the GWP of the system. For napiergrass, the opposite was true − 50% less water improved the GHGI index (i.e., made it more negative), particularly if the C accrual in the deep soil profile was considered. Even if only the surface layer of soil C accrual was accounted for, because the net C sink was not different and the yield was reduced, the deficit irrigation provided a net benefit to the GWP of the system.
Table 2

Greenhouse gas index (GHGI) of the baseline and three alternative cropping scenarios


Sugarcane 100%

Sugarcane 50%

Napiergrassa 100%

Napiergrassa 50%

Surface layer



− 0.10

− 0.23

Deep profile



− 0.27

− 0.77

aFor napiergrass, GHGI were calculated from net C balances that included both surface layer and deep profile soil C accrual


Moving away from intensive sugar production, the napiergrass systems demonstrated the greatest potential for climate change mitigation, whereas sugarcane scenarios that maintain conventional harvest practices are large sources of C to the atmosphere. Reduced emissions from sugarcane systems will rely on removal of burn related emissions and through the reduction of fertilizer use or improvement of nutrient use efficiency in these systems [40]. In the Pacific Basin, future scenarios for energy or fuel production on former sugarcane land should concentrate on ratoon-harvested crops that maintain yields under zero-tillage management for long intervals between kill harvest (if they are required at all). For napiergrass on Maui, deficit irrigation from the conventional sugarcane practice provided the greatest net climate change mitigation and reduced water usage should be an important aspect of planning a sustainable, diversified agricultural landscape under pressure for improved water resource management.



black carbon




carbon equivalent




carbon dioxide


emission factor


Energy Independence and Security Act


greenhouse gas


greenhouse gas index


global warming potential


Hawaii Clean Energy Initiative


Hawaiian Commercial & Sugar


nitrous oxide


Authors’ contributions

MP, MNM, JRK, and SEC contributed to the conception of the problem, theorizing and designing experiments (i.e., “ideas”), and/or interpreting data and drafting the article or revising it critically for important intellectual content (i.e., “writing”). MP and MNM were responsible for data acquisition, data analysis, coding, and/or using the instruments/models (i.e., “work”). MNM, JRK, and SEC were responsible for developing resources, and guaranteeing the overall integrity of the work as a whole before and after publication (i.e., “stewardship”). All authors read and approved the final manuscript [41].


We thank Andrew Hashimoto for his leadership and support of this research; Richard Ogoshi and Adel Youkhana for work on the field trial; Mae Nakahata, Hawaiian Commercial & Sugar and the field crews; Jason Drogowski, Spencer Mawhar, Daniel Richardson, and Michael Ross for assistance in the field; and Jon Wells, Nancy Parker, and Christine Glazer for assistance with laboratory and data analysis.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Data will be made available via ScholarSpace (, an open-access, digital institutional repository for the University of Hawaii Manoa community.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


This work was supported by the Department of Energy [Award Number DE-FG36-08GO88037], the Office of Naval Research [Grants N00014-12-1-0496 and N00014-16-1-2221], and the United States Department of Agriculture-Agricultural Research Service [Interagency Reimbursable Agreement #60-3098-5-006, with funds provided by Office of Naval Research]. This work was further supported by the USDA National Institute of Food and Agriculture, Hatch project [Project HAW01130-H], managed by the College of Tropical Agriculture and Human Resources. These funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Authors’ Affiliations

Dept. of Natural Resources and Environmental Management, University of Hawaii at Manoa, 1910 East-West Rd., Honolulu, HI 96822, USA
Texas A & M AgriLife Research, Blackland Research and Extension Center, 720 E. Blackland Rd., Temple, TX 76502, USA
USDA, Agricultural Research Service, Grassland Soil and Water Research Laboratory, 808 E. Blackland Rd., Temple, TX 76502, USA


  1. Dominguez-Faus R, Powers SE, Burken JG, Alvarez PJ. The water footprint of biofuels: a drink or drive issue? Environ Sci Technol. 2009;43(9):3005–10.View ArticleGoogle Scholar
  2. Gelfand I, Sahajpal R, Zhang X, Izaurralde RC, Gross KL, Robertson GP. Sustainable bioenergy production from marginal lands in the US Midwest. Nature. 2013;493(7433):514–7.View ArticleGoogle Scholar
  3. Lisboa CC, Butterbach-Bahl K, Mauder M, Kiese R. Bioethanol production from sugarcane and emissions of greenhouse gases—known and unknowns: GHG emissions and bioethanol from sugarcane. GCB Bioenergy. 2011;3(4):277–92.View ArticleGoogle Scholar
  4. Robertson GP, Hamilton SK, Del Grosso SJ, Parton WJ. The biogeochemistry of bioenergy landscapes: carbon, nitrogen, and water considerations. Ecol Appl. 2011;21(4):1055–67.View ArticleGoogle Scholar
  5. Samson R, Mani S, Boddey R, Sokhansanj S, Quesada D, Urquiaga S, et al. The potential of C4 perennial grasses for developing a global bioheat industry. Crit Rev Plant Sci. 2005;24(5–6):461–95.View ArticleGoogle Scholar
  6. Drewer J, Finch JW, Lloyd CR, Baggs EM, Skiba U. How do soil emissions of N2O, CH4 and CO2 from perennial bioenergy crops differ from arable annual crops? GCB Bioenergy. 2012;4(4):408–19.View ArticleGoogle Scholar
  7. Crutzen PJ, Mosier AR, Smith KA, Winiwarter W. N2O release from agro-biofuel production negates global warming reduction by replacing fossil fuels. Atmos Chem Phys. 2008;8(2):389–95.View ArticleGoogle Scholar
  8. Walter A, Dolzan P, Quilodrán O, de Oliveira JG, da Silva C, Piacente F, et al. Sustainability assessment of bio-ethanol production in Brazil considering land use change, GHG emissions and socio-economic aspects. Energy Policy. 2011;39(10):5703–16.View ArticleGoogle Scholar
  9. Berndes G. Bioenergy and water—the implications of large-scale bioenergy production for water use and supply. Glob Environ Change. 2002;12(4):253–71.View ArticleGoogle Scholar
  10. Pimentel D. Food for thought: a review of the role of energy in current and evolving agriculture. Crit Rev Plant Sci. 2011;30(1–2):35–44.View ArticleGoogle Scholar
  11. Macedo IC, Seabra JEA, Silva JEAR. Green house gases emissions in the production and use of ethanol from sugarcane in Brazil: the 2005/2006 averages and a prediction for 2020. Biomass Bioenergy. 2008;32(7):582–95.View ArticleGoogle Scholar
  12. Galdos MV, Cerri CC, Lal R, Bernoux M, Feigl B, Cerri CEP. Net greenhouse gas fluxes in Brazilian ethanol production systems. GCB Bioenergy. 2010;2(1):37–44.View ArticleGoogle Scholar
  13. Mosier AR, Wassmann R, Verchot L, King J, Palm CA. Methane and nitrogen oxide fluxes in tropical agricultural soils: sources, sinks and mechanisms. In: Wassmann R, Vlek PLG, eds. Tropical agriculture in transition—opportunities for mitigating greenhouse gas emissions? Dordrecht: Springer; 2004. p. 11–49. Accessed 28 Jun 2016.View ArticleGoogle Scholar
  14. Soil Survey Staff, Natural Resources Conservation Service, United States Department of Agriculture. Web Soil Survey. Accessed 11 Nov 2011.
  15. Giambelluca TW, Chen Q, Frazier AG, Price JP, Chen Y-L, Chu P-S, et al. Online rainfall atlas of Hawai‘i. Bull Am Meteorol Soc. 2013;94(3):313–6.View ArticleGoogle Scholar
  16. Pawlowski MN, Crow SE, Meki MN, Kiniry JR, Taylor AD, Ogoshi R, et al. Field-based estimates of global warming potential in bioenergy systems of hawaii: crop choice and deficit irrigation. PLoS ONE. 2017;12(1):e0168510.View ArticleGoogle Scholar
  17. Kinoshita CM, Zhou J. Siting evaluation for biomass-ethanol production in Hawaii. University of Hawaii Manoa; 1999. Report No.: XXE-8-1709.Google Scholar
  18. Smith DM, Inman-Bamber NG, Thorburn PJ. Growth and function of the sugarcane root system. Field Crops Res. 2005;92(2–3):169–83.View ArticleGoogle Scholar
  19. Anderson-Teixeira KJ, Snyder PK, Twine TE, Cuadra SV, Costa MH, DeLucia EH. Climate-regulation services of natural and agricultural ecoregions of the Americas. Nat Clim Change. 2012;2(3):177–81.View ArticleGoogle Scholar
  20. Lal R. Carbon emission from farm operations. Environ Int. 2004;30(7):981–90.View ArticleGoogle Scholar
  21. Six, Jastrow JD. Organic matter turnover. In: Encyclopedia of Soil Science. New York: Marcel Dekker, Inc.; 2002. p. 936–42.Google Scholar
  22. West TO, Marland G. A synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: comparing tillage practices in the United States. Agric Ecosyst Environ. 2002;91(1–3):217–32.View ArticleGoogle Scholar
  23. Sanhueza, Eugenio. Potential emissions of Kyoto and Non-Kyoto climate active compounds in the production of sugarcane ethanol. Interciencia. 2009;
  24. Del Grosso SJ, White JW, Wilson G, Vandenberg B, Karlen DL, Follett RF, et al. Introducing the GRACEnet/REAP data contribution, discovery, and retrieval system. J Environ Qual. 2013;42(4):1274.View ArticleGoogle Scholar
  25. Kroon PS, Hensen A, van den Bulk WCM, Jongejan PAC, Vermeulen AT. The importance of reducing the systematic error due to non-linearity in N2O flux measurements by static chambers. Nutr Cycl Agroecosyst. 2008;82(2):175–86.View ArticleGoogle Scholar
  26. Smith KE, Runion GB, Prior SA, Rogers HH, Torbert HA. Effects of elevated CO2 and agricultural management on flux of greenhouse gases from soil. Soil Sci. 2010;175(7):349–56.View ArticleGoogle Scholar
  27. Gauder M, Butterbach-Bahl K, Graeff-Hönninger S, Claupein W, Wiegel R. Soil-derived trace gas fluxes from different energy crops—results from a field experiment in Southwest Germany. GCB Bioenergy. 2012;4(3):289–301.View ArticleGoogle Scholar
  28. Meki MN, Kiniry JR, Behrman KD, Pawlowski MN, Crow SE. The Role of Simulation Models in Monitoring Soil Organic Carbon Storage and Greenhouse Gas Mitigation Potential in Bioenergy Cropping Systems. In: Esteves V, editor. CO2 Sequestration and Valorization. InTech; 2014. Accessed 23 Jun 2016.
  29. Wendt JW, Hauser S. An equivalent soil mass procedure for monitoring soil organic carbon in multiple soil layers. Eur J Soil Sci. 2013;64(1):58–65.View ArticleGoogle Scholar
  30. Valencia-Gica RB, Yost RS, Porter GS, Pattnaik R. Retention, accumulation, and movement of phosphorus in a mollisol soil irrigated with dairy effluent in a tropical environment. Soil Sci. 2010;175(10):500–10.View ArticleGoogle Scholar
  31. Kinoshita CM, Ishimura DM, Jakeway L, Osgood R. Production of biomass for electricity generation on the island of Oahu. Hawaii Natural Energy Institute; 1995.Google Scholar
  32. Sumiyoshi Y, Crow SE, Litton CM, Deenik JL, Taylor AD, Turano B, et al. Belowground impacts of perennial grass cultivation for sustainable biofuel feedstock production in the tropics. GCB Bioenergy. 2016.
  33. Richter GM, Agostini F, Redmile-Gordon M, White R, Goulding KWT. Sequestration of C in soils under Miscanthus can be marginal and is affected by genotype-specific root distribution. Agric Ecosyst Environ. 2015;200:169–77.View ArticleGoogle Scholar
  34. Tiemann LK, Grandy AS. Mechanisms of soil carbon accrual and storage in bioenergy cropping systems. GCB Bioenergy. 2015;7(2):161–74.View ArticleGoogle Scholar
  35. Six J, Conant RT, Paul EA, Paustian K. Stabilization mechanisms of soil organic matter: implications for C-saturation of soils. Plant Soil. 2002;241(2):155–76.View ArticleGoogle Scholar
  36. Lal R. Soil carbon management and climate change. Carbon Manag. 2013;4(4):439–62.View ArticleGoogle Scholar
  37. de Richter D, Billings SA. “One physical system”: Tansley’s ecosystem as Earth’s critical zone. New Phytol. 2015;206(3):900–12.View ArticleGoogle Scholar
  38. Cheng W, Parton WJ, Gonzalez-Meler MA, Phillips R, Asao S, McNickle GG, et al. Synthesis and modeling perspectives of rhizosphere priming. New Phytol. 2014;201(1):31–44.View ArticleGoogle Scholar
  39. Knoll JE, Anderson WF, Strickland TC, Hubbard RK, Malik R. Low-input production of biomass from perennial grasses in the coastal plain of Georgia, USA. BioEnergy Res. 2012;5(1):206–14.View ArticleGoogle Scholar
  40. Otto R, Castro SAQ, Mariano E, Castro SGQ, Franco HCJ, Trivelin PCO. Nitrogen use efficiency for sugarcane-biofuel production: what is next? BioEnergy Res. 2016;9:1–18.View ArticleGoogle Scholar
  41. Clement TP. Authorship matrix: a rational approach to quantify individual contributions and responsibilities in multi-author scientific articles. Sci Eng Ethics. 2014;20(2):345–61.View ArticleGoogle Scholar
  42. Cerri CC. Dinaˆmica da mate ́ria orgaˆnica do solo no agrossis- tema cana-de-ac ̧u ́car. Dissertation, Universidade de Sa ̃o Paulo, Piracicaba, Brazil; 1986.Google Scholar


© The Author(s) 2018