ASU Scholarship Showcase

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

A relatively unexplored issue in cybersecurity science and engineering is whether there exist intrinsic patterns of cyberattacks. Conventional wisdom favors absence of such patterns due to the overwhelming complexity of the modern cyberspace. Surprisingly, through a detailed analysis of an extensive data set that records the time-dependent frequencies of attacks over a relatively wide range of consecutive IP addresses, we successfully uncover intrinsic spatiotemporal patterns underlying cyberattacks, where the term “spatio” refers to the IP address space. In particular, we focus on analyzing macroscopic properties of the attack traffic flows and identify two main patterns with distinct spatiotemporal characteristics: deterministic and stochastic. Strikingly, there are very few sets of major attackers committing almost all the attacks, since their attack “fingerprints” and target selection scheme can be unequivocally identified according to the very limited number of unique spatiotemporal characteristics, each of which only exists on a consecutive IP region and differs significantly from the others. We utilize a number of quantitative measures, including the flux-fluctuation law, the Markov state transition probability matrix, and predictability measures, to characterize the attack patterns in a comprehensive manner. A general finding is that the attack patterns possess high degrees of predictability, potentially paving the way to anticipating and, consequently, mitigating or even preventing large-scale cyberattacks using macroscopic approaches.

Supply-demand processes take place on a large variety of real-world networked systems ranging from power grids and the internet to social networking and urban systems. In a modern infrastructure, supply-demand systems are constantly expanding, leading to constant increase in load requirement for resources and consequently, to problems such as low efficiency, resource scarcity, and partial system failures. Under certain conditions global catastrophe on the scale of the whole system can occur through the dynamical process of cascading failures. We investigate optimization and resilience of time-varying supply-demand systems by constructing network models of such systems, where resources are transported from the supplier sites to users through various links. Here by optimization we mean minimization of the maximum load on links, and system resilience can be characterized using the cascading failure size of users who fail to connect with suppliers.

We consider two representative classes of supply schemes: load driven supply and fix fraction supply. Our findings are: (1) optimized systems are more robust since relatively smaller cascading failures occur when triggered by external perturbation to the links; (2) a large fraction of links can be free of load if resources are directed to transport through the shortest paths; (3) redundant links in the performance of the system can help to reroute the traffic but may undesirably transmit and enlarge the failure size of the system; (4) the patterns of cascading failures depend strongly upon the capacity of links; (5) the specific location of the trigger determines the specific route of cascading failure, but has little effect on the final cascading size; (6) system expansion typically reduces the efficiency; and (7) when the locations of the suppliers are optimized over a long expanding period, fewer suppliers are required. These results hold for heterogeneous networks in general, providing insights into designing optimal and resilient complex supply-demand systems that expand constantly in time.

The heterocyclic indole-alkaloid scytonemin is a sunscreen found exclusively among cyanobacteria. An 18-gene cluster is responsible for scytonemin production in Nostoc punctiforme ATCC 29133. The upstream genes scyABCDEF in the cluster are proposed to be responsible for scytonemin biosynthesis from aromatic amino acid substrates. In vitro studies of ScyA, ScyB, and ScyC proved that these enzymes indeed catalyze initial pathway reactions. Here we characterize the role of ScyD, ScyE, and ScyF, which were logically predicted to be responsible for late biosynthetic steps, in the biological context of N. punctiforme. In-frame deletion mutants of each were constructed (ΔscyD, ΔscyE, and ΔscyF) and their phenotypes studied. Expectedly, ΔscyE presents a scytoneminless phenotype, but no accumulation of the predicted intermediaries. Surprisingly, ΔscyD retains scytonemin production, implying that it is not required for biosynthesis. Indeed, scyD presents an interesting evolutionary paradox: it likely originated in a duplication event from scyE, and unlike other genes in the operon, it has not been subjected to purifying selection. This would suggest that it is a pseudogene, and yet scyD is highly conserved in the scytonemin operon of cyanobacteria. ΔscyF also retains scytonemin production, albeit exhibiting a reduction of the production yield compared with the wild-type. This indicates that ScyF is not essential but may play an adjuvant role for scytonemin synthesis. Altogether, our findings suggest that these downstream genes are not responsible, as expected, for the late steps of scytonemin synthesis and we must look for those functions elsewhere. These findings are particularly important for biotechnological production of this sunscreen through heterologous expression of its genes in more tractable organisms.

Most previous works on complete synchronization of chaotic oscillators focused on the one-channel interaction scheme where the oscillators are coupled through only one variable or a symmetric set of variables. Using the standard framework of master-stability function (MSF), we investigate the emergence of complex synchronization behaviors under all possible configurations of two-channel coupling, which include, for example, all possible cross coupling schemes among the dynamical variables. Utilizing the classic Rössler and Lorenz oscillators, we find a rich variety of synchronization phenomena not present in any previously extensively studied, single-channel coupling configurations. For example, in many cases two coupling channels can enhance or even generate synchronization where there is only weak or no synchronization under only one coupling channel, which has been verified in a coupled neuron system. There are also cases where the oscillators are originally synchronized under one coupling channel, but an additional synchronizable coupling channel can, however, destroy synchronization. Direct numerical simulations of actual synchronization dynamics verify the MSF-based predictions. Our extensive computation and heuristic analysis provide an atlas for synchronization of chaotic oscillators coupled through two channels, which can be used as a systematic reference to facilitate further research in this area.

N₂ fixation and ammonia oxidation (AO) are the two most important processes in the nitrogen (N) cycle of biological soil crusts (BSCs). We studied the short-term response of acetylene reduction assay (ARA) rates, an indicator of potential N₂ fixation, and AO rates to temperature (T, -5°C to 35°C) in BSC of different successional stages along the BSC ecological succession and geographic origin (hot Chihuahuan and cooler Great Basin deserts). ARA in all BSCs increased with T until saturation occurred between 15 and 20°C, and declined at 30–35°C. Culture studies using cyanobacteria isolated from these crusts indicated that the saturating effect was traceable to their inability to grow well diazotrophically within the high temperature range. Below saturation, temperature response was exponential, with Q₁₀ significantly different in the two areas (~ 5 for Great Basin BSCs; 2–3 for Chihuahuan BSCs), but similar between the two successional stages. However, in contrast to ARA, AO showed a steady increase to 30–35°C in Great Basin, and Chihuhuan BSCs showed no inhibition at any tested temperature. The T response of AO also differed significantly between Great Basin (Q₁₀ of 4.5–4.8) and Chihuahuan (Q₁₀ of 2.4–2.6) BSCs, but not between successional stages. Response of ARA rates to T did not differ from that of AO in either desert. Thus, while both processes scaled to T in unison until 20°C, they separated to an increasing degree at higher temperature. As future warming is likely to occur in the regions where BSCs are often the dominant living cover, this predicted decoupling is expected to result in higher proportion of nitrates in soil relative to ammonium. As nitrate is more easily lost as leachate or to be reduced to gaseous forms, this could mean a depletion of soil N over large landscapes globally.

The Quadrangles Av-11 and Av-12 on Vesta are located at the northern rim of the giant Rheasilvia south polar impact basin. The primary geologic units in Av-11 and Av-12 include material from the Rheasilvia impact basin formation, smooth material and different types of impact crater structures (such as bimodal craters, dark and bright crater ray material and dark ejecta material). Av-11 and Av-12 exhibit almost the full range of mass wasting features observed on Vesta, such as slump blocks, spur-and-gully morphologies and landslides within craters. Processes of collapse, slope instability and seismically triggered events force material to slump down crater walls or scarps and produce landslides or rotational slump blocks. The spur-and-gully morphology that is known to form on Mars is also observed on Vesta; however, on Vesta this morphology formed under dry conditions.

A variety of geologic landforms and features are observed within quadrangle Av-13 Tuccia in the southern hemisphere of Vesta. The quadrangle covers parts of the highland Vestalia Terra as well as the floors of the large Rheasilvia and Veneneia impact basins, which results in a substantial elevation difference of more than 40 km between the northern and the southern portions of the quadrangle. Measurements of crater size–frequency distributions within and surrounding the Rheasilvia basin indicate that gravity-driven mass wasting in the interior of the basin has been important, and that the basin has a more ancient formation age than would be expected from the crater density on the basin floor alone. Subsequent to its formation, Rheasilvia was superimposed by several mid-sized impact craters. The most prominent craters are Tuccia, Eusebia, Vibidia, Galeria, and Antonia, whose geology and formation ages are investigated in detail in this work. These impact structures provide a variety of morphologies indicating different sorts of subsequent impact-related or gravity-driven mass wasting processes. Understanding the geologic history of the relatively young craters in the Rheasilvia basin is important in order to understand the even more degraded craters in other regions of Vesta.

In this paper we present a time-stratigraphic scheme and geologic time scale for the protoplanet Vesta, based on global geologic mapping and other analyses of NASA Dawn spacecraft data, complemented by insights gained from laboratory studies of howardite–eucrite–diogenite (HED) meteorites and geophysical modeling. On the basis of prominent impact structures and their associated deposits, we propose a time scale for Vesta that consists of four geologic time periods: Pre-Veneneian, Veneneian, Rheasilvian, and Marcian. The Pre-Veneneian Period covers the time from the formation of Vesta up to the Veneneia impact event, from 4.6 Ga to >2.1 Ga (using the asteroid flux-derived chronology system) or from 4.6 Ga to 3.7 Ga (under the lunar-derived chronology system). The Veneneian Period covers the time span between the Veneneia and Rheasilvia impact events, from >2.1 to 1 Ga (asteroid flux-derived chronology) or from 3.7 to 3.5 Ga (lunar-derived chronology), respectively. The Rheasilvian Period covers the time span between the Rheasilvia and Marcia impact events, and the Marcian Period covers the time between the Marcia impact event until the present. The age of the Marcia impact is still uncertain, but our current best estimates from crater counts of the ejecta blanket suggest an age between ∼120 and 390 Ma, depending upon choice of chronology system used. Regardless, the Marcia impact represents the youngest major geologic event on Vesta. Our proposed four-period geologic time scale for Vesta is, to a first order, comparable to those developed for other airless terrestrial bodies.

We report on a preliminary global geologic map of Vesta, based on data from the Dawn spacecraft’s High-Altitude Mapping Orbit (HAMO) and informed by Low-Altitude Mapping Orbit (LAMO) data. This map is part of an iterative mapping effort; the geologic map has been refined with each improvement in resolution. Vesta has a heavily-cratered surface, with large craters evident in numerous locations. The south pole is dominated by an impact structure identified before Dawn’s arrival. Two large impact structures have been resolved: the younger, larger Rheasilvia structure, and the older, more degraded Veneneia structure. The surface is also characterized by a system of deep, globe-girdling equatorial troughs and ridges, as well as an older system of troughs and ridges to the north. Troughs and ridges are also evident cutting across, and spiraling arcuately from, the Rheasilvia central mound.

However, no volcanic features have been unequivocally identified. Vesta can be divided very broadly into three terrains: heavily-cratered terrain; ridge-and-trough terrain (equatorial and northern); and terrain associated with the Rheasilvia crater. Localized features include bright and dark material and ejecta (some defined specifically by color); lobate deposits; and mass-wasting materials. No obvious volcanic features are evident. Stratigraphy of Vesta’s geologic units suggests a history in which formation of a primary crust was followed by the formation of impact craters, including Veneneia and the associated Saturnalia Fossae unit. Formation of Rheasilvia followed, along with associated structural deformation that shaped the Divalia Fossae ridge-and-trough unit at the equator. Subsequent impacts and mass wasting events subdued impact craters, rims and portions of ridge-and-trough sets, and formed slumps and landslides, especially within crater floors and along crater rims and scarps. Subsequent to the formation of Rheasilvia, discontinuous low-albedo deposits formed or were emplaced; these lie stratigraphically above the equatorial ridges that likely were formed by Rheasilvia. The last features to be formed were craters with bright rays and other surface mantling deposits.

Executed progressively throughout data acquisition, the iterative mapping process provided the team with geologic proto-units in a timely manner. However, interpretation of the resulting map was hampered by the necessity to provide the team with a standard nomenclature and symbology early in the process. With regard to mapping and interpreting units, the mapping process was hindered by the lack of calibrated mineralogic information. Topography and shadow played an important role in discriminating features and terrains, especially in the early stages of data acquisition.