Upland and Riparian Surface Soil Processes in an Urban Creek with Native and Non-Native Vegetation after Fire

Wildfires can pose environmental challenges in urban watersheds by altering the physical and chemical properties of soil. Further, invasive plant species in urban riparian systems may exacerbate changes in geomorphological and soil processes after fires. This research focuses on the 2018 Del Cerro fire, which burned upland and riparian areas surrounding Alvarado Creek, a tributary to the San Diego River in California. The study site has dense and highly flammable non-native vegetation cover (primarily Arundo donax) localized in the stream banks and has primarily native vegetation on the hillslopes. We estimated the post-fire organic matter and particle distributions for six time points during water years 2019 and 2020 at two soil depths, 0–15 cm and 15–30 cm, in upland and riparian areas. We observed some of the largest decreases in organic matter and particle-size distribution after the first post-fire rainfall event and a general return to initial conditions over time. Seasonal soil patterns were related to rainfall and variability in vegetation distribution. The riparian soils had higher variability in organic matter content and particle-size distributions, which was attributed to the presence of Arundo donax. The particle-size distributions were different between upland and riparian soils, where the riparian soils were more poorly graded. Overall, the greatest change occurred in the medium sands, while the fine sands appeared to be impacted the longest, which is a result of decreased vegetation that stabilized the soils. This research provides a better understanding of upland and riparian soil processes in an urban and Mediterranean system that was disturbed by non-native vegetation and fire.


Introduction
Wildfires are increasing in frequency and magnitude under a warming climate, impacting natural resources, infrastructure, and millions of people every year [1]. Simultaneously, human encroachment into fire-prone areas has increased the potential for ignition, as well as risks and damages to human communities, and consequently costs to society [2]. In particular, the number of fires under 5 km 2 in the urban riparian environment has also increased in southern California since 2002 [3]. Landscape conversion due to the expansion of the human population and establishments in southern California (United States) has fragmented chaparral ecosystems and intensified the Wildland-Urban Interface (WUI), leading to the potential for fire ignition [4,5]. Urbanization contributes to larger dry weather baseflow due to impermeable land cover ("urban drool"), which alters native riparian vegetation density and community structure [6]. Higher nutrient loads and more frequent flash floods encourage infestations and the rapid settlement of opportunistic and alien vegetation species [7][8][9][10]. The high density of non-native vegetation biomass and plant infestations can alter riparian hydrological and geomorphological regimes [11][12][13] and fire patterns [7,8,11].
Disturbances such as drought, human ignition sources, vegetation type conversion, and fuel accumulation have transformed many WUI riparian areas into corridors for fire movement [8,14,15], which has environmental and economic implications. However, research is needed to improve the prediction of sediment transport in urban Mediterranean systems, where geomorphic and hydrologic regimes are heavily altered by human activities, which, in return, has consequences for fire regimes and frequency [8,11,13]. The feedback from anthropogenic disturbance to fire regime changes result in altered soil processes and the potential to impact drinking water and recreational bodies of water [6,16]. In Mediterranean systems, the cyclical nature of sediment supply and storage is dependent on the vegetation density and type, which acts as dams that store sediment prior to fire [17]. After a fire, organic matter (OM) from vegetation condenses in the soil horizon and forms a water-repellent layer, which can vary with respect to burn severity [18,19]. This hydrophobic layer facilitates sediment mobilization during precipitation and contributes to landscape development [20]. Loss of vegetation, reduced infiltration [21], soil cohesion [22], and soil water repellency [23] alter surface runoff and erosion potential [24][25][26] as well as suspended sediment discharge and sedimentation [27,28]. These processes are influenced by characteristics such as slope, surface roughness, vegetation, and rainfall intensity [29]. Elevated quantities of sediment, coupled with flashier post-fire hydrological regimes, often induce significant geomorphic responses such as aggradation, incision, bank widening, channel narrowing, and braiding (i.e., [30,31]). The natural mitigation of these processes is the reestablishment of native vegetation cover [32].
Previous studies have observed the different impacts of fire between upland and riparian regions [14,33]. In Mediterranean climates, chaparral-vegetated hillslopes transition into riparian areas where the canopy and soil moisture are different [33,34]. After a fire, native riparian species often exhibit rapid recovery in biomass in comparison to upland chaparral, generally reinstating the geomorphic stability of the channel area [33]. However, non-native vegetation disrupts traditional hydrologic and geomorphic processes. In contrast to areas with native vegetation, watersheds that are overgrown with non-native vegetation can contribute continuous fuel for a fire, playing a key role in the invasive grassfire cycle [8,13,15]. Originally introduced to California for erosion control and building materials [8], Arundo donax is now ubiquitous across all major coastal watersheds in southern California [13]. Arundo donax thrives in urban channels and waterways, growing in monotypic stands that are both highly prone to fire [8,13] and facilitate lateral bank erosion through undercutting and collapse [11]. Furthermore, Arundo donax typically recovers more quickly from fire than native vegetation [35].
While there are studies that investigate the effects of fire on the geomorphology of riparian areas with invasive plants (i.e., [11,12]), there is little research that compares upland and riparian soil characteristics under these conditions. Similarly, there are studies that investigate the physical and chemical changes of soil by fire under typical Mediterranean ecosystem conditions [36][37][38], but do not consider the alterations of soil processes driven by human influence. Thus, this research focuses on the 2018 Del Cerro fire, which burned upland and riparian areas surrounding Alvarado Creek, an urban tributary in San Diego, California to understand surface soil processes. We estimate the post-fire organic matter and soil distributions for six time points during water years 2019 and 2020 at two soil depths in upland (primarily native vegetation) and riparian (primarily non-native vegetation) areas.

