Estimating Net Photosynthesis of Biological Soil Crusts in the Atacama Using Hyperspectral Remote Sensing

: Biological soil crusts (BSC) encompassing green algae, cyanobacteria, lichens, bryophytes, heterotrophic bacteria and microfungi are keystone species in arid environments because of their role in nitrogen- and carbon-ﬁxation, weathering and soil stabilization, all depending on the photosynthesis of the BSC. Despite their importance, little is known about the BSCs of the Atacama Desert, although especially crustose chlorolichens account for a large proportion of biomass in the arid coastal zone, where photosynthesis is mainly limited due to low water availability. Here, we present the ﬁrst hyperspectral reﬂectance data for the most wide-spread BSC species of the southern Atacama Desert. Combining laboratory and ﬁeld measurements, we establish transfer functions that allow us to estimate net photosynthesis rates for the most common BSC species. We found that spectral differences among species are high, and differences between the background soil and the BSC at inactive stages are low. Additionally, we found that the water absorption feature at 1420 nm is a more robust indicator for photosynthetic activity than the chlorophyll absorption bands. Therefore, we conclude that common vegetation indices must be taken with care to analyze the photosynthesis of BSC with multispectral data.


Introduction
In arid regions, biological soil crusts (BSCs) are keystone organism communities consisting of cyanobacteria, green algae, bryophytes, heterotrophic bacteria, microfungi and lichens [1]. Many communities encompass species of two coalescing groups: the heterotrophic part is represented by fungi and the phototrophic component by cyanobacteria and/or eukaryotic green algae [2]. The heterotrophic part is responsible for the structure and water supply of the organism and protects the phototrophic organism against ultraviolet radiation by the production of pigments [3]. The phototrophic component provides carbohydrates from photosynthesis. In the case of the absence of water, the phototrophic component is in a dormant state. Under natural conditions where periods of desiccation and rehydration alter frequently, photosynthesis starts usually within seconds after water becomes available [4].
BSCs dominate large areas worldwide that are ruled by harsh abiotic conditions including both hot and cold deserts, where they form the most productive microbial communities [5]. Beside soil protection against erosion and dust trapping by concatenating organic and inorganic material at 1.
To describe for the first time the hyperspectral reflectance signal of BSCs in the Atacama Desert under different water availability conditions. 2.
To test the suitability of hyperspectral remote sensing data for the estimation of net photosynthesis (NP) of BSCs. 3.
To test whether a robust transfer function can be established between NP and hyperspectral images acquired under field conditions, which allows mapping NP across larger scales.

Area of Investigation
The study area is part of the Pan de Azúcar National Park, which is located in the southern part of the Atacama Desert between 25 • 53 and 26 • 15 S and 70 • 29 and 70 • 40 W along the Pacific coast in Chile ( Figure 1). Local terrain is dominated by the steep mountain ridge close to the coast reaching altitudes up to 850 m a.s.l. Heading inland, the terrain slightly decreases to altitudes between 400 and 700 m a.s.l.
Annual rainfall is very low (below 10 mm, [28]), but a high variability leads to episodic and catastrophic precipitation events [29]. Aside from those episodic rain events, the local ecosystems are fed by water from fog and dew [30][31][32]. Mean daily temperatures are between 13 • C and 20 • C during austral winter (July) and summer (January), respectively [28]. Relative humidity under clear sky conditions is between 80 and 85% [28].
The harsh environment leads to a specialized vegetation, which is dominated by cacti (e.g., Eulychnia saint-pieana), species from the genus Euphorbia and the landscape-dominating BSCs. For an overview of the higher vegetation in the park, see [28,33]. Species of BSC change along the humidity gradient from the coastal ridge (high influence of fog) to the hinterland (more arid and less fog). At the coastal ridge, a coverage of 30% was observed mainly dominated by the four species Acarospora cf. gypsi-deserti, Caloplaca santessoniana ad int. [34], Placidium cf. velebiticum and Rinodina sp., all of which are chlorolichens.

Sampling of BSCs
The two most abundant lichen species Acarospora cf. gypsi-deserti and Caloplaca santessoniana were sampled by pressing a sterile 9-cm Petri dish 1 cm into the soil. Excess soil was removed with the Petri dish lid. Samples were air-dried in the field immediately after collection for 2-3 days, until no further condensation occurred. The dry and sealed crust samples were preserved at −20 • C until further processing. For this study, the samples were slowly defrosted under air-tight conditions before they were used for the laboratory analysis.

