Next Article in Journal
Probabilistic Hotspot Prediction Model Based on Bayesian Inference Using Precipitation, Relative Dry Spells, ENSO and IOD
Next Article in Special Issue
PM10 Resuspension of Road Dust in Different Types of Parking Lots: Emissions, Chemical Characterisation and Ecotoxicity
Previous Article in Journal
Features of Multiannual Air Temperature Variability in Poland (1951–2021)
Previous Article in Special Issue
Biofiltration of Waste Gas Containing Cyclohexanol, Cyclohexanone and Butanol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GHG Global Emission Prediction of Synthetic N Fertilizers Using Expectile Regression Techniques

by
Kaoutar Benghzial
1,†,
Hind Raki
1,2,*,†,
Sami Bamansour
1,†,
Mouad Elhamdi
1,2,†,
Yahya Aalaila
1,2,† and
Diego H. Peluffo-Ordóñez
1,2,†
1
SDAS Research Group, Ben Guerir 43150, Morocco
2
Institute of Science, Technology & Innovation (IST&I), Modeling Simulation and Data Analysis, Mohammed VI Polytechnic University, Benguerir 43150, Morocco
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Atmosphere 2023, 14(2), 283; https://doi.org/10.3390/atmos14020283
Submission received: 19 December 2022 / Revised: 23 January 2023 / Accepted: 26 January 2023 / Published: 31 January 2023
(This article belongs to the Special Issue Feature Papers in Air Quality)

Abstract

:
Agriculture accounts for a large percentage of nitrous oxide ( N 2 O ) emissions, mainly due to the misapplication of nitrogen-based fertilizers, leading to an increase in the greenhouse gas (GHG) footprint. These emissions are of a direct nature, released straight into the atmosphere through nitrification and denitrification, or of an indirect nature, mainly through nitrate leaching, runoff, and N 2 O volatilization processes. N 2 O emissions are largely ascribed to the agricultural sector, which represents a threat to sustainability and food production, subsequent to the radical contribution to climate change. In this connection, it is crucial to unveil the relationship between synthetic N fertilizer global use and N 2 O emissions. To this end, we worked on a dataset drawn from a recent study, which estimates direct and indirect N 2 O emissions according to each country, by the Intergovernmental Panel on Climate Change (IPCC) guidelines. Machine learning tools are considered great explainable techniques when dealing with air quality problems. Hence, our work focuses on expectile regression (ER) based-approaches to predict N 2 O emissions based on N fertilizer use. In contrast to classical linear regression (LR), this method allows for heteroscedasticity and omits a parametric specification of the underlying distribution. ER provides a complete picture of the target variable’s distribution, especially when the tails are of interest, or in dealing with heavy-tailed distributions. In this work, we applied expectile regression and the kernel expectile regression estimator (KERE) to predict direct and indirect N 2 O emissions. The results outline both the flexibility and competitiveness of ER-based techniques in regard to the state-of-the-art regression approaches.

1. Introduction

The food and drink industry is the largest producing sector globally, and due to the increased consumer demand for processed food products, it has led to consequential impacts on health and the environment [1]. Biotic and abiotic components of the environment are targeted by air pollution, which is considered one of our era’s greatest scourges. Every substance, solid, liquid, or gas, if being produced in higher concentrations while reducing the quality of our environment, is defined as a pollutant [2]. According to the World Health Organization (WHO), 99% of humans are breathing air that exceeds WHO guideline limits and contains high levels of pollutants, while low- and middle-income countries are subject to the highest exposures. Air quality is closely linked to the Earth’s climate and ecosystem, and is known to be the single largest environmental health risk factor globally. Many of the drivers of air pollution are also sources of greenhouse gas (GHG) emissions [3].
Particulate matter (PM) can be formed directly in the atmosphere by physicochemical reactions between pollutants already present in the atmosphere or can be directly emitted from anthropogenic activities and natural sources to the atmosphere. The United States Environmental Protection Agency defined PM as a term for particles, whose penetration depends on their diminutive size, ranging from particles with diameters of 10 μm (μm) or smaller, called P M 10 , and extremely fine particles with diameters that are generally 2.5 μm (μm) and smaller P M 2.5 [2]. PM is most likely condensed in cities and industrialized areas, as it is geographically shown in Figure 1. PM concentrations levels are represented by color grading, where the intense green represents the highest mean values of these concentrations (μg/m 3 ) The increase in atmospheric P M 2.5 concentration, air movement patterns, and exposure of populations, result in health and economic effects. Food system emissions alone account for about 22.4% of global mortality due to degraded air quality and 1.4% of global crop production losses [4]. A recent study [5] in the United States estimated that 4300 cases annually of premature mortality happen due to maize production. In fact, higher mortality rates were observed within the top five maize-producing states (Iowa, Illinois, Nebraska, Minnesota, and Indiana). Moreover, increased concentrations of P M 2.5 are driven by emissions of ammonia NH3, which result from nitrogen (N) fertilizer use [5].
Industrial facilities, such as power stations, refineries, petrochemicals, chemical and fertilizer industries, and metallurgical, and other industrial plants, are major sources of pollutants emissions. GHG emissions from the agricultural sector increased by 10.1% from 1990 to 2018 and accounted for 9.9% of total US greenhouse gas emissions [6]. In fact, agricultural N 2 O emissions are projected to continue to rise [7]. Agricultural crop production, including farms and the supply chains that produce the chemical and energy inputs, contribute majorly to the emissions of GHGs, which include carbon dioxide ( C O 2 ), nitrous oxide ( N 2 O ), methane ( C H 4 ), and black carbon [4]. Nitrous oxide being one of the most impacting GHG, was chosen in this study. It is estimated that N 2 O emissions in the US account for approximately 75% of total emissions. The truth is that the increased value brought about by nitrogen-based fertilizer applications is outweighed by the expenses of environmental nitrogen pollution, such as the eutrophication of rivers, loss of biodiversity, global warming, and stratospheric ozone depletion, even though N is a limiting component for agricultural production [8].
With the growing rate of big data evolution and its complexity, various prediction methods based on machine learning technologies have been developed for air quality problems [9,10]. Multiple linear regression (MLR) is one of the most popular tools capable of incorporating complex non-linear relationships between the concentration of air pollutants and meteorological variables [11].
This work is based on Greenhouse Gas Emissions from Global Production and the use of Synthetic Nitrogen Fertilizers in the Agriculture dataset, from the Figshare repository, from the year 2018. The authors also used this dataset in a recently published work [12], where they estimated GHG emissions due to synthetic N fertilizer manufacture, transportation, and field use in agricultural systems. Most studies have tackled the GHG emission problems; while integrating ML tools basically focus on C O 2 or C H 4 emissions [13], very few papers are based on ( N 2 O ) emissions. In fact, this gas is 300 times more harmful to the climate than ( C O 2 ) and steadily increases in the atmosphere, with agriculture being the largest contributor, and nitrogen the most used synthetic fertilizer [14].
In this study, we propose two expectile-based regression approaches, namely, expectile regression (ER) and the kernel expectile regression estimator (KERE). Due to their flexibility in application, heavy-tailed distributions and outliers are of interest. In this context, and based on the fact that only a few countries are considered agricultural-producing countries, there is a concentration of information in the tail of the distribution. We used expectile-based regression to take advantage of the parameterized nature, which allows for modeling different aspects of the distribution rather than the simple mean.
The rest of this manuscript is structured as follows: Section 2 briefly outlines the works and studies related to our research. Section 3 and Section 4 investigate and analyze the dataset. Thereof, Section 5 assembles the results and the discussion. Finally, Section 6 portrays concluding remarks and states future works.

2. Background and Related Works

In the last decades, climate changes are becoming more frequent with longer shifts in temperature and weather patterns [15,16]. According to the United Nations reports, human activities are the main drivers of climate change, such as the agricultural production system, which significantly contributes to GHG emissions, such as nitrous oxide N 2 O [17,18].