Study Area
Alvarado Creek is one of the nine tributaries encompassed within the San Diego River watershed located in San Diego, California ( Figure 1). The climate in San Diego is Mediterranean with consistent temperatures throughout the year. Precipitation mainly occurs during late November through March but is variable from year to year. The Alvarado Creek watershed consists of developed urban land comprised of residential, industrial, and institutional land uses. Sources of surface runoff to the watershed include Lake Murray, urban runoff, and nearby steep areas through concrete channels and culverts. The study reach consisted of native plants such as Artemisia californica and Muhlenbergia rigens in the uplands and dense invasive plant species such as Arundo donax and Washingtonia spp. in the riparian areas [12]. Prior to the fire, non-native riparian vegetation constricted the stream channel and contributed to flooding during precipitation events.
varado Creek watershed consists of developed urban land comprised of residential, industrial, and institutional land uses. Sources of surface runoff to the watershed include Lake Murray, urban runoff, and nearby steep areas through concrete channels and culverts. The study reach consisted of native plants such as Artemisia californica and Muhlenbergia rigens in the uplands and dense invasive plant species such as Arundo donax and Washingtonia spp. in the riparian areas [12]. Prior to the fire, non-native riparian vegetation constricted the stream channel and contributed to flooding during precipitation events.
In June of 2018, the Del Cerro fire consumed 38 acres of vegetation in the Alvarado Creek watershed (Figure 1). The soil burn severity was higher in the riparian area compared to the low to moderate soil burn severity in the uplands [3,12]. Following the fire, Arundo donax thrived and matured rapidly after several months, while the native riparian vegetation had a slower regrowth [12].

Precipitation
Daily precipitation data were acquired from MesoWest [39]. Station DW3640 is part of the Citizen Weather Program in San Diego, California, which is near Alvarado Creek. Station DW3640 was assumed to be representative of Alvarado Creek due to its proximity (less than 2 km away) and minimal elevation difference. We identified storm events that had total rainfall greater than 10 mm following the fire and estimated the maximum 15min rainfall intensity (I15) for each storm event between October 2018 to October 2020. The I15 for each storm was associated with the National Oceanic and Atmospheric Association Atlas (NOAA) ATLAS 14-point precipitation frequency to estimate the magnitude of precipitation events.

