Coupling Exergy with the Emission of Greenhouse Gases in Bioenergy: A Case Study Using Biochar

: The loss in the quality of energy throughout any process can be assessed by the thermodynamics magnitude related to its entropic performance—the exergy. This indicator has been suggested as an environmental index, as an alternative to life cycle assessment (LCA), which is a classic tool for this purpose. This study assesses the potential of coupling the life cycle approach and exergy in a bioenergy supply chain environmental performance characterization, examining two scenarios in the sugarcane agroindustry. The ﬁrst one, the reference scenario, is a classical production, and the second includes the returning of a portion of residual biomass from the plant, in the form of biochar, to agricultural soil. The use of biochar engendered an increase in sugarcane productivities and a reduction of nitrous oxide emissions. These changes resulted in scenarios 1 and 2, reducing the exergy destroyed from 390 to 355 MJ/MJ ethanol (9.0%) and decreasing the greenhouse gases emissions (GHG) from 11.8 to 11.0 g CO 2-equivalent /MJ ethanol (6.8%). The latter represents an improvement in the use of carbon. A sensitivity analysis showed that the effect of changing productivity was quite signiﬁcant: The exergy showed a sensitivity of − 0.49, and in total emissions, this ﬁgure was slightly lower, at − 0.41. By changing the emissions of N 2 O in the soil, the sensitivity of exergy was almost null, and the total emissions were 0.077.


Introduction
Since the beginning of the 20th century, the world's population has increased more than four-fold, increasing from 1.7 to 7 billion inhabitants. At the same time, there was a 100-fold increase in industrial production and a 50-fold increase in fossil fuel consumption, reaching 10,700 MTOE (million tons of oil equivalent) in 2016 [1]. More than one decade before, Graedel and Allenby [2] already pointed out that such consumption, either by the finiteness of natural resources or due to the consequential damage, has shown to be not sustainable in relation to the environment. The old industrial model, focused on production (and productivity), has somehow met such demands, although in an immediate manner. A new model has emerged, which is economically and socially feasible, as well as environmentally sustainable.
In order to allow for the development of this new industrial model, it is necessary to quantify the environmental impacts related to any processes, enabling a comparison of possible corrections, in which exergy can be of use. The use of the exergetic method, which, according to Kotas [3], is a relatively new alternative thermodynamic analysis based on the concept of exergy, is loosely defined as a universal measure of the work potential or quality of different forms of energy in relation