2.1. Nitrogen Based-Fertilizers Use and Climate Change

The increasing usage of nitrogen-based fertilizers has a significant influence on the worldwide agricultural soils’ rapid emissions of N 2 O [19]. It represents a major contributing factor to the current rise in the global average temperature and a significant impact on agricultural productivity [20]. On the other hand, in the agricultural sector, there is significant N 2 O discharge, mostly as a result of the usage of nitrogen-based fertilizers as a crop productivity booster [20], or through their manufacturing processes [21]. According to the Food and Agriculture Organization(FAO), the usage of synthetic N-based fertilizers is expected to grow to 50% by 2050 [12]. In 2019, N 2 O levels in the atmosphere have grown by more than 20% from 270 parts per billion in 1750 to 332 parts per billion [14]. Global concerns surround the steadily rising levels of N 2 O in the atmosphere [21], as a strong greenhouse gas with a long lifetime (around 121 years) [7]. This gas is responsible for rising temperatures altering the patterns of the world’s weather, which can cause severe abiotic phenomena, such as droughts and heavy rains.

2.1.1. Nitrogen-Based Fertilizers

Nitrogen supply and crop demand must be synchronized in order to maintain optimal plant growth and reduce environmental losses, with regard to other outside factors [22]. The 4R nutrient stewardship community suggests that nitrogen-based fertilizers must reduce GHG emissions while minimizing any harm to the air quality. The 4Rs concept involves applying fertilizers using the right source, at the right rate, at the right time, and at the right place for each crop [23]. The main sources of nitrogen-based fertilizers are typically divided into ammonium-based ( N H 4 + ) and nitrate-based ( N O 3 ) [24]: ammonium nitrate (AN) ( 33 , 5 % N ) ; Urea ( 46 % N ) ; ammonium sulfate (AS) ( 21 % N ) ; calcium ammonium nitrate (CAN) ( 27 % N ) [25]. Other nitrogen–phosphate (NP) and nitrogen–phosphate–potassium (NPK) fertilizers, which are widely used, such as mono-ammonium phosphate (MAP) with 11 % of N, di-ammonium phosphate (DAP) with 18 % of N and ammonium sulfate ( 21 % of N) [22].

2.1.2. Nitrous Oxide: A Side Effect of Nitrogen-Based Fertilizers Use

N is one of the major macro-nutrients other than P and K, which are crucial for plant growth and development. Under climate-change conditions, enhancing agricultural output and quality while minimizing environmental losses has become an arduous challenge [26]; hence, leading to the usage of synthetic N fertilizers, in order to improve crop production by providing the plant with the necessary nutrient. N is particularly a limiting factor for crop growth, it has an important role within plant cells, which ensures better results in terms of yield and quality as well [8]. Having multiple synergistic interactions with other nutrients, and in accordance with Liebig’s law of the minimum, nitrogen is the first probable nutrient deficiency to occur [27]. The agricultural sector is once again the most impacted by the temperature rise brought on by N 2 O emissions. Research has indicated that it is responsible for 15% of greenhouse gas emissions [20]. According to FAO estimates for the year 2019, nitrogen-based fertilizer field application contributed to 8.3% of world GHG emissions, while fertilizer manufacturing accounted for 0.7% [12]. In other words, the fertilizer supply chain (from the manufacturing to the field application) distributes the whole nitrogen cycle and, thereby, efforts must be taken to lessen the quantity of N 2 O emitted into the atmosphere [28].

2.1.3. Nitrogen Dynamics in Soil and N 2 O Emissions Pathways

Soils naturally release N 2 O due to the microbiological processes that are a part of the nitrogen cycle [29]. As shown in Figure 2, nitrogen-based fertilizers when applied to the soil, some of it is taken up by the plant, while the rest is converted by soil microorganisms or lost through volatilization/leaching. The greatest sources of indirect N 2 O emissions come from agricultural N O 3 leaching and runoff, which account for around 30% of the nitrogen lost from agricultural soils [28]. On the other hand, denitrification and nitrification through microorganisms involve the direct emissions pathway of N 2 O . Under aerobic conditions, the conversion of N H 4 + to N O 3 occurs under the name of nitrification, while the denitrification process takes place involving the conversion of N O 3 to N 2 , with N 2 O acting as an intermediate product [30].

2.2. Machine Learning Tools for Air Quality Predictions: Nitrous Oxide Emissions

When it comes to air quality prediction, tools, such as artificial intelligence (AI) and machine learning (ML), have gained remarkable attention due to their potential to help in GHG emission predictions. Table 1 summarizes various recent and convenient works that integrated ML tools for N 2 O emissions.

3. Dataset

The dataset was acquired from the Figshare repository, entitled Greenhouse Gas Emissions from Global Production and Use of Nitrogen-based Fertilizer in Agriculture. This dataset contains information used to analyze the impact of synthetic N fertilizer production, transportation, and use on global anthropogenic greenhouse gas emissions. The main sources are Agrifootprint 6.0 (AFP6), Food and Agriculture Organization Corporate Statistical Database (FAOSTAT), and International Fertilizer Association (IFASTAT).

3.1. Data Description

The data were derived from a recent study [12], where they estimated GHG emissions per capita terms for the year 2018, due to multiple synthetic N fertilizer parameters, namely, manufacture, transportation, and field use in agricultural systems. This field-level dataset available on N 2 O soil emissions was estimated based on global N 2 O direct emission factors (EFs). From the literature, EFs for indirect N 2 O soil emissions and N fertilizer manufacturing and transportation were extracted. The original dataset is organized in a .xlsx format file, with eight tables. We extracted data from tables 2 and 5 as shown in Figure 3.
To calculate the amount of global GHG emissions generated through the processes of synthetic N fertilizer production, transportation, and agricultural use, they followed the IPCC guidelines for national GHG inventories. They accordingly calculated the emissions at the country level on the basis of the activity data they collected and the appropriate EFs. The IPCC guidelines divide the agricultural N 2 O source into three categories: direct emissions from agricultural land, emissions from animal waste management systems, and indirect emissions associated with volatilized/leached N that are removed in biomass or otherwise exported from agricultural land. Each is estimated to contribute by one-third, while indirect emissions account for approximately two-thirds.

3.2. Data Pre-Processing

In this section, we present the transformations and aggregations applied to the data, in order to explore the relationship between N synthetic fertilizer sources, and N 2 O direct and indirect emissions. This will allow us to set up the necessary structure to achieve the objective of this work.
  • Scaling and missing values: Features with percentage values are calculated to indicate the amount of each nitrogen source used. This will allow for adequate and informative exploratory analysis in later stages. In addition, entries with missing values were removed.
  • Standardizing: Standardization is crucial in feature extraction; especially when features vary in large units or when reported in different measures. To this end, the data are standardized using the following formula: ( X μ ) / σ , where X is the explanatory variable, μ and σ are the mean and standard deviations, respectively. In addition, the total amount of nitrogen used per country was removed due to the direct correlation with the N sources. Moreover, the corresponding features of indirect emissions generated by leaching/run-off and volatilization were added to indicate indirect emissions generally.
  • Exploratory data analysis: To explore the data, two initial investigations are performed. First, the linear regression approach is applied to standardized data to extract feature importance. Second, we investigate the most important features by a distribution graph. In Figure 4, it is apparent that most of the features exude initial signs of heavy-tailed behaviors, which are due to the disparities in agricultural activities (per country).

4. Data Analysis

Generally, ML techniques can be categorized into supervised, unsupervised, semi-supervised, and reinforcement learning. Supervised ML approaches deal with a particular case of problems where each data sample is paired to a label. In particular, regression-type approaches generate an underlying function that provides a real value for each data sample. In this work, the goal is to explore the relationship between N 2 O gas emissions (direct and indirect), the quantity of applied nitrogen, and synthetic fertilizer sources. To this end, we present the following workflow, based on expectile regression models as our learning approach.

