David M. Blair

Preimpact porosity controls the gravity signature of lunar craters

We model the formation of lunar complex craters and investigate the effect of preimpact porosity on their gravity signatures. We find that while preimpact target porosities less than ~7% produce negative residual Bouguer anomalies (BAs), porosities greater than ~7% produce positive anomalies whose magnitude is greater for impacted surfaces with higher initial porosity. Negative anomalies result from pore space creation due to fracturing and dilatant bulking, and positive anomalies result from destruction of pore space due to shock wave compression. The central BA of craters larger than ~215 km in diameter, however, are invariably positive because of an underlying central mantle uplift. We conclude that the striking differences between the gravity signatures of craters on the Earth and Moon are the result of the higher average porosity and variable porosity of the lunar crust.

Formation of the Orientale lunar multiring basin

titOn the origin of Orientale basinle Orientale basin is a major impact crater on the Moon, which is hard to see from Earth because it is right on the western edge of the lunar nearside. Relatively undisturbed by later events, Orientale serves as a prototype for understanding large impact craters throughout the solar system. Zuber et al. used the Gravity Recovery and Interior Laboratory (GRAIL) mission to map the gravitational field around the crater in great detail by flying the twin spacecraft as little as 2 km above the surface. Johnson et al. performed a sophisticated computer simulation of the impact and its subsequent evolution, designed to match the data from GRAIL. Together, these studies reveal how major impacts affect the lunar surface and will aid our understanding of other impacts on rocky planets and moons. Science, this issue pp. 438 and 441 Simulations of the formation of the Orientale basin on the Moon reveal the origin of its multiple rings Multiring basins, large impact craters characterized by multiple concentric topographic rings, dominate the stratigraphy, tectonics, and crustal structure of the Moon. Using a hydrocode, we simulated the formation of the Orientale multiring basin, producing a subsurface structure consistent with high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) spacecraft. The simulated impact produced a transient crater, ~390 kilometers in diameter, that was not maintained because of subsequent gravitational collapse. Our simulations indicate that the flow of warm weak material at depth was crucial to the formation of the basin’s outer rings, which are large normal faults that formed at different times during the collapse stage. The key parameters controlling ring location and spacing are impactor diameter and lunar thermal gradients.

Computer-Aided Discovery: Toward Scientific Insight Generation with Machine Support

The process of scientific discovery is traditionally assumed to be entirely executed by humans. This article highlights how increasing data volumes and human cognitive limits are challenging this traditional assumption. Relevant examples are found in observational astronomy and geoscience, disciplines that are undergoing transformation due to growing networks of space-based and ground-based sensors. The authors outline how intelligent systems for computer-aided discovery can routinely complement and integrate human scientists in the insight generation loop in scalable ways for next-generation science. The pragmatics of model-based computer-aided discovery systems go beyond feature detection in empirical data to answer fundamental questions, such as how empirical detections fit into hypothesized models and model variants to ease the scientists work of placing large ensembles of detections into a theoretical context. The authors demonstrate successful applications of this paradigm in several areas, including ionospheric studies, volcanics, astronomy, and planetary landing site identification for spacecraft and robotic missions.

Clustering Throughput Optimization on the GPU

Large datasets in astronomy and geoscience often require clustering and visualizations of phenomena at different densities and scales in order to generate scientific insight. We examine the problem of maximizing clustering throughput for concurrent dataset clustering in spatial dimensions. We introduce a novel hybrid approach that uses GPUs in conjunction with multicore CPUs for algorithmic throughput optimizations. The key idea is to exploit the fast memory on the GPU for index searches and optimize I/O transfers in such a way that the low-bandwidth host-GPU bottleneck does not have a significant negative performance impact. To achieve this, we derive two distinct GPU kernels that exploit grid-based indexing schemes to improve clustering performance. To obviate limited GPU memory and enable large dataset clustering, our method is complemented by an efficient batching scheme for transfers between the host and GPU accelerator. This scheme is robust with respect to both sparse and dense data distributions and intelligently avoids buffer overflows that would otherwise degrade performance, all while minimizing the number of data transfers between the host and GPU. We evaluate our approaches on ionospheric total electron content datasets as well as intermediate-redshift galaxies from the Sloan Digital Sky Survey. Our hybrid approach yields a speedup of up to 50x over the sequential implementation on one of the experimental scenarios, which is respectable for I/O intensive clustering.

Exploiting Variant-Based Parallelism for Data Mining of Space Weather Phenomena

