Estimation of the desert dust balance and its relationship with environmental factors in the southern Baja California Peninsula

The emission, transport, and deposition of desert dust are influenced by environmental factors evaluated mainly as a function of precipitation, air temperature, and wind. The present study estimates the desert dust balance and its relationship with environmental factors for the period 1981–2020 in the Baja California Sur region, Mexico. Monthly data on wind, air temperature, precipitation, dust emission and deposition from in situ measurements, reanalysis, satellite estimates and a numerical parameterization scheme are used. First, the in-situ time series were reconstructed with a non-linear principal component analysis based on an autoassociative neural network. We then evaluated the performance of the satellite estimates and the reanalysis data using the matching technique, after which we quantified the dust emission in three source zones. To estimate the state of the climate, we obtained the annual cycle of monthly climatologies, and yearly averages. The effect of climate on the dust budget is estimated by means of non-parametric associations between the variables of interest. Finally, we estimate climate trends using Pettitt, Modified Man-Kendall and Theil-Sen tests. The results show that the data sets used represent the mean and variability of in situ measurements. There are three seasons with transitional periods for precipitation and air temperature. The wind shows a longitudinal west–east gradient. Its maximum intensity is from April to September. Northern zone was dominated by highest emissions. Dry deposition mechanism was the most representative. Emissions and dry deposition are controlled by wind intensity and precipitation, while precipitation controls wet deposition. The Pettitt test showed that abrupt changes in the time series coincide with periods where positive sunspot anomalies converge with negative PDO and SOI anomalies. Finally, trend analyses identified a temporal broadening of the dry season, warmer winters, and summers with increases of ~ 0.02 °C, predominantly increasing trends in dust emissions and wind intensity, and a progressive decrease in dust deposition mechanisms.


Introduction
The analysis of climatic conditions in different regions of the world has been established mainly in terms of air temperature, precipitation, and wind (e.g., Aliaga et al. 2017;Almazroui et al. 2020;Jardim et al. 2021). Consequently, the Intergovernmental Panel on Climate Change (IPCC) estimates, in its different scenarios, increases in the probability of adverse events such as droughts, forest fires, intensification of tropical cyclones, and floods . These increases are mainly reflected in gradual increases in mean annual air temperature and precipitation by the end of the Communicated by H. Babaie.
1 3 twentieth century (IPCC 2015). Proof of this is the increasing trends of rainfall that have recently been reported in some regions of the world such as Central Africa (Almazroui et al. 2020), the northwestern United States with 3% increases in rainfall's extremes (Nazarian et al. 2022), and South America (southern Brazil, Uruguay and northeastern Argentina), where the main contributors are natural climate variability, variations in the composition of the atmosphere, and radiative forcing associated with anthropogenic activities (de Barros et al. 2016).
One of the elements resulting from natural and anthropogenic processes, which influence the climate system, is the dust concentration emitted from deserts. These impact the Earth's energy balance through interactions with radiation, clouds, atmospheric chemistry, cryosphere, and biogeochemistry (Kok et al. 2022). For mineral dust particles to reach impact in this way they must complete a cycle, this one divided into three phases: emission, transport, and deposition (Shao et al. 2011). Emission is the result of wind interaction with the earth's surface. While transport depends exclusively on particle size, anthropogenic activities (such as land use change), and wind intensity. Deposition, on the other hand, is established by two mechanisms: dry and wet deposition, which are mainly associated with fluctuations in wind intensity and rainfall, respectively (Bergametti & Forêt, 2014).
Once desert dust is deposited to the ocean could influence on marine primary productivity (Mahowald et al. 2014;Muñoz-Barbosa et al. 2020;Wang et al. 2017), and cause cooling at the ocean Surface (Evan et al. 2009;Martínez Avellaneda 2010;Safaierad et al. 2020). Examples of this are the effects of Saharan intrusions in the Caribbean Sea and the Gulf of Mexico, reported by Martínez Avellaneda (2010). Studies on desert dust trends indicate increases in dust emission from different regions of the world that are congruent with variations in precipitation, base soil and surface wind speed projected by IPCC (Ji et al. 2018;Okin & Reheis 2002;Pu & Ginoux 2017).
The northwestern region of Mexico has a big land extension surrounded by the Gulf of California and the Pacific Ocean. It is divided in two states: Baja California and Baja California Sur. The last one has been characterized as a relevant source of desert dust (Morales-Acuña et al. 2019b), and due to its climatic characteristics, land use and management, it is vulnerable to the effects of climate change, as evidenced by the intensification of desertification processes (Salinas-Zavala et al. 2017;SEMARNAT 2016;Troyo-Diéguez et al. 2014  S are the North, Central, and South regions that act as dust sources, according to Martínez-Flores et al. (2013) These scenarios allow inferring that the desert dust inventory in this region could be increasing, driving impacts on public health, aquatic ecosystems, and climate. For this reason, this work seeks to estimate the climate trends in Baja California Sur and its association with the inventory of desert dust (from January 1981 to December 2020) emitted from point sources distributed in this region.

Study area
The southern region of the Baja California Peninsula (RSBCP) (Fig. 1), located between 22°-24° N and 115°-109° W, has a mountain range with a maximum elevation of 2160 m and is surrounded by the Pacific Ocean to the east and the Gulf of California to the west (Fig. 1). The pluviometric regime for this region establishes two subregions whose differences are framed by the presence of the Sierras and the plains of the region, such that the maximum average annual rainfall is found in the south (310 mm year −1 ) and the minimum in the north (120 mm year −1 ) (Salinas-Zavala et al. 1990;Troyo-Diéguez et al. 2014). The mean annual pattern is distributed in three seasons; two rainy seasons (July-October, November-February) and one dry season (March-June) (Salinas-Zavala et al. 1990).
The RSBCP during the summer the temperature can exceed 40 °C, while in the winter, it oscillates between 5 and 12 °C, excluding the Sierra de la Laguna region . The wind pattern during the summer establishes predominantly westerly direction with intensities < 3 m s −1 (Morales-Acuña et al. 2019a, b); in the other climatic seasons, the wind direction obeys the Pacific wind pattern (Morales-Acuña et al. 2019b).
The arid and semi-arid characteristics of the region allow the identification of extensive areas that Martínez-Flores et al. (2013) classified as a dust source, located mainly to the north of this region, and distributed southward along the coastal area adjacent to the Pacific Ocean. Regarding the predominant vegetation, 83% of the vegetation cover in this region consists of tropical dry forest and xerophytic scrub. Additionally, the zonal vegetation in coastal areas corresponds to sclerophyllous chaparral and succulent scrub along the Thermo-Mediterranean and Inframediterranean belts, respectively, and to different types of thermo-tropical-desert vegetation in the rest of the peninsula (Peinado et al. 2008). In this region of Mexico, it has been determined that the most predominant soil types are Regosol (39.99%) and Leptosol (29.71%). In addition, Arenosol participates in this classification with 2.12% (INEGI 2010).

Land-based rainfall and air temperature measurements
Daily mean precipitation and air temperature data for the period 1981 to 2015 from 10 weather stations distributed as shown in Fig. 1 were obtained from the CLImate COMputing (CLICOM) project available at http:// clicom-mex. cicese. mx/. CLICOM, a climatological data management software system developed by the United Nations, uses a database of surface weather stations from the National Meteorological System of Mexico (S.M.N, 2023), which also performs quality control of the information provided in their data portals. The location of the weather stations used in this study is shown in Table 1.

CHIRPS (Climate Hazards Group InfraRed Rainfall with
Station data) is an estimated rainfall database created and managed by the U.S. Geological Survey (USGS) and the University of California, Santa Barbara (UCSB). The information is distributed in NetCDF4 format with different temporal resolutions, from daily to decades. The spatial resolution of these data is 0.05°. CHIRPS are based on data from the Climate Hazards Rainfall Climatology-(CHPClim), Geostationary Thermal Infrared (IR) sensor, Tropical Rainfall Measuring Mission (TRMM) satellite mission, NOAA Climate Forecast System (CFSv2) atmospheric rainfall field model, and in situ rainfall observations from various sources including U.S. and regional meteorological services (Funk et al. 2015).

Reanalysis products
Monthly data (1981-2020) at 0.1° resolution for the zonal and meridional wind components, and air temperature, were obtained from ERA5-Land (https:// cds. clima te. coper nicus. eu). ERA5-Land is a reanalysis data set that provides a consistent view of the evolution of ground-based variables over several decades with improved resolution compared to ERA5. ERA5-Land has been produced by reproducing the ground-based component of the ECMWF ERA5 climate reanalysis. The reanalysis combines model data with observations from around the world into a globally consistent and complete data set. The approximate monthly time series of area-averaged dry and wet dust deposition for the marine region of the La Paz Bay (LPB) for the period 2003-2020 were obtained from https:// giova nni. gsfc. nasa. gov/; these global data are the product of a MERRA-2 model with a 0.5 × 0.6 • of spatial resolution.

Pre-processing of databases
An algorithm developed in the MATLAB R2022a programming language was used to search for missing data in the time series obtained from the CLICOM database. The series with missing data were reconstructed using a Non Linear Principal Component Analysis (NLPCA) with inverse network architecture and a total of 3000 iterations based on an autoassociative neural network proposed by Scholz et al. (2008). This technique is available at http:// www. nlpca. org/ matlab. htm and has been successfully used in the reconstruction of hydrological time series (Canchala et al. 2020;Miró et al. 2017). In this study, standardized spatial anomalies were initially calculated; subsequently, the NLPCA analysis was applied, obtaining reconstructed spatial anomaly series for each analyzed variable. Finally, the reconstructed series were obtained by multiplying the reconstructed standardized spatial anomaly series by the standard deviation of the original data, and the arithmetic mean of these was summed, as shown in Eq. 1.
where S rc is the reconstructed series, x rc is the reconstructed standardized spatial anomaly series. i y x i , are the standard deviation and mean of the in situ data provided by CLICOM, respectively. To evaluate the model performance and forecasting ability at both training and test time stamps, we used two statistical methods: root means square error (RMSE) and Spearman correlation coefficients ( r Spearman ) , applying Eqs.
(2) -(3), respectively. (1) where S i and S rci are the range of the measurement taken on the i-th element of the series, S y S rc represent the average ranges, respectively. The data from the different sources (reanalysis and satellite) were spatially trimmed, considering the RSBCP location, and the matching pixels with this region. So that a set of geospatial information was obtained and concatenated into three-dimensional matrices: longitude, latitude coordinates, and time. After selecting matching pixels with the study area, the erroneous data in the matrices were replaced by non-existent data (NaN). With the classified monthly arrays, multiyear series were obtained. For each variable, maps of the multi-annual monthly cycle and the annual mean were presented to visualize each variable's spatial distribution in the RSBCP. Climatology of cumulative monthly rainfall and total annual rainfall for 1981-2020 was obtained.

Rainfall and air temperature data validation
The evaluation of CHIRPS and ERA5-Land products was carried out from 1981 to 2015, using data from 10 weather stations located at different points in the RSBCP. First, grid values of rainfall and air temperature estimates for each rainfall station coordinate were extracted. Subsequently, comparative metrics were evaluated (Table 2). Within the comparative metrics, Spearman's correlation coefficient (ρ), mean absolute error (MAE), and bias (B) were considered (See, for example, Alemu & Bawoke 2019; Dinku et al. 2018;Urrea et al. 2016).
D is the difference between two samples ranges, and N is the number of ranges. C i is the amount of rainfall or air temperature derived from weather stations. S i represents the satellite-derived rainfall estimate or ERA5-Land air temperature, as appropriate. n is the total number of data.

Quantification of dust emission
The numerical parameterization scheme for dust emission E mgm −2 s −1 as a function of soil characteristics and grain size Table 2 Comparison metrics were used to evaluate monthly air temperature and rainfall data from CHIRPS and ERA5-Land for 1981-2020 Shao et al. (1993), modified by Nickovic et al. (2001) and adjusted for Baja California Peninsula by Morales-Acuña et al. (2019b), will be used to calculate dust emission from three locations acting as uniformly distributed dust sources as shown in Table 3.
where C(mgm −5 s 2 ) is an empirical dimensional constant (C = 1.7195mgm −5 s 2 , Morales-Acuña et al. 2019b). t * and * ( ms −1 ) are the frictional threshold velocity ( t * = 0.25ms −1 Ginoux 2017; Ginoux et al. 2012;Nickling & Gillies 1993;Westphal et al. 1987) and the wind friction velocity which was calculated for conditions of unstable atmospheric stability as proposed by Morales-Acuña et al. (2019b). is a dimensionless factor referring to dust productivity (Eq. 5) that depends on vegetation, the fraction of the area susceptible to erosion and types of soil, and the relative contribution of each of these factors.
where is the vegetation factor defined as: M is the mask of the desert as a function of vegetation, and N is the number of pixels in the grid of the mask. In the present work α = 1 (Nickovic et al. 2001), because our study area is represented as a single pixel for each point source; vegetation can be neglected because this area has been considered to be sand desert with scarce vegetation (Morales-Acuña, 2015). Tegen & Fung (1994) assume an erodible fraction k (Eq. 5) as a function of particle type. In our calculations we specify k = 0.17 . The soil factor and relative contribution k (Eq. 5) will be considered according to the predominant granulometry in the study areas and to the particle size whose time-life in the atmosphere facilitates long-distance transport, i.e., the relative contribution of sands and clays, so that k = 0.25 (Table 2. Nickovic et al. 2001). After calculating each value in Eq. 5, the dust productivity factor in this work is set as = 0.0425.

Associations of dust emissions with hydrometeorological variables
Once the monthly arrays of each of the hydrometeorological variables were obtained, the nearest neighbor technique was applied to extract from each matrix the monthly time series of the matching pixel with the geographic location of the dust source area in which the emission was calculated. Subsequently, statistical associations were made with Spearman's coefficient as shown in the following section.

Temporary trends
The nature of time series can have abrupt changes or breaks in which the behavior can change to increase or decrease from this. Considering the above, the different variables analyzed in this study will be evaluated using Pettitt's test (Pettitt 1979). This test approximates a sequence of random variables belonging to non-parametric methods. It is used to identify a change point in a time series and can be described as follows: Once the breakpoint of a time series is detected, the data set is divided into two intervals by classifying the samples coming from the same population into before x 1 , x 2 , ...x t and after x t+1 , x t+2 , ...x T the breaking point. The U t,T test statistic is given by (Eq. 7): The most significant change point is obtained when K T = max | | U t,T | | , the significance level is determined in our work for < 0.05 , in such a case the null hypothesis is rejected, this value is calculated by the expression = exp −6K 2 T T 3 +T 2 . If there is a significant change point, the time series is divided into two parts at the location of the change point.
Once the change points of the time series have been obtained, the trend of the time series of the different variables will be evaluated. First, the Modified Mann-Kendall non-parametric test proposed by Hamed & Rao (1998) was applied for each pixel of the matrix belonging to the study area. The Mann-Kendall test obtains the S-statistic (Eq. 8), considering the length n of a time series data and the ele- where, If the data are identical and have an independent distribution, with a mean value of S equal to zero, the variance will be given by (Eq. 9): The statistical value of the Mann-Kendall test is obtained by the following expression: The hypothesis test is with H0 indicating no trend in the time series while H1 indicates a trend. H0 is rejected with a confidence level of 95% if |z| > 1.96 . Suppose the coefficients resulting from this test have autocorrelation. In that case, the trend values may be over or underestimated (Hamed & Rao 1998), so the modified Mann Kendall test considers the structure of the significant autocorrelation coefficients in the time series data. This is achieved by modifying the variance V used to calculate the z-value of the original Mann-Kendall method such that: where r i is the significant autocorrelation coefficient at the i-th order of the data and n * is represented as an effective number of observations to account for the autocorrelation in the data. Sen's slope estimator was used to calculate the actual slope of trends in cumulative rainfall time series data. If a linear trend exists, the magnitude of the monotonic trend in hydrologic time series can be quantified using Sen (1968) nonparametric slope estimator using Eq. 12.
represents the slope mean values between data measurements x j and x i at time steps jei con j > i , respectively. The positive value of indicates an increasing trend while the negative value of indicates a decreasing trend. The sign of reflects the direction of the trend of the data, while its value indicates the slope of the trend. The advantage of this method is that it limits the influence of missing values or outliers on the slope compared to linear regression.

Reconstruction and validation of databases
The quality control for the CLICOM time series identified that only the air temperature series had discontinuities. These were reconstructed at 95% confidence level, a correlation coefficient of 0.92, and a mean square error of 4.6 × 10 −5 , on average for ten selected stations.
For precipitation, the validation metrics allowed us to identify (with a 95% confidence level) that the CHIRPS-v2.0 data slightly underestimate the rainfall values from the weather stations, with a mean bias of 0.95 (Fig. 2a). This behavior is mainly observed at the La Paz station. On the other hand, the monthly variability of the CHIRPS-v2.0 data adequately represents that of the rain gauges, such that Spearman's correlation coefficient has a mean value of 0.8.
Comparisons between the ERA5-Land air temperature databases with those provided by CLICOM (Fig. 3) identified an average correlation above 0.9 (95% confidence) and an average bias of 1.005. At two stations, the ERA5-Land temperature values slightly overestimate those in situ stations. The average MAE obtained (1.54 °C) shows that the absolute difference between the two measurements leads to minor errors.

Climatology of the hydrometeorological variables
The climatology of the accumulated monthly rainfall (Fig. 4) allowed us to identify the existence of two rainy seasons, one dry season and two transition months explained below: (1) High rainfall (> 110 mm month −1 ), comprised from July to September, (2) low rainfall (50-100 mm month −1 ), during December and January, (3) Dry, with rainfall less than 50 mm month −1 from March to June. February and October are transition months between one season and the other. On the other hand, total annual rainfall presented maximum values in the south (> 300 mm month −1 ) and minimum values (< 200 mm month −1 ) in the center and north of the region. The monthly climatology of air temperature (Fig. 5) identified July-September as the months with the highest temperatures (> 28 °C), excluding the highland region located in the south of the study area. The lowest temperatures (< 20 °C) occurred during December-March. From April to June, temperatures are in the range 22-27 °C. November seems to be a transition month.
The mean annual temperature in the study area presents a west-east longitudinal gradient with lower temperatures (< 21 °C) in the southern and northeastern regions. The highest temperatures (> 24 °C) occurred in the coastal region adjacent to the CG.
The analysis of the wind fields for the annual cycle (Fig. 6) established a longitudinal gradient in wind intensity, with higher intensities (> 3 ms −1 ) on the Pacific Ocean coasts and lower intensities (< 1.5 ms −1 ) in the coastal region of the GC. On the other hand, it is observed (Fig. 6) that in the areas of high topography, the wind intensity has lower values concerning the plains regions. The predominance of maximum values for wind intensity occurs from April to September, being the June-August quarter the one with the highest intensities.
From April to September, the predominant wind direction is from the west, with a characteristic behavior in which the air masses that encounter the mountain massif will re-direct themselves taking a northeasterly direction and, in some cases, bordering the mountain. From November to February, the predominant wind is from the north, while October, March, and May, the predominant direction is from the northwest. The mean annual wind field (Fig. 6, left panel) allows us to identify a predominant direction from the northwest and intensities of 1 to 4 ms −1 distributed spatially along a longitudinal west-east gradient. To the south of the study region, the lowest intensities are found just in the mountain zone.

Desert dust budget
The monthly dust climatology emitted from the three source regions in the (Fig. 7) allow identifying a unimodal behavior with maxima in May and minima in December, January, and November, for the North, Central, and South region, respectively. Dust emission rates from these three regions ranged from 0.54 y 466.42 mgm −2 s −1 . The ). The boxplot diagrams show that the highest dispersion of the data is had during summer, and the lowest is observed in winter.
The mechanisms of desert dust deposition on LPB during the period 2003-2020 (Fig. 8) allowed use to establish that dry deposition has the maximum values concerning wet deposition in such a way that it exceeds an average of 0.1661 mgm −2 s −1 to wet deposition. Referring to the annual cycle of the monthly climatology of dry deposition (Fig. 8), a unimodal behavior can be observed whose maxima seem to be associated with the seasons of maximum dust emission obtained in the three studied localities. The On the other hand, the wet deposition of desert dust in the LPB (Fig. 8) allowed identifying a unimodal behavior with maxima during the July to October seasons; this season is associated with the arrival of tropical cyclones with heavy rainfall as reported by Cavazos et al. (2008). Therefore, the minimum (0.0001 mgm −2 day −1 ) and maximum (0.006 mgm −2 day −1 ), values for this desert dust deposition mechanism are established in May and July, respectively.

Associations of desert dust with air temperature, wind, and rainfall
The leading associations between air temperature (at 2 m) and rainfall were found in the northern region and the minimum associations in the central region (Table 4). In addition, it can be observed that dust emission is directly associated with wind and inversely with rainfall. Air temperature presented low associations; for the case of the northern region, the associations were not statistically significant.
The dry deposition mechanism was directly related to wind and inversely related to rainfall, while wet deposition was inversely related to wind, with rainfall dominating. Air temperature presented statistically significant associations only with the wet deposition mechanism (Table 5).

Temporal trends of air temperature, rainfall, and desert dust budget
After evaluating the breakpoints in each rainfall time series, monthly air temperature, wind, and dust emission at 95% confidence level, we found the year 2000 as this breakpoint in all these variables. Both dry and wet deposition mechanisms, the breaking point was 2011.
The trends for the annual rainfall cycle (Figs. 9 and 10) allow identifying a decreasing behavior with values between 0.05 and 0.3 mmyear −1 (at 95% confidence level). The high rainfall season contains the most significant decreases, with September being the most representative month in this behavior. On the other hand, the low rainfall season shows an average decrease of 0.2 mmyear −1 in almost every region.
Months with the most significant temperature increases in almost the entire region were March and November, with the former showing the highest trends (> 0.025 °C year −1 ). A characteristic feature is during the months with lower temperature values, the trends are positive (Fig. 5, December-March). The same behavior occurs during the summer months (Fig. 5, June-August).   The trend analysis for the annual cycle of wind intensity (Fig. 11) allowed identifying increases from 0.008 a 0.02 ms −1 year −1 . Increases in wind intensity were predominant in regions of low orography (< 200 m). In July, this increase could be observed in most the region; in the rest of the months, these increases were mainly located in the north and south. In December, the increases occurred in the central region adjacent to the Pacific Ocean.
The analysis of the temporal trends in the time series of dust emission in the three selected zones in the RSBCP (Table 6) for the 1981-1999 and 2001-2020 intervals shows that dust emission has a predominance of positive trends in the Northern region and a predominance of negative trends in the Central and Southern regions. Increases occurred mainly in May, July, October, and November, while decreases occurred in March, April, June, and August.
The trend analysis for desert dust deposition in LPB during 2003-2020, at 95% confidence level before and after the breakpoint (Table 7), allowed establishing an average decrease of −7 × 10 −4 mgm −2 day −1 year −1 for the wet deposition mechanism and −0.01mgm −2 day −1 year −1 for dry deposition. Fig. 11 Temporal trends of the annual cycle of wind intensity for 1981-2020, at 95% confidence level

Discussion
The inverse NLPCA approach applied to the time series with missing data from the CLICOM network allowed us to identify that the extraction of the nonlinear components directly from incomplete data and their subsequent nonlinear mapping reconstructed the time series, maintaining the variability and values of the original series, verifying the above with the comparative metrics proposed by Scholz et al. (2008). The comparison metrics obtained between CHIRPS-v2.0 data and CLICOM network rain gauges used in this work (Fig. 2)  The spatial distribution of rainfall establishes a latitudinal gradient from south to north, due to the pluviometric regime reported by Salinas-Zavala et al. (1990) and Troyo-Diéguez et al. (2014), who explains these differences due geographical differences. The rainfall climatology allowed classify two rainy seasons and one dry season that coincide with that reported by Salinas-Zavala et al. (1990), but the present work determined two transition months: February and October. The rainfall season and maximum rainfall are in the mountain regions located south of the study area. This season coincides with the tropical cyclone season and the significantly increasing linear trends of heavy rainfall in mountain sites reported by Cavazos et al. (2008). This rainfall contributes 40% of the total annual rainfall for the Baja California Peninsula (Breña-Naranjo et al. 2015). The climatology for air temperature found in the present study (Fig. 5) allowed the identification of warm and cold seasons, consisting of the months of summer (June-September) and winter (December-March), congruent with Troyo-Diéguez et al. (2014) andINEGI (2011). In addition, the month of November is proposed as a transition period between maximum and minimum temperatures.
The longitudinal west-east gradient for the wind intensity found in this study (Fig. 7) can be attributed to the influence of northern winds in the eastern coastal region of the study area, as reported by Morales-Acuña et al. (2019a, b) and Castro & Martínez (2010). Also, the predominance in the wind direction during the summer season (Fig. 8) agrees with the west-east direction reported by Morales-Acuña et al. (2019a, b). The redirection and behavior that air masses have when encountering the mountain located south of the region from April to September (Fig. 8) could be attributed to Froude number values < 1 (Nappo 2013), such that, when this condition is met, atmospheric flows tend to go around obstacles found in their trajectory. The predominant southwesterly direction in the RSBCP during the winter months (November-February) responds to the seasonal winter pattern and the displacement of the North Pacific high-pressure center, which in this season is directed towards the equator (Badan-Dangon et al. 1991;Morales-Acuña et al. 2019b;Romero-Centeno et al. 2007).
The analysis and quantification of emission and deposition (Figs. 7 and 8), showed that in the three source regions studied, there is a unimodal behavior in the emission rates where the maximum emissions are during spring-summer and minimum in winter, coinciding with those reported by Morales-Acuña et al. (2019b), for the Vizcaíno desert. If we consider the annual cycle fluctuations of the hydrometeorological variables analyzed, it can be suggested that this behavior of maximum emission rates responds to favorable conditions that facilitate the atmospheric flows of desert dust from each of the regions (North, Central, and South) to the LPB, mainly during March to June, due to the dry season, average temperatures, and winds with predominantly easterly direction.
Regarding deposition, although the wet deposition mechanism plays an essential role in the supply of nutrients and toxic substances to terrestrial and aquatic environments (Pan et al. 2017). In this study, it was determined that the primary mechanism for carrying out desert dust intrusions resulting from emissions and dispersion from the analyzed source areas to the LPB is through dry deposition, being 0.1661 mgm −2 s −1 higher than wet deposition, this behavior may be due to the climatic characteristics in the RSBCP as a semidesert area where most of the year the accumulated rainfall is less than 500 mmmonth −1 (Fig. 4), and the wind shows optimal conditions of intensity, 75% of the year (Fig. 6), which would favor erosion and subsequent dust deposition. In addition, Zhao et al. (2003) suggest that dry deposition is a dominant dust removal process near the source areas (as is the case of the LPB), and the removal of dust particles by rainfall turns out to be the primary process in the transport pathway in regions distant from the source. Comparing the dust fluxes emitted from each source with the atmospheric depositions in the LPB, we find very marked differences suggesting that not everything emitted from the source regions tends to be deposited in the LPB. This behavior is comparable with the balances of desert dust fluxes obtained from models. In many world regions, dust fluxes are not in balance because the transport models generally do not conserve mass due to the time interpolation techniques that they usually use (Ginoux et al. 2004). Additionally, if we compare the normalized values by area obtained in the present study with the emission and deposition of dust in other regions of the world (Table. 8); it is inferred that: RSBCP and LPB, are not very relevant to the global desert dust budget, but a very significant local contribution for the region, as is the case in Japan and its coasts (Osada et al. 2014).
The analysis of the breaking point in the monthly climatology of time series identified an abrupt change in the year 2000 for the variable's rainfall, air temperature, wind, and dust emission, while for the mechanisms of deposition of airborne dust in the LPB, the change in the series occurred in 2011. By repeating this analysis to the complete series (at 95% confidence level), it was possible to determine that in both cases, the change took place precisely in December of the corresponding year. Taking into account that fluctuations in the time series of atmospheric, hydrological, and oceanographic variables generally depend on atmospheric phenomena such as El Niño Southern Oscillation (ENSO), the Pacific Decadal Oscillation (PDO) (Alvarez-Olguin and Escalante-Sandoval 2017; Arriaga-Ramírez & Cavazos 2010;Llanes-Cárdenas et al. 2018) and solar activity, monthly time series of the indices, Southern Oscillation (SOI, https:// www. ncdc. noaa. gov/ telec onnec tions/ enso/ soi), PDO (https:// www. ncdc. noaa. gov/ telec onnec tions/ pdo/) and sunspots (Sunspot, https:// wwwbis. sidc. be/ silso/ month lyssn plot were obtained. These time series were smoothed with the Lanczos cosine filter (Duchon 1979) low pass, with a 12-month window (Fig. 12). It was found that in the periods corresponding to the breakpoints obtained in the present work, the solar activity presents positive anomalies (Fig. 12, insets). In addition, it could be identified that during this period of positive sunspot anomalies, the PDO and SOI present inverse anomalies between 1 and-2 standard deviations (Fig. 12), out of phase -3 months at 95% confidence level.
The progressive decreases in rainfall found in the present study (Fig. 9) suggest a temporal extension of the dry season, thus affecting the pluviometric regime of the region and establishing a unimodal behavior in which there would be a rainy season (June-September) and the rest of the dry season. These decreases are congruent with those reported by Blunden & Arndt (2016) and Morales-Acuña et al. (2019b) for this region of Mexico, even so, in an exercise for the detection of the trend between 1960 and 2010 carried out by Núñez-González (2020), in 6 stations located in the RSBCP determined the existence of no annual trend, negative trends for winter-summer and favorable trends for autumn. The increases in air temperature (Fig. 10) established during the winter and summer months (which in turn represent the seasons with the most remarkable differences), increases of ~ 0.02 °C that could affect the economic activities of the region, characterized by focusing part of its economy on agricultural activities with cyclical crops during the autumn-winter, spring-summer seasons. In addition, in terms of desert dust emissions, these increases in air temperature would favor fluxes to the atmosphere due to the lack of soil moisture resulting from soil evaporation processes. The temperature increases found in the present study are in agreement with those reported by the IPCC (2019), which estimates that global temperature will increase between 0.1 and 0.3 °C per decade, while seasons in some regions are experiencing warming; in addition, results reported by Blunden & Arndt (2016)  The predominance of increasing trends in dust emission for the northern region of the present study is congruent with the increases in this variable in the Vizcaíno desert reported by Morales-Acuña et al. (2019b). Furthermore, the months in which increases in dust emission occurred (May, July, October, and November) coincide with the increasing trends in wind intensity, air temperature, and decreases in rainfall obtained in the present study. In addition, the months with decreasing trends coincide with the predominance of no statistically significant trend for the climatic variables evaluated in the dust source regions. Finally, the decreases in the wet and dry deposition mechanisms obtained in this work are congruent and respond to the decreasing trends in rainfall and increases in wind intensity, which regulate these mechanisms.
Finally, considering that the direct and indirect radiative effects of terrigenous particles affect climate from regional to global scales (Myhre et al. 2014), specifically "direct effect" by absorbing and scattering solar radiation, desert dust particles reduce the amount of energy reaching the Earth's surface while enhancing the greenhouse effect by absorbing and emitting longwave radiation (Spyrou 2018). The intensification of rainfall in the southern mountainous region found in the present study from July to September (Fig. 4) could be related not only to the presence of tropical cyclones, which in turn have reduced their activity in the North Atlantic due to factors such as dust intrusion from the Sahara (Strong et al. 2018), but also the change in radiative forcing and its impact on local climate parameters such as precipitation and air temperature (Liu et al. 2014;Rap et al. 2013). Specifically for precipitation, studies have shown that desert dust inputs can affect regional water resources by modulating rainfall distribution (Jha et al. 2021), such that these dust inputs sometimes enhance heavy rainfall events and suppress light rainfall events (Alizadeh-Choobari 2018; Li et al. 2011). Consistent with the negative trends in precipitation found in the present work, strong atmospheric inputs of desert dust in turn reduce precipitation by drying out the soil and producing more dust, providing a possible feedback loop to further reduce precipitation (Rosenfeld et al. 2001). On the other hand, the positive air temperature trends obtained in the present work (Fig. 5) could be related to the reduction of incoming solar radiation at the surface (Slingo et al. 2006), he increase in daily heating rates  (Chen et al. 2017) and the warming of the middle troposphere and cooling of the lower troposphere (Mohalfi et al. 1998), resulting from atmospheric intrusions of desert dust. In addition, the sixth leading cause of mortality in the state of Baja California Sur is acute respiratory infections, with a rate of 1.9 predicted by the Ministry of Health in 2018 for 2050, and it is expected that one of the possible causes of this public health problem is the inhalation of dust particles with diameters less than or equal to 10 microns, as has been reported in southern Europe (Stafoggia et al. 2016), Madrid (Diaz-Jiménez et al. 2012, Cuba (Venero-Fernández 2016), among others. Also for larger sizes of dust particles, it has been found that they can potentially damage external organs, irritating eyes and ears (Zhang et al. 2016). This raises the alarm in the study region for the development of studies and action plans focused on the effects of atmospheric dust intrusion on public health.

Conclusions
The time series reconstruction techniques applied to the CLICOM database and validation of the reanalysis and remote sensing products used in this work showed statistically significant metrics that guaranteed the use of the data obtained for the study region after the pre-processing to which they were subjected.
Although the analysis of climate in the RSBCP (1981-2020) temporally and spatially classifies rainfall, air temperature, and wind within the annual cycle, trend analyses determine that this classification is affected by increases in air temperature and wind intensity and decreases in rainfall. It causes the wind (predominant direction is towards the east and southeast) to transport moisture away from the peninsular region, and in turn, it generates more extended droughts, warmer winters, and hotter summers.
Since the emission of desert dust is driven mainly by the wind intensity and rainfall, the primary mechanism of deposition of this terrigenous material on the LPB is dry deposition. It is deduced that the fluctuations and trends of the hydrometeorological variables affect the production of desert dust, thus finding a predominance of positive trends in the region of lower rainfall (North) and negative trends in the opposite case (South). In the months where decreasing trends in rainfall, increases in air temperature, and wind intensity converge, dust emissions from the different zones were increasing. However, dust deposition over the LPB has a decreasing trend because it is controlled mainly by fluctuations in wind intensity, so it is suggested that due to increases in wind intensity and the predominance in the direction towards the east and southeast, most of the dust emitted from the peninsula is deposited in the southern Gulf of California. Finally, the found convergence between the representative indexes (solar activity, the PDO and the SOI, the last two out of phase by -3 months), and the temporal breakpoints (different time series analyzed) represents particular interest to carry out studies on the implications of climate indexes on the trends and temporal fluctuations of hydrometeorological and environmental variables.