4.1. Expectile-Based Regression

4.1.1. Linear Expectile Regression

Linear expectile regression (ER) was first introduced by Newey and Powell in risk measurement [36]. This approach can be defined as the generalization of conditional expectation to model the relationship between a dependent variable and the covariates [37]. Multiple studies were conducted to explore ER performance, particularly when dealing with heavy-tailed distribution [38]. Although the ER provides a complete picture of the data [39], its statistical properties are under-explored in contrast with other methods, such as linear regression and quantile regression [40,41].
Given Y a random variable, the expectile of level τ denoted by μ τ is defined as follows:
μ τ ( Y ) = min τ R E ( ϕ τ ( Y τ ) ) ,
where
ϕ τ ( u ) = τ u 2 , u < 0 ( 1 τ ) u 2 u 0 .
ϕ τ is the asymmetric least square (ALS) loss function that assigns weights τ and 1 τ to positive and negative deviations, respectively. Figure 5 provides example curves of ϕ τ i and expectile values μ τ i with respect to different expectile levels τ i , respectively.
Let us suppose that we have n samples ( y i , x i ) , where x i = ( 1 , x i , 1 , , x i , p ) T are the covariates. The expectiles defined in Equation (1) are used to set up the expectile regression, which assumes a linear model of the following form:
μ τ ( y i | x i ) = x i T β τ
The estimated coefficients β ^ τ can be obtained by minimizing the empirical loss function:
1 n i = 1 n ϕ τ ( y i x i T β τ ) .
To solve the optimization problem in Equation (3), we suggest the following Algorithm 1 based on using iterative reweighted least squares (IRLS) [42]:   
Algorithm 1: ALS for estimating ER coefficients.
Input: Measured dataset ( x i , y i ) i = 1 m .
  1. Initialize β τ , 0 .
  2. Use the empirical loss function in Equation (3);
  3. Update the coefficient using the algorithm of IRLS [42],
Output: The coefficients’ estimates β ^ .
   The linear expectile regression has shown great performance in contrast with classical approaches of regression [36]. However, we may encounter more complex datasets for which ER might be too restrictive in terms of errors. To this end, researchers have developed more flexible methods namely: expectile regression with boosting (ER-boosting) [43] and non-parametric estimator of conditional expectiles based on local linear polynomials with a one-dimensional covariate [44].

4.1.2. Kernel Expectile Regression Estimator (KERE)

In this work, we adopt a recent flexible method introduced in a modern study [45], based on exploiting the properties of the reproducing kernel Hilbert spaces (RKHS) [46]. Let H K denote a Hilbert space generated by a predefined kernel K. Given n samples ( y i , x i ) , the kernel expectile regression estimator is derived from the following optimization problem:
( f ^ n ( x ) , α ^ 0 ) = arg min f H K , α 0 R i = 1 n ϕ τ ( y i α 0 f ( x i ) ) + λ f H K 2 ,
where f spans the Hilbert space, f H K 2 is the norm of f in H K , α 0 is the intercept and λ is the regularization parameter.
Although Equation (4) lies in an infinite dimensional space, the dimension of this formulation is reduced by using the representer theorem and the reproducing property [47]. Thus, the optimization parameter f in Equation (4) and its RKHS norm are expressed as follows:
f ( x ) = k = 1 n α k K ( x k , x ) , f H K 2 = i = 1 n j = 1 n α i α j K ( x i , x j ) .
where K is the kernel function and α k R .
To this end, Equation (4) can be reformulated as follows:
{ α ^ i } = argmin α i R i = 1 n ϕ τ y i α 0 j = 1 n α j K ( x i , x j ) + λ i = 1 n j = 1 n α i α j K ( x i , x j ) .
A compact formulation of Equation (6) using matrix notations is introduced as follows:
α ^ = argmin α R n + 1 i = 1 n ϕ τ ( y i K i α ) + λ α T K α ,
where
K i = ( 1 , K ( x i , x 1 ) , , K ( x i , x n ) ) K 0 = 0 0 n × 1 0 n × 1 K
K = K ( x 1 , x 1 ) K ( x 1 , x n ) K ( x 2 , x 1 ) K ( x 2 , x n ) K ( x n , x 1 ) K ( x n , x n )
The proposed algorithm to solve Equation (7) relies on using maximization–minimization (MM) approaches [48]. The key idea is to find a surrogate function through the Taylor expansion that majorizes the objective function. Optimizing this surrogate function will either improve the value of the objective function or leave it unchanged.
Using the MM approach to solve Equation (7) yields the following formulation for iteratively updating α = α 0 , α 1 , , α n .
α ( t + 1 ) = α ( t ) + K u 1 λ K 0 α ( t ) + 1 2 i 1 n ϕ τ ( r i ( t ) ) K i ,
where K u 1 is the inverse matrix of K u defined as follows:
K u = m a x ( 1 τ , τ ) n 1 T K K 1 K K + λ m a x ( 1 τ , τ ) K .
Algorithm 2 summarizes the steps to reach the KERE estimates y ^ of the output y .
Algorithm 2: Kernel expectile regression estimator.
Input: Dataset ( x i , y i ) i = 1 m , kernel function K ( . , . ) , tolerance ϵ , maximum iterations i m a x .
 1: Calculate the kernel matrix K = ( K ( x i , x j ) ) i , j .
 2: Initialize α ( 0 ) , r i ( 0 ) and t 0 .
 3: While r i ( t ) ϵ and t i m a x :
   ● Calculate updated residue r i ( t ) = y i K i α ( t )
   ● Update α ( t ) based on Equation (10).
   ● t t + 1
 4: Calculate the output estimator y ^ i = α 0 ( t ) + j = 1 n α j ( t ) K x j , x i
Output: The vector of estimates y ^ .

4.2. Experimental Setup

In this work, we evaluated some of the expectile-based approaches, namely, expectile regression (ER) [36] and kernel expectile regression estimator (KERE) [45] on the GHG emission dataset. As detailed in Section 4, both ER and KERE depend on the chosen expectile level w. To depict the performance relative to the expectile level, we construct multiple models of ER and KERE using multiple expectile levels spanning the following values:
τ 0.01 , 0.05 , 0.1 , 0.2 , 0.25 , 0.5 , 0.7 , 0.75 , 0.8 , 0.95 .
In addition, KERE models require the kernel function to be selected. Although various kernels are available for use, we choose the well-known radial basis family (RBF) kernel defined in Equation (11).
K x , y = e 1 σ 2 | | x y | | 2 ,
where σ stands for the bandwidth.
To select the best hyperparameters for KERE models, we perform two dimensional 8-fold cross-validation to select the optimal hyperparameters ( λ , σ ) , where λ and σ stand for the regularization parameter and the RBF kernel bandwidth, respectively. Moreover, the maximum number of iterations i m a x and the tolerance value ϵ is fixed to 4000 and 10 6 , respectively.
Furthermore, we select the KERE model corresponding to the expectile level of interest being w = 0.7 to conduct a benchmark comparison with state-of-the-art regression models. Table 2 summarizes the selected regression methods to be compared with KERE, as well as their respective characteristics. We also conduct an 8-fold cross-validation to tune each model’s hyperparameter.
The evaluation process is two-fold. First, we compare ER and KERE models using a customized error metric, namely, mean absolute deviation (MAD) (w) defined below. This type of error reflects the model fit with an emphasis on the tails by assigning weights τ and 1 − τ to positive and negative deviations, respectively. Second, a subset of the KERE models is selected to be compared with state-of-the-art regression approaches using mean absolute error (MAE) and root-mean-square error (RMSE).
  • Mean absolute deviation MAD (w):
    MAD ( w ) = 1 n i = 1 n ϕ w y y ^
  • Mean absolute error (MAE):
    MAE = 1 n i = 1 n | y y ^ |
  • Root-mean-square error (RMSE):
    RMSE = 1 n i = 1 n 1 n y y ^ 2

4.3. Benchmark Methods