Laboratory Analysis
CO 2 gas exchange measurements were conducted according to [35] with the two most abundant crustose chlorolichen species Acarospora cf. gypsi-deserti and Caloplaca santessoniana with soil adhering to their rhizines. Before measurements, the intact lichen samples underwent a reactivation procedure of one day exposure at 8 • C in the dark. Afterward, the samples were placed in the gas exchange cuvette and sprayed with sterile, filtered water to activate their metabolism 24 h prior to measurement. Ahead of the measurements, full water saturation was achieved by submerging the samples in water for ten minutes. Excessive water and droplets were carefully shaken from the sample before measurements. CO 2 gas exchange measurements were conducted under controlled laboratory conditions using a mini-cuvette system (GFS 3000, Walz Company, Effeltrich, Germany). The response of net photosynthesis (NP) to water content was determined for the lichens. Complete desiccation cycles (from the water-saturated phase to the air-dried status) were measured under saturating light (800 mE) and ambient CO 2 at 24 • C. Samples were weighted between each measurement cycle, and the water content was calculated as mm precipitation equivalent. Sample dry weight was determined after 3 days in a drying oven (Heraeus Instruments T6P, Thermo Fischer Scientific Inc., Waltham, MA, USA) at 60 • C. The CO 2 gas exchange of the samples was related to the lichens' surface and chlorophyll content, the latter determined after [36].

Hyperspectral Measurements
Hyperspectral reflectance data have been acquired with two different instruments. In the field, a Specim ENIR hyperspectral camera (Specim, Spectral Imaging Ltd., Oulu, Finland) has been used, which acquires images consisting of 256 bands sensitive to electromagnetic radiation between 600 nm and 1600 nm. The spectral configuration allows continuously acquiring reflectance values with full-width-half-max values between 4.08 nm (red part of the electromagnetic radiation) and 4.17 in the near-infrared (NIR). In early summer in December 2017, hyperspectral samples were taken under sunny conditions close to noon to minimize bi-directional effects due to varying illumination geometries at acquisition time. The camera is a line scanner, which was mounted on a rotary stage. Rotation speed was automatically adjusted to the integration time so that a full spatial coverage of the target without overlapping pixels was achieved. To convert the raw values into reflectances, a Spectralon white standard (reflectance ≈ 95%) was arranged in the images. Sensor noise was removed by acquiring images under dark conditions immediately before and after each scan. Therefore, the lens of the camera was closed by a screw-cap. Both dark current images were averaged column-wise since noise in the instrument changes across the scanner array. Reflectance values were then calculated: Here, R λ,c is the reflectance in band λ and column c. counts λ,c are the raw count values of the sensor, and d λ,c and w λ are the averaged dark current and white standard values, respectively.
All locally-dominating species of BSCs were sampled and arranged in rows. From each species, at least 5 replicates have been sampled ( Figure 2). After acquisition of the first image under dry conditions, the samples and the surrounding soil were watered. In total, seven images were acquired over a 20-min time period until the samples were visually dry, again.
The second instrument was a Spectrometer (HandySpec Field, Tec5 AG, Oberursel, Germany), which was used in the laboratory with an artificial light source. We could not use the hyperspectral camera in the laboratory because the focus range and size of the instrument did not allow properly measuring the samples in the cuvette. The spectrometer is sensitive to radiation between 305 nm and 1705 nm at a spectral resolution of 1 nm and has two channels crossing at roughly 1050 nm. To convert the raw values into reflectances, the same approach as described above has been used, and Equation (1) has been applied. At time intervals of five minutes, the samples were removed from the GFS 3000 to weigh them and to acquire hyperspectral reflectance values. At each time interval, three independent measurements were taken and averaged to reduce the effects of small movements of the sample and the measuring device. Unfortunately, the sample could not be removed from the cuvette due to fragility. This causes the radiation to interact with the plastic of the cuvette and/or the background may partly influence the absolute reflectance values measured by the spectrometer.

