Microbial and geo-archaeological records reveal the growth rate, origin and composition of desert rock surface communities
- 1Zuckerberg Institute for Water Research, Jacob Blaustein Institutes for Desert Research, Ben-Gurion University of the Negev, Sede Boqer Campus 8499000, Israel
- 2Southern Region, Israel Antiquities Authority, Omer 84965, Israel
- 3Soil and Water Research Infrastructure and Institute of Soil Biology, Biology Centre CAS, České Budějovice, Czechia
- These authors contributed equally to this work.
Correspondence: Roey Angel (firstname.lastname@example.org) and Osnat Gillor (email@example.com)
Biological rock crusts (BRCs) are ubiquitous features of rock surfaces in drylands composed of slow-growing microbial assemblages. BRC presence is often correlated with rock weathering, soiling effect or mitigating geomorphic processes. However, their development rate is still unknown. In this work, we characterised and dated BRCs in an arid environment, under natural conditions, by integrating archaeological, microbiological and geological methods. To this end, we sampled rocks from a well-documented Byzantine archaeological site and the surrounding area located in the central Negev, Israel. The archaeological site, which is dated to the fourth to seventh centuries CE, was constructed from two lithologies, limestone and chalk. BRC started developing on the rocks after being carved, and its age should match that of the site. Using stable carbon and oxygen isotope ratios, we confirmed the biogenic nature of the crusts. The BRC samples showed mild differences in the microbial community assemblages between the site and its surroundings, irrespective of lithology, confirming the dominance of aeolian inoculation sources. All BRCs were dominated by Actinobacteria, Cyanobacteria and Proteobacteria. We further measured the BRC thickness on 1700-year-old building stone blocks and determined it to be 0.1–0.6 mm thick. Therefore, a BRC growth rate was estimated, for the first time, to be 0.06–0.35 mm kyr−1. Our dating method was then validated on a similar archaeological site located ca. 20 km away, giving comparable values. We propose that BRC growth rates could be used as an affordable yet robust dating tool in archaeological sites in arid environments.
In arid and hyper-arid environments where abiotic processes are considered the primary contributor to landform formation, barren rock surfaces, free of vegetation, are a ubiquitous feature (Owen et al., 2011). These surfaces are exposed to multiple stress factors, such as lack of water, high radiation and extreme temperature fluctuations and therefore represent the edge of biotic existence on Earth (Viles, 2008). Such rock surfaces serve as a habitat for microorganisms by providing colonisation strata on the outside (epilith; Pointing and Belnap, 2012), in the inner pores (endolith; Büdel and Wessels, 1991; Friedmann and Kibler, 1980) or underneath the rock (hypolith; Wierzchos et al., 2013). The microorganisms colonising rocks form a hardy biomineral rock coating known as a biological rock crust (BRC; Gorbushina, 2007), which is common in most arid and hyper-arid regions worldwide (Gorbushina, 2007; Lebre et al., 2017; Pointing and Belnap, 2012). BRC structure, when developed on carbonate-rich minerals (e.g. limestone), often consists of a hardpan-laminated structure composed of masses of micritic to microsparitic carbonate layers interbedded with microbial coatings. Such laminated structures were previously reported by Alonso-Zarza and Wright (2010) as terrestrial calcretes. From a microbial perspective, BRCs are typically colonised by Actinobacteria, Cyanobacteria, and other phototrophic and heterotrophic bacteria, yet they lack, or nearly lack, archaea, fungi and algae (Lang-Yona et al., 2018). BRC inoculum has been proposed to originate from settled dust (Viles, 2008) or the surrounding soil (Makhalanyane et al., 2015). Wieler et al. (2019) noted that the microbial composition of the BRC was distinct from the soil communities but shared some of the communities found in the settled dust. Local contribution of aeolian material (Crouvi et al., 2008) may also affect the microbial communities found in the settled dust. The influence of biota on landscapes is a topic that has remained mostly unexplored (Dietrich and Perron, 2006; Viles, 2019). BRCs have been noted to play a crucial role in the functioning of arid and hyper-arid environments because of the limited activity of plants and soil (Pointing and Belnap, 2012). A range of geobiological roles have been suggested for BRCs, including rock weathering (Garcia-Pichel, 2006; Warscheid and Braams, 2000), soiling effect (Viles and Gorbushina, 2003), deterioration of architectural heritage sites (Cabello-Briones and Viles, 2017) and mitigating geomorphic processes (Mcllroy de la Rosa et al., 2014). BRCs have also gained attention from astrobiologists, who point out that such communities, when located in deserts, may act as analogues for possible life on Mars (Corenblit et al., 2019). Yet, despite the suggested link between rock surface morphologies and microbial activity, one basic but important question remains unanswered: how quickly do BRCs develop after a rock surface has been exposed? Illustrating the origin, composition and growth rates of BRCs may shed light on the timescales on which natural processes (e.g. desertification, land degradation) take place in arid regions. Currently, there is no available information about the development rates of epi- and endolithic lithic biofilms on rock surfaces in deserts. Formation rates of sub-aerial BRCs have been noted to be very slow (Pointing and Belnap, 2012) but have not been backed by numbers. The goal of this study was to provide a first quantitative estimation of the growth rate of lithic communities on rock surfaces in arid regions under natural conditions. Conventional dating methods like radiocarbon cannot be applied to these rock features due to the lack of datable material. Cosmogenic radionuclide methods may be useful by yielding maximum‐limiting surface exposure estimates. However, they require large amounts of sample material, which is destructive and, therefore, unsuitable when dating archaeological sites or monuments. Lichenometry, which has been used extensively as a chronological tool in Arctic settings, is inapplicable for studying hot-desert geomorphology since deserts lack the moisture needed for lichen development (Dorn, 2009; Mcllroy de la Rosa et al., 2014). To overcome these obstacles, we used a well-preserved archaeological site built of stones that exhibit developed BRCs, marking the upper limit of the time frame needed for such processes to occur. The site was constructed from chalk and limestones, the dominant rock types in the Negev and in deserts in general (Mabbutt, 1977). Furthermore, because the nearby hillslopes were the source for the building stones, the features and composition of the BRCs can be compared to those under natural conditions. Owing to local environmental parameters, we hypothesised that the BRC origin, composition and growth rates would differ between stones composing the archaeological buildings and natural stones in the adjacent slopes. To test our hypothesis, we combined field observations and geological, archaeological and molecular microbiology characterisation to estimate the development rate and compositional trajectory of BRCs in arid regions.
To estimate the BRC development rate, we applied a quantitative analysis comparing microbial communities from chalk and rocky limestone slopes to their equivalents in a dated archaeological site, the Byzantine city of Shivta (Fig. 1). Evaluating the origin and composition of the BRCs, the study involves archaeological, geological and microbiological methods applied to collected BRC-coated rocks.
2.1 Study site
Shivta is an exceptionally well-preserved village of the Byzantine period (fourth–seventh centuries CE) that continued to be partially occupied into the early Islamic period (eighth to ninth centuries CE; Avni, 2014). Shivta is located in the south-western regions of the Negev, Israel, at an elevation of 350 m a.s.l. (30.88∘ N, 34.63∘ E; WGS 84 grid; Fig. 1), where the environmental settings have maintained arid conditions since the Holocene (10 000 years ago), and characterised by an average annual precipitation of 90 mm yr−1 and an aridity index (P PET) of 0.05, on the borderline between arid and hyper-arid ecosystems (Amit et al., 2011). The site and its vicinity are part of a rocky terrain, underlined predominantly by carbonate rock slopes. Limestone bedrock outcrops of the Turonian Nezer Formation and soft chalk outcrops of the Santonian Menuha Formation surround the city (Fig. 1). BRCs are common on the atmospherically exposed parts of many of the rocks. The site covers an area of 0.8 km2 and contains the remains of three churches, a central watchtower and associated structures, two sizable public water reservoirs, three wine presses, a large inn of the Byzantine period, and numerous private buildings. The main occupation of the residents in the Byzantine period appears to have been agriculture, mainly viticulture for the export of wine, and road services for Christian pilgrims travelling to and from Mt Sinai (in modern-day Egypt). Thanks to these activities, the village appears to have been very prosperous in the early and middle Byzantine periods (fourth–mid-sixth centuries CE) when most of the village was constructed. Excavations in the site in 2016 conducted by a team led by Guy Bar Oz and Yotam Tepper on behalf of the University of Haifa have confirmed the date of the construction of private houses in the south-eastern part of the town (Tepper et al., 2018). Their excavations extended to the bedrock under several houses over which the structures were constructed. The location of the village in an arid setting in the central Negev Highlands forced the inhabitants to construct buildings using locally quarried stone with minimum use of wood, which would otherwise have had to be imported from great distances. Therefore, the upper floors were supported with stone arches and wood was sparsely used for doorframes, doors and shelves. In lieu of wood, many installations, such as animal troughs, were carved from local stone. The walls of the houses are well-preserved, offering researchers a unique view of buildings constructed over 1700 years ago that stand to heights of two and even three stories high. The lower courses of the walls on the ground floor were made using rough limestone blocks taken from the Nezer Formation, while the upper walls were constructed of a lighter chalkstone, taken from the Menuha Formation that is more easily worked into blocks. The spaces between the heavy stones of the lower courses were sealed with mortar, and the interior walls were often covered with a base of mud plaster covered with white lime plaster (Fig. 2).
2.2 Field sampling
A total of 24 rock samples were collected from the Byzantine city of Shivta (30.88∘ N, 34.63∘ E; WGS 84 grid; the samples were named ShivtaSite 1–12) and its surroundings (30.87∘ N, 34.62∘ E; WGS 84 grid; the samples were named ShivtaSlope 1–12). Twelve rock samples were collected from the limey Nezer Formation including six samples from the archaeological site and six from the nearby natural slope. The same sampling procedure was applied for the chalky Menuha Formation. All rock samples taken from the slopes and site walls included whole rocks collected directly in the field and were not subsampled in the lab. To avoid the slope aspect effect that may lead to a different moisture regime, all samples were retrieved from south-facing slopes or walls and were collected during January 2015.
2.3 Geological analysis
The geological methods used in this study were applied to estimate the potential role of the substrate. The methods were based on direct field observations and characterisation of the subjected lithologies using thin sections, XRD analyses, total effective porosity and stable isotope analysis. Petrographic thin sections, 30 µm thick, were prepared for each lithology to test the main components in both the rock crust and the host rocks examined under a light microscope (Zeiss, Oberkochen, Germany). XRD analyses for bulk mineralogical components (Sandler et al., 2015) were conducted separately on the rock crust and the host rock. Three replicates were collected from each lithology. Powdered samples were scanned by a Panalytical X'Pert3 Powder diffractometer equipped with a PIXcel detector. The scanning range was 3–70∘ 2θ; the step size was 0.013∘; speed was 70.1 s per step. Total effective porosity (ϕ) (Scherer, 1999) was measured using a gas permeameter and porosimeter device (Core Laboratories, Amsterdam, the Netherlands) on 12 rock core cylinder samples (radius 18.5 mm and height 26.5 mm). The value of chalk total effective porosity was cited from Schütz et al. (2012) and is reported in Table 1.
a The chalk was dated to the Menuha Formation. b The limestone was dated to the Nezer Formation.
For the stable isotope δ13C and δ18O analyses, 1–2 mg of rock surface powder was obtained using a microdrill (Dremel, Racine, WI, USA) along a cross section of the rock crust and its host rock. Six profile measurements of δ13C and δ18O were performed on the chalk and limestone samples. Measurements (in duplicate) of δ18O H2O and δ13C DIC were performed on a gas source isotope ratio mass spectrometer (GS-IRMS; Thermo Fisher Scientific, Waltham, MA, USA) coupled to a GasBench II interface (Thermo Fisher Scientific) after CO2 equilibration and CO2 extraction by acidification for δ18O H2O and δ13C DIC, respectively. The samples were calibrated against internal laboratory standards: Vienna Standard Mean Ocean Water (VSMOW) and carbonate standard NBS19. δ13C and δ18O values were also referenced relative to the Vienna Pee Dee Belemnite (VPDB) standard as previously described (Uemura et al., 2016) with an SD of 0.1 ‰. All values are reported in per mil (‰).
2.4 Method validation
Method validation was conducted at the Nitzana archaeological site using the same procedure as for the Shivta site. The ancient village of Nitzana (Nessana), ca. 270 m a.s.l., is situated in the western part of the Negev (Fig. S1 in the Supplement), ca. 20 km from the Shivta site. Excavation at the site revealed that it was first established during the Hellenistic period and continued to be occupied through the Roman, Byzantine and early Islamic periods (Langgut et al., 2020, and references therein). The excavation at the Nitzana site exposed domestic quarters, two churches, a fortress and an archive of papyri from the sixth to seventh centuries CE. The papyri make Nitzana the best-documented community of all Byzantine–early Islamic sites in the Negev. The lower courses of the walls at the site were made using limestone blocks taken from the Lower Yeter Formation, while the upper walls were constructed of a lighter chalkstone, taken from the Upper Yeter Formation (Fig. S1).
2.5 DNA extraction, PCR amplification and sequencing
For DNA extraction from rocks, the surface (ca. 100 cm2) was scraped off using a wood rasp (66–67 HRC hardness; Dieter Schmid, Berlin, Germany) that was cleaned with 70 % technical-grade ethanol before each sampling. DNA was then extracted using 0.4 g of a sample using an Exgene Soil DNA extraction kit (GeneAll, Seoul, South Korea) according to the manufacturer's instructions. A 466 bp fragment of the SSU rRNA gene was amplified using the universal bacterial primers 341F (CCTAYGGGRBGCASCAG) and 806R (GGACTACNNGGGTATCTAAT) flanking the V3 and V4 region (Klindworth et al., 2013). Library construction and sequencing were performed at the DNA Services Facility, the University of Illinois at Chicago (USA), using a MiSeq sequencer (Illumina, San Diego, CA, USA) in the 2×250 cycle conﬁguration (V2 reagent kit). The raw SRA ﬁles were deposited into EMBL ENA SRA database and can be found under study accession PRJNA381355 (https://www.ebi.ac.uk/ena/browser/view/PRJNA381355, last access: 2 June 2021).
2.6 Sequence processing and analysis of bacterial communities
Paired reads generated by the MiSeq platform were quality filtered and clustered into OTUs using the UPARSE pipeline (Edgar, 2013), with modifications. Contig assembly was performed using the fastq_mergepairs command. Then, contigs were dereplicated with the derep_fulllength command, and singleton sequences were removed. OTU centroids were then determined with the cluster_otus command (set at 3 % radius). Abundances of OTUs were determined by mapping the ﬁltered contigs (before dereplication, including singletons) to the OTU centroids using the usearch_global command (set at 0.97 % identity). Following these steps, a total of ca. 1.4 million reads remained. OTU representatives were classified using mothur's implementation of a naïve Bayesian sequence classiﬁer (Schloss et al., 2009; Wang et al., 2007) against the SILVA 119 SSU NR 99 database (Quast et al., 2013). All downstream analyses were performed in R V3.4.4 (R Core Team, 2020). Data handling and manipulation were performed using the package phyloseq (McMurdie and Holmes, 2013). For alpha-diversity analysis, all samples were subsampled (rarefied) to the minimum sample size using bootstrap subsampling with 1000 iterations to account for library size differences, while for beta-diversity analysis library size normalisation was performed using GMPR (Chen et al., 2018). The inverse Simpson and Shannon's H diversity indices and the Berger–Parker dominance index were calculated using the function EstimateR in the vegan package (Oksanen et al., 2018) and tested using ANOVA in the stats package, followed by a post hoc estimated marginal means test (Searle et al., 1980) from the emmeans package (Lenth et al., 2020). Variance partitioning and testing were performed using the PERMANOVA (McArdle and Anderson, 2001) function adonis in vegan using Horn–Morisita distances (Horn, 1966). Differences in phyla composition between the sample type were tested using the non-parametric aligned rank transform ANOVA (ART ANOVA; Wobbrock et al., 2011) in the package ARTool and FDR corrected using the Benjamini–Hochberg method (Ferreira and Zwinderman, 2006) (p.adjust in function stats). Detection of differentially abundant OTUs was performed using a beta-binomial regression model (Martin et al., 2019) in the package corncob. Plots were generated using the package ggplot2 (Wickham, 2016). The scripts for reproducing the microbial analysis can be found at https://github.com/roey-angel/BRC_growth_microbiome (last access: 2 June 2021).
3.1 Geotechnical properties
BRCs were observed on all rock surfaces (i.e. limestone and chalk) both in the Byzantine site of Shivta and adjacent slopes. To evaluate the geological differences between BRCs, we performed geological characterisation of 24 limestone (Nezer Formation) and chalk (Menuha Formation) samples. Table 1 depicts the geotechnical parameters tested, including mineralogy, total effective porosity and bulk density. Several large quarries correlated with the Byzantine period (i.e. the Byzantine rule over Syria–Palestine 390–636 CE) are found within the chalk and limestone slopes of the Menuha and Nezer formations in the northern Negev (Fig. 1).
The limestone rocks constructing the Byzantine site were carved from the Turonian Nezer Formation; a well-bedded (30–50 cm), bio-micritic, fine to coarse crystalline. Being well-bedded, limestone is easily useable as a building material for different purposes (walls and ceilings). These blocks were only lightly polished since they were used mainly as the foundation of walls, leaving their original BRC intact. Consequently, most rocks were irregular in shape and sealed with mortar covered with white lime plaster (Fig. 2).
The chalk blocks were carved from the Santonian Menuha Formation: a massive soft white chalk found to be applicable as a building material thanks to the ease of polishing it. New rock surfaces were exposed in the carving process while building Shivta, leaving no remnants of the original BRCs. Soft chalk blocks were used mainly for building the upper parts of the walls (Fig. 2). The absence of plaster between the chalk blocks, unlike between the limestone blocks, indicates the intense polishing of these building blocks in order to fit them together. Therefore, the carved chalk blocks could be considered succession planes for new microbial colonisers.
3.2 Thickness measurements
Microscopic examinations of thin sections (30 µm thick) prepared from 24 limestone and chalk blocks showed that, as expected, the BRCs were restricted to the atmospherically exposed parts of the rocks (Wieler et al., 2019). The BRCs were characterised by a hardpan-laminated structure composed of masses of micritic to microsparitic carbonate layers interbedded with microbial coatings covering the lime and chalk host rocks. BRC thickness significantly differed between those BRCs found on limestone compared to those found on chalk building blocks, even when located on a single wall in the archaeological site (Table 1, Fig. 2). Limestone BRC thickness at both the Byzantine site and along the rock outcrops located in the limestone quarries at the adjacent slopes ranged between 1.4–4.3 mm (3.2±0.26 and 2.6±0.38; Fig. 2d, e). In contrast, chalk BRC thickness differed between the blocks in the Byzantine site and the natural chalk slopes. The chalk BRC thickness at the Byzantine site ranged between 0.1–0.6 mm (0.29±0.07; Fig. 2b), while at the adjacent slopes, it ranged between 0.5–1.8 mm (1.28±0.18; Fig. 2c). Dating of the Shivta site is mainly associated with the early–middle Byzantine period (fourth–sixth centuries CE) (1700–1500 years ago) as was recently confirmed by Guy Bar Oz and Yotam Tepper on behalf of the University of Haifa (Tepper et al., 2018). Dividing the observed thickness of the BRC from the chalk building blocks by the age of the Byzantine Shivta site (1700–1500 years), located under long-term stable arid conditions and free of anthropogenic effects, suggests a growth rate of 0.06–0.35 mm kyr−1 (0.17±0.04).
3.3 Method validation
We validated our BRC dating method by conducting a parallel thickness measurement for BRCs sampled from a different archaeological site at Nizzana (see “Materials and methods”). We used an identical procedure to that used in the Shivta site. All samples were retrieved from south-facing walls and were processed similarly. We measured rock crust thickness of 0.1–0.2 mm (1.67±0.01) on the carved limestone and chalk building blocks. Dividing the thickness by the age of the Nitzana site, we found the BRC growth rates to be consistent with the one calculated at the Shivta site (0.1 mm kyr−1).
3.4 Isotopic composition
The biogenic nature of the crusts was confirmed using a cross-section analysis of the stable carbon and oxygen isotope ratios in the crust and host rock. For the limestone sample, a shift was found between the δ13C values for the BRC (0–2 mm) and the host rock (2–5 mm) layers, with BRC values ranging between −4 ‰ and −5 ‰ VPDB and host-rock values between 0 ‰ and 1 ‰ VPDB (Fig. 3a). However, for the chalk, δ13C values for the BRC ranged between −0.2 ‰ and −1.9 ‰ VPDB, and in the host rock, they ranged between 0.1 ‰ and −2.7 ‰ VPDB. These results were consistent for both the slope and archaeological site samples. The limestone δ13C values are typical indicators of carbon isotope exchange of primary marine CaCO3 (abundant in the bedrock) with CO2 released by microbial respiration (i.e. of carbon originating from photosynthesis) with the subsequent precipitation of pedogenic calcrete (Brlek and Glumac, 2014). The differences in δ13C values between the chalk and limestone are suggested to result from the BRC thickness. Thicker BRCs, as observed in the limestone samples, should produce more biogenic activity than for the chalk. The limestone δ18O values ranged between −3 ‰ in the BRC and −7.3 ‰ VPDB in the host rock. This decrease in δ18O values in the host rock is suggested to result from meteoric water substitution of the marine limestone (Alonso-Zarza and Tanner, 2006). The chalk δ18O values ranged between −3.3 ‰ and −5.4 ‰ VPDB for both the BRC and the host rock; these negative δ18O values were consistent with values of other Santonian chalk (Clarke and Jenkyns, 1999; Liu, 2009).
3.5 Bacterial diversity of the BRCs in the Byzantine Shivta site and the nearby slopes
To elucidate the identity of the bacterial communities on the BRCs, we performed multiplexed barcoded amplicon sequencing of the small subunit ribosomal RNA gene (SSU rRNA gene). As expected, we found simple BRC communities (i.e. low richness and low diversity) in all BRC samples. Comparing the BRCs collected from the Byzantine site to the nearby slopes from the two primary lithologies (the limestone and chalk) showed no statistically significant differences in the observed number of OTUs or their inverse Simpson or Shannon indices, with all samples averaging 312±11 OTUs (Fig. 4a, Table S1). However, the city samples showed a slightly higher dominance of the most abundant OTU compared to the slope samples (Berger–Parker) values of 0.24±0.03 vs. 0.14±0.01; p=0.0017). The bacterial community composition of the BRCs was very similar in both lithologies and sample sources and was heavily dominated by members of the phylum Actinobacteria, followed by Cyanobacteria, Proteobacteria (mainly Alphaproteobacteria from the orders Sphingomonadales and Rhodospirillales) and Bacteroidetes (Fig. 4c). A variance partitioning analysis did, however, show a small difference between the samples from the city and the ones from the slopes, explaining 8.1 % of the variance in the Morisita–Horn distance matrix (p=0.023). In contrast, neither the lithology nor the interaction between the lithology and sample source correlated with differences in bacterial community composition (Table S2). The differences in the communities between the samples from the city and the ones from the slopes were not detected when comparing the dominance of each phylum between locations using ART ANOVA testing (Fig. 4d). However, using a differential abundance test, 64 individual OTUs were differentially more abundant in the city samples while 57 were differentially more abundant in the slope samples (from a total of 694 OTUs; Fig. S2). These OTUs came from all the dominant phyla and belonged to orders that are commonly found in hot deserts such as Rubrobacterales and Solirubrobacterales (Actinobacteria), Cytophagales (Bacteroidetes), Thermomicrobiales (Chloroflexi), Chroococcidiopsidales (Cyanobacteria), Rhodospirillales (Alphaproteobacteria), and Burkholderiales (Betaproteobacteria) (Table S3; Angel and Conrad, 2013; Makhalanyane et al., 2015). However, there does not seem to be a discernible taxonomic pattern for the differentially abundant OTUs.
BRCs are a common and important feature of atmospherically exposed rock surfaces in drylands around the globe, yet reliable growth estimations based on field data are rare. Only a handful of studies have tried to estimate BRC growth rates, among them Lange (1990) measured radial growth rates of epilithic lichens in the arid Negev, Israel, and reported an average growth rate of 0.371 mm yr−1. Another study conducted by Krumbein and Jens (1981) found that desert varnish, a typical rock crust in desert regions, shows black fungi microcolonies after 8 weeks of isolation in the lab. Liu and Broecker (2000) noted the accumulation rates of desert varnish in natural conditions to be much slower, in the range of <0.001 to 0.04 mm kyr−1, and to rarely reach a thickness exceeding 200 µm.
In this study, we observed a 0.1 to 0.6 mm thick biological rock crust coverage on chalk building blocks of the Byzantine site (Fig. 2b) dated to the fourth to seventh centuries CE (Tepper et al., 2018). The building blocks experienced long-term stable arid conditions and were free from direct anthropogenic influences. Thanks to the processing method of the bricks used for construction, leaving BRC-free surfaces exposed to the elements, the BRC is confined by the dating of the archaeological site (Tepper et al., 2018). This provides us with an estimated growth rate of 0.06–0.35 mm kyr−1 for BRCs of primarily bacterial origin. Besides, the presence of a 1.4–4.3 mm thick BRC in all limestone rocks (Fig. 2d, e) and a 0.5–1.8 mm thick BRC on the chalk rocks collected from the slopes (Fig. 2c) provide a plausible maximum for BRC growth under these conditions. Calibrating local climate conditions and mineral composition is needed for BRC growth rates when applying the methods of this paper to different conditions. The presence of similarly thick crusts on the rocks of a different site which is dated to the same period but is located 20 km away and is constructed of slightly different rocks, stemming from a different formation, provides validation for our methods.
Applying stable isotopes, we noted a consistent trend where more negative δ13C values were found in the limestone crust compared to in their host rock (Fig. 3a), reflecting a biogenic agency in the rock crust production, regardless of the substrate type. The oxygen isotope profile of the chalk samples collected from the archaeological site showed no difference between the bedrock and the crust (Fig. 3b) Thus, it indicates the early stages of BRC formation. The limestone samples showed more negative δ13C values in the BRC compared to the host rock. Yet, the δ18O values were more negative at the host rock compared to the BRC, which may indicate sub-aerial exposure of the BRC and substitution of the marine limestone.
From a microbiological perspective, all the samples studied here showed very similar microbial community composition, irrespective of lithology or BRC thickness. This similarity in composition demonstrates the indifference of microorganisms to the type of attachment surface in this case and that the community probably changes very little after establishing. However, the differences in BRC thickness between the chalk and limestone sampled from the slopes could indicate that the latter can better support BRC growth.
Differences in the dominance of the most abundant OTU as well as some differences in the relative abundance of about 17 % of the OTUs were observed between BRCs from the city and slope samples. The OTUs belonged to taxa that are commonly found in hyper-arid deserts. However, the differences in abundance were relatively minor. Moreover, they lacked a clear taxonomic pattern, with different OTUs belonging to the same taxa ending up as differentially either more or less abundant. This could indicate a strong influence of stochastic processes (founder effect) over environmental filtering at the OTU level. Therefore, our findings indicate that the rock surfaces at the archaeological site act as succession planes for BRC development. Moreover, similarities in the microbial composition between rocks that were in contact with the ground to those that were detached from the ground (chalk BRCs found in Shivta) indicate a major role for aeolian processes (e.g. settled dust particles) in determining the community composition of BRCs in deserts as has been previously reported (Wieler et al., 2019).
Assisted recovery of lithobiontic communities has not been conducted in natural settings, and most research has focused on the regeneration of soil biocrusts (Velasco Ayuso et al., 2017). Artificial cultivation of soil biocrusts (Zhang et al., 2018) has suggested that inoculation of cyanobacteria and algal communities enhance the recovery of biocrusts. Testing the soiling impact on rock surfaces has indicated that bacterial colonisation plays an important part in the development of fungal biofilms (Viles and Gorbushina, 2003). In fact, cyanobacteria and fungi are believed to be critical components of soil and rock biocrusts (Gorbushina, 2007; Weber et al., 2016). In arid BRCs, fungi are scarce (Lang-Yona et al., 2018), but our results suggest that not only Cyanobacteria but also other taxa like Alphaproteobacteria and Actinobacteria play a key role in BRC development (Fig. 4c). Similarly, the dominance of Alphaproteobacteria and Actinobacteria in early soil biocrusts formation has been observed in arid regions (Ji et al., 2017; Pepe-Ranney et al., 2016).
Unlike more humid dryland soils, in hyper-arid soils and rock biocrusts, Actinobacteria remain the dominant phylum even in later stages of development (Holmes et al., 2000; Angel and Conrad, 2013; Idris et al., 2017). Kuske et al. (2012) identified a deep branching subclass within the Actinobacteria phylum abundant worldwide, the Rubrobacteridae, a taxon extremely resistant to desiccation and UV stress (Holmes et al., 2000; Rainey et al., 2005). The presence of Rubrobacteridae suggests that they may be involved in shaping the rock surface structure during biocrust formation (Mummey et al., 2006; Rainey et al., 2005).
Our analyses suggest that both rocks experience the same regional-scale environmental factors (Fig. 4b) in shaping the BRC composition. This demonstrates the ecological filtering effect of the rock surfaces, which imposes unique abiotic challenges for the inhabiting microbes and implies that dust particles are the main potential source of colonisation for the microbial communities (Wieler et al., 2019). We thus suggest that local-scale environmental parameters play a major role in shaping the microbial taxa that colonise fresh rock surfaces in arid regions. This also suggests that the BRCs cannot be regarded as passive deposits of microbial cells but should rather be seen as a filter selecting for a specific subset of adapted microbes that can persist and form a biofilm under these harsh conditions.
This study provides the first estimate to an unfathomed question on the rate of BRC growth rate under a natural setting. The growth rate observed here, 0.06–0.35 mm kyr−1, validates the extremely slow nature of such a succession process. Using a well-dated archaeological site in the arid Negev, Israel, we demonstrate the possibility of using such human-made artefacts to document and confine long microbial developmental processes that are otherwise too slow to monitor. We confirm the biological origin of the crusts through analysis of stable carbon and oxygen isotope ratios. We further show that local-scale environmental parameters are found to shape the microbial taxa that colonise the fresh rock surfaces, including Actinobacteria, Cyanobacteria and Proteobacteria. This also suggests that the BRCs cannot be regarded as passive deposits of microbial cells but should rather be seen as a filter selecting for a specific subset of adapted microbes that can persist and form a biofilm under these harsh conditions. Conversely, once a growth rate has been established for a region, it could be used, by itself, to date the age of atmospherically exposed archaeological artefacts.
The code is available at https://doi.org/10.5281/zenodo.4747587 (Angel, 2021a) and https://github.com/roey-angel/BRC_growth_microbiome (last access: 2 June 2021).
The data are available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA381355 (Angel, 2021b).
The supplement related to this article is available online at: https://doi.org/10.5194/bg-18-3331-2021-supplement.
NW, TEG, OG and RA conceptualised the study; NW performed the fieldwork and lab work; NW and RA analysed the data; NW, TEG, OG and RA wrote the paper. OG and RA contributed equally to this study.
The authors declare that they have no conflict of interest.
Roey Angel was supported by BC CAS, ISB and SoWa RI (MEYS; projects LM2015075 and EF16_013/0001782 – SoWa Ecosystems Research).
This paper was edited by Denise Akob and reviewed by two anonymous referees.
Alonso-Zarza, A. M. and Tanner, L. H. (Eds.): Paleoenvironmental record and applications of calcretes and palustrine carbonates (Vol. 416), Geological Society of America, 2006. a
Alonso-Zarza, A. M. and Wright, V. P.: Chapter 5 Calcretes, in: Developments in Sedimentology, edited by: Alonso-Zarza, A. M. and Tanner, L. H., vol. 61 of Carbonates in Continental Settings: Facies, Environments, and Processes, 225–267, Elsevier, Amsterdam, the Netherlands, Oxford, UK, https://doi.org/10.1016/S0070-4571(09)06105-6, 2010. a
Amit, R., Simhai, O., Ayalon, A., Enzel, Y., Matmon, A., Crouvi, O., Porat, N., and McDonald, E.: Transition from Arid to Hyper-Arid Environment in the Southern Levant Deserts as Recorded by Early Pleistocene Cummulic Aridisols, Quaternary Sci. Rev., 30, 312–323, https://doi.org/10.1016/j.quascirev.2010.11.007, 2011. a
Angel, R.: Rock metagenome (ID 381355) – BioProject – NCBI, available at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA381355, last access: 2 June 2021b. a
Angel, R. and Conrad, R.: Elucidating the Microbial Resuscitation Cascade in Biological Soil Crusts Following a Simulated Rain Event, Environ. Microbiol., 15, 2799–2815, https://doi.org/10.1111/1462-2920.12140, 2013. a, b
Avni, G.: The Byzantine-Islamic Transition in Palestine: An Archaeological Approach, Oxford University press, Oxford, UK, 2014. a
Brlek, M. and Glumac, B.: Stable Isotopic (δ13C and δ18O) Signatures of Biogenic Calcretes Marking Discontinuity Surfaces: A Case Study from Upper Cretaceous Carbonates of Central Dalmatia and Eastern Istria, Croatia, Facies, 60, 773–788, https://doi.org/10.1007/s10347-014-0403-7, 2014. a
Büdel, B. and Wessels, D. C. J.: Rock Inhabiting Blue-Green Algae/Cyanobacteria from Hot Arid Regions, Algological Studies/Archiv für Hydrobiologie, Supplement Volumes, 385–398, 1991. a
Cabello-Briones, C. and Viles, H. A.: Evaluating the Effects of Open Shelters on Limestone Deterioration at Archaeological Sites in Different Climatic Locations, Int. J. Archit. Herit., 11, 816–828, https://doi.org/10.1080/15583058.2017.1300710, 2017. a
Chen, L., Reeve, J., Zhang, L., Huang, S., Wang, X., and Chen, J.: GMPR: A Robust Normalization Method for Zero-Inflated Count Data with Application to Microbiome Sequencing Data, PeerJ, 6, e4600, https://doi.org/10.7717/peerj.4600, 2018. a
Clarke, L. J. and Jenkyns, H. C.: New Oxygen Isotope Evidence for Long-Term Cretaceous Climatic Change in the Southern Hemisphere, Geology, 27, 699–702, https://doi.org/10.1130/0091-7613(1999)027<0699:NOIEFL>2.3.CO;2, 1999. a
Corenblit, D., Darrozes, J., Julien, F., Otto, T., Roussel, E., Steiger, J., and Viles, H.: The Search for a Signature of Life on Mars: A Biogeomorphological Approach, Astrobiology, 19, 1279–1291, https://doi.org/10.1089/ast.2018.1969, 2019. a
Crouvi, O., Amit, R., Enzel, Y., Porat, N., and Sandler, A.: Sand Dunes as a Major Proximal Dust Source for Late Pleistocene Loess in the Negev Desert, Israel, Quaternary Res., 70, 275–282, https://doi.org/10.1016/j.yqres.2008.04.011, 2008. a
Dorn, R. I.: Desert Rock Coatings, in: Geomorphology of Desert Environments, edited by: Parsons, A. J. and Abrahams, A. D., 153–186, Springer Netherlands, Dordrecht, https://doi.org/10.1007/978-1-4020-5719-9_7, 2009. a
Holmes, A. J., Bowyer, J., Holley, M. P., O'Donoghue, M., Montgomery, M., and Gillings, M. R.: Diverse, yet-to-Be-Cultured Members of the Rubrobacter Subdivision of the Actinobacteria Are Widespread in Australian Arid Soils, FEMS Microbiol. Ecol., 33, 111–120, https://doi.org/10.1111/j.1574-6941.2000.tb00733.x, 2000. a, b
Idris, H., Goodfellow, M., Sanderson, R., Asenjo, J. A., and Bull, A. T.: Actinobacterial Rare Biospheres and Dark Matter Revealed in Habitats of the Chilean Atacama Desert, Sci. Rep.-UK, 7, 8373, https://doi.org/10.1038/s41598-017-08937-4, 2017. a
Ji, M., Greening, C., Vanwonterghem, I., Carere, C. R., Bay, S. K., Steen, J. A., Montgomery, K., Lines, T., Beardall, J., van Dorst, J., Snape, I., Stott, M. B., Hugenholtz, P., and Ferrari, B. C.: Atmospheric Trace Gases Support Primary Production in Antarctic Desert Surface Soil, Nature, 552, 400–403, https://doi.org/10.1038/nature25014, 2017. a
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., and Glöckner, F. O.: Evaluation of General 16S Ribosomal RNA Gene PCR Primers for Classical and Next-Generation Sequencing-Based Diversity Studies, Nucleic Acids Res., 41, e1–e1, https://doi.org/10.1093/nar/gks808, 2013. a
Krumbein, W. E. and Jens, K.: Biogenic Rock Varnishes of the Negev Desert (Israel) an Ecological Study of Iron and Manganese Transformation by Cyanobacteria and Fungi, Oecologia, 50, 25–38, https://doi.org/10.1007/BF00378791, 1981. a
Kuske, C. R., Yeager, C. M., Johnson, S., Ticknor, L. O., and Belnap, J.: Response and Resilience of Soil Biocrust Bacterial Communities to Chronic Physical Disturbance in Arid Shrublands, ISME J., 6, 886–897, https://doi.org/10.1038/ismej.2011.153, 2012. a
Lang-Yona, N., Maier, S., Macholdt, D. S., Müller-Germann, I., Yordanova, P., Rodriguez-Caballero, E., Jochum, K. P., Al-Amri, A., Andreae, M. O., Fröhlich-Nowoisky, J., and Weber, B.: Insights into Microbial Involvement in Desert Varnish Formation Retrieved from Metagenomic Analysis, Env. Microbiol. Rep., 10, 264–271, https://doi.org/10.1111/1758-2229.12634, 2018. a, b
Lange, O. L.: Twenty-Three Years of Growth Measurements on the Crustose Lichen Caloplaca Aurantia in the Central Negev Desert, Israel J. Bot., 39, 383–394, https://doi.org/10.1080/0021213X.1990.10677162, 1990. a
Langgut, D., Tepper, Y., Benzaquen, M., Erickson-Gini, T., and Bar-Oz, G.: Environment and Horticulture in the Byzantine Negev Desert, Israel: Sustainability, Prosperity and Enigmatic Decline, Quaternary Int., https://doi.org/10.1016/j.quaint.2020.08.056, online first, 2020. a
Lenth, R., Buerkner, P., Herve, M., Love, J., Riebl, H., and Singmann, H.: Emmeans: Estimated Marginal Means, Aka Least-Squares Means, available at: https://cran.r-project.org/web/packages/emmeans/index.html (last access: 2 June 2021), 2020. a
Liu, K.: Oxygen and Carbon Isotope Analysis of the Mooreville Chalk and Late Santonian-Early Campanian Sea Level and Sea Surface Temperature Changes, Northeastern Gulf of Mexico, U.S.A., Cretaceous Res., 30, 980–990, https://doi.org/10.1016/j.cretres.2009.02.008, 2009. a
Liu, T. and Broecker, W. S.: How Fast Does Rock Varnish Grow?, Geology, 28, 183–186, https://doi.org/10.1130/0091-7613(2000)28<183:HFDRVG>2.0.CO;2, 2000. a
Mabbutt, J. A.: Desert Landforms, Australian National University Press, Canberra, Australia, 1977. a
Makhalanyane, T. P., Valverde, A., Gunnigle, E., Frossard, A., Ramond, J.-B., and Cowan, D. A.: Microbial Ecology of Hot Desert Edaphic Systems, FEMS Microbiol. Rev., 39, 203–221, https://doi.org/10.1093/femsre/fuu011, 2015. a, b
McArdle, B. H. and Anderson, M. J.: Fitting Multivariate Models to Community Data: A Comment on Distance-Based Redundancy Analysis, Ecology, 82, 290–297, https://doi.org/10.1890/0012-9658(2001)082[0290:FMMTCD]2.0.CO;2, 2001. a
Mcllroy de la Rosa, J. P., Warke, P. A., and Smith, B. J.: The Effects of Lichen Cover upon the Rate of Solutional Weathering of Limestone, Geomorphology, 220, 81–92, https://doi.org/10.1016/j.geomorph.2014.05.030, 2014. a, b
McMurdie, P. J. and Holmes, S.: Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data, PLOS ONE, 8, e61217, https://doi.org/10.1371/journal.pone.0061217, 2013. a
Mummey, D., Holben, W., Six, J., and Stahl, P.: Spatial Stratification of Soil Bacterial Populations in Aggregates of Diverse Soils, Microb. Ecol., 51, 404–411, https://doi.org/10.1007/s00248-006-9020-5, 2006. a
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., Minchin, P. R., O'Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. H., Szoecs, E., and Wagner, H.: Vegan: Community Ecology Package, available at: https://cran.r-project.org/web/packages/vegan/index.html (last access: 2 June 2021), 2018. a
Owen, J. J., Amundson, R., Dietrich, W. E., Nishiizumi, K., Sutter, B., and Chong, G.: The Sensitivity of Hillslope Bedrock Erosion to Precipitation, Earth Surf. Proc. Land., 36, 117–135, https://doi.org/10.1002/esp.2083, 2011. a
Pepe-Ranney, C., Koechli, C., Potrafka, R., Andam, C., Eggleston, E., Garcia-Pichel, F., and Buckley, D. H.: Non-Cyanobacterial Diazotrophs Mediate Dinitrogen Fixation in Biological Soil Crusts during Early Crust Formation, ISME J., 10, 287–298, https://doi.org/10.1038/ismej.2015.106, 2016. a
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., Peplies, J., and Glockner, F. O.: The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools, Nucleic Acids Res., 41, D590–D596, https://doi.org/10.1093/nar/gks1219, 2013. a
Rainey, F. A., Ray, K., Ferreira, M., Gatz, B. Z., Nobre, M. F., Bagaley, D., Rash, B. A., Park, M.-J., Earl, A. M., Shank, N. C., Small, A. M., Henk, M. C., Battista, J. R., Kämpfer, P., and da Costa, M. S.: Extensive Diversity of Ionizing-Radiation-Resistant Bacteria Recovered from Sonoran Desert Soil and Description of Nine New Species of the Genus Deinococcus Obtained from a Single Soil Sample, Appl. Environ. Microb., 71, 5225–5235, https://doi.org/10.1128/AEM.71.9.5225-5235.2005, 2005. a, b
Sandler, A., Meunier, A., and Velde, B.: Mineralogical and Chemical Variability of Mountain Red/Brown Mediterranean Soils, Geoderma, 239–240, 156–167, https://doi.org/10.1016/j.geoderma.2014.10.008, 2015. a
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., Lesniewski, R. A., Oakley, B. B., Parks, D. H., Robinson, C. J., Sahl, J. W., Stres, B., Thallinger, G. G., Van Horn, D. J., and Weber, C. F.: Introducing Mothur: Open Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities, Appl. Environ. Microb., 75, 7537–7541, https://doi.org/10.1128/AEM.01541-09, 2009. a
Schütz, F., Norden, B., and Förster, A.: Thermal Properties of Sediments in Southern Israel: A Comprehensive Data Set for Heat Flow and Geothermal Energy Studies, Basin Res., 24, 357–376, https://doi.org/10.1111/j.1365-2117.2011.00529.x, 2012. a
Searle, S. R., Speed, F. M., and Milliken, G. A.: Population Marginal Means in the Linear Model: An Alternative to Least Squares Means, Am. Stat., 34, 216–221, https://doi.org/10.1080/00031305.1980.10483031, 1980. a
Sneh, A., Avni, Y., Bartov, Y., Zilberman, E., Braun, M., Lasman, N., and Weinberger, R.: The geological map of Israel, 1 : 50,000. Sheet 18-III: Shivta, available at: https://www.gov.il/he/departments/general/shivta-map (last access: 2 June 2021), 2011. a
Tepper, Y., Erickson-Gini, T., Farhi, Y., and Bar-Oz, G.: Probing the Byzantine/Early Islamic Transition in the Negev: The Renewed Shivta Excavations, 2015–2016, Tel Aviv, 45, 120–152, https://doi.org/10.1080/03344355.2018.1412058, 2018. a, b, c, d
Uemura, R., Nakamoto, M., Asami, R., Mishima, S., Gibo, M., Masaka, K., Jin-Ping, C., Wu, C.-C., Chang, Y.-W., and Shen, C.-C.: Precise Oxygen and Hydrogen Isotope Determination in Nanoliter Quantities of Speleothem Inclusion Water by Cavity Ring-down Spectroscopic Techniques, Geochim. Cosmochim. Ac., 172, 159–176, https://doi.org/10.1016/j.gca.2015.09.017, 2016. a
Velasco Ayuso, S., Giraldo Silva, A., Nelson, C., Barger, N. N., and Garcia-Pichel, F.: Microbial Nursery Production of High-Quality Biological Soil Crust Biomass for Restoration of Degraded Dryland Soils, Appl. Environ. Microb., 83, e02179-16, https://doi.org/10.1128/AEM.02179-16, 2017. a
Viles, H. A. and Gorbushina, A. A.: Soiling and Microbial Colonisation on Urban Roadside Limestone: A Three Year Study in Oxford, England, Build. Environ., 38, 1217–1224, https://doi.org/10.1016/S0360-1323(03)00078-7, 2003. a, b
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R.: Naive Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy, Appl. Environ. Microb., 73, 5261–5267, https://doi.org/10.1128/AEM.00062-07, 2007. a
Weber, B., Büdel, B., and Belnap, J. (Eds.): Biological Soil Crusts: An Organizing Principle in Drylands, Springer, New York, NY, 1st Edn., 2016. a
Wickham, H.: Ggplot2: Elegant Graphics for Data Analysis, Springer-Verlag, New York, USA, 2016. a
Wieler, N., Ginat, H., Gillor, O., and Angel, R.: The origin and role of biological rock crusts in rocky desert weathering, Biogeosciences, 16, 1133–1145, https://doi.org/10.5194/bg-16-1133-2019, 2019. a, b, c, d
Wierzchos, J., Davila, A. F., Artieda, O., Cámara-Gallego, B., de los Ríos, A., Nealson, K. H., Valea, S., Teresa García-González, M., and Ascaso, C.: Ignimbrite as a Substrate for Endolithic Life in the Hyper-Arid Atacama Desert: Implications for the Search for Life on Mars, Icarus, 224, 334–346, https://doi.org/10.1016/j.icarus.2012.06.009, 2013. a
Wobbrock, J. O., Findlater, L., Gergle, D., and Higgins, J. J.: The Aligned Rank Transform for Nonparametric Factorial Analyses Using Only Anova Procedures, in: Proceedings of the SIGCHI Conference on Human Factors in Computing Systems, CHI '11, 143–146, Association for Computing Machinery, New York, NY, USA, https://doi.org/10.1145/1978942.1978963, 2011. a
Zhang, Y., Duan, P., Zhang, P., and Li, M.: Variations in Cyanobacterial and Algal Communities and Soil Characteristics under Biocrust Development under Similar Environmental Conditions, Plant Soil, 429, 241–251, https://doi.org/10.1007/s11104-017-3443-2, 2018. a