In order to assess the performance of the kernel expectile regression estimator on the proposed dataset, we compare it to twelve other benchmark regression approaches. As detailed in Table 2, we use support vector regression, lasso, light gradient boosting machine, random forest, K-neighbor, extra trees, AdaBoost, gradient boosting, decision tree, Huber, multilayer perceptron, and ridge regressors. Table 2 summarizes the techniques considered as well as their hyperparameters to be tuned using K-fold cross-validation.

4.4. Computational Software

The computational software for this study was written using both RStudio and Python. RStudio was used for comparing kernel expectile regression estimator and linear expectile regression. Both approaches were implemented using “KERE” and “Expectreg” libraries [49], respectively. On the other hand, Python was used to compute the comparison with regression benchmark approaches detailed in Table 2. The benchmark comparison was conducted using the PyCaret library (version 3.0.0rc4), specifically PyCaret.regression module.
All of the aforementioned regression techniques were computed using an 8-fold cross-validation to tune the corresponding hyperparameters. The advantage of the PyCaret library in Python is the agility of its classes, particularly Compare_models in setting up the proper framework for a fair comparison.

5. Results and Discussion

Firstly, we report the results of two expectile-based approaches, namely kernel expectile regression estimator (KERE) and expectile regression (ER). The two methods were applied to predict both N 2 O direct and indirect emissions. First, KERE and ER are compared on both the training and testing phases. We evaluate the models primarily using mean absolute deviation error (MAD), which varies with respect to the expectile level. In addition, we report the mean absolute error (MAE), root-mean-square error (RMSE), and R². Second, we chose the KERE model corresponding to the expectile level w = 0.7 for comparison with state-of-the-art regression approaches using R², MAE, and RMSE.
We report the results for both the training and testing of KERE and ER regarding the direct emissions in Table 3 and Table 4, respectively. It is noticeable that KERE outperforms ER in all reported metrics, namely MAD, RMSE, and MAE. This is because KERE is able to depict non-linear behavior utilizing the kernel trick. In addition, ER approach reports an increasing MAD error as the expectile levels increase, whereas KERE stays relatively stable as the expectile level increases. This is reflected by the mean absolute deviation of ER (0.8) being 0.308 compared to KERE (0.7) being 0.041. In addition, it is apparent that R² values drop significantly between training and testing for both KERE and ER which highlights the failure to explain the direct emissions variance.
One KERE model corresponding to expectile level w = 0.7 is selected from Table 4 to be compared with the benchmark approaches summarized in Table 1. Table 5 summarizes the performance of the benchmark regressors in contrast with KERE models, where R 2 , MAE and RMSE are reported.
It is shown that the KERE R-squared values are slightly low but significantly better than the rest of the regression approaches (18%). We can argue that the usage of fewer fertilizers implies less factors implicated in the agricultural processes. For example, in small countries where agriculture is not the main activity, the source of N-based fertilizers does not necessarily explain N 2 O emissions. However, when N-based fertilizers are applied in bigger quantities, the source explains more about the direct emissions.
Secondly, we report the performance analysis of KERE and ER with respect to indirect emissions. Table 6 and Table 7 summarize the performance evaluation of both approaches in the training and testing phases. The MAD, RMSE, and MAE metrics are reported with respect to the various expectile levels. Similar to the models reported previously, KERE models significantly outperform ER models. The MAD evaluation metric of KERE and ER on the testing set is 0.046 and 2.762 corresponding to expectile level w = 0.7 . Whereas KERE and ER report an MSE metric of 0.391 and 3.010 with respect to the same expectile level.
Similar to the previously mentioned results, the KERE model corresponding to the expectile level w = 0.7 is selected to be compared with the benchmark regression methods by reporting MAE, RMSE, and R². As outlined in Table 8, the selected KERE model performed significantly better, especially with regard to the RMSE metric evaluation.
Reporting the training and testing results for both direct and indirect emissions outlined the out-performance of kernel expectile regressions estimator approach to its counterpart Expectile Regression, especially with regard to MAD evaluation error. Furthermore, the KERE model corresponding to expectile level w = 0.7 was selected to focus on the data closer to the tail of the data. The latter corresponds to the range of medium to large countries, showing that such a model, in addition to being flexible, performs better than all other considered benchmark regression approaches.
The KERE technique is indeed an explainable approach, allowing us to explore the relationship between synthetic N fertilizers use and global N 2 O emissions. Reducing N rates is not the main factor for reducing GHG emissions, mainly nitrous oxide. The adoption of lower N rates underestimates N 2 O emissions [50], as there are many factors at stake, such as the management of fertilizer applications while enhancing N use efficiency (NUE). N-based fertilizers vary depending on the N-form they contain, either ammonium N H 4 + or nitrate N O 3 . In fact, N H 4 + is the starting point by which the soil microorganisms perform the nitrification process to form N O 3 , from which other soil microorganisms convert it to N 2 gas through the process of denitrification, while emitting the N 2 O gas during the whole process and this is the direct pathway of the N 2 O emissions. Referring to Figure 2 which highlights N pathways, and based on a modern study [51], it appears that there is a positive correlation between soil moisture content and cumulative N 2 O emissions. When water content is high in the soil, it was suggested that the required conditions for denitrification are met, leading to higher N O 3 concentrations in the soils, providing N substrate for the production of N 2 O . In fact, floods and rain may have an impact on GHG emissions, increasing precipitation may enhance soil N 2 O emissions [52].
A first on-farm study [50] was conducted to report N 2 O response to multiple fertilizer rates on production-scale fields. They observed linear and nonlinear increases in N 2 O depending on the study locations and year. However, the nonlinear exponential response models best represented the overall N 2 O response to N fertilizer across all site years. A more recent study [53] also demonstrated the nonlinear relationship between the application of nitrogen-based fertilizers and N 2 O emissions, explaining how it changes depending on the meteorological circumstances, while the correlation between the N 2 O emissions and the N-fertilizer rate used remains unclear. In another study, it has been proven that monitoring nitrogen application alone is not capable of stimulating N 2 O emissions as much as the combination of nitrogen addition and rainfall reduction [54]. Nonetheless, nitrogen fertilization is an external factor; other management practices generate N 2 O as a side effect of their applications, such as irrigation or tillage practices or even the crop type used. Cropping systems have an impact on soil quality and soil GHG emissions. In a recent study [55], the importance of combining winter cover crop cultivation for single cropping systems with reduced N fertilizer application was investigated. This work supports our results; the cropping system and N rate application impact the GHG emissions and N 2 O direct emissions. In comparison with traditional cotton cultivation, where N 2 O and other GHG emissions were increased, it appears that both cover cropping with reduced N helped to mitigate soil GHG emissions. Furthermore, the geographical disposition may have an impact as well. For example, countries, such as China, the USA, Canada, India, and Brazil are known for their enormous agricultural production, hence, diverse soil characteristics, diverse climates, and various crop types are greatly responsible for the increase of GHG emissions. From another point of view, and in the case of China, rice represents the most economically important crop, where rice paddies are an excellent environment for the biological activities of nitrification and denitrification processes which have been accelerated especially in flooded soils, leading to enormous N 2 O production as described by [56]. Other N-fertilization management techniques have an impact on N 2 O emissions, with an increase in N fertilizer’s use, adopting the 4R Nutrient Stewardship has a significant potential to lower N 2 O emissions [57].
The results also show that the source of N synthetic fertilizers does not contribute as much to indirect emissions of N 2 0 . When it comes to indirect emissions, research has proven that the leaching and runoff of nitrate from the application of synthetic N fertilizers is a substantial indirect source of N 2 O emissions from groundwater. This indicates that indirect emissions are not related to the source of N fertilizer itself, but depend mainly on the N content in agricultural soils that could be lost through leaching/Runoff [58].

6. Conclusions and Future Works