This paper studies a form of parallelism termed variant-based parallelism, which exploits commonalities and reuse among variant computations in order to improve multithreading scalability. The problem is motivated by space weather studies that aim to identify changes in the Earths ionosphere caused by auroral activity, tsunamis, and earthquakes. Today it is common to execute cluster algorithm variants with different parameters in order to determine which ones best explain phenomena in empirical data. We propose a novel approach and a set of optimizations to maximize throughput in such clustering algorithms. This is achieved by executing multiple clustering algorithm variants in parallel and developing efficient approaches to concurrently cluster data and maximize the reuse of results from completed variants. We present evaluations on real-world space weather datasets with up to 5 million ionospheric total electron content data points as well as synthetic datasets with up to a million data points. Results show a 1101% performance improvement due to indexing tailored for variant-based clustering, and a 2209% performance improvement when applying all of our proposed optimizations. Our optimizations enable new approaches in computer-aided discovery and could enable the short run times required for early warning systems for natural hazards.

Computer aided detection of transient inflation events at Alaskan volcanoes using GPS measurements from 2005–2015

Abstract Analysis of transient deformation events in time series data observed via networks of continuous Global Positioning System (GPS) ground stations provide insight into the magmatic and tectonic processes that drive volcanic activity. Typical analyses of spatial positions originating from each station require careful tuning of algorithmic parameters and selection of time and spatial regions of interest to observe possible transient events. This iterative, manual process is tedious when attempting to make new discoveries and does not easily scale with the number of stations. Addressing this challenge, we introduce a novel approach based on a computer-aided discovery system that facilitates the discovery of such potential transient events. The advantages of this approach are demonstrated by actual detections of transient deformation events at volcanoes selected from the Alaska Volcano Observatory database using data recorded by GPS stations from the Plate Boundary Observatory network. Our technique successfully reproduces the analysis of a transient signal detected in the first half of 2008 at Akutan volcano and is also directly applicable to 3 additional volcanoes in Alaska, with the new detection of 2 previously unnoticed inflation events: in early 2011 at Westdahl and in early 2013 at Shishaldin. This study also discusses the benefits of our computer-aided discovery approach for volcanology in general. Advantages include the rapid analysis on multi-scale resolutions of transient deformation events at a large number of sites of interest and the capability to enhance reusability and reproducibility in volcano studies.

Evidence of large empty lava tubes on the Moon using GRAIL gravity

NASAs GRAIL mission employed twin spacecraft in polar orbits around the Moon to measure the lunar gravity field at unprecedentedly high accuracy and resolution. The low spacecraft altitude in the extended mission enables the detection of small-scale surface or subsurface features. We analyzed these data for evidence of empty lava tubes beneath the lunar maria. We developed two methods, gradiometry and cross-correlation, to isolate the target signal of long, narrow, sinuous mass deficits from a host of other features present in the GRAIL data. Here, we report the discovery of several strong candidates that are either extensions of known lunar rilles, collocated with the recently discovered “skylight” caverns, or underlying otherwise unremarkable surfaces. Owing to the spacecraft polar orbits, our techniques are most sensitive to east-west trending near-surface structures and empty lava tubes with minimum widths of several kilometers, heights of hundreds of meters, and lengths of tens of kilometers.

Optimizing Parallel Clustering Throughput in Shared Memory

This article studies the optimization of parallel clustering throughput in the context of variant-based parallelism, which exploits commonalities and reuse among variant computations for multithreading scalability. This direction is motivated by challenging scientific applications where scientists have to execute multiple runs of clustering algorithms with different parameters to determine which ones best explain phenomena observed in empirical data. To make this process more efficient, we propose a novel set of optimizations to maximize the throughput of Density-Based Spatial Clustering of Applications with Noise (DBSCAN), a frequently used algorithm for scientific data mining in astronomy, geoscience, and many other fields. Our approach executes multiple algorithm variants in parallel, computes clusters concurrently, and leverages heuristics to maximize the reuse of results from completed variants. As scientific datasets continue to grow, maximizing clustering throughput with our techniques may accelerate the search and identification of natural phenomena of interest with computational support, i.e., Computer-Aided Discovery. We present evaluations on a whole spectrum of datasets, such as geoscience data on space weather phenomena, astronomical data from the Sloan Digital Sky Survey on intermediate-redshift galaxies, as well as synthetic datasets to characterize performance properties. Selected results show a 1,115 percent performance improvement due to indexing tailored for variant-based clustering, and a 2,209 percent performance improvement when applying all of our proposed optimizations.