Materials and Methods
The methodology includes the characterization of a case study in the bioenergy sector, namely, sugarcane-based ethanol production, coupled with a thermodynamic variable that represents the exergy extraction (CEENE-cumulative exergy extraction from the natural environment), with GHG emissions.
Exergy and emissions, as well as the relationship between them, were treated in two scenarios. Scenario 1, which is the reference system (RS), describes a sugarcane, typical for ethanol in the southwestern region of Brazil, calculated from the inventory reported by Macedo et al. [11], whose data were adjusted and/or updated, as presented in Table 1, with appropriate references. Scenario 2 describes the changes due to the reusing of the industrial waste filtercake, which, after drying and pyrolysis, is transformed into biochar. These changes affect some items of the inventory, detailed in Table 1, thus altering both exergy consumption and GHG emissions.
Exergy is divided into four parts: potential, which is associated with a potential difference (X po ) that can be gravitational or even magnetic; kinetics (X ki ), which is associated with motion (speed); Physics (X ph ), which occurs as a result of differences in temperature and pressure, with respect to the environment and chemistry (X ch ), which are due to differences in the chemical composition between the initial and final states [3].
Exergy requires a reference level for its final state, which is considered dead when modified matter is under an equilibrium with the reference environment. As for physical exergy, both references, temperature and pressure, are described, with a zero index in the form of T 0 and P 0 , standardized in 298 K and 1 bar, respectively [8]. Similarly, enthalpy (H 0 ) and entropy (S 0 ) refer to this same condition of T 0 and P 0 . Thus, the physical exergy can be calculated according to Equation (1).
where X: physical exergy (J); H(J): enthalpy of the flow in the process under T, P conditions; H 0 (J): enthalpies of the flow under the T 0 , P 0 reference conditions; S (JK −1 ): entropy of the flow in the process under conditions T, P; S 0 (JK −1 ): entropy of the flow in the process under T 0 , P 0 conditions; T 0 (K): reference temperature; and P 0 : reference pressure (100 KPa). Concerning chemical exergy, the level or reference state for most elements is where they present themselves in the most oxidized or chlorinated form. When in balance with the atmosphere, hydrosphere and crust, these elements, under such conditions, are taken as the reference state, with their state, concentration and composition defined, that is, the exergy value of a chemical compound is determined by comparing it with the standard condition of this compound in the terrestrial biosphere [12][13][14]. The general formulation proposed by Szargut [12] is presented in Equations (2) and (3).
where the exergy part of the molar component i, b ch,i is calculated by Equation (3).
where b 0 ch,i : molarpart of i component exergy under the T 0 , P 0 (Jmol −1 ) reference conditions; R: the ideal gas constant (Jmol −1 K −1 ); x i : molar composition of i component; and ϕ i : fugacity coefficient of i component.
According to Maes and Passel [14], there are two possible approaches for including solar irradiation in the exergy calculation. The first one is to include it only indirectly, only accounting for the biomass provided to the industrial process [15], and a second one, more comprehensive and inclusive approach, accounting for all solar irradiation necessary to produce biomass, which is the expanded concept of the cumulative exergy consumption (CExC) of Szargut et al. [8] and the entitled CEENE methodology (cumulative exergy extraction from the natural environment), assessing all kinds of resources that are excluded from the natural ecosystem, proposed by Dewulf et al. [9]. The value of CEENE, normalized here by megajoules of the ethanol unit, is described in Equation (4).
where CEENE: the process (MJ/MJ ethanol); ∑ExC i (in): sum of the accumulated exergy of n inputs (MJ/MJ ethanol); and ∑ExC i (out): sum of the accumulated exergy of n outputs (MJ/MJ ethanol). GHG emissions were also normalized by megajoules of ethanol, and this was accounted for in the general Equation (5).
where ε GHG : greenhouse gas emissions of the process (g CO 2-equivalent /MJ ethanol); and ε i : sum of the emissions of n inputs/products. The ratio of these quantities, which is between CEENE and GHG emissions, in the considered scenarios, is already presented in Equation (6).
where κ j : CEENE and emissions ratio in scenario j (in MJ/kg CO 2-equivalent ); CEENE j : Cumulative Exergy Extraction from the Natural Environment in scenario j (MJ/MJ ethanol); and ε GHG-j : GHG emissions in scenario j (kg CO 2-equivalent /MJ ethanol). In general, the use of this same ratio, when altering n possible input variables, can be presented as a dependent Y function, as shown in Equation (7).
where Y: variable that represents the rate of CEENE and GHG emissions (MJ/g CO 2-equivalent ); X 1 , X 2 . . . X n : input variables; CEENE: (MJ/MJ ethanol); and ε: GHG emissions (g CO 2-equivalent /MJ ethanol). When addressing the sensitivity issue of variables, in the study, it was briefly indicated that, when comparing the output dependent variable behavior, it is usual to change the input variable value from 0 to 200% of its base or reference value. The behavior between the input and output variation is what makes it possible to quantify the sensitivity of this variable, both in relation to the other variable and the other points of the scanning [16,17].
Regarding sensitivity, when there is an explicit algebraic relationship, describing the relationship between an independent X i variable and a dependent Y variable, this analysis is easy to carry out, and in this case, the Φ i sensitivity coefficient for an independent variable in the study can be calculated on the basis of an X i variable. For small changes in the input parameters, the partial derivative can be approximated as a finite difference, where non-linearity is neglected, and the partial derivative can be approximated in numerical form as Equation (8).
where Φ i : sensitivity in X i ; % ∆Y: percentage variation of Y magnitude; and % ∆X i : percentage variation of X i magnitude.