The food supply chain is a salient contributor to global emissions of air pollutants. These major air pollutant compounds, such as N 2 O , are emitted by different stages of the food system, from food production, processing, packaging, transport, retail, consumption, and disposal. This work used an explainable approach in order to explore the relationship between synthetic N fertilizers and global N 2 O emissions.
According to our results and the findings within the literature, we conclude, two major points:
  • Using the kernel expectile regression estimator approach is highly suitable when dealing with air quality related data.
  • Integrating additional external factors is highly recommended, for more accurate results and in order to interpret our outputs.
Based upon the results from the KERE method, the source of N-based fertilizers is not capable of explaining N 2 O emissions on their own. On the one hand, the disparity in the degree of N 2 O emissions between countries, where many national GHG inventories employ the IPCC’s first-tier method, which uses a set of global default emission factors (EF) to estimate the N 2 O emissions depending on varying geography. On the other hand, N-fertilizer sources and quantities are not the only major contributors in nitrous oxide emissions, in fact, a significant portion of greenhouse gas emissions is now caused by farming activities either directly or indirectly, leading to a decay in air quality. Agricultural practices are of great importance; hence, we propose that reducing the usage of synthetic N fertilizers is not the most suitable solution to reduce N 2 O emissions. However, there is a need to monitor nitrogen-based fertilizer usage, while keeping an eye on the whole farming management practices, such as conventional tillage, which is a major contributor in N 2 O emissions, and no-tillage can decrease the N 2 O emissions, in the presence of different fertilizer treatments. Furthermore, our research suggests the necessity for a more specified and detailed estimation of N 2 O emissions, based on each country for a longer period of time. In order to establish a more effective study, we also encourage agribusinesses, organizations, and farmers to empower data availability. Bridging the gap between research and industry is a challenging step. Gaining insight from the data is very important; hence, increasing the usage of explainable machine learning tools is of great significance toward integrating interpretability for air quality models.  

Author Contributions

