key: cord-0073798-sjpivy9a authors: de Souza, Amaury; Aristone, Flavio; Abreu, Marcel Carvalho; de Oliveira-Júnior, José Francisco; Fernandes, Widinei Alves; Pobocikova, Ivana title: Spatio-temporal variations of tropospheric nitrogen dioxide in South Mato Grosso based on remote sensing by satellite date: 2022-01-20 journal: Meteorol Atmos Phys DOI: 10.1007/s00703-021-00855-5 sha: d6c1eae6cb213474808d1dff010c50df6a95f3ed doc_id: 73798 cord_uid: sjpivy9a The study evaluates the characteristics of tropospheric nitrogen dioxide (NO(2)) concentrations from 2005 to 2020 in eight cities in the State of Mato Grosso do Sul (MS), Midwestern Brazil, using data available from the Ozone Monitoring Instrument (OMI). The average concentration varies from 2981 × 10(15) molecules/cm(2) in Campo Grande, where commercial, industrial and vehicular activities take place, to 2906 × 10(15) molecules/cm(2) in Corumbá and 3035 × 10(15) molecules/cm(2) in Porto Murtinho, regions of livestock and biomass burning. The results, based on the Mann–Kendall (MK) test, show a significant increase (p < 0.05) in the NO(2) column levels in the region. For each of the eight cities studied, a significant seasonal cycle of NO(2) columns was determined. The maximum value of NO(2) concentration was observed in the dry period, from July to September, while the minimum value was registered in the rainy period, from October to March. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00703-021-00855-5. Nitrogen dioxide (NO 2 ) is a reactive and short-lived gas, with natural and anthropogenic sources, are harmful contaminants, as they are associated with adverse effects on human health, especially with regard to respiratory diseases, and the environment (Souza et al. 2018a, b) . The main sources of NO 2 are fossil fuel combustion, biomass burning, soil emissions and lightning (Hagenbjörket et al. 2017) . NO 2 is a toxic air pollutant in the high-concentration condition and plays an important role in tropospheric chemistry as a precursor of tropospheric ozone and secondary aerosols Fan et al. 2021; Castellanos and Boersma 2012; Grajales and Baquedo-Bernal 2014; Fantozzi et al. 2015; Griffin et al. 2018; Xue et al. 2020; Nakada and Urban 2020) . Observations of the space-time variations of NO 2 form the basis for understanding the spatial distributions and temporal trends of NO 2 . Nonetheless, the limitation on insufficient spatial and temporal coverage of ground-level monitoring networks is the main challenge to understand the NO 2 dynamics in the atmosphere. Many techniques and methods have been used successfully in monitoring atmospheric NO 2 based on in situ surface measurements, remote sensing from satellite sensors and ground-based instruments (Xiao et al. 2018; Drosoglou et al. 2018) . Although in situ measurements and remote sensing of ground-based instruments show high accuracy and precision, their usefulness in determining the spatiotemporal distributions of trace gases is limited due to their sparse spatial and temporal coverage. Space-based measurements provide information on large-scale NO 2 distributions and in areas where in situ and ground-based systems cannot be easily deployed (Lamsal et al. 2010) . Measurement of atmospheric NO 2 is ordinarily obtained from in situ data of air quality stations (AQS) (Fan et al. 2020 ). However, AQS lack space-time coverage, particularly in countries like Brazil, where measurements are restricted to metropolitan regions (MR)- (Zeri et al. 2011 (Zeri et al. , 2016 . On the other hand, the data obtained via Remote Sensing provide data with larger spatiotemporal coverage in relation to AQS and without failures or discontinuities in their time series. Satellite observations provide NO 2 data from the total tropospheric column over an area, while AQS measures locally (Richter et al. 2005 ). The first NO 2 monitoring in the troposphere started in 1995 using the Global Ozone Monitoring Experiment (GOME) on the European Remote Sensing Satellite-2 Responsible Editor: Sang-Woo Kim. (ERS-2). Ozone Monitoring Instrument (OMI) started aboard the Aura satellite in 2002, and in 2004 the new GOME-2 began to collect data. The Tropospheric Monitoring Instrument (TROPOMI) aboard the Precursor Copernicus Sentinel-5 satellite started in 2017. The OMI measures atmospheric NO 2 daily since late 2004 with spatial resolution of 13 × 24 km in the nadir visualization, increasing the size to 24 × 135 km for widest viewing angles. Several studies have been performed using satellite observation data of tropospheric NO 2 columns (Richter et al. 2005; Castellanos and Boersma 2012; Schneider and Van Der 2012; Lelieveld et al. 2015; Khokhar et al. 2015; Irie et al. 2016; Cai et al. 2018; Yanfang et al. 2019; Yavaşlı 2020; Zheng et al. 2018) , and they show that its concentrations are dependent on economic, industrial, and other human-controlled activities but not correlated to emission policies to improve air quality. Several inventories on both global and regional scales of NO 2 concentrations normally use regression models for trending analysis. Nevertheless, verification of tropospheric NO 2 trends has not yet been carried out for the State of MS, Brazil. In this paper, the tropospheric concentration of NO 2 is analyzed for the period 2005-2020 throughout the available OMI data. It has been highlighted of the characteristics of the spatial distribution of NO 2 using the Kendall seasonal test for each pixel. The State of MS is located in the Midwest Region of Brazil and encompasses approximately 358,159 km 2 , see Fig. 1a . It is known for its agricultural activities, the main economic products of MS being soybean and beef cattle production. The State topography is represented in Fig. 1b and shows elevation ranging from 24 to 1100 m (Teodoro et al. 2016) . The average annual temperature varies from 20 to 26 °C, and the mean annual rainfall ranges from 1000 to 1900 mm. The climate presents well-defined dry seasons from April to September, when the higher rainfall records are observed in the southern portion of the state. The rainy season goes from October to March, the northern region of the State receives higher rainfall records compared to the southern region. The Köppen's classification divides the climate diversity of MS into several regions: (1) "Aw" in the Southeast and Northern parts of the state; (2) "Am" for the central region; (3) "Af" in the Southwest; and (4) "Cfa" in the Southern, see Fig. 1c . In the Southwest portion of MS, south of the Pantanal region, which lies between − 21° and − 22° latitude, the climate is characterized as tropical forest ("Af"), with rainfalls distributed throughout the year. The central portion of the State is predominantly characterized by tropical monsoon climate ("Am"), with a small dry season during winter. In the northern part, in a small portion of the central region, and the Southeastern region of the state, the climate is characterized as savannah ("Aw"), which tends to present dry winters and rainy summers. Only in the Southern of the state, the climate is humid during all seasons ("Cfa"), with hot summers and average temperatures above 22 °C (Souza et al. 2012; Alvares et al. 2013) . The biome diversity of MS is represented in Fig. 1d . It includes areas of the Atlantic Forest, the Cerrado, and the Pantanal with 14%, 61%, and 25% of the state-encompassed area. The Atlantic Forest is particularly an important biome due to its abundant biological diversity. It has gained substantial interest as a conservation area since this biome has been considerably reduced in the past. The Brazilian Cerrado is a vast tropical savanna ecoregion widely known for its native habitats and rich biodiversity. It represents the second largest biome in South America, after the Amazon. The Cerrado of South Mato Grosso is located in two hydrographic regions of Brazil, represented by the rivers Parana and Paraguay. The Pantanal is the world's largest inland flooded grasslands and savannas wetland region. It is home to rich wildlife and worldwide known for its unique biome. The region is considered delicate regarding biodiversity due to environmental degradation and damage (Teodoro et al. 2016) . Eight cities of South Mato Grosso have been studied in this paper. These cities are, the State capital Campo Grande, Chapadão do Sul, Corumbá, Coxim, Dourados, Ponta Porã, Porto Murtinho, and Três Lagoas, all represented in Fig. 1 . They have been chosen accordingly the availability of meteorological and NO 2 concentration data, expressed in 10 15 molecules/cm 2 , for the period from 2005 to 2020. The meteorological data series have been collected from the Brazilian National Water Agency (ANA) database, from Hidroweb (www. snirh. gov. br/ hidro web/ publi co/ medic oes_ histo ricas_ abas. jsf), and from the Brazilian National Institute of Meteorology (INMET) (www. inmet. gov. br), all accessed in in January 2021. The data include time series of maximum and minimum air temperature, i.e., T max and T min , respectively, relative humidity (RH) in %, and precipitation (Prec) in mm. A summary of the meteorological conditions for the studied period is presented in Table 1 . The eight cities have different characteristics related to geographic position, population and main economic activities, as follows: (1) Coxim, located in the Cerrado region of influence, with small patches of open fields. It presents itself in its different physiognomies and in enclaves with the Seasonal Forest. The planted pasture is expressive in the central portion of the municipality. Its economy is strongly linked to fishing, culture, livestock and agriculture tourism; (2) Três Lagoas, the predominant vegetation is the Cerrado with strips of Atlantic Forest. Its economy is strongly linked to fishing, livestock and pulp tourism. It is the third most populous city in MS State, with a population of around 123,000 inhabitants; (3) Campo Grande, the State capital, is located on the watershed of the Paraná and Paraguay River basins. The Guarani Aquifer passes under the city, being the capital of the State holding the largest percentage of the Aquifer within the Brazilian territory, the predominant vegetation being the Cerrado. Its economy is strongly linked to tourism and agriculture. It is the third largest and most developed urban center in the Midwest Region of Brazil, and the most populous city in MS and the 19th most populous in Brazil; (4) Chapadão do Sul, located in the Cerrado region The nutritious sludge layer that remains in the soil after flooding allows for the development of a rich flora. In areas where floods dominate, but which are dry during the winter, vegetations, such as carandá and paratudal palms, occur. Its economy is quite diversified, highlighting the activities of mining, fishing and tourism, and agriculture. The highest rates of fires in MS occur in the region of Corumbá. The population of Corumbá is 112,000 inhabitants. (7) Porto Murtinho, in the Pantanal biome, in a transition area with the Cerrado, distributed almost equitably, the coverage is typical of the Pantanal (Cerrado Estépico). There is also cultivated pasture and some crops, the most representative activity is agriculture and was a great producer of mate. The population of Porto Murtinho is around 15,000 inhabitants. (8) Dourados, its natural vegetation is clear fields, also having large parts of Cerrados and large patches of tropical forests, since it is located in area Atlantic Forest area but in a transition area Atlantic Forest-Cerrado. Its economy is mainly based in agriculture and planted pasture. It is the most populous city in the interior of MS State, with a population of approximately 228,000 inhabitants. For the analysis, data were collected in grid format with spatial resolution of 1° × 1° and daily temporal resolution referring to areas in the regions of Mato Grosso do Sul aboard the Aura satellite in 2004. The main objective of the OMI is to obtain global measurements at high spatial and spectral resolution of a series of trace gases in the troposphere and stratosphere . It can distinguish between types of aerosols, such as smoke, dust and sulfates (Duncan et al. 2013 ). According to Torres et al. (2013) , this is done through the OMAERUV algorithm that uses pre-computed reflectance from the upper atmosphere for a set of 21 aerosol models composed of three types of aerosols (dust, carbonaceous aerosols, sulfate-based aerosols). As described on the project's own page, the OMI is a nadir imaging spectrograph that measures solar radiation backscattered by the Earth's atmosphere and surface over the entire wavelength range from 270 to 500 nm with a spectral resolution of about 0.5 nm. The ratio of the calibrated radiance and irradiance is equal to the Earth's atmospheric bidirectional scattering distribution function (BSDF). Light entering the telescope is depolarized using a scrambler and then split into two channels: the UV channel with a wavelength range 270-380 nm and the VIS channel with a wavelength range 350-500 nm ). These factors are important to obtain the concentration of trace gases and to allow the IMO to monitor tropospheric pollution phenomena, such as biomass burning and industrial pollution, as the tropospheric pollution record is essential to study the human impact on the Earth's atmosphere and climate (Dobber 2006). The preliminary analyses consist of establishing position and dispersion statistics of notched box-plot for the NO 2 concentration. The medians obtained for the listed cities in MS on a monthly scale are then compared. The trend analysis has been applied to the monthly data. The Mann-Kendall (MK) non-parametric test (Mann 1945; Kendall 1975) has been applied to define the trend. The linear regression analysis to verify the proportion of increase or decrease of the response variable is based on the slope of the characteristic function. The Pettitt test (1979) has been applied to verify the point of interruption of the line (beginning of the trend). The estimated slope of the trend is calculated using a generalized version of Sen's slope estimator (Sen 1968 ). The MK is the trend verification test recommended by the World Meteorological Organization (WMO), widely used in hydrological studies (Burn et al. 2011; Sa'adi et al. 2019) . It tests whether the series presents a downward or upward trend. The null hypothesis means that there is no trend. For a time series, the Mann-Kendal (S) statistic is defined by: where x i is the i-th observation and x k is the k-th observation immediately after the i-th observation, n is the number of elements in the time series. The variance is determined by the following relationship: where m is the number of tied groups, t i is the number of loop lengths equal to i. The MK test requires that the sample data are serially independent, i.e., absence of autocorrelation. Correlated sample data in series affect the test ability to correctly assess the significance of the trend. To eliminate any effect of serial correlations in the MK test, autocorrelation analysis has been performed. A correlogram gives a summary of correlations at different periods of time. It consists of the graphic analysis of the autocorrelation function for sequential values of delay k = 0, 1, 2, …, n. In the occurrence of verified autocorrelation, the MK test was applied with the modification proposed by Yue and Yang (2004) . This modification involves a variance correction to address the issue of serial correlation in the trend analysis, maintaining the significance of the test (Sa'adi et al. 2019) . Regression analysis provides the trend rate using the slope of the resulting function. The slope value indicates the magnitude of the trend while the sign (±) indicates an increase or decrease behavior of the tested variable. However, the slope must be confirmed by the Student's t test for statistical significance, i.e., whether it differs from zero or not. The Pettitt test (1979) uses the Mann-Whitney U (t, M) statistic to verify whether two samples belong to the same population. Thus, the test places the data series in two samples: x 1 , x 2 , …, x t and x t+1 , x t+2 , …, x n , to calculate U (t, M) as follows (Kengni et al. 2019 ): with: Pettitt statistics account for the number of exceedances in which the value of the first sample exceeds that of the following sample. The null hypothesis of the test is that there is no breaking point in the data series. The level of significance is adopted 5%, as assumed by other authors (Uliana et al. 2015 ). Principal component analysis (PCA) was performed to explain the variance structure of the data through linear combinations of the original variables. The variables used were: NO 2 concentration (10 15 molecules/cm 2 ), maximum (T max ) and minimum (T min ) temperature of the air (°C), (RH) relative humidity (%), and (Prec) precipitation (mm) for the MS cities. PCA allows the use of variables correlated with each other, to form components capable of explaining the variance of the data. The main components were extracted using the correlation matrix between the variables (Shi and Harrison 1997) . PCA was used after calculating the Kaiser-Meyer-Olkin (KMO) index. Bartlett's sphericity tests were first performed to confirm the adequacy of PCA to the NO 2 concentration and meteorological data. The KMO test verifies the correlation measure of the independent variables. The value of the test varies from 0 to 1, whereas values below 0.5 indicate that the application of PCA is inappropriate. The Bartlett's test of sphericity evaluates whether the correlation matrix is an identity matrix, which would indicate that there is no data correlation, indicating that the factorial model is inappropriate. The position and statistic dispersion of the NO 2 tropospheric concentration data are presented in Table 2 . The box-plot displaying the monthly distribution of the tropospheric NO 2 concentration in the MS studied cities between 2005 and 2020 is exhibited in Fig. 2 . The city of Campo Grande present as average of NO 2 concentration 2.98 × 10 15 molecules/ cm 2 , and the median 2.87 × 10 15 molecules/cm 2 . Ponta Porã, which is in the south of the State and borders Paraguay, registered higher values: 3.038 × 10 15 and 2.935 × 10 15 molecules/cm 2 as average and median, respectively. Porto Murtinho, located in the west of the State of MS and on the border with Bolivia, registered the second-highest mean 3.035 × 10 15 molecules/cm 2 and median 2.814 × 10 15 molecules/cm 2 . The distribution of data points over the 15 years shows that, in addition to each city being uniquely identifiable, Ponta Porã and Porto Murtinho exhibit some heavy tail properties. The resulting parameters of asymmetry and kurtosis (Skew, kurtosis) for these cities are the largest in comparison to the other studied cities. NO 2 concentrations in the State of MS indicated less seasonal amplitude among the annual averages of 0.20 × 10 15 molecules/cm 2 , with higher values in July, August and September, and constant values between January and June (Fig. 2) . From October to December, there is a decrease in NO 2 concentration, probably due to the excessive use of power plants and domestic heating, which increase the NO 2 levels during the winter. Another point is the burning of fields in the region to clean the pastures (Caúla et al. 2015; Oliveira Júnior et al. 2020) . The notch box-plot shows differences between some of the cities in their medians, since the notches do not overlap. Ponta Porã presents higher median, especially from January to August, while Dourados shows increase of the median from September to December. Both cities are located in the south portion of the MS, at the boundary between the Cerrado and the Atlantic Forest biomes, where the crops and pasture are expressive economic activities (Tomei et al. 2020 ). All the cities displayed a delay of one or two months in relation to the highest concentrations of NO 2 , from July and October relative to the driest months, from June to September (Teodoro et al. 2016; Abreu et al. 2020) . Cities in the north-northwestern portions of the State of MS present lower medians. It is the case of Coxim, located in the transition between the Pantanal and the Cerrado, and Chapadão do Sul, in the Cerrado. These cities have lower crop density and lower fire foci (Oliveira-Júnior et al. 2020), which may contribute for lower concentrations of NO 2 . The fires are associated with the cleaning of areas for agriculture and contribute to increase the concentrations of NO 2 . In contrast, the monthly behavior demonstrates similarity between the concentration of NO 2 in the cities throughout the year, with a higher NO 2 concentration between July and November, reaching the maximum values in August and (Feng et al. 2001; Garrett et al. 2010; Zhao et al. 2020) . The anthropic action of using fire to clean areas for planting agricultural commodities in the rainy season (Silva Junior et al. 2020) occurs during the dry period, mainly in the months of August and September. This period precedes the rainy season, which begins in October (Teodoro et al. 2016; Abreu et al. 2020) . The burning of biomass has been reported as one of the responsible factors to increase the NO 2 concentrations in some periods of the year (Ossohou et al. 2019) . The effect of seasonality is therefore related to anthropogenic human activities, not exclusively natural events. The analysis of the NO 2 levels data series showed autocorrelation and trends for the studied period, as detected by the modified MK test. The cities of Dourados and Porto Murtinho showed decrease in NO 2 concentrations, Coxim presented increase instead, considering all the monthly data exhibited in Table 3 . Calculation of Sen's inclinations of significant trends indicate changes up to ± 0.001 × 10 15 molecules/cm 2 . These trends changed approximately 1% on average. In terms of magnitude, the cities of Dourados and Porto Murtinho decreased by 0.084 × 10 15 and 0.064 × 10 15 molecules/cm 2 respectively, while in Coxim, there was an increase of 0.095 × 10 15 molecules/cm 2 . These results are expressed in Fig. 3 . Additional relevant aspects about the concentrations of NO 2 in the studied period are the peaks and valleys in several areas in the MS. In general, the effects of the 2008 global economic crisis can be observed during the period 2009-2011 in most of the assessed cities. Therefore, variations in the measurements during this period are associated with depletion in industrial production. The blocking measures due to COVID-19 in 2020 affected industries, transports, and other human activities in Brazil and many other countries. Therefore, emissions of air pollutants during the pandemic have been drastically reduced in several countries Wang et al. 2021) . In Brazil, the impact of COVID-19 containment actions in terms of air pollutants is still being evaluated. In this study, reductions of the averaged NO 2 concentrations have been observed during 2020 compared to the period from 2005 to 2019 in only two cities, Dourados and Três Lagoas, in the order of 0.02 × 10 15 and 0.03 × 10 15 molecules/cm 2 . In the other studied cities, higher concentrations of NO 2 have been observed during 2020, with an average of 0.01 × 10 15 to 0.22 × 10 15 molecules/cm 2 , see Table 4 . Therefore, although for some months of 2020, the concentration of NO 2 was lower, in general it was higher compared to the period from 2005 to 2019. These results suggest the general non-commitment of the containment values during the pandemic in Brazil, which may have contributed to the excessive number of deaths, since the relationship between COVID-19 and its mortality with the concentration of NO 2 in the atmosphere has been evidenced. The concentration of NO 2 as a precursor to the formation of secondary particulate material may favor COVID-19 and make the respiratory system more susceptible to this infection (Mele et al. 2021) . Figure 4 shows the monthly dispersion of NO 2 concentration as a function of time in years, and the regression curves where most of the cases can be observed. The MK test for the monthly data series shows some trends for specific months. In addition, during the dry period in MS, which is from April to October, some relevant trends have been detected. Increase trends are observed for Corumbá in July and for Coxim in July, August, and December. These cities are located in the driest area of South Mato Grosso With the exception of Chapadão do Sul, the decreasing trends in the NO 2 concentrations during April and May may be associated with the characteristics of precipitation in the southern region of the state, in which there are higher total precipitates (Abreu et al. 2020 ). There is a noticeable relationship connecting the environmental impact of agricultural expansion and fires to the NO 2 concentrations. Although few trends have been detected, they occurred in months of area cleaning activities with the use of fires that result in large emissions of gasses, include NO x , which contribute to increase the greenhouse gas emissions and global warming (Oliveira- . Table 5 shows the MK test results, Tau, Sen's slope, and the break point, i.e., the year of change by the Pettitt test validating a significant trend in MK test. The changes indicated an increase of NO 2 concentration between 0.03 × 10 15 and 0.070 × 10 15 molecules/cm 2 , especially in drier months for cities in the south of the State of MS. The decrease of values varied between 0.02 × 10 15 and 0.04 × 10 15 molecules/ cm 2 . The break points in the monthly series occurred from 2010 to 2014, i.e., from the first to the third quartile of the data series. Overall, the eigenvectors are elevated for all variables, indicating that PC1 is an overview of the relationship between meteorological variables and NO 2 concentration, in which the higher the RH and Prec, the lower the NO 2 concentrations. The role of temperature is ambiguous, indicating that in cities where these variables had positive/ negative eigenvectors, the highest temperatures are associated with the highest/lowest concentrations of NO 2 . Higher scores indicate places with higher concentrations of NO 2 , lower precipitations and relative humidity. These results also reaffirm the seasonal factor of NO 2 concentrations, which are mainly related to humid climatic conditions (Hua et al. 2021 ). The scores are higher from July through August, the driest months of the year. The higher temperatures and RH in the atmosphere, which occur during summers in MS, are typically beneficial for the removal of NO 2 through Photolysis (Feng et al. 2001; Dong et al. 2009 ). Nevertheless, at some cities (Dourados, Corumbá, Ponta Porã and Chapadão do Sul, the lower temperatures would be allied with lower concentrations of NO 2 . PC2 also originated higher eigenvectors for almost all variables in all cities varying its pattern among cities, i.e., presenting positive and negative eigenvector values. In general, this component indicates the driest months with the highest scores and the most humid months with the lowest scores. Its eigenvectors for NO 2 is always positive. Thus, PC2 indicates differences in the dynamics of NO 2 throughout the atmosphere of every region of MS. Generally, moisture and precipitation are NO 2 dispersing agents (Cai et al. 2018) , as well as high temperatures cause NO 2 photolysis, reducing its concentration (Feng et al. 2001; Dong et al. 2009 ). The correlation of meteorological elements with the concentration of atmospheric NO 2 can be modeled using the multivariate data technique, represented by the eigenvectors of Table 6 , avoiding multi-collinearity in regression models that would include several factors correlated not merely with their response of variable NO 2 , but also to each other through weather variables. Therefore, it is possible to estimate meteorological patterns and to establish scores for each month through meteorological variables, by recognizing higher concentrations of atmospheric NO 2 . Both the satellite and the terrestrial NO 2 data sets of measurements for MS-Brazil from 2005 to 2020 present basic levels during the rainy season, which occurs from October and March, and maximum levels during the dry season from April to August. The study shows high relative humidity in the MS from October to March, which contributes to the removal of pollutants such as NO 2 . Temperature is also an important factor during a few hours due to the photochemical effect that produce OH radicals, which reacts with NO 2 causing it to sink in terms of HNO 3 . Climatic seasonality is, therefore, an influencing factor in the atmospheric specifications of NO 2 . Additionally, biomass-burning activities, especially during the driest period of the year, contribute to the increase of atmospheric pressure followed by an increase of NO 2 concentration. Burning of areas for cleaning and subsequent implantation of agricultural crops is notable in South Mato Grosso. Even during the COVID-19 pandemic, the levels of pollutant were similar to those of previous years, with the exception of months in which there were more important restrictions of commercial and industrial activities. The largest tropospheric and total NO 2 concentrations occurred in industrialized and highly agglomerated regions. The highest Table 6 Parameters of the analysis of main components and the eigenvectors of the first two main components PC1 and PC2 = principal component 1 and 2, respectively, KMO = Kaiser-Meyer-Olkin parameter, NO 2 = concentration of NO 2 (10 15 molecules/cm 2 ), T max = maximum temperature (°C) and T min = minimum temperature (°C), RH = relative humidity (%) and Prec = precipitation (mm) Cities NO 2 rules were also observed in biomass burning regions, the trend is not explained as a whole. They reflect differences in trends of variations due to anthropogenic connections of NO x in the regions. It has been investigated stationarity in the regional NO 2 column time series and some trends have been found. It is especially noteworthy that trends of increasing NO 2 concentration occur in drier months of the year at some cities. The results demonstrate that the NO 2 of the tropospheric column increased by around 3% and the breaking point of the series varied from 2010 to 2014. This study successfully produced the analysis of the vertical column of tropospheric NO 2 with values that are smooth and consistent in terms of space and time. A multivariate approach using eigenvectors to predict the NO 2 rules has been applied. The main components within 89% variance of the explained data and the eigenvectors are coherent to climatic seasonality and the mandatory application of NO 2 . Satellite data can be used to assess spatial and temporal variations of NO 2 . The integration of satellite inclusion with terrestrial data provides an overview of NO 2 in the atmosphere, which is an essential tool for quality managers and politicians to understand the spatial distribution of air quantities to define air quality policy and mitigation plan for South Mato Grosso. The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00703-021-00855-5. Analysis of monthly and annual rainfall variability using linear models in the state of South Mato Grosso Köppen's climate classification map for Brazil Trends and Variability in Extreme Rainfall Events in British Columbia Spatial-temporal variations in NO 2 and PM 2.5 over the Chengdu-Chongqing Economic Zone in China during 2005-2015 based on satellite remote sensing Reductions in nitrogen oxides over Europe driven by environ-mental policy and economic recession Overview of fire foci causes and locations in Brazil based on meteorological satellite data from Ozone monitoring instrument calibration Influence of precipitation on air quality in several cities of China MAX-DOAS NO 2 observations over Guangzhou, China; ground-based and satellite comparisons The observed response of Ozone Monitoring Instrument (OMI) NO 2 columns to NO x emission controls on power plants in the United States A comprehensive analysis of the spatio-temporal variation of urban air pollution in China during The role of primary emission and transboundary transport in the air quality changes during and after the COVID-19 lockdown in China Spatio-temporal variations of ozone and nitrogen dioxide concentrations under urban trees and in a nearby open area Chemical composition of precipitation in Beijing area Assessing the relative contributions of transport efficiency and scavenging to seasonal variability in Arctic aerosol Inference of surface concentrations of nitrogen dioxide (NO 2 ) in Colombia from tropospheric columns of the ozone measurement instrument (OMI) High resolution mapping of nitrogen dioxide with TROPOMI: first results and validation over the Canadian oil sands The spatial variation of O 3 , NO, NO 2 and NO x and the relation between them in two Swedish cities Quantitative estimation of meteorological impacts and the COVID-19 lockdown reductions on NO 2 and PM 2.5 over the Beijing area using Generalized Additive Models (GAM) Turnaround of tropospheric nitrogen di-oxide pollution trends in China Toward understanding the variation of air quality based on a comprehensive analysis in hebei province under the influence of COVID-19 lockdown Rank correlation methods Rainfall variability on the southern slope of the Bambouto mountain (West-Cameroon) and impact on the crop cultivation calendar Detection of trends and seasonal variation in tropospheric nitrogen dioxide over Pakistan Indirect validation of tropospheric nitrogen dioxide retrieved from the OMI satellite instrument: insight into the seasonal variation of nitrogen oxides at northern Midlatitudes Abrupt recent trend changes in atmospheric nitrogen dioxide over the Middle East The ozone monitoring instrument Nonparametric tests against trend NO 2 levels as a contributing factor to COVID-19 deaths: the first empirical estimate of threshold values COVID-19 pandemic: impacts on the air quality during the partial lockdown in São Paulo state Fire foci related to rainfall and biomes of the state of South Mato Grosso Trends and seasonal variability of atmospheric NO 2 and HNO 3 concentrations across three major African biomes inferred from long-term series of ground-based and satellite measurements A non-parametric approach to the change-point problem Increase in tropospheric nitrogen dioxide over China observed from space Trends analysis of rainfall and rainfall extremes in Sarawak, Malaysia using modified Mann-Kendall test A global single-sensor analysis of 2002-2011 tropospheric nitrogen dioxide trends observed from space Estimates of the regression coefficient based on Kendall's tau Regression modelling of hourly NO x and NO 2 concentrations in urban air in London Persistent fire foci in all biomes undermine the Paris Agreement in Brazil Análise de Agrupamento da Precipitação e da Temperatura no South Mato Grosso Estudo da associação entre óxidos de nitrogênio e concentração de ozônio com parâmetros meteorológicos Análise da Relação Entre O 3 , NO e NO 2 Usando Técnicas de Regressão Linear Múltipla Cluster analysis applied to the spatial and temporal variability of monthly rainfall in South Mato Grosso State Assessing the relationship between sugarcane expansion and human development at the municipal level: a case study of South Mato Grosso Improvements to the OMI near-UV aerosol algorithm using A-train CALIOP and AIRS observations Modeling of tropospheric NO 2 column over different climatic zones and land use/land cover types in South Asia Análise de tendência em séries históricas de vazão e precipitação: uso de teste estatístico não paramétrico Impacts of COVID-19 lockdown, Spring Festival and meteorology on the NO 2 variations in early 2020 over China based on in-situ observations, satellite retrievals and model simulations NO 2 , O 3 , and CO) in the inland basin city of Chengdu Spatio-temporal variations in NO 2 and SO 2 over Shanghai and Chongming Eco-Island measured by Ozone Monitoring Instrument (OMI) during 2008-2017 Spatiotemporal variations of tropospheric column nitrogen dioxide over Jing-Jin-Ji during the past decade Spatio-temporal Variations of tropospheric nitrogen dioxide in Turkey based on satellite remote sensing The Mann-Kendall test modified by effective sample size to detect trend in serially correlated hydrological series Spatiotemporal analysis of particulate matter, sulfur dioxide and carbon monoxide concentrations over the city of Rio de Janeiro, Brazil Assessment of the variability of pollutants concentration over the metropolitan area of São Paulo, Brazil, using the wavelet transform Impacts of the meteorological condition versus emissions reduction on the PM2.5 concentration over Beijing-Tianjin-Hebei during the COVID-19 lockdown Impact of precipitation with different intensity on PM 2.5 over typical regions of China Spatial and temporal distribution of NO 2 and SO 2 in Inner Mongolia urban agglomeration obtained from satellite remote sensing and ground observations Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations The authors would like to thank the universities The authors declare no conflicts of interest.