Bioenergy Systems, GHG and Exergy
The simplified flowchart of the process in the sugar ethanol sector, which was the basis of this study, is presented in Figure 1. Within the limits of this system, composed of both agricultural and industrial parts, are catalogued as both emissions and exergy flows.
It is highlighted that the main exergy flow entering the system is sunlight, and the main outputs are ethanol and electricity, which are exported to the plant. The latter was produced from steam from bagasse and straw. All emissions throughout the process were catalogued, from the capture of carbon in photosynthesis, turning it into sugarcane and straw, until its return to the atmosphere during the cycle.
The limits of the system under study are the limits of the agro-industry itself, from agricultural biomass production to industrial processing. In sucroenergy, the vinasse, ash and filtercake that were returned to the soil, without any treatment, have not been accounted for. The exergy outcomes of CEENE and the associated GHG emissions, by ethanol megajoules, will be grouped at the same levels, as proposed by Szargut et al. [8], as follows:  Level 1: Direct inputs, considered beyond the Sun, fuel consumption and the electric energy occasionally taken from an external source;  Level 2: Indirect energy inputs, where it is added to the first level of the energy/emissions associated with the production of other inputs, either by the crop or through the manufacturing process;  Level 3: The energy/emissions, related to the production and maintenance of equipment and installations, added to the previous levels. In this study, level 3 was used.
Emissions during biofuel production, with which this study is concerned, at the levels already mentioned, as well as the carbon cycle in sucroenergy, are represented in a summarized way in Figure 2. At the entrance, the carbon photosynthesis capture is represented, and at the exit, both the biotic and nonbiotic carbon are represented. Non-biotic carbon-related emissions, as well as those that occur during the use of soil, represented by nitrogenous fertilizer's N2O, are the agents considered in this comparative LCA, that is, the biotic cycle is considered to be neutral.
Exergy inputs are grouped in the same way as the emissions, and their calculations were carried out in accordance with the aforementioned CEENE Equation (4). In Figure 3, the exergy inputs/products to the sucroenergy sector are briefly presented, both the high (energy drinks) and low exergy. In the case of fertilizers, accumulated exergy is calculated by the sum of the exergy consumed during production and by the chemical exergy present in the fertilizer itself (input), as well as the exergy of nitrous oxide emissions. The main exergy output of the bioenergy system is the ethanol itself, but there is also a large contribution of bagasse biomass and surplus straw to the production of electricity since the internal usage is irreversible in the process.  The exergy outcomes of CEENE and the associated GHG emissions, by ethanol megajoules, will be grouped at the same levels, as proposed by Szargut et al. [8], as follows: • Level 1: Direct inputs, considered beyond the Sun, fuel consumption and the electric energy occasionally taken from an external source; • Level 2: Indirect energy inputs, where it is added to the first level of the energy/emissions associated with the production of other inputs, either by the crop or through the manufacturing process; • Level 3: The energy/emissions, related to the production and maintenance of equipment and installations, added to the previous levels. In this study, level 3 was used.
Emissions during biofuel production, with which this study is concerned, at the levels already mentioned, as well as the carbon cycle in sucroenergy, are represented in a summarized way in Figure 2. At the entrance, the carbon photosynthesis capture is represented, and at the exit, both the biotic and non-biotic carbon are represented. Non-biotic carbon-related emissions, as well as those that occur during the use of soil, represented by nitrogenous fertilizer's N 2 O, are the agents considered in this comparative LCA, that is, the biotic cycle is considered to be neutral.
Exergy inputs are grouped in the same way as the emissions, and their calculations were carried out in accordance with the aforementioned CEENE Equation (4). In Figure 3, the exergy inputs/products to the sucroenergy sector are briefly presented, both the high (energy drinks) and low exergy. In the case of fertilizers, accumulated exergy is calculated by the sum of the exergy consumed during production and by the chemical exergy present in the fertilizer itself (input), as well as the exergy of nitrous oxide emissions. The main exergy output of the bioenergy system is the ethanol itself, but there is also a large contribution of bagasse biomass and surplus straw to the production of electricity since the internal usage is irreversible in the process. The exergy outcomes of CEENE and the associated GHG emissions, by ethanol megajoules, will be grouped at the same levels, as proposed by Szargut et al. [8], as follows:  Level 1: Direct inputs, considered beyond the Sun, fuel consumption and the electric energy occasionally taken from an external source;  Level 2: Indirect energy inputs, where it is added to the first level of the energy/emissions associated with the production of other inputs, either by the crop or through the manufacturing process;  Level 3: The energy/emissions, related to the production and maintenance of equipment and installations, added to the previous levels. In this study, level 3 was used.
Emissions during biofuel production, with which this study is concerned, at the levels already mentioned, as well as the carbon cycle in sucroenergy, are represented in a summarized way in Figure 2. At the entrance, the carbon photosynthesis capture is represented, and at the exit, both the biotic and nonbiotic carbon are represented. Non-biotic carbon-related emissions, as well as those that occur during the use of soil, represented by nitrogenous fertilizer's N2O, are the agents considered in this comparative LCA, that is, the biotic cycle is considered to be neutral.
Exergy inputs are grouped in the same way as the emissions, and their calculations were carried out in accordance with the aforementioned CEENE Equation (4). In Figure 3, the exergy inputs/products to the sucroenergy sector are briefly presented, both the high (energy drinks) and low exergy. In the case of fertilizers, accumulated exergy is calculated by the sum of the exergy consumed during production and by the chemical exergy present in the fertilizer itself (input), as well as the exergy of nitrous oxide emissions. The main exergy output of the bioenergy system is the ethanol itself, but there is also a large contribution of bagasse biomass and surplus straw to the production of electricity since the internal usage is irreversible in the process.   All the calculation presented in the course of this work were carried out, in vector form, in the numerical software, Scilab [18], from the values listed in Table 1. From the dataset of the table from Macedo et al. [11], the values of the consumed energies and the associated emissions are quantified integrally in fossil form, where the values of the energy consumption and constants in the form of LHV (Lower Heat Value) are approximated as the consumption of exergy, considered as energy inputs (fuels), and these values are close.  All the calculation presented in the course of this work were carried out, in vector form, in the numerical software, Scilab [18], from the values listed in Table 1. From the dataset of the table from Macedo et al. [11], the values of the consumed energies and the associated emissions are quantified integrally in fossil form, where the values of the energy consumption and constants in the form of LHV (Lower Heat Value) are approximated as the consumption of exergy, considered as energy inputs (fuels), and these values are close.  The production of the biochar occurs after the collection of the filtercake residue inside the distillery when it has been dried, and pyrolysis is transformed into biochar. As already mentioned, originally (scenario 1), this residue returns to the agricultural soil without any treatment and, unlike what occurs in scenario 2, already in the form of biochar. As shown in Figure 4, inside the distillery, the collection of the filtercake occurs during the separation of impurities from the sugarcane juice through the process known as juice clarification. The production of the biochar occurs after the collection of the filtercake residue inside the distillery when it has been dried, and pyrolysis is transformed into biochar. As already mentioned, originally (scenario 1), this residue returns to the agricultural soil without any treatment and, unlike what occurs in scenario 2, already in the form of biochar. As shown in Figure 4, inside the distillery, the collection of the filtercake occurs during the separation of impurities from the sugarcane juice through the process known as juice clarification. According to Jeffery et al. [28], reporting an increase in crop yield, and through a statistical metaanalysis by compiling results from several studies, it was possible to assess the relationship between biochar and productivity. The outcomes showed an important benefit, with an average increase of 10%.
Similarly, Cayuelaa et al. [29], also using the meta-analysis methodology, obtained a mean value in the reduction of nitrous oxide emissions of 23% for low concentrations of biochar in the soil.
Then, for scenario 2, the calculation of the emissions and the energies, present in all processes in the agro-industrial system, will be considered:  An increase in productivity of 10%;  A decrease in nitrous oxide emissions of 23%.  According to Jeffery et al. [28], reporting an increase in crop yield, and through a statistical meta-analysis by compiling results from several studies, it was possible to assess the relationship between biochar and productivity. The outcomes showed an important benefit, with an average increase of 10%.