Conceptualization, K.B., H.R. and D.H.P.-O.; methodology, Y.A., S.B. and M.E.; software, Y.A., S.B. and M.E.; validation, Y.A.; formal analysis, H.R. and K.B.; investigation, Y.A. and D.H.P.-O.; resources, H.R. and D.H.P.-O.; data curation, H.R.; writing—original draft preparation, H.R. and K.B.; writing—review and editing, H.R., D.H.P.-O. and Y.A.; visualization, Y.A. and M.E.; supervision, D.H.P.-O.; project administration, D.H.P.-O.; funding acquisition, D.H.P.-O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work is supported by the SDAS Research Group accessed on 29 January 2023 (https://sdas-group.com). Modeling Simulation and Data Analysis (MSDA) Research Program accessed on 29 January 2023 (https://msda.um6p.ma/home).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GHGgreenhouse gas
KEREkernel expectile regression estimator
ERexpectile regression
MLmachine learning
Nnitrogen
MADmean absolute deviation
MAEmean absolute error
RMSEroot-mean-square error

References

  1. Watson, N.J.; Bowler, A.L.; Rady, A.; Fisher, O.J.; Simeone, A.; Escrig, J.; Woolley, E.; Adedeji, A.A. Intelligent sensors for sustainable food and drink manufacturing. Front. Sustain. Food Syst. 2021, 5, 408. [Google Scholar] [CrossRef]
  2. Manisalidis, I.; Stavropoulou, E.; Stavropoulos, A.; Bezirtzoglou, E. Environmental and health impacts of air pollution: A review. Front. Public Health 2020, 8, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. World Health Organization. Report-WHO Air Quality Database. 2022. Available online: https://www.who.int/data/gho/data/themes/air-pollution/who-air-quality-database (accessed on 15 December 2022).
  4. Crippa, M.; Solazzo, E.; Guizzardi, D.; Van Dingenen, R.; Leip, A. Air pollutant emissions from global food systems are responsible for environmental impacts, crop losses and mortality. Nat. Food 2022, 3, 942–956. [Google Scholar] [CrossRef]
  5. Hill, J.; Goodkind, A.; Tessum, C.; Thakrar, S.; Tilman, D.; Polasky, S.; Smith, T.; Hunt, N.; Mullins, K.; Clark, M.; et al. Air-quality-related health damages of maize. Nat. Sustain. 2019, 2, 397–403. [Google Scholar] [CrossRef]
  6. Alhashim, R.; Anandhi, A. Global Warming and Toxicity Impacts: Peanuts in Georgia, USA Using Life Cycle Assessment. Sustainability 2022, 14, 3671. [Google Scholar] [CrossRef]
  7. Tian, L.; Cai, Y.; Akiyama, H. A review of indirect N2O emission factors from agricultural nitrogen leaching and runoff to update of the default IPCC values. Environ. Pollut. 2019, 245, 300–306. [Google Scholar] [CrossRef]
  8. Rütting, T.; Aronsson, H.; Delin, S. Efficient use of nitrogen in agriculture. Nutr. Cycl. Agroecosyst. 2018, 110, 1–5. [Google Scholar] [CrossRef] [Green Version]
  9. Rosero-Montalvo, P.D.; Caraguay-Procel, J.A.; Jaramillo, E.D.; Michilena-Calderón, J.M.; Umaquinga-Criollo, A.C.; Mediavilla-Valverde, M.; Ruiz, M.A.; Beltrán, L.A.; Peluffo, D.H. Air quality monitoring intelligent system using machine learning techniques. In Proceedings of the 2018 International Conference on Information Systems and Computer Science (INCISCOS), Quito, Ecuador, 13–15 November 2018; pp. 75–80. [Google Scholar]
  10. Rosero-Montalvo, P.D.; López-Batista, V.F.; Arciniega-Rocha, R.; Peluffo-Ordóñez, D.H. Air Pollution Monitoring Using WSN Nodes with Machine Learning Techniques: A Case Study. Log. J. IGPL 2022, 30, 599–610. [Google Scholar] [CrossRef]
  11. Bekkar, A.; Hssina, B.; Douzi, S.; Douzi, K. Air-pollution prediction in smart city, deep learning approach. J. Big Data 2021, 8, 1–21. [Google Scholar] [CrossRef]
  12. Menegat, S.; Ledo, A.; Tirado, R. Greenhouse gas emissions from global production and use of nitrogen synthetic fertilisers in agriculture. Sci. Rep. 2022, 12, 14490. [Google Scholar] [CrossRef]
  13. Sánchez-Pozo, N.N.; Trilles-Oliver, S.; Solé-Ribalta, A.; Lorente-Leyva, L.L.; Mayorca-Torres, D.; Peluffo-Ordóñez, D.H. Algorithms Air Quality Estimation: A Comparative Study of Stochastic and Heuristic Predictive Models. In Proceedings of the International Conference on Hybrid Artificial Intelligence Systems, Bilbao, Spain, 22–24 September 2021; pp. 293–304. [Google Scholar]
  14. Tian, H.; Xu, R.; Canadell, J.G.; Thompson, R.L.; Winiwarter, W.; Suntharalingam, P.; Davidson, E.A.; Ciais, P.; Jackson, R.B.; Janssens-Maenhout, G.; et al. A comprehensive quantification of global nitrous oxide sources and sinks. Nature 2020, 586, 248–256. [Google Scholar] [CrossRef] [PubMed]
  15. Ummenhofer, C.C.; Meehl, G.A. Extreme weather and climate events with ecological relevance: A review. Philos. Trans. R. Soc. B Biol. Sci. 2017, 372, 20160135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Elahi, E.; Khalid, Z.; Tauni, M.Z.; Zhang, H.; Lirong, X. Extreme weather events risk to crop-production and the adaptation of innovative management strategies to mitigate the risk: A retrospective survey of rural Punjab, Pakistan. Technovation 2022, 117, 102255. [Google Scholar] [CrossRef]
  17. Tubiello, F.N.; Rosenzweig, C.; Conchedda, G.; Karl, K.; Gütschow, J.; Xueyao, P.; Obli-Laryea, G.; Wanner, N.; Qiu, S.Y.; De Barros, J.; et al. Greenhouse gas emissions from food systems: Building the evidence base. Environ. Res. Lett. 2021, 16, 065007. [Google Scholar] [CrossRef]
  18. Elahi, E.; Khalid, Z. Estimating smart energy inputs packages using hybrid optimisation technique to mitigate environmental emissions of commercial fish farms. Appl. Energy 2022, 326, 119602. [Google Scholar] [CrossRef]
  19. Wang, J.; Smith, P.; Hergoualc’h, K.; Zou, J. Direct N2O emissions from global tea plantations and mitigation potential by climate-smart practices. Resour. Conserv. Recycl. 2022, 185, 106501. [Google Scholar] [CrossRef]
  20. Malhi, G.S.; Kaur, M.; Kaushik, P. Impact of climate change on agriculture and its mitigation strategies: A review. Sustainability 2021, 13, 1318. [Google Scholar] [CrossRef]
  21. Chai, R.; Ye, X.; Ma, C.; Wang, Q.; Tu, R.; Zhang, L.; Gao, H. Greenhouse gas emissions from synthetic nitrogen manufacture and fertilization for main upland crops in China. Carbon Balance Manag. 2019, 14, 20. [Google Scholar] [CrossRef] [Green Version]
  22. Fageria, N.K.; Baligar, V. Enhancing nitrogen use efficiency in crop plants. Adv. Agron. 2005, 88, 97–185. [Google Scholar]
  23. Bruulsema, T. Nutrient Stewardship: Taking 4R Further. Crop. Soils 2022, 55, 34–40. [Google Scholar] [CrossRef]
  24. Cao, P.; Lu, C.; Yu, Z. Historical nitrogen fertilizer use in agricultural ecosystems of the contiguous United States during 1850–2015: Application rate, timing, and fertilizer types. Earth Syst. Sci. Data 2018, 10, 969–984. [Google Scholar] [CrossRef]
  25. Demirer, T.; Röck-Okuyucu, B.; Özer, I. Effect of different types and doses of nitrogen fertilizers on yield and quality characteristics of mushrooms (Agaricus bisporus (Lange) Sing) cultivated on wheat straw compost. J. Agric. Rural Dev. Trop. Subtrop. (JARTS) 2005, 106, 71–77. [Google Scholar]
  26. Mridha, K.; Hasan, S.M.A. Artificial Intelligence (AI) for Agricultural Sector. In Proceedings of the 2021 International Conference on Control, Automation, Power and Signal Processing (CAPS), Jabalpur, India, 10–12 December 2021; pp. 1–6. [Google Scholar]
  27. Lal, R.; Stewart, B.A. Soil Nitrogen Uses and Environmental Impacts; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  28. Ramzan, S.; Rasool, T.; Bhat, R.A.; Ahmad, P.; Ashraf, I.; Rashid, N.; Mir, I.A. Agricultural soils a trigger to nitrous oxide: A persuasive greenhouse gas and its management. Environ. Monit. Assess. 2020, 192, 436. [Google Scholar] [CrossRef] [PubMed]
  29. Carbonell-Bojollo, R.M.; Veroz-González, Ó.; González-Sánchez, E.J.; Ordóñez-Fernández, R.; Moreno-García, M.; Repullo-Ruibérriz De Torres, M.A. Soil Management, Irrigation, and Fertilisation Strategies for N2O Emissions Mitigation in Mediterranean Agricultural Systems. Agronomy 2022, 12, 1349. [Google Scholar] [CrossRef]
  30. He, T.; Xie, D.; Ni, J.; Li, Z.; Li, Z. Nitrous oxide produced directly from ammonium, nitrate and nitrite during nitrification and denitrification. J. Hazard. Mater. 2020, 388, 122114. [Google Scholar] [CrossRef]
  31. Marzadri, A.; Amatulli, G.; Tonina, D.; Bellin, A.; Shen, L.Q.; Allen, G.H.; Raymond, P.A. Global riverine nitrous oxide emissions: The role of small streams and large rivers. Sci. Total Environ. 2021, 776, 145148. [Google Scholar] [CrossRef]
  32. Adjuik, T.A.; Davis, S.C. Machine Learning Approach to Simulate Soil CO2 Fluxes under Cropping Systems. Agronomy 2022, 12, 197. [Google Scholar] [CrossRef]
  33. Pan, B.; Lam, S.K.; Wang, E.; Mosier, A.; Chen, D. New approach for predicting nitrification and its fraction of N2O emissions in global terrestrial ecosystems. Environ. Res. Lett. 2021, 16, 034053. [Google Scholar] [CrossRef]
  34. Bastos, L.M.; Rice, C.W.; Tomlinson, P.J.; Mengel, D. Untangling soil-weather drivers of daily N2O emissions and fertilizer management mitigation strategies in no-till corn. Soil Sci. Soc. Am. J. 2021, 85, 1437–1447. [Google Scholar] [CrossRef]
  35. Saha, D.; Basso, B.; Robertson, G.P. Machine learning improves predictions of agricultural nitrous oxide (N2O) emissions from intensively managed cropping systems. Environ. Res. Lett. 2021, 16, 024004. [Google Scholar] [CrossRef]
  36. Newey, W.K.; Powell, J.L. Asymmetric least squares estimation and testing. Econom. J. Econom. Soc. 1987, 55, 819–847. [Google Scholar] [CrossRef]
  37. Waltrup, L.S.; Sobotka, F.; Kneib, T.; Kauermann, G. Expectile and quantile regression—David and Goliath? Stat. Model. 2015, 15, 433–456. [Google Scholar] [CrossRef]
  38. Sigman, K. A primer on heavy-tailed distributions. Queueing Syst. 1999, 33, 261. [Google Scholar] [CrossRef]
  39. Barry, A.D.; Charpentier, A.; Oualkacha, K. Quantile and Expectile Regression for random effects model. HAL 2016, 2016, hal-01421752. [Google Scholar]
  40. Seber, G.A.; Lee, A.J. Linear Regression Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  41. Sobotka, F.; Kauermann, G.; Schulze Waltrup, L.; Kneib, T. On confidence intervals for semiparametric expectile regression. Stat. Comput. 2013, 23, 135–148. [Google Scholar] [CrossRef] [Green Version]
  42. Holland, P.W.; Welsch, R.E. Robust regression using iteratively reweighted least-squares. Commun. Stat. Theory Methods 1977, 6, 813–827. [Google Scholar] [CrossRef]
  43. Yang, Y.; Zou, H. Nonparametric multiple expectile regression via ER-Boost. J. Stat. Comput. Simul. 2015, 85, 1442–1458. [Google Scholar] [CrossRef]
  44. Yao, Q.; Tong, H. Asymmetric least squares regression estimation: A nonparametric approach. J. Nonparametric Stat. 1996, 6, 273–292. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, Y.; Zhang, T.; Zou, H. Flexible expectile regression in reproducing kernel Hilbert spaces. Technometrics 2018, 60, 26–35. [Google Scholar] [CrossRef]
  46. Berlinet, A.; Thomas-Agnan, C. Reproducing Kernel Hilbert Spaces in Probability and Statistics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  47. Wahba, G. Spline Models for Observational Data; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 1990. [Google Scholar]
  48. Lange, K. The MM Algorithm. In Numerical Analysis for Statisticians; Springer: New York, NY, USA, 2010; pp. 189–221. [Google Scholar]
  49. Ooms, J. Magick: Advanced Graphics and Image-Processing in R; R Package Version 2.7.3; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar]
  50. Hoben, J.; Gehl, R.; Millar, N.; Grace, P.; Robertson, G. Nonlinear nitrous oxide (N2O) response to nitrogen fertilizer in on-farm corn crops of the US Midwest. Glob. Chang. Biol. 2011, 17, 1140–1152. [Google Scholar] [CrossRef]
  51. Mosongo, P.S.; Pelster, D.E.; Li, X.; Gaudel, G.; Wang, Y.; Chen, S.; Li, W.; Mburu, D.; Hu, C. Greenhouse Gas Emissions Response to Fertilizer Application and Soil Moisture in Dry Agricultural Uplands of Central Kenya. Atmosphere 2022, 13, 463. [Google Scholar] [CrossRef]
  52. Miller, L.T.; Griffis, T.J.; Erickson, M.D.; Turner, P.A.; Deventer, M.J.; Chen, Z.; Yu, Z.; Venterea, R.T.; Baker, J.M.; Frie, A.L. Response of Nitrous Oxide Emissions to Individual Rain Events and Future Changes in Precipitation. J. Environ. Qual. 2022, 3, 312–324. [Google Scholar] [CrossRef] [PubMed]
  53. Song, X.; Liu, M.; Ju, X.; Gao, B.; Su, F.; Chen, X.; Rees, R.M. Nitrous oxide emissions increase exponentially when optimum nitrogen fertilizer rates are exceeded in the North China Plain. Environ. Sci. Technol. 2018, 52, 12504–12513. [Google Scholar] [CrossRef] [PubMed]
  54. Geng, S.; Chen, Z.; Han, S.; Wang, F.; Zhang, J. Rainfall reduction amplifies the stimulatory effect of nitrogen addition on N2O emissions from a temperate forest soil. Sci. Rep. 2017, 7, 43329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Sun, G.; Zhang, Z.; Xiong, S.; Guo, X.; Han, Y.; Wang, G.; Feng, L.; Lei, Y.; Li, X.; Yang, B.; et al. Mitigating greenhouse gas emissions and ammonia volatilization from cotton fields by integrating cover crops with reduced use of nitrogen fertilizer. Agric. Ecosyst. Environ. 2022, 332, 107946. [Google Scholar] [CrossRef]
  56. Yao, Y.; Zeng, K.; Song, Y. Biological nitrification inhibitor for reducing N2O and NH3 emissions simultaneously under root zone fertilization in a Chinese rice field. Environ. Pollut. 2020, 264, 114821. [Google Scholar] [CrossRef] [PubMed]
  57. Ezui, G.; Haugen-Kozyra, K.; Heaney, D.; Nirjan, L.; Graham, C.; Njoroge, S.; Zingore, S.; Bruulsema, T. Can 4R Practices Limit the Nitrous Oxide Emissions from Increasing Fertilizer Use in Sub-Sahara Africa? Fertilizer Canada. 2022. 4rsolution.org. Available online: https://ureaknowhow.com/wp-content/uploads/2022/04/2022-Fertilizer-Canada-Can-4R-Practices-Limit-the-Nitrous-Oxide-Emissions-from-Increasing-Fertilizer-Use-of-Sub-Saharan-Africa.pdf (accessed on 18 December 2022).
  58. Zhao, L.; Fadong, L.; Zhang, Q.; Wang, J.; Peifang, L.; Tian, C.; Li, X. Influence of land use and change in the proportion of electron donors required for denitrification on N2O in groundwater. Environ. Sci. Pollut. Res. Int. 2021, 28, 17684–17696. [Google Scholar]
