1. Introduction
Atmospheric deposition (AD) of nutrients to water bodies often includes quantifying the amounts of nitrogen and phosphorous deposited directly from the atmosphere to the water surface [
1,
2,
3,
4,
5,
6,
7]. AD of nitrogen can occur through dry (includes gaseous), dust, and bulk (precipitation) deposition. Unlike nitrogen, phosphorus does not have a stable gaseous phase, but is often associated with small particles < 2.5 µm, and larger dust particles that behave as gases with limited settlement from gravity, while phosphorous does not have a gas phase, dry, dust, and bulk (precipitation); therefore, AD deposition associated with particulates occurs [
8]. Many human activities, including combustion, mining, agriculture, and transportation, increase particulate dust generation, and therefore atmospheric deposition [
8]. In recent history, dust generation has increased due to anthropogenic activities and climate change [
9,
10], with one study estimating that over the last 100 years, 40% more dust has been put into the atmosphere due to anthropogenic activity [
11]. Some anthropogenic activities, like agriculture, not only increase dust production [
12], but can increase the nutrient content of the dust, especially during planting and fertilization [
13]. Human diversions of water have combined with climate change to cause lakes around the world to dry up, exposing lakebeds that further increase mineral dust generation [
14,
15,
16]. In addition to nutrients, atmospheric dust can also contribute other materials such as heavy metals, organic pollutants, and plastics [
17,
18,
19,
20,
21], and effect lake geochemistry [
22].
The way particulates settle from the atmosphere depends on their size and weight. Larger particulates settle faster than smaller particulates [
23]. Hinds and Zhu [
23] state that particle size is the most important parameter for characterizing aerosol behavior and state that particles with equivalent diameters of 0.1, 1.0, 10, and 100 µm settle in perfectly calm air at 8.8 × 10
−7, 3.7 × 10
−5, and 3.1 × 10
−3, and 2.5 × 10
−1 m per second (m/s), respectively. This means that for particles smaller than about 10 µm, a light breeze can keep the particle from settling. Particulates, such as photochemical smog (mostly nitrogen particles), smokes, and fine dust less than 2.5 µm, do not settle from the atmosphere, but are essentially kept aloft by Brownian motion and wind currents [
23]. Particles in the poor-settling 0.2–10 µm range are often enriched in phosphorus, leading to higher nutrient deposition when that fraction is involved [
24].
For discussion and analysis purposes we break AD in this paper into three categories: dust, dry, and precipitation. The sum of these three processes is bulk or total nutrient AD. We define dust AD as larger particles that settle, dry AD as smaller particles that are less apt to settle, such as gases, fumes, and smokes that are deposited by contact with the wet lake surface, and precipitation AD as nutrients deposited by a precipitation event, which includes any dry, dust, gases or smaller particles in the atmosphere at the time of the precipitation event. Nutrients related to dry AD and the smaller particles associated with precipitation AD can travel long distances, are well mixed, and can be difficult to attribute to a specific source. Dust particles (>2.5 µm) settle out of the atmosphere, and, subsequently, often do not travel as far, making source attribution possible. Dust AD is often associated with high-wind events or mechanical activities that mobilize and transport dust from the surface.
Utah Lake is a large but shallow eutrophic freshwater lake located in the center of Utah county. The lake has large population centers on the north and east sides that potentially increase nutrient pollution loadings to the lake. While the south and west sides of the lake are not yet densely populated, they are being developed and will likely continue to grow in the future. The lake is classified as hyper-eutrophic by the Larsen–Mercier Trophic State model, but is only moderately eutrophic by the Carlson Trophic State Index [
25]. The lake has a surface area of roughly 375 km
2 (~92,665 acres), but an average depth of less than 3 m (~10 ft) during non-drought years, and even shallower during drought years [
26]. The shallow depths create a surface-area-to-volume ratio that makes the lake more susceptible to nutrient deposition from the atmosphere than deeper lakes [
26,
27,
28]. While nutrients like nitrogen and phosphorus are often the limiting factor in aquatic plant growth [
29], Utah Lake plant growth may be limited by high turbidity, causing reduced light penetration [
25]. This is supported by recent studies showing algal blooms, as measured by remote sensing methods, showing no trends over the last 40 years while both population and anthropogenic nutrient loads have increased [
30,
31,
32].
Recently, interest in the nutrient loadings to Utah Lake has grown significantly as the Utah Division of Water Quality (DWQ) looks to reduce wastewater treatment plant phosphorus discharge limits to 0.1 mg/L [
26]. It is estimated that, for Utah Lake, these reductions could cost up to 1 billion dollars to meet the limits; however, it is not certain that the algal blooms are directly related to nutrient loadings from wastewater treatment plants because of other nutrient sources for the lake and other limiting factors to algal growth, so these limits may not significantly affect algal blooms [
26]. Using satellite data, researchers showed that the algal blooms have slightly decreased over the last 40 years despite the likely increase in nutrient loadings during that time [
30,
33,
34]. Abu-Hmeidan et al. [
35] have shown that phosphorus levels in ancient Utah Lake sediments are not statistically different from current sediments, suggesting that the historically recent wastewater plants are not causing nutrient increases in the lakebed [
35], and Randall et al. [
26] presented evidence that water column phosphorous concentrations are mediated by sediment recycling rather than external inputs. The lake has been shown to retain 90 percent of influent phosphorus loadings, likely stored in the sediments, further supporting a system controlled by sediments [
25].
Several past studies, including research performed by DWQ, have looked at nutrient loadings to the lake through natural water inflows and wastewater treatment effluents and have quantified nutrient loads [
25]. However, the studies failed to consider atmospheric deposition or geologic conditions, such as sediments and soils, as a significant source of nutrients [
36]. More recent studies have shown that AD is a major contributor of both phosphorus and nitrogen to waterbodies, and even suggest that atmospheric deposition alone could keep Utah Lake eutrophic [
11,
27,
36,
37,
38,
39]. In Utah Lake, via research evaluation, nutrients from sediments were found to be a significant source, and sediment sources alone, without external loadings, provide as much as 19 mg/L of phosphorus to laboratory water column samples [
26,
35]. The Utah Lake watershed contains major phosphorus deposits in several geologic formations, which result in high-phosphorous concentrations in sediments and soils and contributes to the large potential load from this source [
26,
35,
40,
41].
There has been published research focused on characterizing different nutrient AD processes to Utah Lake. Brown et al. [
42] measured precipitation AD and estimated that approximately 100 tons/year of total phosphorus falls on Utah Lake. Unpublished works by Brahney [
8] and Carling for the Utah Lake Science Panel have focused on dust AD only and estimate between 4 and 13 tons/year of total phosphorus. Barrus et al. [
36] used lakeside measurements to estimate total AD rates and estimated between 130 and 260 tons/year of total phosphorus through deposition [
36]. These studies all use different measurement methods that preferentially measure AD from different processes.
Research by Carling et al. [
43] looked at sources and composition of atmospheric deposition of dust on the Wasatch Front, the region bordering Utah Lake, finding distant desert playas in southwest Utah to be the main source of deposition [
43,
44], with dust being transported by southern winds funneled by area geography. However, these studies focused on larger wind events that average only 4.3–4.7 events per year, used deposition sampler designs that discriminate against local sources, and never tested the samples against local areas. While Carling et al. [
43] found distant desert playas to be the main source of dust on the Wasatch Front, other works have shown proximal sources to be an underrepresented source in AD studies [
45]. We extend this work by testing both local and distant sources, and, specifically, evaluating dust deposited near the shore of Utah Lake.
This paper aims to provide further understanding of the sources and composition of atmospheric deposition of dust or solids to Utah Lake. We hypothesize that smaller, local storms transport local soil near the lake in smaller quantities, but much more frequently. With elemental fingerprinting of area soils, we show that local sources provide most of the dust AD to Utah Lake, even though dust is a minor portion of the total nutrient load to Utah lake from AD. Through the sampling program that generated the material, our study was not designed to quantify dust deposition to the lake; the resulting data support the conclusion that dust is not a major source of AD. This is based on the large number of samples required to obtain enough dust for measurements, which indicates that the mass of dust deposited to Utah Lake is small. This manuscript is based on the Master’s thesis by Telfer [
46].
2. Materials and Methods
2.1. Utah Lake Deposition Samples
The geologic materials in the deposition samples were obtained from the sample network and design reported by Barrus et al. [
36]. This earlier work, and the samples we analyzed here, were designed to quantify total AD to the lake, not to collect dust. For this dust attribution study, we used samples taken at five sites surrounding the lake. The sampler design, seen in
Figure 1, consisted of a dry deposition bucket with 3–4 L of water and an empty wet deposition bucket for precipitation with a lid moving to cover one bucket at a time. The dry bucket collects any particulates that contact the water surface, and the wet bucket collects precipitation. Details on sampler placement, design, and operation are provided in [
36]. The operation, maintenance, and sampling procedures were performed in accordance with the design of previous studies [
27,
36].
We collected samples weekly, when possible, from May to October of 2022. Once in the laboratory, we decanted 300–400 mL of liquid collected in the field for nutrient testing, with the remainder stored and later used for this study. After decanting the remaining collected solids in the remaining sample, some solids were likely lost. The wet deposition samples (precipitation events) rarely contained more than 300–400 mL of total liquid. Because this left over sample liquid was not stored, we did not use wet deposition samples for this study, only solids collected in the dry sample buckets. At times, the samplers did not open for precipitation events, and thus, the dry buckets did contain some deposition from precipitation events.
The individual samples did not contain much dust, so we combined individual samples to obtain enough dust to analyze based on season. We combined the weekly samples from May through July into a spring sample and the weekly samples from August through October into a fall sample for each sample location.
Once the samples were combined, we filtered the water using a 0.45 µm cellulose acetate vacuum filter to separate particulates from water to obtain the dust (i.e., solids) for analysis. Traditionally, 0.45 µm is considered the boundary between dissolved and non-dissolved solids in water, so this filtration step separated all suspended solids including colloids [
36].
2.2. Prevailing Winds
While smaller dust particles can be transported significant distances at relatively low wind speeds, larger particles require higher wind speeds. In addition, higher wind speeds are required to suspend dust from potential source areas. We expect a significant portion of the dust may be suspended because of mechanical disturbance such as vehicle travel, though we have no data to support this statement, so we discuss only the impacts of wind speed.
We obtained wind data from the Utah Climate data center for 6 different weather stations around Utah Lake [
47]. The data consists of 10 min average data for wind speed and direction from 1 January 2017 to 7 July 2022. The data set contains 1.6 million total records.
Figure 2 shows five different wind roses from three weather stations located around Utah Lake, selected to demonstrate different wind regimes. The color scale for all the plots is the same. The length of each segment on the roses indicates how frequently the wind blows in that direction (expressed as a percentage of time over the study period), while the color represents the wind speed. The panels, from left to right, show wind roses based on all winds, and winds over 15, 20, 25, and 30 km/h. For panels 2–5, the segment length is scaled to represent percentages only for winds in that speed category rather than the entire study duration. Winds at the Tintic weather station (top row) show that the prevailing wind direction, considering all winds, is from the northeast (left panel); however, for winds over 15 km/h (other panels), it is from the south/southwest. Close examination of the left panel shows these higher winds as small segments, but difficult to see. The middle row, Lincon Point weather station, shows low-speed prevailing winds almost exclusively from the southeast, while higher speed winds are almost exclusively from the southwest. The bottom row, Lindon weather station, shows prevailing winds without a single dominant direction, though nearly all winds over 15 km/h are from the northwest.
These plots show that wind directions for winds over 15, 20, 25, and 30 km/h are similar.
Figure 3 shows a map with wind roses for each of the six weather stations. The stations are marked with a small green diamond, with the wind rose for all winds and for winds over 15 km/h located above and below for each diamond, respectively. The color scale is the same as that for
Figure 2.
Figure 3 shows that, south of the lake, strong winds are predominantly from the south/southwest, while on the northern portion of the lake, strong winds are predominantly from the northwest. These wind data will help explain the attribution results later in the analysis, and were used to help select potential source sample locations.
2.3. Source and Deposition Samples
Figure 4 shows the location of the five deposition sample locations in light purple, and the 17 potential source sample locations in light red. The deposition sample locations are: Pump Station (northwest corner of lake), Orem (northeast region), Bird Islan (in the lake), Lake Shore (southcentral), and Mosida (southwest). The Lake Shore location is isolated from the Mosida by West Mountain. The 17 potential source locations are, with name and number: 1 = Lake Shore, 2 = Eagle Mountain, 3 = 5-Mile Pass, 4 = Chimney Rock Pass, 5 = Goshen WMA, 6 = Elberta, 7 = Mouth of Spanish Fork Canyon, 8 = Cherry Creek, 9 = Fumarole Butte, 10 = Sunstone Knoll, 11 = Sevier Dry Lake, 12 = Cricket Mountain, 13 = White Hills near Sevier Lake, 14 = Mid-Sevier Lake near road, 15 = Highway 6 South of Delta, 16 = Burraston Ponds, and 17 = Miners Canyon. There is a deposition sample location and a potential source location both named “Lake Shore”, but they are in different locations.
Previous work by Goodman et al. [
44] and Carling et al. [
43] sampled surface soils from southwestern playa sources near Sevier Dry Lake. We collected our samples representing potential dust source areas (
Figure 4) in a larger area. We collected three samples as close as possible to the locations reported in Goodman et al. [
44], which are called Sevier Dry Lake, Sunstone Knoll, and Fumarole Butte (
Figure 4). We collected additional surface soil samples of potential dust source areas in the open areas between Sevier Dry Lake and Utah Lake (
Figure 4), from potential sources areas surrounding the lake, and to the west of the Lake (
Figure 4).
We collected soils samples following the EPA AP-42 method, which was used by previous researchers studying surface soils as dust sources [
48,
49]. We scraped the surface soils with a shovel in an upward direction to a depth of ~2.5 cm (1 inch) or the diameter of the largest particle, whichever was less. As recommended by the EPA method, we collected all the particles in the area and did not avoid large particles. The EPA recommends collecting 10 pounds per site, but we only collected 1–2 pounds to allow for easier storage and transport.
2.4. Sample Preparation and Analysis
We first sieved the potential source soil samples with a 250-µm sieve for 5 min to remove the large particles unlikely to be transported by wind. This sieve size was chosen based on previous particle size analysis performed by Dansie et al. [
50], Chepil [
51], and Péwé [
52], who recommended this size to eliminate particles too large to travel long distances. We did not treat the soil samples to remove organics if the organic material could pass the sieve, as this indicates it was small enough to be transported with the dust. When we collected the Sevier Dry Lake sample it was wet and was allowed to dry over time; this caused the sample to dry into hard clumps that were not easily broken. We broke these clumps apart using a soil flail grinder. Six of the source dust samples were collected when damp and dried in a drying chamber at 45 degrees Celsius overnight. This drying procedure caused some small clumping of the samples; we broke these clumps gently using a clean spatula before sieving.
The previous study filled buckets with 3–4 L of water and collected the water from the buckets weekly to collect dry atmospheric deposition nutrients (and solids) [
27,
36]. They transported each weekly sample to the laboratory and decanted 300–400 mL of liquid, ~7.5–15% of the sample from the top, and used this liquid to measure nutrient concentrations. Decanting the small sample for nutrient analysis left most of the collected solids in the remaining sample, though some solids were likely lost in the decanting. After decanting, the samples were stored for later analysis. We used these samples to represent deposition to Utah Lake. We only used samples collected from the dry deposition buckets, shown in
Figure 1, that had been stored after analysis by the previous study.
We had to discard a few samples that had algae growing in the archived samples, as these would have clogged the filter. We used the remaining archived samples, i.e., the liquid remaining after decanting 300–400 mL for nutrient analysis, to perform source attribution analysis and to quantify mass deposition. For source attribution, this mass loss from either decanting or the few excluded samples is not an issue. Estimates of the mass deposition rates are only approximate as some samples and solids were excluded.
To separate the solids deposited as dust from the water samples, we filtered the sample liquid through 0.45-µm cellulose acetate vacuum filters. We selected this filter type to separate any particles too large to be dissolved and so that the filter material would not affect the nitric acid digestion and ICP-OES analysis. We explored other separation techniques and filter types with little success. The last dust sample collected on 10 October at the Pump Station and had heavy contamination from the seal material on our sampler’s lid. We processed this sample, labeled “Pump Station Contamination” or PS contamination, separately from the other Pump Station samples, rather than including it in the fall sample. However, after analysis, the samples results were not significantly different.
We selected the cellulose acetate filters because they have extremely low-ash content. After filtering the samples to separate the particulate matter, we placed the filters and accompanying solids in a muffle furnace at 500 degrees Celsius for 2.5 h to remove the filter from the samples before digestion [
53].
Once the samples were either filtered and reduced to ash (dust samples) or sieved (soil source samples), we digested the resulting samples using a nitric acid microwave digestion method. This digestion method uses 69.9% nitric acid at 170 degrees Celsius for 20 min with 50 min of warmup and cooldown time. As the temperature is above boiling, the digesting is performed in a pressure vessel.
We analyzed the digested samples for element concentrations using a Thermo-Scientific iCAP7400 Radial ICP-OES for 25 elements Al, As, B, Ba, Ca, Cd, Co, Cr, Cu, Fe, K, Mg, Mn, Mo, Na, Ni, P, Pb, S, Se, Si, Sr, Ti, V, and Zn. The results showed that cadmium and selenium were below detection limits for all samples, except for one sample taken at Goshen WMA that had 2 mg/kg selenium. Because of these low concentrations, we removed these elements from the analysis. Thus, we measured 23 elements for each sample. We treated these data as vectors in a 23-dimension sample space. Treating these data as vectors allowed us to compare similarity between samples and potential sources. We computed detection limits according to the EPA method detection limit procedure using the BYU Environmental Analytic Lab soil standard.
2.5. Data Analysis
We used four different methods, two statistical and two visual, to compare the elemental spectral signatures of the dust samples collected near the shore of Utah Lake with soil samples collected in suspected source areas. We used these metrics to determine which source samples were most similar to the deposition samples. The analysis methods were spectral angle (SA) [
54], coefficient of determination (R
2), and Multidimensional Scaling (MDS) to visualize these two distance measurements [
55]. We used multidimensional radar plots to visualize similarities and differences between the samples and potential sources. We computed SA, R
2, and the radar plots using python. We used SAS JMP Pro version 16 to compute and visualize MDS plots results. The R
2 value is similar to a Euclidean distance as it is the squared pair-wise difference between two samples normalized by the difference of the sample with the average. While the Euclidean distance is the square route of the squared distances, R
2 does not apply the square root, but can be thought of as a distance measure in a higher dimensional space.
Figure 5 presents the difference between the two distance measures we used, showing one sample in two-dimensional space (orange) with two potential sources shown in grey and blue for Source A and Source B, respectively. This is analogous to using only 2 of the 23 metals for attribution. The left pane shows that the distance between the two source measurements, and the sample in Euclidian measured using Euclidean distance. In this panel, Source A is much closer to the sample than Source B. In the right panel, we show the SA distance between the sources and samples. In this plot, spectral angle B is much closer to the sample that spectral angle A, indicating that the sample is most likely from Source B. In these plots, if you dilute the sample values with a 1:1 dilution, you obtain Source B, which is closer in chemical composition to the sample than Source A.
SA, rather than measuring the Euclidean distance between two measurement vectors, measures the angle between two measurement vectors in multi-dimensional space [
56]:
where
and
are the source and measurement data vectors, each with
a of length
c, which for our data is a value of 23. SA indicates how well the shape of the two series match with an angle of 0, indicating vectors pointing in identical directions. SA is not influenced by differences in magnitude. Spectral angles are the angles created between two vectors in multidimensional space. We created vectors for our samples by assigning each analyzed value a point in an n-dimensional space, where n is the number of values, or 23 for our data. We assume each metal concentration is a separate dimension. The angle between these vectors is then measured, small spectral angles indicate similar samples. Using spectral angles allows us to compare datasets that have similar ratios between the different elements so that concentration or magnitude differences do not interfere. This technique was used to find similarities between geologic samples in previous studies [
57]. We used a spectral angle script found in the Hydroerr Python package [
58].
MDS is a parameter reduction technique that attempts to preserve pairwise distances from higher dimensions in a lower dimension. We compute a complete pair-wise distance matrix with values for the distance between every set of points. The MDS plot tries to distribute samples such that the Euclidean distance in two-dimensional space between each set of points is proportional to the distance in the higher dimensional space.
MDS allows visualization of similarity between samples, generally using a two-dimensional plot [
54]. An MDS plot is interpreted by comparing the proximity of samples in the reduced MDS space where similar samples, those close in higher dimensional space, are placed close together, and dissimilar samples are spread apart. This technique has been used to interpret similarities in geologic samples [
59,
60]. We evaluated the MDS transform using a Shepard Diagram, which plots the pairwise higher-dimensional distance against the MDS distance for each pair of points. Then, we computed an R
2 value of the correlation of the MDS distance with the actual distance for both metrics.
We computed MDS plots using both a pair-wise distance matrix, computed using Euclidean distance, and a pair-wise distance matrix, computed using the SA metric. We then used the MDS algorithm to plot these points in two-dimensional space while retaining the proportional distance between each pair of points. We computed distance matrices using phyton with the Hydroerr Package for SA and the Skikitlearn Package for R2. We computed the MDS transform values using SAS JMP Pro V16.
Visualizing high-dimensional data is difficult. To show how similar and dissimilar samples are in 23 dimensions, we created radar plots of our data. Radar plots allow for the visualization of the similarities and differences between the samples via the shape of the plot. One feature is that the shapes presented in the plots are not influenced by concentration or magnitude of the data. We used the raw uncorrected ICP data for the radar plots because the dilution-corrected results span a large difference in magnitudes that are difficult to visualize on a single plot without normalizing each measurement. The uncorrected data can be thought of as data normalized for concentrations of each element as we dilute the samples to be within the operating range of the ICP. For these plots, we decreased the uncorrected calcium concentrations by a factor of 12 to fit the radar diagram scale. Soils in this study are calcareous with very high calcium values. With radar plots, magnitude is insignificant; the ratios between elements, and the shape they create, are what are visualized and analyzed. However, if one element is significantly larger than the others, the resulting plot does not show variations in the less concentrated elements well.
4. Discussion
As discussed in
Section 2.1, the samples we used for this analysis were not designed to quantify the amount of solids in atmospheric deposition to Utah Lake, rather they were designed to measure nutrient deposition [
36]. To obtain a rough estimate of the mass deposition rate, we used weekly samples collected from the previous study collected over 6 months from May to October 2022. We discarded some samples (
Section 2.1).
For source attribution we needed a large enough mass of the deposited solids for ICP analysis, about 0.03–0.04 g (i.e., 30–40 mg). Solids deposition rates at Utah Lake are so low that we needed to combine the individual weekly samples from any given deposition location to obtain enough dust to analyze. We combined the weekly samples into two seasonal samples, with samples from May through July combined into a spring sample and samples from August through October combined into a fall sample. This was performed for each sample location. The resulting composite samples consisted of the solids from about 12 weekly samples. At the Bird Island deposition location, we combined all the weekly samples into a single composite sample because the total mass was too low for analysis if split. One sample at the Pump Station deposition location was visibly contaminated so we did not combine this sample with the others but analyzed it separately.
As discussed in
Section 2.1, we filtered the samples on 0.45 µm filters, then dried the filters to measure the mass of solids in each sample; 0.45 µm is considered the dividing line between dissolved and suspended solids, so this preparation method captured all the suspended solids. The filters were then reduced to ash and dissolved for analysis.
The resulting data (
Table 1) show that over the 6-month collection period, the total mass collected at any given location ranged from 50 mg to 220 mg at Bird Island and Lake Shore, respectively. The buckets have a surface area of 0.0615 m
2 [
27,
36] and Utah Lake has a surface area of approximately 40,000 ha (95,000 acres), which gives a total deposition mass over the 6-month period of about 325 to 1400 kg [
42]. This is a unit deposition rate of 0.8 g/m
2 to 3.6 g/m
2 over the 6-month sample period for Bird Island and Lake Shore, respectively.
Brahney [
8] reported dust deposition in Provo from 28 to 49 g/m
2/yr, significantly higher than our values. We attribute this to our shorter sampling period and the fact that we excluded any precipitation samples. Brown et al. [
42] found that about 50% of phosphorous deposition is related to precipitation. We did not capture any dust deposited during the winter months and we also excluded some samples.
While the samples we analyzed were not collected to measure the mass of dust falling on the lake, we did collect and weigh all solids greater than 0.45 µm in the archived samples (
Table 1). These results indicate that the mass of solids greater than 0.45 µm deposited on Utah Lake is low, as are the corresponding estimated mass deposition rates.
Our estimates of mass deposition are not truly representative of the amount falling on the lake, but provide good order-of-magnitude data. The few discarded samples and decanted liquid are unlikely to have changed the mass estimates significantly. Members of the Utah Lake Science Panel have argued that the main source of nutrient loading to Utah Lake is dust [
8]. They assumed that the dust had nutrient levels similar to those measured in Utah Lake sediments [
26,
35] and showed that if these concentrations were correct, then it would require unreasonable dust deposition to match the nutrient loadings calculated by Barrus et al. [
36]. Our data show that neither dust is the primary vehicle for atmospheric deposition of phosphorous to Utah Lake, or that the small particles, less than 10–20 µm, can have significantly higher phosphorous concentrations. This is supported by the literature, which shows that phosphorous readily sorbs onto small particles.
This level of mass deposition strongly suggests that dust is not the primary source of nutrient loading to Utah Lake, as estimated phosphorous nutrient loads range from as low as 150 to 200 Mg per year [
27,
36,
42]. These mass dust measurements support the research by Brahney [
8], who estimated about 15 Mg/yr of phosphorus when considering solid particulates alone, whereas the Utah Lake Science Panel estimated 30 Mg/yr [
62]. Brown et al. [
42] measured about 120 Mg/yr from precipitation. Barrus et al. [
36] measured about 230 Mg/yr from all sources. This implies that deposition from dust is about 6.5% to 13% of total phosphorus deposition, precipitation is about 52%, and dry or contact deposition ranges from 35% to 41.5% of the total.
This is likely due to dust deposition consisting of larger settleable particulates compared to bulk deposition, which includes dry or contact processes, containing particulates that do not settle. Bulk deposition is able to wash the PM2.5 fraction during precipitation events, and these small particles contain higher phosphorus levels [
24] than larger particles. The period of this study was during high-dust production for the area; in the months before and after the study, the area is often wet or snowy, preventing local dust generation. While we have no data, we do not expect significant solid mass deposition during the winter and spring months, i.e., November through April.
The dust weights provide insight into spatial dust deposition across the lake. The Orem and Mosida samplers are on opposite sides of the lake and collected nearly identical amounts of dust. The MDS plot in
Figure 8 shows that dust collected at the Orem and Mosida samplers are similar and plot close to each other in MDS space. The similar weight and composition of the two samples suggests that deposition does not change significantly across the lake. The Bird Island sampler collected far less dust. Several Bird Island samples showed zero or near zero dust collected for the week, while a single sample in the last week of collection more than doubled the total dust collected at this site. This might be caused by Bird Island’s proximity to West Mountain. It is possible that the mountain is channeling the winds in a way that misses Bird Island while still crossing from Mosida to Orem. However, as shown by Barrus et al. [
36], nutrient deposition at Bird Island is not statistically different from deposition at the other stations. Our data show that dust deposition is different. We attribute this to the fact that dust deposition is likely less than 15% of total deposition. Also, large particles make up the bulk of the mass for dust deposition, but have a lower associated nutrient concentration [
24]. This supports the conclusion that nutrient deposition does not fall off across the lake, while dust deposition does.
While this study focused on Utah Lake dust attribution, the methods and approach presented here can be applied to many situations involving source attribution for particulate matter such as dust or sediments. Published dust source attribution studies often use methods such as strontium isotopes [
43,
63], atmospheric transport modeling [
16], advanced spectroscopy methods [
64,
65], X-ray diffraction [
65], and X-ray fluorescence [
66]. In their introduction, Zhang et al. [
66] provide a broad overview of methods used for dust attribution that includes additional methods. Our approach, using ICP analysis, generally requires a simpler sample preparation and analysis approach that can be cost effective compared to these more advanced geochemical methods. We showed that this approach is sensitive and clearly differentiates soil particles from different sources. The MDS plots clearly group locations from similar geologic formations together and match our understanding of atmospheric transport. The methods presented in this study can be used for dust, soil, and sediment attribution in an easy, cost-effective manner.
5. Conclusions
Our data show that most of the dust deposition on Utah Lake comes from local fields south and west of the lake. Both MDS distances and Spectral Angle analysis show that the Cherry Creek source location is most similar to all the deposition locations, closely followed by the Chimney Rock Pass source location. The more distant desert playa sources likely still contribute to deposition, but contribute less than local areas. The two Lake Shore samplers, both source and deposition, are less similar to the other local sources and show the strongest similarity to the desert playas. This is likely due to the sampler’s proximity to West Mountain, which blocks local transport from Cherry Creek, but not long-range transport from the desert playas. While these two sample locations are closer to the playa sources than the other deposition locations, they do not plot closely and are most similar to each other.
Dust deposition is likely a minor source of the nutrient loading onto Utah Lake. This was shown by the research of Brahney [
8], who evaluated dust deposition in the literature and found small nutrient loads compared to Brown et al. [
42] and Barrus et al. [
36] who measured deposition using different methods at locations around the lake. We see evidence that dust is a minor contributor to nutrient loads when the dust weights at Bird Island are lower than the shore locations, but nutrient loadings at Bird Island are not statistically different from shore locations [
36]. Large dust particles cannot be the main source of nitrogen and phosphorus when small amounts of dust were collected, but nutrients were still present. We attribute this to much higher nutrient concentrations on small, PM10 and PM2.5 size particles [
24].
Dust deposition sources around Utah Lake appear to be relatively similar with similar deposition rates, as evidenced by the similarity between the Orem and Mosida deposition samples, which are most similar to the Cherry Creek and Chimney Rock Pass source locations. The near identical dust weights and similar source compositions suggest that the dust concentrations are not changing across the lake. The Bird Island sampler is located just south of the path between Orem and Mosida and collected similar dust compositions but far less dust. This may be due to Bird Island’s proximity to West Mountain, which may affect wind patterns and shield the island to some degree.
While this study focuses on source attribution for Utah Lake, we present methods and data that can be used for dust attribution in any setting. ICP data are relatively inexpensive and easy to obtain. They are precise and provide good fingerprints for attribution. The use of MDS and Radar plots simplifies attribution analysis, both visually showing both proximity in high-dimensional space. MDS plots provide analysts with clear metrics for attribution. We demonstrated MDS plots with two distance metrics, Euclidean distance and SA, but MDS plots can be used for any other error metric between two measurements with many parameters, such as R2. Data from other sources can also be used in MDS plots, as long as the data are normalized so that the various parameters have similar magnitudes.
We present this study to both further the understanding of dust sources to Utah Lake and provide a method for soil and dust attribution. We encourage others to use and extend these tools. We showed that this approach is sensitive and clearly differentiates soil particles from different sources. The methods presented in this study can be used for dust, soil, and sediment attribution in an easy, cost-effective manner.