Accounting for Emissions and Exergy
Similarly, Cayuelaa et al. [29], also using the meta-analysis methodology, obtained a mean value in the reduction of nitrous oxide emissions of 23% for low concentrations of biochar in the soil.
Then, for scenario 2, the calculation of the emissions and the energies, present in all processes in the agro-industrial system, will be considered:

•
An increase in productivity of 10%; • A decrease in nitrous oxide emissions of 23%. Table 2 presents all the inventoried processes in the agro-industrial system under study. For the final accounting of CEENE and emissions, all items are not necessarily considered, as detailed below. Accounting for the emissions, we initially considered the CO 2-equivalent values related to each ton of cane produced. Later, we considered the productivity of sugarcane suitable for each hectare planted, and finally, we considered the production and exergy of ethanol, leading to the emissions in megajoules functional unit of ethanol produced (g CO 2-equivalent /MJ ethanol), as shown in Equation (9).

Accounting for Emissions and Exergy
where ε GHG : GHG emissions in grams of CO 2-equivalent per megajoule of ethanol produced; Prod_cane: productivity of sugarcane per hectare (Mg/ha); Prod_ethanol: ethanol productivity per hectare (kg/ha); and Ex ethanol : exergy of hydrated ethanol (MJ/kg). The life cycle approach to the entire system was used in this study, with GHG emissions in kilograms of CO 2-equivalent per megajoules of ethanol (ε), which are included in As an example, will further detail how was obtain the item ε 15 during the use of nitrogen fertilizer, by inserting emissions per ton of sugarcane in Equation (10). The aggregation of the emissions per megajoule of ethanol occurred in two ways: • Scenario 1, without the use of biochar, as shown in Equation (11); • Scenario 2, with the use of biochar, as shown in Equation (12).
CEENE initially accounted for each ton of cane produced, then the productivity of sugarcane for each hectare planted was considered, and finally, the productivity of ethanol was considered, which resulted in the irreversibilities in megajoules per megajoule of ethanol produced (MJ/MJ ethanol), as shown in Equation (13).
where CEENE: Cumulative Exergy Extraction from the Natural Environment per megajoule of ethanol produced (MJ/MJ); Prod_cane: productivity of sugarcane in tons of sugarcane per hectare (Mg/ha); Prod_ethanol: ethanol productivity per hectare (kg/ha); and Ex ethanol : exergy of hydrated ethanol (MJ/kg) As in the above-mentioned case of the emissions of nitrogen fertilization (ε 15 ) in Equation (10), we also individually calculated the exergy consumed (Ex 15 ) during the use of nitrogen fertilizer, initially per ton of sugarcane, as in Equation (14), and then, we transported this value to Equation (13). This item, as well as several other ones, has the characteristic of needing to add the energy consumed in the production of compost to the chemical exergy of compost because these two resources are being degraded.
Ex (15) where m N / SCT: nitrogen fertilizer mass per sugarcane ton (kg/SCT); EN: exergy consumed in the production of nitrogen fertilizer (MJ/kg); and X N : chemical exergy of nitrogen fertilizer (MJ/kg).
Being an exergy input, the exergy of the sun (Ex 1 ), whose value must also be transferred to Equation (13), shall be calculated as shown in Equation (15). X sun = (10, 000)(365)(0.9)(AR) (15) where Xsun: absorbed solar exergy per hectare in the cycle of one year (MJ/ha); 10,000 m 2 = 1 ha; 365 days a year; AR: average annual radiation for one day in the southeast region of Brazil (16 MJ/m 2 day); and "0.9": exergy quality factor for sunlight. In accordance with Equations (16) and (17) and inserting data from the exergy values in Table 2, the total extraction of exergy of the system was calculated for both scenarios, where this extraction (irreversibilities) is in megajoule units per kilogram of ethanol.

The Biochar Changing Scenarios
With the use of biochar, both the productivity of sugarcane and the nitrous oxide are altered, and the impacts of these variations are initially analyzed punctually and, later, more generally, when the same variables are scanned. Figure 5 presents the exergy (CEENE) and GHG emissions behavior of scenarios 1 and 2. Exergy is calculated according to Equations (16) and (17), in MJ/MJ ethanol, and GHG emissions, in g CO 2-equivalent /MJ ethanol, is calculated according to Equations (11) and (12). The increase in productivity had the greatest contribution to the decrease of exergy extraction, and the influence of the nitrous oxide emissions is much less significant. In addition to the two causes are GHG emissions due to increased productivity and lower emissions of nitrous oxide itself, where increased productivity is more sensitive.
Energies 2019, 12 FOR PEER REVIEW 10 With the use of biochar, both the productivity of sugarcane and the nitrous oxide are altered, and the impacts of these variations are initially analyzed punctually and, later, more generally, when the same variables are scanned. Figure 5 presents the exergy (CEENE) and GHG emissions behavior of scenarios 1 and 2. Exergy is calculated according to Equations (16) and (17), in MJ / MJ ethanol, and GHG emissions, in g CO2-equivalent/MJ ethanol, is calculated according to Equations (11) and (12). The increase in productivity had the greatest contribution to the decrease of exergy extraction, and the influence of the nitrous oxide emissions is much less significant. In addition to the two causes are GHG emissions due to increased productivity and lower emissions of nitrous oxide itself, where increased productivity is more sensitive.

Scanning Production and N2O Emissions
From the inventory database inserted in the Scilab numerical software, Figure 7 presents individually, through two functions, the relationship between CEENE and GHG emissions quantities. In the first function, the productivity of sugarcane varied between 50 and 150% (the other values are the same as those in RS, the reference system). The second function had the same percentages when varying the emissions of nitrous oxide. This figure shows a great sensitivity of  Figure 6 represents the values alteration concerning exergy-emissions rates for scenarios 1 and 2, according to Equation (6), in MJ/g CO 2-equivalent . It is observed that the exergy-emissions relationship remained relatively constant due to the simultaneous occurrence of increased productivity and the decrease in the nitrous oxide emissions with the use of biochar.
Energies 2019, 12 FOR PEER REVIEW 10 With the use of biochar, both the productivity of sugarcane and the nitrous oxide are altered, and the impacts of these variations are initially analyzed punctually and, later, more generally, when the same variables are scanned. Figure 5 presents the exergy (CEENE) and GHG emissions behavior of scenarios 1 and 2. Exergy is calculated according to Equations (16) and (17), in MJ / MJ ethanol, and GHG emissions, in g CO2-equivalent/MJ ethanol, is calculated according to Equations (11) and (12). The increase in productivity had the greatest contribution to the decrease of exergy extraction, and the influence of the nitrous oxide emissions is much less significant. In addition to the two causes are GHG emissions due to increased productivity and lower emissions of nitrous oxide itself, where increased productivity is more sensitive.

Scanning Production and N2O Emissions
From the inventory database inserted in the Scilab numerical software, Figure 7 presents individually, through two functions, the relationship between CEENE and GHG emissions quantities. In the first function, the productivity of sugarcane varied between 50 and 150% (the other values are the same as those in RS, the reference system). The second function had the same percentages when varying the emissions of nitrous oxide. This figure shows a great sensitivity of

