The Baja California Peninsula, a Signiﬁcant Source of Dust in Northwest Mexico

: Despite their impacts on ecosystems, climate, and human health, atmospheric emissions of mineral dust from deserts have been scarcely studied. This work estimated dust emission ﬂux (E) between 1979 and 2014 from two desert regions in the Baja California Peninsula (BCP) using a modiﬁed dust parameterization scheme. Subsequently, we evaluated the processes controlling the variability of E at intra- and interannual scales. During the period 1979–2014 peak E were generally recorded in summer (San Felipe) and spring (Vizcaino), and the lowest emissions occurred in autumn (San Felipe) and winter (Vizca í no). Intra- and interannual variability in E was associated with ﬂuctuations in wind speed and direction, precipitation, and soil moisture, which, in turn, were controlled by the seasonal displacement of the North Paciﬁc high-pressure center. Key drivers of the interannual variability of E are strong El Niño Southern Oscillation (ENSO) events. These climatic events and the hydrometeorological variables mentioned above played a major role in the onset and occurrence of dust events, with the highest annual emissions at Vizca í no. Besides, a lag of 19 months (San Felipe) and 21 months (Vizcaino) was recorded between the occurrence of relevant E and ENSO events, apparently in response to the e ﬀ ect of this climatic event on precipitation. The climate variability of E in both desert regions was evidenced by the positive trends associated with increases in wind speed and air temperature, and with decreases in precipitation and soil moisture. Finally, our ﬁndings suggest that the BCP should be considered as a signiﬁcant source of dust for the regional inventory of particulate matter emissions from the Earth’s surface.


Introduction
Emissions of mineral dust (MD) from surface soil have been described as resulting from three main mechanisms: (1) Aerodynamic lift, (2) saltation bombardment (sandblasting), and (3) disintegration of aggregates (auto-abrasion) [1,2]. These mechanisms are part of the wind erosion processes and are subject to weather conditions as well as to the properties, characteristics, and uses of land [1,3].

Climate
In general, the BCP is an arid area with mean annual precipitation of ~200 mm distributed across the latitudinal range and a bimodal pattern characterized by peaks in summer (in the south) and winter (in the north) [27]. San Felipe is one of the hottest and driest deserts in North America. Air summer temperature exceeds 45 °C and the total mean annual rainfall is less than 200 mm year −1 . In Vizcaino desert, the maximum temperature is 35 °C, with the rainy season between July and October, and minimum precipitation between December and February [27].
Given the importance of moisture and dust transport in the desert, studies on the circulation of wind in the BCP have focused mainly on the summer. These show that surface wind responds to: (1) The thermal contrast between the Pacific Ocean (PO) and the GC [21,26,28,29], (2) orography of the BCP [21,26,[30][31][32], the difference in surface air temperature between the margins of the BCP (peninsular PO margin and GC) [21], and the high-pressure belt in the northern hemisphere [24].

Dust Emission Flux from San Felipe and Vizcaino
The dust emission flux, E (mg m s ), as a function of the characteristics of soil and particle size distribution proposed by Shao et al. [33] and modified by Nickovic et al. [14] was used: For u * ≥ u * , where (mg m s ) is an empirical dimensional constant and was derived from the slope of the linear regression applied to the data of monthly atmospheric dust fluxes from January to December 2011 provided by the project CONACyT-N° 166897 (Further reference on the sampling design and data collection, refer to Muñoz-Barbosa et al. [20]). Using that information, =

Climate
In general, the BCP is an arid area with mean annual precipitation of~200 mm distributed across the latitudinal range and a bimodal pattern characterized by peaks in summer (in the south) and winter (in the north) [27]. San Felipe is one of the hottest and driest deserts in North America. Air summer temperature exceeds 45 • C and the total mean annual rainfall is less than 200 mm year −1 . In Vizcaino desert, the maximum temperature is 35 • C, with the rainy season between July and October, and minimum precipitation between December and February [27].
Given the importance of moisture and dust transport in the desert, studies on the circulation of wind in the BCP have focused mainly on the summer. These show that surface wind responds to: (1) The thermal contrast between the Pacific Ocean (PO) and the GC [21,26,28,29], (2) orography of the BCP [21,26,[30][31][32], the difference in surface air temperature between the margins of the BCP (peninsular PO margin and GC) [21], and the high-pressure belt in the northern hemisphere [24].

Dust Emission Flux from San Felipe and Vizcaino
The dust emission flux, E mg m −2 s −1 , as a function of the characteristics of soil and particle size distribution proposed by Shao et al. [33] and modified by Nickovic et al. [14] was used: For u * ≥ u t * , where C (mg m −5 s 2 ) is an empirical dimensional constant and was derived from the slope of the linear regression applied to the data of monthly atmospheric dust fluxes from January to December 2011 provided by the project CONACyT-N • 166897 (Further reference on the sampling design and data collection, refer to Muñoz-Barbosa et al. [20]). Using that information, C = 1.7195 mg m −5 s 2 . u t * and u * (m s −1 ) are threshold friction velocity (u t * = 0.25 m s −1 [5][6][7][8]) and friction velocity. δ is a dimensionless factor referring to dust productivity (Equation (2)) that depends on vegetation, Atmosphere 2019, 10, 582 4 of 18 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 ( Table 1 in Nickovic et al. [14]), because our study area is represented as a single pixel for each zone. Vegetation can be neglected because this area has been considered to be sand desert with scarce vegetation [26]. Tegen and Fung [34] assume a erodible fraction γ k (Equation (2)) as a function of particle type. In our calculations we specify γ k = 0.17. The soil factor and relative contribution β k (Equation (2)) 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. [14]). After calculating each value in Equation (2), the dust productivity factor in this work is set as δ = 0.0425.
where u(z) is wind velocity from the North American Regional Reanalysis (NARR) database (ftp: //ftp.cdc.noaa.gov/Datasets/NARR), κ = 0.41 is the Von Kármán's constant, L is Monin-Obukhov length. L < 0 is refered to unstable atmosphere conditions. z N = 0.02 m represents the height used for calculating the friction velocity and z 0 = 0.0002 m is the ground roughness for sand deserts reported by Wieringa [37]. Stability function, Ψ N , is determinated for L < = (unstable): The α 1 value used in Equation (7) is the same as in Son et al. [38] (α 1 = 28). The correspondence between the Monin-Obukhov length and the stability classes proposed by Pasquill [39] allowed us to assume that in our case 1 L = −0.036. From Equation (5) where γ. represents the heat flux under unstable conditions [36]; g is the gravity; T, c p and ρ a are the temperature, specific heat and density for air respectively; and Q s is the sensible heat flux. Under conditions not too far from neutrality, the log-linear profile of the wind near the surface is used for the calculation of u * [40]. For strong and unstable conditions, the profile proposed by Lonin and Linero [36] can be used. Both profiles for atmospheric unstable conditions lead to a strong wind velocity gradient near the surface, with a dimished velocity flux due to convective mixed, usualy stronger than mechanical mix [36].

Association between Hydrometeorological Variables and Dust Emission Flux
The data on soil volumetric moisture (SVM) and the zonal and meridional components of wind were obtained from the NARR database and monthly averages were calculated. The spatial resolution of the grid is 0.3 • (~32 km), with 349 × 277 cells for 29 atmospheric pressure levels and a temporal resolution of three hours. This work only considered the first level of this grid (1000 hPa ≈ 100 m). Monthly averages of precipitation (Precp) and air temperature (T) were obtained from the University of Delaware (https://www.esrl.noaa.gov/psd/data/.../data.UDel_AirT_Precip.htm). Likewise, in this database, Precp and T data were extracted from the pressure level closest to the surface for the grid points corresponding to San Felipe and Vizcaino coordinates ( Figure 1). With the grid points selected (one for each study area) monthly series were derived for each hydrometeorological variable for the period from 1979 to 2014. Subsequently, monthly climatologies were calculated from monthly time series, which were displayed graphically to illustrate the intra-annual behavior of each hydrometeorological variable.

Intra-Annual Variability
Intra-annual variability of hydrometeorological variables and dust emission fluxes were determined from an analysis of power spectral density (PSD), using the periodogram method (refer to [40]). In addition, and considering that wind is the primary driver of the dust emission process, the variance ellipses proposed by Emery and Thompson [41] were elaborated to conduct an in-depth analysis of the relevant aspects of their variability.
In order to describe and quantify the associations between dust emission flux and hydrometeorological variables in each study area, an analysis of mathematical association (principal component analysis, PCA, in its numeric representation), was performed following the criteria of Santamaria-Del-Angel et al. [42].

Interannual Variability
An annual series for the period 1979-2014 was obtained based on monthly time series for each hydrometeorological variable. Subsequently, the normality of these annual series was contrasted with the Shapiro and Wilk test [43]. Based on the results obtained, the non-parametric techniques (modified Mann-Kendall and Theil-Sen) used in works, such as Hamed and Rao [44], Dahmen and Hall [45] and Sen [46], were applied to estimate the trend and slope of each time series. Then, to explore the influence of hydrometeorological variables on MD, we analyzed the associations between wind speed, temperature, soil volumetric moisture, precipitation, and dust emission flux obtained in each of the study areas for the period 1979-2014.

Effects of ENSO on the Frequency and Intensity of Dust Events
To identify the effects of ENSO events on dust production, we selected the monthly time series for dust emission flux that coincided with strong ENSO events. Subsequently, standardized anomalies of dust emission flux (Z i , [47]) were calculated with Equation (9) and then depicted graphically on the Southern Oscillation Index (SOI) (https://www.ncdc.noaa.gov/teleconnections/enso/indicators/soi/).
where X i represents individual data for each variable, X is the average, and SD X is the standard deviation.
In order to calculate the inter-annual variability in the time series of dust emission flux, a low-pass Lanczos filter with a 12-month window was applied. Finally, the standardized anomalies of the filtered series were calculated and, from these, cross-correlations were performed [48], with the SOI for each of the study regions.

Intra-Annual Variability
The monthly climatologies of dust emission ( Figure 2a Figure 2c) display a unimodal behavior in both regions, with a peak in August (28.6 • C in San Felipe and 26.50 • C in Vizcaíno). In addition, higher temperatures were recorded in San Felipe relative to Vizcaíno between April and October.
Atmosphere 2018, 9, x FOR PEER REVIEW 6 of 17 filtered series were calculated and, from these, cross-correlations were performed [48], with the SOI for each of the study regions.

Intra-Annual Variability
The The climatologies of the annual cycle of soil moisture (Figure 2d) reveal the same variability in the two regions, with a minimum in June and a peak in January, the highest values predominated in Vizcaíno. Finally, precipitation ( Figure 2e) showed a bimodal behavior in both regions, with peaks in August (1.75 mm month −1 ) and February (1.66 mm month −1 ) in San Felipe, and in December (2.11 mm month −1 ) and September (1.64 mm month −1 ) in Vizcaíno.
It is evident (Figure 2) that the peaks of dust emission flux are consistent with the peaks of wind speed and the minima of soil moisture and precipitation. The mean values for the variables in different seasons (Table 1) show peak dust emission flux in summer (17.41 mg m −2 d −1 ) and spring (281.11 mg m −2 d −1 ) for San Felipe and Vizcaino, respectively. Table 1. Seasonal climatologies of dust emission flux (E), soil volumetric moisture (SVM), temperature (T), wind speed (v) and precipitation (Precp) for winter, spring, summer, and autumn for the period 1979-2014. N = 36.

San Felipe
Vizcaino With regard to the meteorological variables (Table 1), optimal conditions for dust emission flux in San Felipe occurred during the summer, when low moisture (7.18%) and precipitation (0.7 mm month −1 ) coincided with peak values of wind speed (2.62 m s −1 ) and temperature (26.9 • C). In Vizcaino, minimum moisture (15.59%) coincided with maximum temperature (24.7 • C) in the summer, although the peak of wind speed (4.42 m s −1 ) occurred in spring.
PSD showed two dominant peaks of energy at 12 and 6 months in dust emission flux and the hydrometeorological variables in the Vizcaino region (Figure 3b,d,f,h,j). San Felipe showed a single peak of energy at 12 months in wind speed, soil moisture, and precipitation (Figure 3c,e,g,i), and statistically significant peaks at six months in soil moisture and precipitation (Figure 3e,g). Dust emission flux showed a low annual signal below the established statistical significance.
A comparison between variance ellipses for wind in the different seasons of the year (Figure 4) shows the highest variability in winter in the two study regions (Figure 4). The standard deviation in the semi-axis of highest variability was 0.61 m s −1 for San Felipe and 0.64 m s −1 for Vizcaino (Table 2).
Average vectors for San Felipe showed a prevailing southeast direction pointing toward the GC, although the eastward direction dominated in summer. In Vizcaino, most vectors also pointed to the southeast, heading to the semi-axis of highest variability (Figure 4). Additionally, the major axis of the variability ellipses is perpendicular and parallel to the BCP in the San Felipe and Vizcaíno regions.
The lowest wind variability occurred during autumn and spring, with standard deviations on the semi-axis of highest variability of 0.38 m s −1 and 0.45 m s −1 for San Felipe and Vizcaino, respectively ( Table 2).
During spring-summer, the axis of highest variability of wind speed showed similar behavior in Vizcaíno and San Felipe (~0.45 m s −1 ). Mean wind speed was~1.49 m s −1 for Vizcaino (Table 2). Wind maps for the BCP region ( Figure 4) show that the dominant wind direction is from the BCP to the GC, except during winter. The highest intensities predominated during spring and summer. Atmosphere 2018, 9, x FOR PEER REVIEW 8 of 17  although the eastward direction dominated in summer. In Vizcaino, most vectors also pointed to the southeast, heading to the semi-axis of highest variability (Figure 4). Additionally, the major axis of the variability ellipses is perpendicular and parallel to the BCP in the San Felipe and Vizcaíno regions. The lowest wind variability occurred during autumn and spring, with standard deviations on the semi-axis of highest variability of 0.38 m s −1 and 0.45 m s −1 for San Felipe and Vizcaino, respectively ( Table 2). During spring-summer, the axis of highest variability of wind speed showed similar behavior in Vizcaíno and San Felipe (~0.45 m s −1 ). Mean wind speed was ~1.49 m s −1 for Vizcaino (Table 2). Wind maps for the BCP region ( Figure 4) show that the dominant wind direction is from the BCP to the GC, except during winter. The highest intensities predominated during spring and summer.  Table 2. Parameters of variance ellipses and average wind speed for the different seasons of the year during the period 1979-2014. Mag and Dir refer to average wind speed and direction, respectively, S u is the axis of highest variability, S v is the axis of lowest variability, and θ is the inclination angle.

Seasons
Zone The association between Spearman's correlation coefficients of Z-transformed variables vs. the significant components calculated in the PCA (Table 3) revealed, through component 1 (which accounts for 37.15% and 46.25% of variance in San Felipe and Vizcaino, respectively), a direct association between dust emission flux and wind, except for the fact that dust emission flux and wind speed are inversely related to precipitation in Vizcaíno. In San Felipe component 2 (that explains 35.31% of variance), soil volumetric moisture and precipitation are directly related, with variables being controlled by temperature. Finally, in Vizcaino component 2 (31.55% of variance), moisture and temperature are inversely related (Table 3).

Interannual Variability, Trends in Dust Emissions Flux and Hydrometeorological Variables
Mean annual values of dust emission flux and hydrometeorological variables during the period 1979-2014 in San Felipe and Vizcaino deserts are shown in Table 4. With regard to dust emission flux, values for Vizcaíno and San Felipe are 170.77 and 24.05 mg m −2 d −1 , respectively (Table 4). The hydrometeorological variables reported higher values for wind speed (3.66 m s −1 ), soil moisture (18.79%) and precipitation (1.01 mm month −1 ) in Vizcaíno, while temperature (20.77 • C) was higher in San Felipe.
When the normality of the dataset is contrasted with the Shapiro and Wilk test [43], the null hypothesis is rejected for N = 36. Hence, the set of hydrometeorological data and dust emission flux is not normally distributed at a significance level of P < 0.05, as the statistics obtained are lower than the critical value provided by the tables (<0.935).

Influence of El Niño on Dust Emission Flux
The cross-correlation between the time series of El Niño and standardized anomalies of dust emission flux showed a lag of 19 months in San Felipe (ρ = 0.4, P < 0.05) and 21 months in Vizcaino (ρ = 0.33, P < 0.05), i.e., there are positive anomalies in dust emission flux after the onset of El Niño indicating a significant increase in dust production in both study areas.
The results obtained in the cross-correlation are shown in Figure 5. These show that when the climate event starts losing intensity, dust emission flux shows positive anomalies ranging between one and three standard deviations in both regions for each event studied. Atmosphere 2018, 9, x FOR PEER REVIEW 12 of 17

Discussion
The differences between the location and geographical features of the study areas suggest, through the analysis of intra-annual variability, that dust emission flux rates are regulated by the synergy between hydrometeorological variables. Dust emissions flux are relevant when low soil moisture, high air temperature, and high wind speeds prevail.
Precipitation levels in the two desert areas during winter, spring, autumn, and summer, as calculated in the present study, are similar to those reported by Hastings and Turner [49].
Because of the extensive plains (elevation <200 m a.s.l.) between 26-28° N, Vizcaíno is affected by winds from the northwest (values in the Pacific off the BCP as reported by Castro and Martinez [50] and Morales-Acuña [26] are above 5 m s in spring). In contrast, the natural barrier (elevation > 1200 m a.s.l.) located in the west margin of San Felipe (Figure 1) minimizes this effect.
The main feature found to affect the annual variability of the wind field in the region is the displacement during winter and summer (frequency ~6 months) of the North Pacific High-Pressure Center (NPHPC), located approximately at 32° N-140° W [51,52]. The findings using spectral analysis suggest that the influence of the NPHPC on dust emission flux and the hydrometeorological variables is intense not only in the Pacific Ocean decreasing toward the coast, as reported by Castro and Martinez [50], but also in Vizcaíno, being weaker in San Felipe, as reported by Salinas-Zavala et al. [24]. The above suggests that the climate conditions in Vizcaino are better suited for promoting dust events relative to San Felipe. Particularly, the seasonal variability of precipitation is attributed mainly to processes at local (upwelling) and global (El Niño) levels in addition to tropical cyclones and the passage of cold fronts [52,53].
The direction of average vectors of variance ellipses in the period of highest variability (winter) (Figure 4) indicates that the variability in wind speed partly responds to the winter pattern and the NPHPC, which in this season moves toward the equator (~25° N-125° W) producing strong pressure gradients and a component of wind blowing in the same direction of displacement (~40° N) [51].

Discussion
The differences between the location and geographical features of the study areas suggest, through the analysis of intra-annual variability, that dust emission flux rates are regulated by the synergy between hydrometeorological variables. Dust emissions flux are relevant when low soil moisture, high air temperature, and high wind speeds prevail.
Precipitation levels in the two desert areas during winter, spring, autumn, and summer, as calculated in the present study, are similar to those reported by Hastings and Turner [49].
Because of the extensive plains (elevation <200 m a.s.l.) between 26-28 • N, Vizcaíno is affected by winds from the northwest (values in the Pacific off the BCP as reported by Castro and Martinez [50] and Morales-Acuña [26] are above 5 m s −1 in spring). In contrast, the natural barrier (elevation > 1200 m a.s.l.) located in the west margin of San Felipe (Figure 1) minimizes this effect.
The main feature found to affect the annual variability of the wind field in the region is the displacement during winter and summer (frequency~6 months) of the North Pacific High-Pressure Center (NPHPC), located approximately at 32 • N-140 • W [51,52]. The findings using spectral analysis suggest that the influence of the NPHPC on dust emission flux and the hydrometeorological variables is intense not only in the Pacific Ocean decreasing toward the coast, as reported by Castro and Martinez [50], but also in Vizcaíno, being weaker in San Felipe, as reported by Salinas-Zavala et al. [24]. The above suggests that the climate conditions in Vizcaino are better suited for promoting dust events relative to San Felipe. Particularly, the seasonal variability of precipitation is attributed mainly to processes at local (upwelling) and global (El Niño) levels in addition to tropical cyclones and the passage of cold fronts [52,53].
The direction of average vectors of variance ellipses in the period of highest variability (winter) (Figure 4) indicates that the variability in wind speed partly responds to the winter pattern and the NPHPC, which in this season moves toward the equator (~25 • N-125 • W) producing strong pressure gradients and a component of wind blowing in the same direction of displacement (~40 • N) [51].
The associations between dust emission flux and the hydrometeorological variables in each study area, using PCA analysis show that dust emission flux is driven mainly by fluctuations in wind speed, as reported by Pu and Ginoux [19] for deserts adjacent to our study area. These associations also highlight the direct relationship between soil moisture and precipitation, both controlled by temperature.
Annual averages of the hydrometeorological variables (Table 4) revealed that Vizcaíno showed higher values than San Felipe, due to its exposure to the northwest winds, the latitudinal and anomalous distribution of precipitation [53], and its proximity to the California Current [52]. Separately, dust emission flux rates are higher in Vizcaíno than in San Felipe. Annual precipitation values are consistent with those reported by Arriaga-Ramírez and Cavazos [53] for the middle of the Baja California Peninsula, north of Sonora, and southwest Arizona.
A comparison between annual emissions in San Felipe and MD fluxes reported by Muñoz-Barbosa et al. [20] and Félix-Bermúdez et al. [54] for the regions of San Felipe and Ensenada Bay (a region adjacent to the San Felipe desert), respectively (Table 7) identified a rising trend, thereby reflecting increases in dust emission flux that is consistent with those obtained by Pu and Ginoux [19] for the BCP. The collate of MD emissions in the most relevant desert areas worldwide (Table 7) with atmospheric inputs of MD from the Vizcaíno and San Felipe deserts evidence that Vizcaíno is an important source of dust emission flux, since its mean annual emission (170.77 mg m −2 d −1 ) is of the same order of magnitude as emissions from the Sahel and Sahara deserts [55,56]. The comparison of seasonal dust emission flux recorded in Vizcaíno and San Felipe (Table 7) with the standardized emissions per unit area reported by Xuan et al. [58] for northern China (Gobi desert) show that peak dust emission flux occurs in spring, as observed in Vizcaíno. Besides, MD emissions per unit area in Vizcaíno and San Felipe are higher than those in the Gobi desert ( Table 7).
The trend analysis of the hydrometeorological variables (Table 5) identified increasing trends in wind speed and temperature. The rising trend in temperature observed in this study is consistent with the one reported by Blunden et al. [59], who reported positive anomalies of air temperature ranging between 1 • C and 5 • C in Vizcaino and San Felipe for the period 1981-2010. In addition, our positive trends confirm the reports of Englehart and Douglas [60] and Pavia et al. [61] regarding a shift to positive trends in the period between 1970 and 2004 in Mexico. Our results are also consistent with those reported by Prein et al. [62] for the southwest United States. The decrease in precipitation observed in this work is consistent with the one reported by Blunden et al. [59] for northern Baja California and differs from the results obtained by Arriaga-Ramírez and Cavazos [53] for the central part of the BCP, northern Sonora, and southwest Arizona, where annual series show no statistically significant trend. The negative trends for soil moisture coincide with a lower amount of humidity caused by the decrease in precipitation and the rising of air temperature.
Regarding trends in dust emission flux, the increases found for the past four decades are consistent with the rise in dust emission flux from spring to autumn reported by Pu and Ginoux [19]. The association between dust emission flux and hydrometeorological variables ( Table 6) is consistent with studies conducted in southeastern United States [19], China, and Mongolia [18], suggesting that increases in dust emission flux in recent decades are regulated by variations in wind speed, temperature, precipitation, and the extension of bare soil. The lag in the occurrence of relevant dust emission flux due to El Niño events (strong and very strong) observed in this work is consistent with Okin and Reheis [16] in the southwestern United States. In our study area, such effects can be explained by the impact of this climate phenomenon on precipitation in the BCP [52,63,64].
Regarding the potential implications of the dust deposition in surface waters of adjacent seas, if the particles contains Fe, Cu, or Cd, these metals may affect primary production [65][66][67], as they play a key role in biogeochemical processes such as photosynthesis [68]. Studies on dissolved Fe suggest that it may be a limiting element for primary production [20]. Félix-Bermúdez et al. [54] suggest that the impact of dissolved Fe on marine ecosystems probably depends on the composition of water in terms of the N:Fe ratio. Accordingly, when water is already enriched with Fe, the atmospheric iron deposited will be an additional source of a non-limiting nutrient and N will probably act as the limiting element for primary production. Otherwise, dissolved Fe will be the limiting nutrient with a key role as a driver of primary production.

Conclusions
This study showed that seasonal and inter-annual variability in dust emission flux in the Vizcaíno and San Felipe desert areas is controlled by fluctuations in the hydrometeorological variables and that these, in turn, may respond to the influence of the NPHPC. Consequently, under optimal hydrometeorological conditions, dust emission flux peaks in summer (San Felipe) and spring (Vizcaíno).
The increases in dust emission flux during the study period  are associated with the rising trends in wind speed and air temperature, and with the decreases in moisture and precipitation that varied in both desert regions.
In addition, the frequency and intensity of dust events in San Felipe and Vizcaino appear to be related to strong and very strong ENSO events, with lags of 19 and 21 months in San Felipe and Vizcaino, respectively.
Finally, given the amount of dust emitted from each desert areas evaluated and the rising trend in these emissions, the BCP should be considered as a significant source of dust at local and regional levels.