Hyperspectral Analysis
The hyperspectral analysis was conducted in R statistical software [37] using the raster- [38] and the hsdar-packages [39]. To remove noise from the raw spectra, a Savitzky-Golay filter has been applied (filter length of 15 and 25 bands for Specim ENIR and HandySpec Field, respectively). To cope with the different spectral configuration between the spectrometer and the hyperspectral camera, the spectra acquired with the spectrometer were resampled to the bands of the hyperspectral camera by considering Gaussian spectral response functions defined by the center and full-width-half-maximum values of the Specim ENIR sensor.
Since the spectra measured in the laboratory may be partly contaminated by background reflection, we compare two different approaches to remove artifacts arising from the measurement setup. The first one is based on normalized ratio indices (NRI), which were calculated for all possible band combinations [40]: Here, R is the reflectance at wavelength λ. It can be assumed that the spectra in the laboratory acquired with an artificial light source may differ systematically from those taken in the field. Therefore, the NRI-values were further normalized using the dry spectra as reference: The second approach to normalize the spectra was performed by applying continuum removal ensuring that only local differences within spectral absorption features were taken into account.
We constructed a segmented continuum line to each spectrum and calculated the band depth values [41]: Here, R λ and CV λ denote the reflectance and continuum line value at wavelength λ. For each identified feature, the integral and the width between lower and upper full-width-half-maximum values have been calculated.
The relationship between water content and NP for BSC is known to be non-linear [42,43]. Therefore, species-specific second degree polynomial functions were fitted to the hyperspectral signal of Caloplaca santessoniana and Acarospora cf. gypsi-deserti and NP as the dependent variable. Independent variables were each normalized ratio index and each variable of the absorption features: Here, x denotes the normalized ratio index or absorption feature variable. Please note that we did not directly fit models to the band depth values and that only univariate models were fitted because the indices are highly correlated. Based on the R 2 -values of the models, the model with the highest proportion of explained variance has been selected.
The normalized ratio index with the highest predictive performance regarding NP under laboratory conditions was calculated from all hyperspectral images acquired in the field. Predictive performance was compared by the R 2 -values of the models. Normalization of the index values has been performed according to the normalization of the spectrometer data by using the image acquired under dry conditions as the reference. Afterward, the species-specific models derived from the laboratory measurements were used to predict NP in the hyperspectral images. For the models with the highest proportion of explained variance derived from continuum removal variables, the same approach has been conducted.
The soil close to the BSCs was used as a control to determine whether the relationships between NP and hyperspectral signal is a consequence of the reflectance properties of the BSCs or the soil beneath the crusts. Therefore, pixels covered only by soil without any crusts have been extracted from the hyperspectral images taken in the field and processed in the same manner as for the pixels covered by BSCs.

Results
Spectra measured in the laboratory differed clearly between dry and wet stages ( Figure 3). Under dry conditions, the spectra of all species were similar to typical soil spectra, because the chlorophyll absorption feature at approximately 680 nm was only marginally developed. If specimens were watered, the spectra changed immediately. The strongest differences could be observed at the chlorophyll absorption bands and at the water absorption bands around 1450 nm. Reflectance values in the NIR increased, whereas those in the red and mid-infrared decreased. The dryer the specimen became, the smaller were the spectral differences compared to the spectra acquired under initial dry conditions. After 4-5 h, the spectra were almost identical to the initial ones.
Spectra were also different between wet and dry stages if hyperspectral data were acquired under field conditions ( Figure 4). Here, the pattern was less pronounced compared to the laboratory spectra. Spectral signatures of the BSC species differed from soil reflectances in terms of shape only in the red part, where the soil reflectance was characterized by a slower increase towards the NIR. Both soil and the BSC spectra had a distinct absorption feature at approximately 1150 nm. Large spectral differences were found among the BSC species. Similar reflectance values were observed for the two wide-spread species Caloplaca santessoniana and Acarospora cf. gypsi-deserti, whereas reflectances of Placidium cf. velebiticum were generally substantially lower. Additionally, no chlorophyll absorption feature could be observed for Placidium cf. velebiticum in the dry stage. Rinodina sp. was characterized by intermediately high reflectances and an intermediately developed chlorophyll absorption feature under dry conditions.    The relationship between NP and spectral reflectances was evaluated for the two most common species Caloplaca santessoniana and Acarospora cf. gypsi-deserti. Here, normalized ratio indices revealed similar patterns in both species ( Figure 5). Close relationships were found for indices calculated from reflectances in the NIR including the red edge and the red part of electromagnetic radiation. Additionally, indices derived from bands at the water absorption feature at 1420 nm and any band from the NIR revealed close relationships. In contrast, band combinations from the NIR and the red-edge were not useful to estimate NP. In general, the R 2 -values were higher for Caloplaca santessoniana than Acarospora cf. gypsi-deserti. Regarding Caloplaca santessoniana, the best index was calculated from reflectance values at 1475 nm and 1433 nm, whereas the best index for Acarospora cf. gypsi-deserti used bands at wavelengths of 1388 nm and 679 nm (Table 1). Consequently, for both species, at least one band was selected from the water absorption feature around 1420 nm. For both species, differences between the best and the indices with consecutively lower explanatory power were small. Four absorption features have been identified by applying continuum removal. The first one was the chlorophyll absorption feature at approximately 680 nm. The second and third feature were smaller and located around 950 nm and 1100 nm. The fourth one around 1420 nm was the large water absorption feature. Variables from all four absorption features have been tested for their usability to predict NP of Acarospora cf. gypsi-deserti and Caloplaca santessoniana ( Table 2). The best variable for both species was the integral of the water absorption feature around 1400 nm. Additionally, the width of the absorption feature around 1400 nm performed well for the prediction of NP in Caloplaca santessoniana, but not in Acarospora cf. gypsi-deserti. All other features including the chlorophyll absorption feature around 680 nm were not useful to estimate NP of both BSC species.  A non-linear relationship between NP and water availability was observed for Acarospora cf. gypsi-deserti (not shown). In the low and high water availability stages, NP was close to zero, and the highest NP was measured at an intermediate water content of approximately 2 × 10 −4 mL m −2 . For normalized ratio indices and continuum removal, the index and the absorption feature derived variables that revealed the closest relationship to NP were tested for the applicability to the spectra acquired under field conditions. Regarding normalized ratio indices, the best index to predict NP of Acarospora cf. gypsi-deserti was derived from bands at 1388 nm and 679 nm and revealed a clearly non-linear relationship with NP (Figure 6a). In the dry and wet stages, the index values were small and large, respectively. The highest values of NP were observed for intermediate index values. If applied to spectra acquired under field conditions, no pattern within the estimates over time could be observed (Figure 6b).  For Caloplaca santessoniana, a close linear correlation was observed between NP and the water content of the sample in the laboratory. The best normalized ratio index for Caloplaca santessoniana was derived from reflectance values at 1475 nm and 1433 nm. In contrast to Acarospora cf. gypsi-deserti, the relationship was almost linear (Figure 6c). If applied to field spectra, a clear pattern was observed (Figure 6d). NP was highest directly after artificial watering of the samples in the field and decreased until the samples were visually dry. However, the highest and lowest predicted NP were one order of magnitude lower compared to the laboratory measurements.
The variable derived from continuum removal with the highest predictive performance was the integral of the feature at 1420 nm for both species. To test if the relationships identified between NP and the hyperspectral signal are a consequence of the BSCs or the soil beneath the crust, the change in the integral of the absorption feature at 1420 nm for pixels covered only by soil was investigated over time (Figure 7). In general, the integral values slightly declined during drying-out of the soil. Regarding Acarospora cf. gypsi-deserti, a non-linear relationship was observed, which was similar to the relationship derived for the normalized ratio index (Figure 8a). If the empirical transfer function was applied to the field spectra, a clear pattern was observed in NP predictions (Figures 8b and 9a). NP was highest in all individuals in the first image after artificial watering and decreased with progressive drying-out. The estimates were in the range of the laboratory measurements.  For Caloplaca santessoniana, the relationship was linear except one measurement causing the regression line to decrease at the higher range of the values (Figure 8c). As for Acarospora cf. gypsi-deserti, a clear pattern in the predictions of NP for the field samples was observed, which were in the range of laboratory measurements (Figures 8d and 9b).

Discussion
We found differences in the spectral signal of different species of BSCs. Irrespective of the water content, spectral properties of Caloplaca santessoniana and Acarospora cf. gypsi-deserti were similar, whereas Rinodina sp. and Placidium cf. velebiticum were characterized by generally lower reflectance values. Water content affected the spectra of all species mainly in the water absorption bands around 1420 nm and the chlorophyll absorption feature at approximately 680 nm. The spectral signal in both absorption features changed within minutes after water became available to the BSC, which confirms the result of [44] for Australia, while NP initialization is faster than the observations reported by [12] from the Arabian Peninsula. In contrast to the results from Australia, the spectral signal was almost identical after the sample became dry both in the laboratory and the field [44]. Since the distinct spectral feature at 1120 nm is present in the soil and all BSC species, it is most likely a consequence of an andesite-rich bedrock [45]. This shows that even spectra from the laboratory samples were partly influenced by the reflection of the soil beneath the crust, which could not be completely removed from the samples.
Robust relationships between water content and NP were found for the two most wide-spread species Acarospora cf. gypsi-deserti and Caloplaca santessoniana in the laboratory. The non-linear, positively-skewed response of Acarospora cf. gypsi-deserti corresponds well to prior measurements of other species, e.g., of the genus Lecidella [42]. The fast and quasi-linear relationship between NP and water content for Caloplaca santessoniana has been observed, e.g., for genus Peltula [46]. Hyperspectral reflectances changed with the water content of the sample, which allowed deriving robust transfer functions between the NP measurements and the spectra acquired in the laboratory. In this respect, the most relevant spectral bands for NP estimation were situated in the water absorption feature at 1420 nm and the chlorophyll absorption bands. Since no direct measurements of chlorophyll activity via, e.g., chlorophyll fluorescence were available to date, it is unknown if water content or chlorophyll activity is the main influencing factor on normalized ratio indices.
The spectra in the laboratory may be contaminated partly by radiation that interacted with the plastic cuvette, because the fragile samples could not be removed for hyperspectral data acquisition. Consequently, two different normalization methods of the spectra were compared: normalized ratio indices and variables based on absorption features derived by continuum removal. Regarding the normalized ratio indices, we found that good predictors of NP for both species were derived from bands in the water absorption feature, the NIR between 1100 nm and 1300 nm and the red part lower than 700 nm. Contrasting patterns between normalized ratio indices and NP between both species were mainly found in the red-edge where NP of Acarospora cf. gypsi-deserti could not be predicted by indices derived from bands close to the shoulder of the red-edge, whereas models for Caloplaca santessoniana revealed a generally good predictability in the same spectral range. This shows that hyperspectral responses to changes in water content were partly species dependent. Since the red-edge is a frequently-used indicator for NP (e.g., via the NDVI and other vegetation indices), care must be taken that the estimates of NP from vegetation indices are validated against data from other sensors such as chlorophyll fluorescence imagers [18].
The variables derived from the absorption features defined by continuum removal had a generally lower predictive performance for NP compared to the normalized ratio indices. The only variable that was related to NP was the integral of the water absorption feature at 1420 nm. Surprisingly, no relationship between the chlorophyll absorption feature and NP was detected for either species after continuum removal. This underlines that common vegetation indices should be treated with care if they are used as predictors for chlorophyll activity of BSCs in space and time.
Irrespective of the applied normalization method, the parameters in the polynomial fits differed highly among both species as a consequence of the species-specific form of the relationship between water content and NP. This clearly shows that no universal empirical transfer function across different BSC species between NP and the hyperspectral reflectance signal can be derived. Instead, species-specific relationships at least for the most wide-spread species must be established for any estimation of NP in space and time based on remotely-sensed data in the future.
We tested whether our empirical transfer functions derived from laboratory measurements can be applied to field data even though two different sensors have been used. We found that normalized ratio indices were not useful in our case because estimates of NP based on field data differed from laboratory measurements partly by more than one magnitude. Especially for Acarospora cf. gypsi-deserti, the empirical transfer functions derived from normalized ratio indices were unable to reproduce any kind of pattern between NP and water content. One reason for the failure of normalized ratio indices for Acarospora cf. gypsi-deserti could be that differences between the samples in the field and the laboratory sample (e.g., thickness of BSC, chlorophyll content) caused the empirical transfer function to be ambivalent, because only reflectance values at two wavelengths were considered. The other factor influencing the transferability of the laboratory-derived functions to field data is the usage of two different sensors in this experiment. This was unfortunately necessary because no line scanner camera could be used in the laboratory since the focus range and size of the instrument did not allow properly measuring the samples in the cuvette. To cope with this limitation, we spectrally integrated the reflectance values of the spectrometer to the bands of the hyperspectral camera assuming Gaussian spectral response functions, because no spectral response function could be provided by the manufacturer. If this assumption is not fully valid, it may partly explain why spectral features in parts of the electromagnetic radiation with high heterogeneity (e.g., the red edge) were not useful to predict NP of BSCs in this study.
Although the predictive performance of the integral at 1420 nm was lower compared to normalized ratio indices in the laboratory spectra, the predictions of those empirical transfer functions based on field data were in astonishingly good agreement with the theoretical expectations [42]. We found that after watering in the field, hyperspectrally-predicted NP for all individuals decreased over time, which was in the range of the laboratory estimates. In the meantime, changes in the hyperspectral signal of the soil were much less pronounced compared to the BSCs. This implies that the relationships between NP and the hyperspectral signal were a consequence of the reflectance properties of the BSCs rather than the soil beneath the crust. However, it has to be noted that NP was measured under controlled laboratory conditions in this study. Consequently, the absolute predicted values based on field data must be considered carefully. Comparing the predictions between the different normalization methods, we conclude that continuum removal is a more robust technique compared to normalized ratio indices, especially under different conditions (e.g., field vs. laboratory) and when different sensors are used. Additionally, the usage of full absorption features reduces the sensitivity of the empirical transfer functions to differences among individuals such as the thickness of the BSCs.
The spatial resolution in our analysis is extremely high, which leads to the question of whether the results can be applied to data featuring lower spatial resolutions such as data acquired by satellites or aircraft. Since our transfer functions have been fitted to laboratory measurements where the entire field of view was covered by BSCs, any kind of transfer to mixed pixels is speculative. Therefore, additional investigations based on field data are required in order to verify our results for large-scale data. Irrespective of the mixed pixel effects, our finding that the most robust spectral indicator for NP of BSCs is located in the water absorption feature around 1420 nm causes severe implications if satellite or airborne hyperspectral data will be used for the investigation of NP. First, the absorption feature at 1420 nm does not only depend on the water content of the BSCs, but also on water vapor content of the atmosphere (e.g., [47]), which could be neglected in the present study because distances between the sensor and the BSCs were small. Consequently, if satellite or airborne data will be used, a proper atmospheric correction is mandatory. Taking the low signal to noise ratio of current hyperspectral sensors into account (e.g., [48]), it must be questioned if present hyperspectral satellite data are appropriate to estimate NP of BSCs in the Atacama Desert. Second, the usage of satellite data is generally restricted since the BSCs are only active for a short time after water becomes available. In the Atacama Desert, the latter strongly depends on fog input, which impedes satellite data in remotely sensing NP of BSCs, because satellite scenes of active BSCs are cloud-covered.

Conclusions
The aims of the present study were three-fold:

1.
We described for the first time the hyperspectral reflectance signal of BSCs in the Atacama Desert under different water availability conditions. Here, we could demonstrate that hyperspectral reflectance signals among wide-spread species of BSCs in the Atacama Desert differed largely, but water content affected the spectra in a similar manner. Changes in water availability immediately influenced the chlorophyll absorption bands in the visible and the water absorption bands in the near-infrared part of the electromagnetic spectrum.

2.
We tested the suitability of hyperspectral remote sensing data for the NP estimation of BSCs and found that the relationship between water content and NP is highly species dependent, urging the need for species-specific empirical transfer functions between NP and hyperspectral reflectance values. In this respect, the species-dependent transfer functions between the size of the water absorption feature at 1420 nm were better predictors than any variable derived from chlorophyll absorption bands.

3.
We tested whether the transfer function derived under laboratory conditions can be applied to hyperspectral images acquired in the field, which allows mapping NP across larger scales. Our results were in astonishingly good agreement with the theoretical expectations if the transfer function relied on the water absorption feature at 1420 nm, suggesting that the spectral patterns between laboratory and field measurements were highly comparable and underlining the general possibility for area-wide predictions in the field. However, the use of the water absorption bands limits the usability of space-and air-borne data in future applications because of the strong water absorption in the atmosphere accompanied by the low signal to noise ratio of current sensors. Therefore, we suggest using drones flying at low elevations above the ground to reduce the influence of the atmosphere on the reflectance values measured at the platform. Such kinds of data can be used to provide area-wide NP estimations of BSCs in the southern part of the Atacama Desert in the future, where BSCs are keystone organisms providing key ecosystem functions such as protection against soil erosion, weathering of nitrogen and phosphorus and dust trapping.