Scanning Production and N 2 O Emissions
From the inventory database inserted in the Scilab numerical software, Figure 7 presents individually, through two functions, the relationship between CEENE and GHG emissions quantities. In the first function, the productivity of sugarcane varied between 50 and 150% (the other values are the same as those in RS, the reference system). The second function had the same percentages when varying the emissions of nitrous oxide. This figure shows a great sensitivity of CEENE and GHG emissions due to the sweep of productivity and the almost null sensitivity in the sweeping of nitrous oxide emissions in soil.
Energies 2019, 12 FOR PEER REVIEW 11 CEENE and GHG emissions due to the sweep of productivity and the almost null sensitivity in the sweeping of nitrous oxide emissions in soil. Using the algorithm in the Scilab numerical programming software, Equations (6), (11) and (16) are applied in an amplified way, when scanning from 50% to 150% around the base value, which is the reference value (X-axis). For example, for sugarcane productivity (Prod_cane), as shown in Equations (9) and (13), this value is 70 tons per hectare, i.e., scanning takes place from 35 and 105 tonnes per hectare, impacting (or not) on the outcomes of the Equation. The notable function behavior around the base value represents a greater or lesser sensitivity of the input.   Figure 9 presents the values of GHG emissions "ε" and the variables between 50% and 150% for the inputs or products: (a) sugarcane productivity; and (b) nitrous oxide emissions. The sensitivity of Using the algorithm in the Scilab numerical programming software, Equations (6), (11) and (16) are applied in an amplified way, when scanning from 50% to 150% around the base value, which is the reference value (X-axis). For example, for sugarcane productivity (Prod_cane), as shown in Equations (9) and (13), this value is 70 tons per hectare, i.e., scanning takes place from 35 and 105 tonnes per hectare, impacting (or not) on the outcomes of the Equation. The notable function behavior around the base value represents a greater or lesser sensitivity of the input. CEENE and GHG emissions due to the sweep of productivity and the almost null sensitivity in the sweeping of nitrous oxide emissions in soil. Using the algorithm in the Scilab numerical programming software, Equations (6), (11) and (16) are applied in an amplified way, when scanning from 50% to 150% around the base value, which is the reference value (X-axis). For example, for sugarcane productivity (Prod_cane), as shown in Equations (9) and (13), this value is 70 tons per hectare, i.e., scanning takes place from 35 and 105 tonnes per hectare, impacting (or not) on the outcomes of the Equation. The notable function behavior around the base value represents a greater or lesser sensitivity of the input.   Figure 9 presents the values of GHG emissions "ε" and the variables between 50% and 150% for the inputs or products: (a) sugarcane productivity; and (b) nitrous oxide emissions. The sensitivity of  Figure 9 presents the values of GHG emissions "ε" and the variables between 50% and 150% for the inputs or products: (a) sugarcane productivity; and (b) nitrous oxide emissions. The sensitivity of cane productivity from the base value is higher, according to Equation (8), at −0.41, and the nitrous oxide emissions are smaller, at 0.07.

The GHG Behavior
Energies 2019, 12 FOR PEER REVIEW 12 cane productivity from the base value is higher, according to Equation (8), at −0.41, and the nitrous oxide emissions are smaller, at 0.07. Figure 9. GHG behavior of the variable cane productivity and N2O emissions (50%-150%).

The Behavior ratio of CEENE: GHG and the Linear Function Result
The alterations of the ratio CEENE: emissions have been treated individually using Scilab numerical programming. The original vector and linear Y function (X1) are presented, in Figure 10, in the same way as the original vector and linear Y function (X2) were presented in Figure 11, and X1 for sugarcane productivity and X2 for nitrous oxide emissions are also presented. These functions presented a fairly linear behavior, which made it possible to represent them with a linear function.
A linear function simply represents the vector outcomes from which they originated. A simple linear regression (RLS), a and b values as well as the standard deviation were addressed with the reglin command [a, b, sig], where mean(Y(Xi)), representing the average and the coefficient of determination (R 2 ), which were also calculated using Scilab. Equation (18), whose outcomes are presented in Figure 10 and Equation (19), is shown in Figure 11. After completing such equations, the aforementioned statistical data are also presented.