Soil Sampling
Soil samples were obtained from two locations along a previously established transect [12,40]. The upland location was located where chaparral vegetation was burned and the riparian was located adjacent to the stream, where Arundo donax was present ( Figure  1). Between October 26, 2018 and October 22, 2020, soils were sampled at an upper (0-15 cm) and lower (15-30 cm) depth at each location (Table 1). Three samples were taken at each location to consider inherent variability. Soils collected on October 26, 2020 represent conditions immediately after the fire, as precipitation is the main driver of geomorphic and hydrologic change after the fire. Thus, measurements were generally before and after the wet seasons or, when feasible, after notable rainfall events.  (Figure 1). The soil burn severity was higher in the riparian area compared to the low to moderate soil burn severity in the uplands [3,12]. Following the fire, Arundo donax thrived and matured rapidly after several months, while the native riparian vegetation had a slower regrowth [12].

Precipitation
Daily precipitation data were acquired from MesoWest [39]. Station DW3640 is part of the Citizen Weather Program in San Diego, California, which is near Alvarado Creek. Station DW3640 was assumed to be representative of Alvarado Creek due to its proximity (less than 2 km away) and minimal elevation difference. We identified storm events that had total rainfall greater than 10 mm following the fire and estimated the maximum 15-min rainfall intensity (I15) for each storm event between October 2018 to October 2020. The I15 for each storm was associated with the National Oceanic and Atmospheric Association Atlas (NOAA) ATLAS 14-point precipitation frequency to estimate the magnitude of precipitation events.

Soil Sampling
Soil samples were obtained from two locations along a previously established transect [12,40]. The upland location was located where chaparral vegetation was burned and the riparian was located adjacent to the stream, where Arundo donax was present ( Figure 1). Between 26 October 2018 and 22 October 2020, soils were sampled at an upper (0-15 cm) and lower (15-30 cm) depth at each location (Table 1). Three samples were taken at each location to consider inherent variability. Soils collected on 26 October 2020 represent conditions immediately after the fire, as precipitation is the main driver of geomorphic and hydrologic change after the fire. Thus, measurements were generally before and after the wet seasons or, when feasible, after notable rainfall events.

Loss on Ignition
Soil organic matter (SOM) content was measured as a loss on ignition (LOI) for two samples per location and depth and in triplicate (twelve total per location). Samples were hand-homogenized, sieved to <2 mm, placed in a weighing dish, and dried overnight (>16 h) in a Yamato DX400 drying oven at 105 • C. Subsequently, 10 g of dried samples were separated into replicates that were placed in aluminum weighing boats of known mass and introduced into a Lindberg Blue Muffle Furnace at 550 • C for four hours [41]. The organic matter (OM) was calculated as: where W dry is the weight of the dry soil and W a f ter ignition is the weight of soil after ignition. Student's unpaired two-tailed t-tests were used to evaluate the statistical significance between sites and depths. If the p-value was less than the significance level of 0.05, we rejected the null hypotheses that the means were equal and the difference was statistically significant.

Soil Classification
We used the ASTM D6913-17 standard test methods for particle-size distribution of soils using sieve analysis [42] and ASTM D4318-17 standard test methods for liquid limit, plastic limit, and plasticity index of soils [43] to classify soil samples for each location and depth in duplicate. Visible foreign objects such as fabrics, shells, or metals were removed from the samples. Soils were pre-treated by removing the organic matter contents and separating the soil particles [44,45]. The samples were submerged into large beakers with distilled water and stirred vigorously. Fine particles were allowed to settle for 10 min, then suspended organic contents were removed using fine strainers. This process was repeated three times to increase the removal of organics. While this process did not completely remove the presence of organics in soil samples, it was minimized to limit the impact on the sieve analysis.
For each location and depth, the soil retained in each sieve was measured. We also expressed the results of the sieve analysis in terms of the percentage of the total weight of soil that passed through different sieves, presented by semi-logarithmic plots of the particle diameters in log scale and the corresponding percent finer. The percentages of gravel, sand, silt, and clay-size particles that were present in the soil were obtained from the particle-size distribution curve. According to the Unified Soil Classification System (USCS) [46], the sizes for gravel, sand, and fines (silt and clay) are defined as >4.75 mm, 4.75 to 0.075 mm, and <0.075 mm, respectively. Fine, medium, and coarse sand are defined as 0.075 to 0.425 mm, 0.425 to 2.0 mm, and 2.0 to 4.75 mm, respectively. To quantify the change from October 2018 with respect to each sampling date, a two-sample Kolmogorov-Smirnov (KS) test was used to compare the distributions. The KS null hypothesis is that distributions were the same at α = 0.05 [47].

Stream-Bed Material
The pebble count method [48] is a quick and simple technique to characterize streambed materials that can be utilized to estimate if disturbances are introducing fine sediments into streams [49]. Following the Wolman (1954) method, we sampled the intermediate axis of 100 pebbles picked from the bed of the stream channel. Using this sample, we established a frequency distribution and estimated the D 50 , which represents the median substrate size, of which 50% of the samples are equal to or smaller than. We utilized this method to monitor the stream-bed material within the channel that intersects the location of the soils sampled ( Figure 1). Samples were available for October 2018, March 2019, and September 2020. We also sampled a reference reach that was unburned and approximately 1.5 km upstream; samples were available for February 2018 and June 2020.