Figure 1. Geographic visualization representing the latest global distribution of annual mean levels of P M 2.5 concentrations (μg/m 3 ) in cities (population weighted) according to the WHO Air Quality Database [3].
Figure 1. Geographic visualization representing the latest global distribution of annual mean levels of P M 2.5 concentrations (μg/m 3 ) in cities (population weighted) according to the WHO Air Quality Database [3].
Atmosphere 14 00283 g001
Figure 2. Direct and indirect pathways of N 2 O emissions from agricultural soil after nitrogen-based fertilizer applications.
Figure 2. Direct and indirect pathways of N 2 O emissions from agricultural soil after nitrogen-based fertilizer applications.
Atmosphere 14 00283 g002
Figure 3. Data extraction and retrieving process from the source dataset.
Figure 3. Data extraction and retrieving process from the source dataset.
Atmosphere 14 00283 g003
Figure 4. Density plots of selected nitrogen fertilizer consumption per tonne (t) per country: (a) density plot of ammonium nitrate (AN), (b) density plot of ammonium sulfate (AS), (c) density plot of calcium ammonium nitrate (CAN), (d) density plot of other NS, (e) density plot of urea, (f) density plot of other NP.
Figure 4. Density plots of selected nitrogen fertilizer consumption per tonne (t) per country: (a) density plot of ammonium nitrate (AN), (b) density plot of ammonium sulfate (AS), (c) density plot of calcium ammonium nitrate (CAN), (d) density plot of other NS, (e) density plot of urea, (f) density plot of other NP.
Atmosphere 14 00283 g004
Figure 5. The right figure shows the expectile loss function for different expectile levels τ , and the left figure shows the estimated expectiles for the lognormal distribution.
Figure 5. The right figure shows the expectile loss function for different expectile levels τ , and the left figure shows the estimated expectiles for the lognormal distribution.
Atmosphere 14 00283 g005
Table 1. Machine learning models used recently for predicting N 2 O and GHG emissions.
Table 1. Machine learning models used recently for predicting N 2 O and GHG emissions.
RefTitleML ModelGHG Emissions
[11]Air-pollution prediction in a smart city, deep learning approachHybrid CNN-LSTMCombining historical and meteorological data of pollutants to predict P M 2.5 concentrations in Beijing, China.
[31]Global riverine nitrous oxide emissions: The role of small streams and large riversData-driven random forest ML (RFML)Integrating a data-driven ML model with a physically based up-scaling model to identify N 2 O global emissions, from streams and rivers.
[32]Machine learning approach to simulate soil  C O 2  Fluxes under cropping systemsRegression models: KNN, SVR, and RF.Including ML models for regression purposes in order to predict soil GHG emissions without the biogeochemical expertise.
[33]A new approach for predicting nitrification and its fraction of  N 2 O   emissions in global terrestrial ecosystemsStochastic gradient boosting (SGB)Integrating an ML-SGB model to estimate and predict the global nitrification rate ( R nit ) and N 2 O emissions from the process of nitrification.
[34]Untangling soil weather drivers of daily  N 2 O   emissions and fertilizer management mitigation strategies in no-till cornConditional inference tree (CIT)Using a ML-CIT model to identify the main soil–weather drivers of daily N 2 O hot moments and fertilizer management options to mitigate them.
[35]Machine learning improves predictions of agricultural nitrous oxide ( N 2 O ) emissions from intensively managed cropping systemsRandom forestCoupling an ML-RF model with a cropping system model to predict daily soil N 2 O emissions.
Table 2. The benchmark regression-based approaches used to predict direct and indirect Nitrogen emissions.
Table 2. The benchmark regression-based approaches used to predict direct and indirect Nitrogen emissions.
TechniqueHyperparameters
Support vector machine (SVM)Kernel function, its parameters, and regularization C.
Lasso regressionRegularization parameter λ .
Light gradient boosting MachineBoosting approach.
Random forest regressorNumber of trees and features.
K-neighbor regressorNumber of neighbors.
Extra trees regressorNumber of trees and features.
AdaBoost regressorLearning rate.
Gradient boosting regressorLearning rate η .
Decision tree regressorMinimum sample leaf, maximum depth, split rule.
MLP regressorNumber of layers
Huber regressorStrength parameter α .
Ridge regressionRegularization parameter λ .
Table 3. Direct emissions. The training MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
Table 3. Direct emissions. The training MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
wMADRMSEMAE R 2
KEREERKEREERKEREERKEREER
0.010.0080.0080.9190.9240.2480.2900.540.16
0.050.0260.0400.7260.8800.2030.2800.920.23
0.10.0880.0790.9380.8500.2730.2760.100.25
0.150.0400.1120.5110.8300.1600.2750.940.27
0.20.0540.1440.5050.8140.1680.2800.840.28
0.250.0460.1740.4170.8000.1510.2870.930.30
0.30.0430.2020.3710.7880.1360.2990.940.30
0.40.0530.2520.3560.7690.1530.3280.920.31
0.50.0440.2900.2970.7620.1360.3640.940.31
0.60.0680.3160.3550.7720.1840.4100.880.31
0.70.0560.3240.3310.8100.1900.4760.900.30
0.750.0430.3200.3310.8480.1850.5240.910.30
0.80.0410.3080.3330.9060.2060.5860.900.29
0.950.0670.2410.8371.4830.7391.0850.690.25
Table 4. Direct emissions. The testing MAD, RMSE and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
Table 4. Direct emissions. The testing MAD, RMSE and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
wMADRMSEMAE R 2
KEREERKEREERKEREERKEREER
0.010.0159.3071.2053.3950.4120.9240.0330.005
0.050.07432.141.1365.9960.4131.2860.1970.002
0.10.14459.8681.2008.2800.4211.5950.0640.001
0.150.20680.0111.0809.7960.4321.7870.1970.0007
0.20.30897.2111.13311.0970.4661.9620.0580.0006
0.250.316114.241.05712.3990.4482.1420.1980.0005
0.30.367130.651.05313.7060.4582.3260.1940.0004
0.40.461157.0041.04816.1980.4652.6960.1820.0003
0.50.541171.651.04018.5280.4873.0660.1870.0002
0.60.663169.591.08020.5650.5113.4060.1120.0002
0.70.663147.291.03222.0970.5243.7070.1870.0001
0.750.691126.8331.04022.4380.5293.8180.1750.0001
0.80.708101.6661.04722.4270.5743.9030.1680
0.950.58925.0941.26622.2291.0104.3930.1200.0001
Table 5. The benchmark regression-based approaches used to predict direct nitrogen emissions.
Table 5. The benchmark regression-based approaches used to predict direct nitrogen emissions.
PhaseTrainingTesting
Technique MAERMSERMAERMSER
Support Vector Regression1.1181.199−13.560.4941.2700.02
Lasso Regression0.4180.736− 0.5050.5931.302−0.024
Light Gradient Boosting Machine0.39980.715−0.4520.5751.2700.025
Random Forest Regressor0.2430.745−0.120.4481.3625−0.1216
K Neighbors Regressor0.2550.702−0.080.4541.2680.028
Extra Trees Regressor0.2530.720−0.050.4451.320−0.053
AdaBoost Regressor0.3050.776−0.410.4891.2470.060
Gradient Boosting Regressor0.4110.727−0.460.5851.288−0.003
Decision Tree Regressor0.2490.748−0.1290.4481.362−0.121
Huber Regressor0.3110.767−1.601.0985.179−15.21
MLP Regressor0.5131.105−22.191.1373.587−6.77
Ridge Regression0.4150.732−0.800.6431.369−0.13
Kernel expectile regression 10.1900.3310.900.5241.0320.187
1 KERE model with respect to expectile level = 0.7.
Table 6. Indirect emissions. The training MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
Table 6. Indirect emissions. The training MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
wMADRMSEMAE R 2
KEREERKEREERKEREERKEREER
0.010.0120.0131.1371.1620.1600.2560.930.18
0.050.0440.0640.9431.0890.1390.2710.980.42
0.10.0930.1180.9511.0170.1550.2830.820.50
0.150.2130.1641.1900.9580.1890.2900.080.53
0.20.0680.2030.5800.9100.1060.2950.990.54
0.250.0660.2350.5110.8720.1030.3000.990.55
0.30.0640.2610.4600.8410.1000.3030.990.55
0.40.1070.2960.5040.8020.1360.3090.920.55
0.50.0560.3110.3360.7890.1040.3150.990.55
0.60.0520.3050.3010.8020.1180.3240.980.55
0.70.0990.2790.4480.8420.1810.3380.880.55
0.750.0490.2570.2850.8750.1540.3480.980.54
0.80.0480.2280.2960.9200.1790.3620.970.53
0.950.0730.0830.9481.1940.5170.5350.60.48
Table 7. Indirect emissions. The testing MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
Table 7. Indirect emissions. The testing MAD, RMSE, and MAE metrics for both KERE and ER are reported with respect to the various expectile levels ranging between 0.01 and 0.99 .
wMADRMSEMAE R 2
KEREERKEREERKEREERKEREER
0.010.00040.00050.0920.1750.1320.0920.1220.004
0.050.0050.0530.1040.2810.1820.1040.1230.887
0.10.0080.3130.1050.6140.2580.1050.3500.933
0.150.0010.6780.0820.9110.3170.0820.8390.943
0.20.0241.0540.1801.1630.3660.1800.1160.945
0.250.0321.4100.2091.3840.4080.2090.1390.945
0.30.0351.7540.2261.5930.4480.2260.1390.944
0.40.0322.3440.2341.9820.5280.2340.3560.941
0.50.0352.7390.2652.3400.6080.2650.1150.937
0.60.0402.8910.3182.6790.6950.3180.1580.929
0.70.0462.7620.3913.0100.7920.3910.3060.917
0.750.0312.5810.3523.1770.8460.3520.1140.908
0.80.0302.3160.3873.3490.9050.3870.1130.897
0.950.0510.8431.0103.8241.2301.0100.2530.821
Table 8. The benchmark regression-based approaches used to predict indirect nitrogen emissions.
Table 8. The benchmark regression-based approaches used to predict indirect nitrogen emissions.
PhaseTrainingTesting
TechniqueMAERMSERMAERMSER
Huber Regressor0.07360.2891−0.28460.2101.5640
Support Vector Regression0.89941.0082NA0.3081.5440.01
AdaBoost Regressor0.0850.312−0.0280.4462.009−0.65
Lasso Regression0.14890.3416−8.4870.2791.5660
Decision Tree Regressor0.08450.3213−0.11500.2271.576−0.02
Extra Trees Regressor0.08350.3196−0.04010.2251.575−0.02
Gradient Boosting Regressor0.14810.3408−8.43460.4491.990−0.62
K Neighbors Regressor0.09030.3217−0.8590.2371.573−0.01
Random Forest Regressor0.13330.53000.14150.2251.574−0.01
Light Gradient Boosting Machine0.22860.7311−338.590.2791.5660
Ridge Regression0.16000.4165−17.720.3111.578−0.02
MLP Regressor0.4150.732−0.800.7332.557−1.68
Kernel Expectile Regression 10.1540.2850.900.5240.980.18
1 KERE model with respect to expectile level = 0.7.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Benghzial, K.; Raki, H.; Bamansour, S.; Elhamdi, M.; Aalaila, Y.; Peluffo-Ordóñez, D.H. GHG Global Emission Prediction of Synthetic N Fertilizers Using Expectile Regression Techniques. Atmosphere 2023, 14, 283. https://doi.org/10.3390/atmos14020283

AMA Style

Benghzial K, Raki H, Bamansour S, Elhamdi M, Aalaila Y, Peluffo-Ordóñez DH. GHG Global Emission Prediction of Synthetic N Fertilizers Using Expectile Regression Techniques. Atmosphere. 2023; 14(2):283. https://doi.org/10.3390/atmos14020283

Chicago/Turabian Style

Benghzial, Kaoutar, Hind Raki, Sami Bamansour, Mouad Elhamdi, Yahya Aalaila, and Diego H. Peluffo-Ordóñez. 2023. "GHG Global Emission Prediction of Synthetic N Fertilizers Using Expectile Regression Techniques" Atmosphere 14, no. 2: 283. https://doi.org/10.3390/atmos14020283

APA Style

Benghzial, K., Raki, H., Bamansour, S., Elhamdi, M., Aalaila, Y., & Peluffo-Ordóñez, D. H. (2023). GHG Global Emission Prediction of Synthetic N Fertilizers Using Expectile Regression Techniques. Atmosphere, 14(2), 283. https://doi.org/10.3390/atmos14020283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop