Evaluation of Atmospheric Downward Longwave Radiation in the Brazilian Pampa Region

: Atmospheric downward longwave radiation ﬂux ( L ↓ ) is a variable that directly inﬂuences the surface net radiation and consequently, weather and climatic conditions. Measurements of L ↓ are scarce, and the use of classical models depending on some atmospheric variables may be an alternative. In this paper, we analyzed L ↓ measured over the Brazilian Pampa biome. This region is located in a humid subtropical climate zone and characterized by well deﬁned seasons and well distributed precipitation. Furthermore, we evaluated the performance of the eleven classical L ↓ models for clear sky with one-year experimental data collected in the Santa Maria experimental site (SMA) over native vegetation and high relative humidity throughout the year. Most of the L ↓ estimations, using the original coefﬁcients, underestimated the experimental data. We performed the local calibration of the L ↓ equations coefﬁcients over an annual period and separated them into different sky cover classiﬁcations: clear sky, partly cloudy sky, and cloudy sky. The calibrations decreased the errors, especially in cloudy sky classiﬁcation. We also proposed the joint calibration between the clear sky emissivity equations and cloud sky correction function to reduce errors and evaluate different sky classiﬁcations. The results found after these calibrations presented better statistical indexes. Additionally, we presented a new empirical model to estimate L ↓ based on multiple regression analysis using water vapor pressure and air temperature. The new equation well represents partial and cloudy sky, even without including the cloud cover parameterization, and was validated with the following ﬁve years in SMA and two years in the Cachoeira do Sul experimental site (CAS). The new equation proposed herein presents a root mean square error ranging from 13 to 21 Wm − 2 and correlation coefﬁcient from 0.68 to 0.83 for different sky cover classiﬁcations. Therefore, we recommend using the novel equation to calculate L ↓ over the Pampa biome under these speciﬁc climatic conditions.


Introduction
Atmospheric downward longwave radiation flux (L↓) is an important component of radiative balance and a relevant variable for meteorological and climate studies [1]. The L↓ is the thermal infrared radiation emitted from the atmosphere mainly by the H 2 O, CO 2 , O 3 molecules, and aerosols, and its maximum energy is in the longer wavelengths at about 10 µm [2]. Due to the wavelength range (4-100 µm) and intense interaction with water vapor in the atmosphere, which varies widely during the day, seasons, and from one region to another, L↓ estimates are highly sensitive to the amount of water in the atmosphere [3].
Theoretically, L↓ is obtained from Stefan-Boltzmann law, which states that all bodies emit energy in the form of radiation. According to this law, the radiation flux depends solely on the emissivity and the fourth power of the body's absolute temperature. However, it is not easy to determine L↓ due to difficulties to obtain the emissivity and temperature of the entire atmosphere (air column above the observer). Recent studies have suggested L↓ measurements to be performed using remote sensing and machine learning [4][5][6][7][8][9]. However, empirical and physical models are classically used to estimate L↓ and widely used worldwide. Physical models are based on the atmospheric radiative transfer phenomena, in which detailed profiles of atmospheric constituents are required, despite often being unavailable [10,11]. On the other hand, empirical models use atmospheric variables derived from surface-level observational data [12][13][14][15][16][17][18][19][20], which are specifically for definite atmospheric conditions and in general, need to be locally calibrated [21,22]. This hypothesis is justified by 90% of the L↓ coming mainly from the lower layers of the atmosphere (800-1600 m above the surface) where the molecules' concentration is high [23,24]. Most L↓ models are conceived under clear sky conditions, in which clear sky effective emissivity (ε e ) depends mainly on water vapor pressure near the surface (e a ) and air temperature (T a ). However, cloud presence increases L↓ reaching the surface [25]. Through numerous methods, various studies have combined clear sky emissivity classical models with cloud sky corrections using algorithms based on cloud cover or solar index to estimate L↓ [12,19,22,[26][27][28][29][30][31][32][33][34]. Flerchinger et al. [26] evaluated the combination of 13 clear sky emissivity models over ten cloud sky correction models with two algorithms for sky cover fraction (or cloud cover fraction) using data from sites across North America and China. Marthwes et al. [28] combined eighteen models for downward longwave radiation based on air emissivity and with six different cloud sky fraction models analyzing a dataset over the tropical forest in eastern Amazonia. These studies demonstrated the significant variability of results throughout the regions and models, but they do not performed a local model calibration. In the present study, we propose a joint calibration between clear sky emissivity and cloud sky correction models to estimate L↓ using a basic structure based on the formulation proposed by Maykurt and Church [35], Jacobs [36], Iziomon et al. [34], Sugita and Brutsaert [37], among others. The analysis was carried out using eleven classic clear sky equations for atmospheric downward longwave radiation and combined with three sky cover fraction functions. This methodology is expected to estimate L↓ more accurately since emissivity equations, in general, are obtained empirically, and it is believed that testing many combinations with sky cover fraction may generate better parameterization for applications in different areas of research, including agriculture, climate and solar energy. Additionally, we evaluated L↓ estimation results in different specific daily sky cover classifications: clear sky, partly cloudy sky, and cloudy sky. Moreover, a new empirical equation for emissivity under clear sky conditions is defined using air temperature and water vapor pressure.
We evaluated all these estimates over a set of experimental data obtained in native vegetation over the Brazilian Pampa biome. The Pampa biome is characteristic of southern Brazil, Argentina and Uruguay, although it is still not well characterized by L↓. In southern Brazil, this biome presents mostly grassland vegetation interspersed with gallery forests [38]. The subtropical climate is characterized by well distributed precipitation throughout the year (i.e., no dry or wet seasons), which favors cattle production and agriculture. Therefore, it is necessary to study L↓ and develop local models that can be utilized for both meteorological and agricultural models.

Site and Instrumentation
The data used in this study were collected in two experimental sites located in the Pampa biome in Southern Brazil that were roughly 60 km away distant from each other ( Figure 1). The Santa Maria site (SMA) is located in an area of native grassland that belongs to the Federal University of Santa Maria (UFSM) in Santa Maria, Rio Grande do Sul State (29 • 43 26.76" S and 53 • 45 34.92" W, altitude of 88 m). This site is commonly used for cattle production. The Cachoeira do Sul site (CAS) is located on a flooded rice paddy field in Cachoeira do Sul, Rio Grande do Sul State (30 • 16 37.59" S and 53 • 8 52.25" W, altitude of 40.5 m). The rice is cultivated from November to April and the surface remains fallow for the rest of the year. The climate in the experimental sites is classified as subtropical humid (Cfa, Köppen climate classification; [39]), with high summer temperatures (~35 • C) and low winter temperatures, which can even reach negative temperatures and with the occurrence of frost. The precipitation for both sites is well distributed year-round (between 83 and 157 mm by month  according to Brazilian National Meteorological Service (INMET) [40], and the daily average of relative humidity is high throughout the year [41]. The data used in this study were collected in two experimental sites located in the Pampa biome in Southern Brazil that were roughly 60 km away distant from each other ( Figure 1). The Santa Maria site (SMA) is located in an area of native grassland that belongs to the Federal University of Santa Maria (UFSM) in Santa Maria, Rio Grande do Sul State (29° 43′ 26.76″ S and 53° 45′ 34.92″ W, altitude of 88 m). This site is commonly used for cattle production. The Cachoeira do Sul site (CAS) is located on a flooded rice paddy field in Cachoeira do Sul, Rio Grande do Sul State (30° 16′ 37.59″ S and 53° 8′ 52.25″ W, altitude of 40.5 m). The rice is cultivated from November to April and the surface remains fallow for the rest of the year. The climate in the experimental sites is classified as subtropical humid (Cfa, Köppen climate classification; [39]), with high summer temperatures (~35 °C) and low winter temperatures, which can even reach negative temperatures and with the occurrence of frost. The precipitation for both sites is well distributed year-round (between 83 and 157 mm by month  according to Brazilian National Meteorological Service (INMET) [40], and the daily average of relative humidity is high throughout the year [41]. In the SMA site, the L↓ and global solar radiation (Rg) were measured by the net radiation sensor (CNR4, Kipp & Zonen, Deft, The Netherlands), while the relative humidity (RH) and air temperature (Ta) were measured using the thermo-hygrometer sensor (HMP155, Vaisala, Finland), all at the height of 3 m. The data period used in this study was from January 2014 to December 2019. More details on the SMA site are described in Zimmer et al. [42] and Rubert et al. [43].
In the CAS site, the L↓ was measured by a pyrgeometer (CG1, Kipp & Zonen, Deft, The Netherlands) and Rg by a pyranometer (LI200S, LI-COR Inc., Lincoln, NE, USA) at 6 m height and RH and Ta by a thermo-hygrometer (CS215-L, Campbell Scientific Inc., Logan, UT, USA) at 3 m height. Data were collected between 20 February 2013 and 19 February 2016. More details on the CAS site are described in Diaz et al. [44] and Souza et al. [45]. In the SMA site, the L↓ and global solar radiation (R g ) were measured by the net radiation sensor (CNR4, Kipp & Zonen, Deft, The Netherlands), while the relative humidity (RH) and air temperature (T a ) were measured using the thermo-hygrometer sensor (HMP155, Vaisala, Finland), all at the height of 3 m. The data period used in this study was from January 2014 to December 2019. More details on the SMA site are described in Zimmer et al. [42] and Rubert et al. [43].
In the CAS site, the L↓ was measured by a pyrgeometer (CG1, Kipp & Zonen, Deft, The Netherlands) and R g by a pyranometer (LI200S, LI-COR Inc., Lincoln, NE, USA) at 6 m height and RH and T a by a thermo-hygrometer (CS215-L, Campbell Scientific Inc., Logan, UT, USA) at 3 m height. Data were collected between 20 February 2013 and 19 February 2016. More details on the CAS site are described in Diaz et al. [44] and Souza et al. [45].
Measurements were processed in 30 min in both sites. Actual vapor pressure (e a , hPa) was obtained according to Allen

RH 100
(1) where: T a is the air temperature ( • C) and RH is the relative humidity (%).

Clearness Index
The daily clearness index of the atmosphere (K t ) over shortwave wavelengths (0.3-0.7 µm) represents the atmospheric transmissivity and is defined as the ratio between the daily integrated shortwave radiation received on surface-level (measured global solar radiation, R g ) and the daily shortwave radiation coming from the sun that is theoretically received at the top of the atmosphere (R 0 ) (extraterrestrial solar radiation) [2].
According to Kuye and Jagtap [47], the daily clearness index is used to classify the daily sky cover as cloudy sky (CD) for K t values between 0 and 0.35, partly cloudy sky (CP) for K t values between 0.35 and 0.65, and clear sky (CS) for K t values between 0.65 and 1. We adopt such classification in the present work.

Evaluated Parameterizations
The longwave radiation emitted from the atmosphere to the surface can be estimated by the Stefan-Boltzmann equation: where: σ is the Stefan-Boltzmann constant (5.67 × 10 −8 Wm −2 K −4 ) and ε e and T a (K) are the effective emissivity and the temperature of the atmosphere. The temperature and emissivity are integrated quantities over an atmospheric column above the ground that emits longwave radiation. However, ε e and T a are typically estimated from ground-based meteorological observations in most models for L↓ based on Equation (2). In general, the methods to estimate the effective emissivity ε e corresponding to the atmospheric emissivity under clear sky conditions are based on the ambient temperature and water vapor pressure (ε CLR (T a , e a )) [12][13][14][15][16][17][18]20]. The clouds increase the amount of L↓ reaching the surface by changing the emissivity and temperature of the atmosphere. Therefore, the emissivity should be able to include the effects of the cloud (i.e., effects of the sky cover), and various equations have attempted to estimate the increased L↓ produced by clouds [27,[35][36][37][48][49][50].
In order to obtain the effective emissivity of the atmosphere for all sky cover conditions, a general form for ε e combines the atmospheric emissivity under clear sky conditions (ε CLR (T a , e a )) with the cloud sky correction (g(K t )). To obtain the cloud sky correction, many models of sky cover fraction have been proposed, which is also called cloud cover fraction [27,48,[51][52][53][54][55]. To simplify, we suggest the cloud sky correction function using the clearness index (K t ). Thus, the effective emissivity can be written as defining: where µ and λ are parameters to be adjusted (hereafter called coefficients), and C(K t ) is the sky cover fraction function (as defined in Marthews et al. [28]). The term g(K t ) in Equation (3) guarantees increased emissivity when clouds are present, representing a correction function that considers sky coverage. Equation (4) is based on the models proposed by Iziomon et al. [34], Jacobs [36], Maykut and Church [35], Sugita and Brutsaert [37] and Choi et al. [56]. Therefore, if µ is zero, Equation (3) represents the clear sky emissivity; in other cases, the influence of cloudiness will be analyzed from the joint calibration of the sky's emissivity equations (ε CLR (T a , e a )) and the sky cover correction function g(K t ). Some equations for ε CLR and its original coefficients are presented in Table 1. It is important to note that each coefficient has different units to make atmospheric emissivity dimensionless. Equations EAN [12,18,19,57] only use e a in their formulation, while EST [12,13] uses T a , and EBT [7], EST [8], EID [9], EPR [30], and EKZ [48] use both meteorological variables. The EBT and EPR models are physically based models. The EBT model is based on the analytic solution of the Schwarzschild equation for the standard atmospheric lapse rates of T a and e a , and the EPR model was based on radiative transfer simulations. The other models were proposed empirically by using observations of longwave radiation for definite atmospheric conditions. Both physical and empirical model parameters and performance are significantly affected by geographical location and local atmospheric conditions and require site-specific validation and parameterization [57]. Table 1. Equations for atmospheric emissivity under clear sky conditions, ε CLR (T a , e a ), and original coefficients. T a is the air temperature (K) and e a is the water vapor pressure (hPa).

Reference
Code -Three parameterizations for estimating the sky cover fraction function C(K t ) are presented in Table 2 and used in this work. Many schemes to estimate sky cover fraction function are present in the literature. Kasten and Czeplak [53], Crawford and Duchon et al. [27], Konzelmann et al. [48], and Jegede et al. [55] have suggested equations that depend on theoretical clear sky estimations, in which the estimation of atmosphere turbidity and pressure-correction air function are necessary (as described in Marthews et al. [28]). Here, we decided to only use parameterizations using K t as a variable. Table 2. Equations for sky cover fraction function C(K t ).

Reference
Code

Proposed Model
A new empirical clear sky emissivity model for atmosphere downward longwave radiation is proposed: 16 (5) where the values of a 1 , a 2 , and a 3 are the parameters to be adjusted. The acronym EAI is used for this proposed equation. This proposal was formulated empirically using mathematical expressions based on classical equations variables ( Table 1). The Equation (5) describes a multiple regression analysis performed between the experimental [L↓/(σT a 4 )] and two atmospheric variables (T a , e a ) for clear sky days. The first term of Equation (5) connects the air temperature and the water vapor pressure. Water vapor in the atmosphere varies widely for small differences in the temperature and pressure conditions. The second term is related to relative humidity, which is represented by the actual vapor pressure. The third and final term is relative to air temperature.

Statistical Indexes and Analysis
L↓ estimations were compared with experimental data using the correlation coefficient (R 2 ), the root mean square error (RMSE), and percentage bias (PBias). These statistical indices are defined by where i is the sampler, L↓ is the simulated atmospheric downward longwave radiation, L exp ↓ is the experimental atmospheric downward longwave radiation, ( L ↓ ) is the average of the simulated values, ( L exp ↓ ) is the average of experimental values, and N is the number of simulated data. The analysis of the correlation coefficient results follows the classification suggested by Cohen [62].
Half-hour experimental atmospheric downward longwave radiation measurements obtained at the SMA site for the entire year of 2014 were used to evaluate the models. The sequence of analyses was: Step 1-evaluation of Equation (2) using Equations (3) and (4) with µ = 0 and ε CLR (T a , e a ) from Table 1 and original coefficients; Step 2-calibration of Equation (2) using Equations (3) and (4) with µ = 0 and ε CLR (T a , e a ) from Table 1 and the new proposed equation (Equation (5)); Step 3-calibration of Equation (2) using Equations (3) and (4) with µ = 0, C(K t ) from Table 2, ε CLR (T a , e a ) from Table 1, and the new proposed equation (EAI-Equation (5)); Step 4-validation of Equation (2) using the best results from Step 2 and 3 for both experimental sites.
Because K t cannot be calculated during the night, the C(K t ) was determined using the daily K t (or daily average). Therefore, in Steps 3 and 4, we used the equations of ε CLR for half-hour measurements combined with one value for C(K t ) for the entire day. The analyses were performed for different daily sky cover classifications, as described in Section 2.2.
The combination of Equations (3) and (4) yielded 48 possible solutions for Equation (2). In the coefficient calibration process, each one of them was adjusted by the non-linear least squares method. We used the MATLAB (Mathworks Inc.) program and lsqcurvefit function from Optimization Toolbox™ to solve the least squares method. Only SMA site data from 2014 were considered in the calibration algorithm and individual analysis for each class of sky cover performed. The performance of the models' combination was evaluated using the best RMSE (near zero).

Evaluation of L↓ Models for SMA Site for 2014
The SMA experimental data for 2014 used to evaluate the L↓ equations presents the following sky cover classification: 27.40% for clear sky days (CS), 46.57% for partly cloudy days (CP), and 21.64% for cloudy days (CD). The clearness index did not show seasonality for SMA (data not shown).

Using Original Coefficients-Step 1
The atmospheric downward longwave radiation estimates were obtained by following Step 1. The statistical results of the eleven equations used to estimate L↓ are shown in Table 3. For entire period (EP) and all sky classifications, the EID presented the best results in all analyzed statistical indexes. The ESW and EIJ equations presented the worst results for R 2 in all periods, while the EKZ showed the worst results for RMSE and PBias. All equations presented better R 2 for CS with values above 0.67. This was expected because the original equations were developed for clear sky conditions. In general, all results showed the considered models to underestimate L↓ (negative PBias) with errors above 20 Wm −2 . The errors increased as the clearness index decreased. According to Marthews et al. [28], the analysis of L↓ over a tropical lowland forest revealed that RMSE values below or equal to 20 Wm −2 are acceptable for L↓ estimates. Generally, both physical and empirical model parameters and performance are significantly affected by geographical location and local atmospheric conditions and require site specific validation and parameterization [56]. Therefore, all L↓ equations of Table 1 need to be calibrated for the Brazilian Pampa regions to increase the accuracy of the models, as also highlighted by Li et al. [29].

Calibrating the Coefficients-Step 2
We calibrated the coefficients (a 1 , a 2 , a 3 ) of the L↓ according to Step 2. The new proposed equation (L↓ EAI) was also calibrated. The results of the models' performances are presented in Table 4 with the respective calibrated parameter values. The calibrated coefficient values are small for some parameters and even large for the same parameterization in others, as also noted by Carmona et al. [19] and Zhu et al. [63]. By comparing the statistical indices of the models from Tables 2 and 3, the PBias was reduced drastically by the calibration (around two order of magnitude) in all sky cover classification, resulting in the overestimation for most of them. In general, RMSE was reduced for all sky cover classifications, from 1% for EID to 85% for EKZ. Nevertheless, the analysis of the different sky cover classifications shows that the RMSE decreased by around 40% when the clearness index decreased, in addition to R 2 slightly decreasing. These results are likely due to the small variability of the downward longwave radiation in experimental data, as shown in the x axis of Figure 1 when the clearness index decreased.
In general, the EST presented the smallest errors of the classic calibrated models, with a slight underestimation of the experimental data (negative Pbias), followed by an EID model with errors around 3% higher, although with overestimation (positive Pbias). The worst results were found for ESW, with errors greater than 30 Wm −2 except for CD sky cover classification (21.18 Wm −2 ).
Guo et al. [64] evaluated six models (which were also evaluated in this study: EBR, ESW, EIJ, EBT, EID, EPR) for 71 globally distributed sites and reported that the worst results were obtained by equations that did not consider e a in their parametrization (ESW and EIJ). Similar results were found by Carmona et al. [19], who evaluated the same models for a site in the sub-humid Pampean region in Argentina. Here, after calibration, ESW also had the worst results, followed by EKZ and EGR, which, however, are functions of e a (EKZ is also a function of T a ).
The dispersion between the experimental and estimated EAI L↓ are shown in Figure 2. The EAI overestimated the lowest experimental L↓ values, especially during night periods ( Figure 3). Principally in CD conditions, the high experimental L↓ values are better represented by EAI. In general, the EAI presented a better performance for all sky cover classifications compared to the classical equations (Table 4). Moreover, this equation showed the smallest variability in the statistical index for different sky cover classifications, indicating that this formulation accurately represents all sky conditions even using the EP calibration. The best performance of the EAI can be due to the equation was obtained over periods of entire days classified as clear sky. On these days, the sky may have presented some moments of lower clearness indices (some clouds in the sky), but that can be masked by the daily average sky cover classification.   Table 4.

Calibrating the Coefficients-Step 3
We calibrated the coefficients (a 1 , a 2 , a 3 , λ, µ) of the L↓ obtained from Equations (2) and (3) with the ε CLR (T a , e a ) from (Table 1) and C(K t ) from ( Table 2). The results of the five best combinations for each sky cover classification are presented in Table 5 with the respective calibrated parameter values. By adding a cloud sky correction, all combinations presented better results than without this correction for all sky cover classifications, where, in general, RMSE decreased by around 5 Wm −2 and R 2 increased by almost 10%. The EST and EID models are present in all five of the best results, especially CQB and CCB sky cover fraction functions. The smallest errors were obtained in CD classification, where EST_CWU presented the best result, although the R 2 was lower in this sky classification. For the other sky cover classifications, in general, the best results were obtained by EST_CCB and EST_CQB.  The results of all combinations using the new proposed equation (EAI) are presented in Table 6. The EAI equation with sky cover fraction function combinations had better performance than the results demonstrated before, with an RMSE below 20 Wm −2 in all sky cover classifications. The RMSE decreased by no more than 20% from the calibrated one without C(K t ) (comparing Table 4 with Table 6). For CD, there was no significant difference (<3%) in the results with C(K t ). As the EAI was formulated considering the daily mean of K t , when implementing the parameterization of cloud cover C(K t ), no significant influence was observed. In general, EAI presented better results than the classical equations for all sky covers (CS, CP, CD) and did not significantly improve the results using the cloud sky correction, indicating that the EAI is recommended to use in this region, even without considering the cloud sky correction. Here, we will evaluate the EAI calibrated over the entire period (Table 6).
To evaluate whether the L↓ (EAI) equation obtained by calibration for SMA adequately describes the downward longwave radiation at different locations in the Pampa biome, we used data from the SMA (the entire years of 2015-2019) and CAS (from 20 February 2013 to 16 April 2016) experimental sites. To compare EAI with a classical model, we will also present EST results, which is a model that showed the best results in all evaluations ( Table 5).
The EAI equation applied to CAS presented similar results to SMA for R 2 and PBias (~1%) and a difference of around 10% regarding the RMSE ( Table 7), showing that the proposed equation efficiently estimates L↓ in the Pampa biome. It is also noted that when we implemented the parameters for cloud sky correction in EAI, the results improved by 20%; the same behavior was found for CAS. When comparing EAI with EST and its original coefficients without the cloud sky correction, the EAI presented the best results. After incorporating the parameterizations in EST, the statistical indices significantly improved except for using CCB. Nevertheless, using the CQB and CWU parameterizations presented the same statistical indices in both sites for both EAI and EST. The SMA and CAS experimental sites are approximately 60 km far from each other, although they have a significantly different land cover. The SMA surface is covered with natural grassland used by cattle, while CAS is in a lowland rice paddy. However, comparing the diurnal cycle of atmospheric variables for the same data period for the SMA and CAS sites, we can see very similar behavior from 1 January 2014 to 16 April 2016 ( Figure 3). The most notable difference is in diurnal RH when SMA presents smaller values than CAS. The experimental L↓ in CAS and SMA are quite similar during the day, but in the night period, SMA is around 10 Wm −2 greater than CAS (primarily in the night and early morning, from 0 to 8 h). The new model presented proved to be valid, given the results reported herein, in which EAI fit well without requiring specific calibration for CAS. Therefore, the EAI equation can be an excellent alternative to represent L↓ regarding similar weather conditions. Thus, the proposed equation (EAI) presented a good performance to estimate L↓ in the Brazilian Pampa biome with values similar to those obtained by Carmona et al. [19], who developed a model for a region with similar vegetation to that analyzed here (Argentinian Pampa) but included a sky cover fraction (cloud fraction) in their formulation for L↓. In Carmona et al. [19], two equations are proposed, one only for clear sky and another for all sky conditions. However, they used the RH as atmospheric variables in emissivity equations, and to calculate the sky cover fraction, it is necessary to estimate the theoretical solar radiation received in clear-sky conditions, depending on the sky conditions such as the turbidity coefficient and atmospheric pressure. Nonetheless, here we chose to analyze only equations using e a and/or T a , and we suggest employing equations for a sky cover fraction that only needs the clearness index calculated using the extraterrestrial solar radiation, making the proposed approach much easier to estimate L↓.

Conclusions
An L↓ estimation analysis using different equations for emissivity and cloud sky correction concerning experimental data measured in the SMA site was presented. First, eleven well known L↓ models were analyzed for SMA to test their accuracy under different sky cover classifications (CS, CP, and CD). The results showed that all models presented high error and underestimated values compared to the experimental data. Thus, it was necessary to calibrate the equation's coefficients to adjust the L↓ models, leading to better statistical indexes.
A general structure to the cloud sky correction was added to the clear sky emissivity models to improve L↓ estimates and determine how they influence cloud presence in the different periods analyzed. Moreover, we did a joint calibration between the clear sky emissivity model and cloud sky correction. The results showed a considerable improvement in the RMSE compared to the models without cloud sky correction function, especially under CD conditions. However, one limitation of this study was calculated for C(Kt) in the daily mean.
A new equation (EAI) was developed to reproduce L↓ data more accurately. The EAI was calibrated for SMA and showed better performance than the other models tested, with an RMSE ranging from 13 to 21 Wm −2 , PBias close to zero, and higher values than the correlation coefficient.
The EAI equation was evaluated for a long-term dataset for SMA and CAS with good results. The model proposed has a simple parameterization using the fundamental meteorological variables water vapor pressure and air temperature, and its formulation does not require cloud sky correction to satisfactorily represent the L↓. Therefore, we recommend using the proposed equation to calculate L↓ over the Pampa biome under these specific climatic conditions and without sky cover distinction.
As a complementary part of this work, the results can be used in the verification and calibration of the Brazilian Global Atmospheric Model-BAM [65], the atmospheric component of the Brazilian Earth System Model (BESM), since most of the BAM evaluation studies have been carried out for the most tropical portion of Brazil [66], with no work using the Pampa region from southern Brazil.