Precipitation
Sixteen storm events occurred over the study period following the Del Cerro fire ( Table 2). The average total rainfall was 38 mm, the average I15 was 3.7 mm/hr, and the maximum recurrence interval (RI) was 2-5 years. There were no significant storms immediately following the fire in June 2018, thus soil sampling on October 28, 2018 was considered a baseline or representative of burned soil conditions. Between the first two samples, 26 October 2018 to 15 February 2019, there were eight storm events ( Table 2, #1-8) with a total of 309 mm of rainfall, and the largest I15 was 8.4 mm/hr, which was a 2-year event. There was a small storm (#9) between the February 2019 and September 2019 samples (<1 year, Table 2). There were four storm events ( Table 2, #10-13) between the September 2019 and February 2020 sampling, with a total rainfall of 199 mm, and the largest I15 was 7.4 mm/hr (RI < 1 year). There were three storm events ( Table 2, #14-16) between February 2020 to June 2020 samples, with a total rainfall of 217 mm and the largest I15 was 7.6 mm/hr. Finally, there were no storms between the June and October samples collected in 2020. Table 2. Precipitation, maximum 15-min intensity (I15), recurrence interval, and total rainfall for significant events during the study period. Soil data were collected in between the dates above and below each horizontal line.

Soil Organic Matter
We evaluated the organic matter present in the upland and riparian sites at two depths ( Figure 2). The average upland and riparian organic matter ranged from 6.04 to 6.56%. Across all the samples available, the average organic matter between the upland and the second year, the largest change in organic matter content observed across both si and depths was in February 2020. The smallest change in SOM for both depths was served on 15 February 2019 for upland areas and on 21 September 2019 for riparian are In the upland areas, the organic matter contents measured in the upper and low depths were statistically different for all sampled dates except 22 October 2020. Me while, the SOM in both depths was similar for all dates except 21 September 2019 in riparian. Over the study period, the upland upper depth SOM was statistically similar February and June of 2020 but different in October 2020. However, the organic mat observed in the riparian upper and bottom layers was statistically different for the th dates following the initial conditions in October 2018 and returned to baseline conditio in June and October 2020.

Initial Soil Classification
We estimated the particle-size distribution of the soil collected from the upper a lower soil depths of the upland and riparian areas ( Figure 3). Due to COVID-19 saf protocols, access to laboratories was limited in 2020 and, thus, sieve analysis was limi and not performed for the upland samples in October 2020 and the riparian samples In the upland areas, the organic matter contents measured in the upper and lower depths were statistically different for all sampled dates except 22 October 2020. Meanwhile, the SOM in both depths was similar for all dates except 21 September 2019 in the riparian. Over the study period, the upland upper depth SOM was statistically similar in February and June of 2020 but different in October 2020. However, the organic matter observed in the riparian upper and bottom layers was statistically different for the three dates following the initial conditions in October 2018 and returned to baseline conditions in June and October 2020.

Initial Soil Classification
We estimated the particle-size distribution of the soil collected from the upper and lower soil depths of the upland and riparian areas (Figure 3). Due to COVID-19 safety protocols, access to laboratories was limited in 2020 and, thus, sieve analysis was limited and not performed for the upland samples in October 2020 and the riparian samples in June 2020 and October 2020. Before the first rainfall event, the percentage of the fine particles in the upland soils was 18%. The liquid limit of fine particles was 41, and the plasticity index was 5; thus, the soil was classified as silty sand [46]. In the riparian area, the percentage of fine particles was 13%. The liquid limit and plasticity index were 45 and 2, respectively; thus, the soil was also classified as silty sand [46]. In some cases, the percentage of fine particles decreased below 12%, indicating that soil classifications were not silty sand after the first storm season (described further below). We also investigated the mass retained from the sieve analysis in general; the percen of mass retained in the upland soils was consistent across the clay and fine and medium sand (<0.075 to 0.425 mm; Figure 4A,C), while the percentage of coarse (4.75 mm) mass retained was lower. There was a tendency towards medium mass retained (0.85 mm) in both depths in the uplands. The upland upper soil had a general increase in mass retention in diameters 0.6 to 4.75 mm and a decrease in diameters <0.075 to 0.425 mm between Oc tober 2018 and February 2019 ( Figure 4A). In the upper soil depth of the uplands, there was over a 530% increase in coarse mass retained from October 2018 to February 2019 and a general return to initial conditions over the following dates. In the upper soils, there was also a 25% and 50% decrease in medium and fine sands, respectively, from October 2018 to February 2019 and a return to initial conditions over the following dates. In the upland lower depth, there was over a 170% and 92% increase in the coarse mass retained from October 2018 to February 2019 and September 2019, respectively, and ~30% decrease in the fine sands for several samples following October 2018. Silt and clay in the upland up per and lower depths also had a significant decrease in the first storm season and re mained below initial conditions for the remainder of the samples.
The upper riparian generally followed a more normal distribution of mass retained ( Figure 4B), where there were higher percentages of mass retained in the fine sands (0.18 to 0.425 mm). The riparian upper soil had a general increase in mass retention in diameters > 0.30 mm and a decrease in sizes < 0.18 mm ( Figure 4B). In the upper riparian, the coarse sand greatly increased from 0.02% in October 2018 to 1.44 to 3.19% and medium sand increased from 6.67% to 14.90 to 18.70% in the three subsequent samples. The fine sand and silt and clay had a continued loss and lower mass retention in all samples after Octo

Particle-Size Distribution
Generally, in the uplands, the percent finer by mass was more well-graded and distributed across all sieve sizes (0.075 to 4.75 mm), where larger percent finer occurred in the larger sieve openings and decreased with smaller sieve sizes ( Figure 3A,C). Both locations generally contained high frequencies of medium and coarse sand particles (0.425 to 4.75 mm). However, the riparian soils were more poorly graded compared to the uplands ( Figure 3B,D), especially for medium and coarse grain sizes for both layers in October 2018 and the bottom layer on 30 June 2020.
In the upland and riparian top and bottom layers, there was a noticeable decrease in the percent finer distributions between October 2018 to February 2019 (Figure 3). The largest difference was observed in the medium sands, 0.425 mm. The upper and lower soils in the uplands had an initial shift or decrease in the frequency of percent finer than by mass across all sand grades during the first wet season ( Figure 3A,C). By the start of the second wet season, observed sediment sizes in the upper upland soils nearly returned to initial conditions, stabilizing in the following soil samples ( Figure 3A). In the lower upland soils, the distributions remain similar across all sampling dates ( Figure 3C). However, the soil type after the first post-fire storm season transitioned to poorly graded sand with silt (fines were between 5-12%). The fine sand in both upper and lower layers did not increase to initial conditions ( Figure 3A,C).
The upper riparian soils generally were poorly graded sand with silt (fines were generally between 5-12%) for all samples. Similar to the uplands, there was an initial shift that occurred during the first wet season in the upper and lower riparian soils ( Figure 3B,D). However, in the upper riparian soils, the particle-size distribution continued to decrease in September 2019, February 2020, and October 2020 ( Figure 3B). In the lower riparian  Figure 3D).
The KS null hypotheses that the distributions were the same for all dates compared to the October 2018 distributions for each soil layer and location, were accepted, thus, all shifts observed were not significant.
We also investigated the mass retained from the sieve analysis in general; the percent of mass retained in the upland soils was consistent across the clay and fine and medium sand (<0.075 to 0.425 mm; Figure 4A,C), while the percentage of coarse (4.75 mm) mass retained was lower. There was a tendency towards medium mass retained (0.85 mm) in both depths in the uplands. The upland upper soil had a general increase in mass retention in diameters 0.6 to 4.75 mm and a decrease in diameters <0.075 to 0.425 mm between October 2018 and February 2019 ( Figure 4A). In the upper soil depth of the uplands, there was over a 530% increase in coarse mass retained from October 2018 to February 2019 and a general return to initial conditions over the following dates. In the upper soils, there was also a 25% and 50% decrease in medium and fine sands, respectively, from October 2018 to February 2019 and a return to initial conditions over the following dates. In the upland lower depth, there was over a 170% and 92% increase in the coarse mass retained from October 2018 to February 2019 and September 2019, respectively, and~30% decrease in the fine sands for several samples following October 2018. Silt and clay in the upland upper and lower depths also had a significant decrease in the first storm season and remained below initial conditions for the remainder of the samples.

Stream-Bed Material
We sampled the stream-bed material adjacent to the upland and riparian soil. T D50 after the fire and prior to the storm season was 115 mm. After the first storm seaso the D50 was 55 mm in March 2019 and returned to 102 mm in September 2020. Meanwhi the reference reach maintained a D50 of 100 and 105 mm for February 2018 and June 202

Organic Matter
Immediately after the fire, the organic matter in the uplands was higher in the upp soils than in deeper soils compared to the riparian area, where no significant differenc were observed (Figure 2). This could be indicative of the buildup of dense and matu chaparral stands in the uplands that accumulated on the surface prior to the fire [50]. Th The upper riparian generally followed a more normal distribution of mass retained ( Figure 4B), where there were higher percentages of mass retained in the fine sands (0.18 to 0.425 mm). The riparian upper soil had a general increase in mass retention in diameters > 0.30 mm and a decrease in sizes < 0.18 mm ( Figure 4B). In the upper riparian, the coarse sand greatly increased from 0.02% in October 2018 to 1.44 to 3.19% and medium sand increased from 6.67% to 14.90 to 18.70% in the three subsequent samples. The fine sand and silt and clay had a continued loss and lower mass retention in all samples after October 2018. The lower riparian showed similar patterns as the upper soil (increase in medium and coarse sands and decrease in fine sand and clay and silt), however, the mass retained across most of the diameters returned to initial conditions in February 2020.

Stream-Bed Material
We sampled the stream-bed material adjacent to the upland and riparian soil. The D 50 after the fire and prior to the storm season was 115 mm. After the first storm season, the D 50 was 55 mm in March 2019 and returned to 102 mm in September 2020. Meanwhile, the reference reach maintained a D 50 of 100 and 105 mm for February 2018 and June 2020.

Organic Matter
Immediately after the fire, the organic matter in the uplands was higher in the upper soils than in deeper soils compared to the riparian area, where no significant differences were observed (Figure 2). This could be indicative of the buildup of dense and mature chaparral stands in the uplands that accumulated on the surface prior to the fire [50]. This was in contrast to the riparian area, which may be influenced by the stream and coarser sand material. Charred vegetation was deposited on the soil surface and subsurface following the fire, which was observed in the October 2018 analysis (Figure 2A). Oguntunde et al. [37] showed that the combination of fire-formed charcoal and charred residue can change the soil physical properties and, thus, affecting post-fire soil processes. Specifically, saturated hydraulic conductivity, total porosity, bulk density, sand fraction, and infiltration rates of soils were altered and attributed to changes in soil structure [37]. Following eight rainfall events in the first post-fire storm season, the February 2019 upland and riparian upper and lower soils had a statistically significant decrease in organic matter, which was attributed to the significant loss of unprotected soil due to fire. In the uplands, the decrease in organic matter between October 2018 and February 2019 was the largest observed. Vegetation mortality during or after a fire can contribute to the washout of organic compounds by throughfall. Additionally, fires can alter soil properties that may induce water repellency and reduce the amount of time water is retained in soils, which has implications for dissolved organic carbon mobilization [38,40,51]. On subsequent dates, the upland soils had an increasing organic matter trend, which likely followed the steady regrowth and deposition of plant matter that is often observed in the post-fire recovery of chaparral vegetation [50,52]. Further, we observed that the organic matter in the lower upland soils was relatively stable compared to the upper uplands and riparian soils ( Figure 2C).
While the vegetation biomass of Arundo donax was greater than native chaparral in the uplands immediately after the fire, the uplands generally had a higher organic matter content. The riparian area was dominated by Arundo donax, which likely impacted the difference in the observed organic matter patterns. While the upland soils appeared to have a more stable increase, there was more variability observed in the riparian organic matter. Both locations were sensitive to rainfall, which was noticeable in February 2020, and may be due to decay of the deposited plant material and incorporation into the SOM pool under wetter conditions. Similar to the uplands, there was a significant decrease between October 2018 and February 2019. However, there was a second and larger decrease in organic matter in the riparian soils between February 2020 and June 2020, which is attributed to the hydrologic and vegetation dynamics. The higher variability in the riparian soils may be an indicator of the sensitivity of invasive vegetation to available water. Grasses are particularly sensitive to rainfall due to shallow roots [53], which may contribute to the rapid growth and senescence of Arundo donax and, consequently, SOM. Additionally, Arundo donax has extremely dense root structures mostly within a depth of 30 cm and the root diameter can range between 1.0 mm to 20.0 mm [11], which can contribute to more mobile soils that can dislodge the organic horizon. This contrasts the deep penetrating or shallow wide-spreading or branching fibrous roots of some southern Californian chaparral [54]. The difference between the root structures may alter the presence of organic matter. Further, we followed the standard LOI technique to separate sediment and residues larger than 2 mm prior to estimating the organic content, which could have contributed to overlooking roots larger than 4.25 mm.

