ASU Scholarship Showcase

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.

We know very little about how soil-borne pollutants such as selenium (Se) can impact pollinators, even though Se has contaminated soils and plants in areas where insect pollination can be critical to the functioning of both agricultural and natural ecosystems. Se can be biotransferred throughout the food web, but few studies have examined its effects on the insects that feed on Se-accumulating plants, particularly pollinators. In laboratory bioassays, we used proboscis extension reflex (PER) and taste perception to determine if the presence of Se affected the gustatory response of honey bee (Apis mellifera L., Hymenoptera: Apidae) foragers. Antennae and proboscises were stimulated with both organic (selenomethionine) and inorganic (selenate) forms of Se that commonly occur in Se-accumulating plants. Methionine was also tested. Each compound was dissolved in 1 M sucrose at 5 concentrations, with sucrose alone as a control. Antennal stimulation with selenomethionine and methionine reduced PER at higher concentrations. Selenate did not reduce gustatory behaviors. Two hours after being fed the treatments, bees were tested for sucrose response threshold. Bees fed selenate responded less to sucrose stimulation. Mortality was higher in bees chronically dosed with selenate compared with a single dose. Selenomethionine did not increase mortality except at the highest concentration. Methionine did not significantly impact survival. Our study has shown that bees fed selenate were less responsive to sucrose, which may lead to a reduction in incoming floral resources needed to support coworkers and larvae in the field. If honey bees forage on nectar containing Se (particularly selenate), reductions in population numbers may occur due to direct toxicity. Given that honey bees are willing to consume food resources containing Se and may not avoid Se compounds in the plant tissues on which they are foraging, they may suffer similar adverse effects as seen in other insect guilds.

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.

Octopamine plays an important role in many behaviors in invertebrates. It acts via binding to G protein coupled receptors located on the plasma membrane of responsive cells. Several distinct subtypes of octopamine receptors have been found in invertebrates, yet little is known about the expression pattern of these different receptor subtypes and how each subtype may contribute to different behaviors. One honey bee (Apis mellifera) octopamine receptor, AmOA1, was recently cloned and characterized. Here we continue to characterize the AmOA1 receptor by investigating its distribution in the honey bee brain. We used two independent antibodies produced against two distinct peptides in the carboxyl-terminus to study the distribution of the AmOA1 receptor in the honey bee brain. We found that both anti-AmOA1 antibodies revealed labeling of cell body clusters throughout the brain and within the following brain neuropils: the antennal lobes; the calyces, pedunculus, vertical (alpha, gamma) and medial (beta) lobes of the mushroom body; the optic lobes; the subesophageal ganglion; and the central complex. Double immunofluorescence staining using anti-GABA and anti-AmOA1 receptor antibodies revealed that a population of inhibitory GABAergic local interneurons in the antennal lobes express the AmOA1 receptor in the cell bodies, axons and their endings in the glomeruli. In the mushroom bodies, AmOA1 receptors are expressed in a subpopulation of inhibitory GABAergic feedback neurons that ends in the visual (outer half of basal ring and collar regions) and olfactory (lip and inner basal ring region) calyx neuropils, as well as in the collar and lip zones of the vertical and medial lobes. The data suggest that one effect of octopamine via AmOA1 in the antennal lobe and mushroom body is to modulate inhibitory neurons.

Neural representations of odors are subject to computations that involve sequentially convergent and divergent anatomical connections across different areas of the brains in both mammals and insects. Furthermore, in both mammals and insects higher order brain areas are connected via feedback connections. In order to understand the transformations and interactions that this connectivity make possible, an ideal experiment would compare neural responses across different, sequential processing levels. Here we present results of recordings from a first order olfactory neuropile – the antennal lobe (AL) – and a higher order multimodal integration and learning center – the mushroom body (MB) – in the honey bee brain. We recorded projection neurons (PN) of the AL and extrinsic neurons (EN) of the MB, which provide the outputs from the two neuropils. Recordings at each level were made in different animals in some experiments and simultaneously in the same animal in others. We presented two odors and their mixture to compare odor response dynamics as well as classification speed and accuracy at each neural processing level. Surprisingly, the EN ensemble significantly starts separating odor stimuli rapidly and before the PN ensemble has reached significant separation. Furthermore the EN ensemble at the MB output reaches a maximum separation of odors between 84–120 ms after odor onset, which is 26 to 133 ms faster than the maximum separation at the AL output ensemble two synapses earlier in processing. It is likely that a subset of very fast PNs, which respond before the ENs, may initiate the rapid EN ensemble response. We suggest therefore that the timing of the EN ensemble activity would allow retroactive integration of its signal into the ongoing computation of the AL via centrifugal feedback.

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.