The Behavior ratio of CEENE: GHG and the Linear Function Result
The alterations of the ratio CEENE: emissions have been treated individually using Scilab numerical programming. The original vector and linear Y function (X 1 ) are presented, in Figure 10, in the same way as the original vector and linear Y function (X 2 ) were presented in Figure 11, and X 1 for sugarcane productivity and X 2 for nitrous oxide emissions are also presented. These functions presented a fairly linear behavior, which made it possible to represent them with a linear function.
A linear function simply represents the vector outcomes from which they originated. A simple linear regression (RLS), a and b values as well as the standard deviation were addressed with the reglin command [a, b, sig], where mean(Y(X i )), representing the average and the coefficient of determination (R 2 ), which were also calculated using Scilab. Equation (18), whose outcomes are presented in Figure 10 and Equation (19), is shown in Figure 11. After completing such equations, the aforementioned statistical data are also presented. where Y(X 1 ): the function that relates the exergy extraction with GHG emissions in accordance with productivity of sugarcane (MJ/g CO 2-equivalent ); and X 1 : productivity of sugarcane (SCT/ha).  Figure 10. Y function behavior for sugarcane productivity (50% to 150%). Figure 11. Y function behavior for nitrogen fertilization (50% to 150%).

Discussion
Even at a very high value, the exergy of the sun is accounted for, being the great input of exergy in the system. The magnitude of this value results in large irreversibilities in the system, but CEENE's behavior is proportional to this input since the emission values are independent of the greater or lesser radiation considered, which makes the CEENE: emissions ratio also proportional. As for the sensitivities, their magnitudes are not altered, but the behavior of the other variables that depend on the solar radiation, whose values are also proportional, is just as important.
As for the behavior shown in Figure 10 and Figure 11, it is important to note that all the variation occurs around the reference value of 33.1 MJ/g CO2-equivalent, which is Y's value at the center of the X axis. The variables that represent the productivity and emissions of N2O, when undergoing changes due to the use of biochar, were specifically considered in these graphs. To a greater or lesser extent, the following variables were sensitive to the decrease of the value of Yi with the increase of Xi: ethanol productivity; the production of electricity; diesel transport/field; and the CO2 generated in the fermentation and yield of sugar cane. In nitrogen fertilization and nitrous oxide emissions, the values of Yi remained practically constant.
For each megajoule of ethanol produced, the increase in productivity had a greater influence on the decrease in the exergy extraction, compared to the decrease in emissions, and Figure 10 presents this relation, which results in better use of resources when increasing productivity.
As for the irreversibilities in the numerator of Equation 7 and the emissions in the denominator, it should be noted that smaller values of this ratio are not necessarily indicative of a smaller number

Discussion
Even at a very high value, the exergy of the sun is accounted for, being the great input of exergy in the system. The magnitude of this value results in large irreversibilities in the system, but CEENE's behavior is proportional to this input since the emission values are independent of the greater or lesser radiation considered, which makes the CEENE: emissions ratio also proportional. As for the sensitivities, their magnitudes are not altered, but the behavior of the other variables that depend on the solar radiation, whose values are also proportional, is just as important.
As for the behavior shown in Figure 10 and Figure 11, it is important to note that all the variation occurs around the reference value of 33.1 MJ/g CO2-equivalent, which is Y's value at the center of the X axis. The variables that represent the productivity and emissions of N2O, when undergoing changes due to the use of biochar, were specifically considered in these graphs. To a greater or lesser extent, the following variables were sensitive to the decrease of the value of Yi with the increase of Xi: ethanol productivity; the production of electricity; diesel transport/field; and the CO2 generated in the fermentation and yield of sugar cane. In nitrogen fertilization and nitrous oxide emissions, the values of Yi remained practically constant.
For each megajoule of ethanol produced, the increase in productivity had a greater influence on the decrease in the exergy extraction, compared to the decrease in emissions, and Figure 10 presents this relation, which results in better use of resources when increasing productivity.
As for the irreversibilities in the numerator of Equation 7 and the emissions in the denominator, it should be noted that smaller values of this ratio are not necessarily indicative of a smaller number Figure 11. Y function behavior for nitrogen fertilization (50% to 150%).