Sediment Transport Patterns
There was an overall decrease in the particle size across all soils (Figure 3). The greatest change occurred in the medium sands, while the fine sands appeared to be impacted the longest (Figure 3). We observed the release of smaller sediment and coarsening of the soil (Figure 4), which is a natural geomorphic response to stabilize the landscape with larger grain sizes [55][56][57]. This also corroborates field observations by Mathews et al. [12], who noted that fine sediments were removed from the uplands, where rainfall intensity surpassed the saturated hydraulic conductivity of the finest layer of the hillslope. Severely heating soils can also alter the particle-size distribution and result in coarser textures due to an increase in silt fraction [36]. Coarsening patterns also corresponded to our observations in the adjacent channel, where there was a removal of the fines, leaving coarser particles in the streambed, and suggesting the role of the fire.
In the uplands, the first post-fire storm season had the largest shift ( Figure 3A), which was likely due to the loss of vegetation. Vegetation acts as dams that retain sediment; when vegetation is disturbed or removed by fire, sediment is released to the stream margins due to gravity through a mechanism known as dry ravel [17]. Further, during rain events, unprotected soils are more mobile and sediment is more susceptible to transport [27]. Finer sediment leaves the soil profile, increasing the coarser particles after the first wet season [57], which we observed in this study (Figure 4). Decreased vegetation cover followed by rainfall was the primary driver of fine sediment removal; expedited and elevated quantities of fine sediment transport were also noted by studies such as Syvitski et al. [30] and Chin et al. [31].
The natural mitigation of fire-driven erosion and the resultant geomorphic change is the reestablishment of native vegetation cover [32]. In Alvarado Creek, very little vegetation recovery occurred on the upland hillslope between June 2018 to the end of March 2019 and corresponded with the largest measured geomorphic response following the first post-fire storm event, which progressively decreased thereafter [12]. This is common in systems with chaparral species, which have deep root systems that can stabilize soil mass and reduce surface soil movement [54]. Thus, the immediate release of sediment, coupled with an armored hillslope due to deep chaparral root systems [54] and large clast grain structure [57], can result in a smaller sediment discharge, despite the flushing of soils by rain events. In Alvarado Creek, the removal of the topsoil and ash deposits after the first storm event was documented in the field by Mathews et al. [12], leaving only the charred layer that was beneath the initial horizon and a coarser grain size distribution. The removal of the smaller grain sizes likely occurred by erosional mechanisms such as rain splash, overland sheet flow, or thin debris flows [58,59], resulting in a coarsening of soils. In our study area, the higher frequency of larger grain sizes may have contributed to the geomorphic stability of the landscape after the first storm event by increasing the resilience to shear stress that could initiate particle transport during future precipitation events [55,56]. This is due to the low burn severity observed in the uplands, leaving the native chaparral roots mostly intact, and fortifying the structural integrity of the hillslope from hydraulic events [57].
The sediment processes in the riparian zone were more dynamic. In the floodplains of southern California, the high tensile strength and dense root structure of Arundo donax can induce bank instability and subsequently bank collapse during severe floods [8,11]. The high variability in soil characteristics and preferential retention of coarse and medium sands in the riparian soils ( Figure 3B,D) were attributed to the fast vegetation regrowth dominated by monocultures of Arundo donax and resultant geomorphic instability [12]. We hypothesize that Arundo donax outcompeted native vegetation and facilitated a coarsening of sediment and a decrease in fine sand, silt, and clay particles. The roughness and rhizome root structures imposed by Arundo donax infestations decreased hydraulic velocities and flooding, which contributed to preferential sediment deposition, armoring of the channel, and lateral bank movement [60,61]. The aggradation of the floodplain, channel incision, and obliteration of step-pool morphology observed by Mathews et al. [12] resulted in a sustained topographic change in the riparian area surrounding Alvarado Creek and supports our observations.

Conclusions
This study investigated soil organic matter and particle-size distribution in Alvarado Creek, an urban Mediterranean system, to understand the post-fire soil patterns in upland hillslopes with native vegetation and riparian areas with non-native vegetation after fire. In the upland hillslopes, with predominantly native chaparral, fine grains comprised the main sediment transported. This research builds upon our knowledge of how the presence of Arundo donax in urban Mediterranean systems can disrupt soil regimes. The presence of monocultures of Arundo donax and rapid regrowth after fire had a significant impact on the retention of coarse sediment, which armored channels to decrease faster flows. This has implications for river restoration activities that are intended to manage non-native vegetation and alleviate fire hazards and channel capacity to reduce flooding. Further, the presence of Arundo donax in stream channels and riparian areas creates a more dynamic system and imposes higher risks to future fire events, which can further exacerbate geomorphology and can create hazardous conditions.