Discussion
Even at a very high value, the exergy of the sun is accounted for, being the great input of exergy in the system. The magnitude of this value results in large irreversibilities in the system, but CEENE's behavior is proportional to this input since the emission values are independent of the greater or lesser radiation considered, which makes the CEENE: emissions ratio also proportional. As for the sensitivities, their magnitudes are not altered, but the behavior of the other variables that depend on the solar radiation, whose values are also proportional, is just as important.
As for the behavior shown in Figures 10 and 11, it is important to note that all the variation occurs around the reference value of 33.1 MJ/g CO 2-equivalent , which is Y's value at the center of the X axis. The variables that represent the productivity and emissions of N 2 O, when undergoing changes due to the use of biochar, were specifically considered in these graphs. To a greater or lesser extent, the following variables were sensitive to the decrease of the value of Yi with the increase of Xi: ethanol productivity; the production of electricity; diesel transport/field; and the CO 2 generated in the fermentation and yield of sugar cane. In nitrogen fertilization and nitrous oxide emissions, the values of Yi remained practically constant.
For each megajoule of ethanol produced, the increase in productivity had a greater influence on the decrease in the exergy extraction, compared to the decrease in emissions, and Figure 10 presents this relation, which results in better use of resources when increasing productivity.
As for the irreversibilities in the numerator of Equation (7) and the emissions in the denominator, it should be noted that smaller values of this ratio are not necessarily indicative of a smaller number of irreversibilities generated, since the increase in the denominator also gives this result. From this, it can be understood that: the improvement of the system, with regard to the aspect of the efficient use of resources, depends on whether variable Xi is related to exergy or the emissions. An increase in the variable Xi, which represents nitrous oxide and ammonia, leads to apparently controversial results since they cause the decrease of the Y value. Due to the numerator of the Equation, the exergy consumption remains practically constant during the sweep from 50%-150% in X due to its low contribution to the exergy impacts.
The Y function presented in Figures 10 and 11 depends on the exergy in the numerator and the emissions in the denominator, as already mentioned. The Equations presented already coupled these quantities as a function of X 1 and X 2 .
The representation of the best use of carbon in the Y function may be the effect of the decrease in the value of its numerator (exergy) and denominator (emissions). Thus, more nitrous oxide represents a poor use of carbon, as well as a low productivity of sugarcane. Then, from Figure 10, the best use of carbon occurs with increasing productivity and, from Figure 11, with decreasing N 2 O emissions.
As previously shown in Figure 9, the emissions value is already in grams of CO 2-equivalent per megajoules of ethanol, and this already represents the efficient use of carbon. Figures 10 and 11 indicate that, from a base value of the reference system, the efficiency increased with the increase of X 1 or decrease of X 2 .

Conclusions
In addition to presenting biochar as an industrial waste that has positive results when returned to the soil, with an increased productivity and reduction of nitrous oxide emissions in the soil, this study showed that, using the system, it is possible to perform a coupling between CEENE (cumulative exergy extraction from the natural environment) and air pollution in the form of greenhouse gases, where the ratio of these quantities is represented by the Y function. In the reference system (RS), the exergy:emissions ratio was calculated, with a value of 33.1 MJ/g CO 2-equivalent , which is changed as the variables present in the system under study are changed. By performing the irreversibility coupling, in the form of CEENE, with GHG emissions, initially through a function Y (X 1 ) and later Y (X 2 ), where X 1 represents the productivity and X 2 , the emissions of nitrous oxide, the behavior of this ratio was almost linear, making it possible, through linear regression, to write an equation in the form, ax + b, which produced a = −0.075, b = 38.3 and R 2 = 0.997 for productivity (Mg SCT/ha year) and a = −3.71, b = 38.2 and R 2 = 0.995 for nitrous oxide emissions (kg N 2 O/ha year).
In order to advise as to the better use of carbon, it is not necessary to calculate the CEENE: emissions ratio, since emissions are already accounted for in grams of CO 2-equivalent per megajoule of ethanol produced, where, in the graph dealing with the emissions behavior from the scans of sugarcane productivity and N 2 O emissions in the soil, it is quite clear that less carbon is released into the atmosphere for each megajoule of ethanol produced by increasing productivity and decreasing N 2 O emissions.
However, with some considerations, it is possible, through the equations that have been made, to achieve the same result using CEENE-emissions coupling. As a general rule, it can be said that, during the scanning of variables with greater sensitivity in the numerator of the Y function, lower values of Y indicate a better use of the carbon, and already for variables with a greater sensitivity in the denominator of the said function, lower values of Y indicate a worse use of carbon. In other words, the question is whether the variable under study is more sensitive to CENNE or emissions. For cane productivity, the CEENE sensitivity is 20% higher than the emissions (CEENE predominates, −0.49 versus −0.41). For emissions, it is almost zero, and for CEENE and emissions, it is 0.077, which predominates.
Thus, in a bioenergetic system, it was possible to correlate different metrics when making a coupling between the CEENE and the GHG emissions, whose ratio in the reference system was 33.1 MJ/g CO 2-equivalent . Besides relating these quantities, this ratio, to a greater or lesser extent, is altered both by the change in the exergy extraction and by the GHG emissions in the system, or even κ j exergy: emissions ratio scenario j (MJ/kg CO 2-equivalent ) Y(X i ) exergy: emissions ratio variable X i (MJ/g CO 2-equivalent ) Φ i coefficient of sensitivity in X i ϕi: fugacity coefficient of the i component