A Unifying Phase Diagram for the Dynamics of Sheared Solids and Granular Materials

We present a simple unifying model that can be used to analyze, within a single framework, different dynamic regimes of shear deformation of brittle, plastic, and granular materials. The basic dynamic regimes seen in the response of both solids and granular materials to slowly increasing loading are scale-invariant behavior with power law statistics, quasi-periodicity of system size events, and persisting long term mode switching between the former two types of response. The model provides universal analytical mean field results on the statistics of failure events in the different regimes and distributed versus localized spatial responses. The results are summarized in a phase diagram spanned by three tuning parameters: dynamic strength change (weakening, neutral or strengthening) during slip events, dissipation of stress transfer (related to the void fraction in granular materials and damaged solids), and the ratio of shear rate over healing rate controlling the regaining of cohesion following failures in brittle solids. The mean field scaling predictions agree with experimental, numerical, and observational data on deformation avalanches of solids, granular materials, and earthquake faults. The model provides additional predictions that should be tested with future observation and simulation data.     

Parity and disparity between two Chama oysters: the reproductive biology of the Indo-Pacific C. pacifica Broderip,invasive to the Mediterranean Sea; and C. savignyi Lamy,indigenous to the Red Sea

The Indo-Pacific oyster Chama pacifica Broderip, 1835 (Mollusca, Bivalvia, Chamidae) is rarely found in the Northern Red Sea reefs of Eilat (Gulf of Aqaba), where it is outnumbered by its indigenous congener, Chama savignyi Lamy, 1921. The influx of Eritrean biota from the Red Sea into the Mediterranean Sea via the Suez Canal has led to the formation of massive Chama oyster beds along the Eastern Mediterranean shore. However, unlike the Northern Red Sea, the Eastern Mediterranean Chama beds are dominated by C. pacifica oysters, whereas C. savignyi is absent from this region. In an attempt to understand this difference in their respective distribution, the reproductive biology of both species was compared. Histological analysis of the male and female gonads, monitored monthly from March 2009 to August 2010 in both regions, revealed a similar reproductive cycle, comprising six stages: onset of gametogenesis, advanced gametogenesis, ripe, ready to spawn, spent and sexual rest. Female gonads demonstrated an additional, seventh stage – restoration, coinciding with inferred spawning of ripe gametes. Both species were found to be dioecious spawners, with a single, annual, temperature-dependent inferred spawning period. Chama pacifica was found to reproduce efficiently in maximal Mediterranean seawater temperatures not experienced by the Northern Red Sea Chama populations. This study demonstrates the high invasive potential of an oyster species despite its rarity in its source region.     

Application of high resolution DEM data to detect rock damage from geomorphic signals along the central San Jacinto Fault

We analyze geomorphic properties extracted from LiDAR and SRTM (Shuttle Radar Topography Mission) data to test whether the damage zone along the central San Jacinto Fault (SJF) zone can be resolved with remotely-sensed data in a quantitative fashion. The SJF is one of the most active faults in southern California, with well expressed geomorphology and a fast slip rate, as seen in the geology and by GPS. We use ArcMap and the TauDEM toolbox to compare several morphometric parameters, including drainage density (Dd), on both sides of the fault, using a 1 km and a 5 km buffer for the LiDAR and SRTM data, respectively. We also analyze the spatial patterns of Dd near the fault, using two different definitions of spatial Dd. The high resolution of the LiDAR data allows us to focus on a single fault, eliminating the effects of parallel nearby faults. From the LiDAR data we find that the highest Dd values occur in areas between two fault strands, followed generally by rocks on the northeast side of the fault, with the lowest Dd values occurring on the southwest side of the fault. The SRTM data shows a band of high Dd values centered on the main fault trace with ~ 1 km width. Our results indicate that there is a strong correlation between drainage density and proximity to the fault, with zones of structural complexity along the fault displaying the highest Dd. We interpret this to largely be an effect of degree of rock damage, as these are areas that are expected to be more damaged, and field observations support this contention. If we are correct, then it appears that the northeast side of the SJF is generally more damaged. South of the trifurcation area there is evidence that the signal is reversed on the larger scale, with more damage on the southwest side of the fault inferred from the SRTM data, possibly caused by extension between the Coyote Creek and Clark faults. The implications of the observed asymmetry could be geological evidence for rupture propagation direction, because a preferred propagation direction is predicted to produce asymmetric damage structure that would be recorded in the volume of rock surrounding a fault.     

Pre- and post-1998 ENSO records of shallow-water octocorals (Alcyonacea) in the Chagos Archipelago

When compared, principal octocorals collected in the Chagos Archipelago before and after the 1998 ENSO shared many common taxa. While a few discontinuities in their biodiversity revealed subtle changes in more persistent genera (Lobophytum, Sarcophyton), some fast-growing “fugitive” genera (e.g. Cespitularia, Efflatounaria, Heteroxenia) disappeared after the ENSO-related coral bleaching. Such transient fugitives might thus be eliminated from soft coral communities on isolated reef systems, possibly in the long term, by events of this nature. The appearance of Carijoa riseii, a species often considered a fouling organism, even an invasive, may well be indicative of reef degradation during the ENSO event. The post-ENSO recovery manifested by this fauna nevertheless gives cause for hope for their survival in the face of climate change.     

Decay of dynamic fracturing based on three-dimensional measurements of clastic-dike geometry

Several fundamental issues of fracture mechanics during the post-dynamic stage are yet not fully understood, including fracture arrest mechanisms, effects of the three-dimensional fracture propagation on fracture aperture and height relations, and the role played by fracture tips on fracture termination. We studied these issues in the seismically active Dead Sea basin, where clastic dikes (>10 m) and numerous smaller dikelets (<1 m) dynamically intruded the late-Pleistocene soft rock of the Lisan Formation. A three-dimensional study of the dikelets shows that they form arrest zones at the tips of the larger clastic dikes. Geometrically, the dikelets are divided into two parts: (1) the main dikelet, in which the aperture profile along the dikelet height is approximately elliptic; and (2) the elongated tip, in which the aperture profile along the tip height is approximately constant. The dikelet aspect ratio is defined as A/H, where A is dikelet aperture and H is dikelet height. A plot of A/H versus H describes power relations with two different zones: (1) Zone A, with a small variation of A/H, between 0.02 and 0.06, for dikelets in height range of 100–700 mm; and (2) Zone B, where the aspect ratio sharply increases to 0.23 in dikelets with heights less than 100 mm. We interpret that during deceleration, when the height of the elongated tip became greater than 1/10 of the dikelet height, inelastic conditions are dominant. Under these conditions, the fracture velocity decreases significantly and the dikelet aspect ratio increases. The present observations and analyses indicate that formation of elongated tips and dike (fracture) segmentation are essential for the decay of the dynamic fracturing.     

The Role of Heterogeneities as a Tuning Parameter of Earthquake Dynamics

We investigate the influence of spatial heterogeneities on various aspects of brittle failure and seismicity in a model of a large strike-slip fault. The model dynamics is governed by realistic boundary conditions consisting of constant velocity motion of regions around the fault, static/kinetic friction laws, creep with depth-dependent coefficients, and 3-D elastic stress transfer. The dynamic rupture is approximated on a continuous time scale using a finite stress propagation velocity (quasidynamic model). The model produces a brittle-ductile transition at a depth of about 12.5 km, realistic hypocenter distributions, and other features of seismicity compatible with observations. Previous work suggested that the range of size scales in the distribution of strength-stress heterogeneities acts as a tuning parameter of the dynamics. Here we test this hypothesis by performing a systematic parameter-space study with different forms of heterogeneities. In particular, we analyze spatial heterogeneities that can be tuned by a single parameter in two distributions: (1) high stress drop barriers in near-vertical directions and (2) spatial heterogeneities with fractal properties and variable fractal dimension. The results indicate that the first form of heterogeneities provides an effective means of tuning the behavior while the second does not. In relatively homogeneous cases, the fault self-organizes to large-scale patches and big events are associated with inward failure of individual patches and sequential failures of different patches. The frequency-size event statistics in such cases are compatible with the characteristic earthquake distribution and large events are quasi-periodic in time. In strongly heterogeneous or near-critical cases, the rupture histories are highly discontinuous and consist of complex migration patterns of slip on the fault. In such cases, the frequency-size and temporal statistics follow approximately power-law relations.on leave from CNRS Rennes, France     

Software tools for accessing the GPS Seamless Archive

The Scripps Orbit and Permanent Array Center (SOPAC) has completed development for the UNAVCO community of first-generation GPS Seamless Archive (GSAC) software. The GSAC is a virtual archive composed of an assembly of agencies and investigators exchanging information about their respective GPS-related data holdings in a well defined, cohesive manner. The superset of this published information is collected and ingested into centralized databases administered currently by two data brokers (Retailers), who make the data available to the public in a seamless manner. There are three user interfaces available: the interactive GSAC Wizard, a command-line Unix-style executable called gsac-client, and a front door HTTP service called the GSAC Retailer Service Interface. Each user interface provides access to the data collections of 6 different GPS archives (GSAC Wholesalers) in North America. Together these archives have published more than 2 million GPS data files pertaining to over 10,000 different geodetic monuments. These datasets are composed in large part of data collected by US scientists and their collaborators over the period 1986 to the present in Western North America and other tectonically active regions around the globe, as well as the holdings of two IGS global data centers. In this article, we describe how the three GSAC user interfaces provide the community a powerful set of tools for seamlessly mining information and collecting data files from a distributed network of GPS archives.The GPS Toolbox is a column dedicated to highlighting algorithms and source code utilized by GPS Engineers and scientists. If you have an interesting program or software package you would like to share with our readers, please pass it along; e-mail it to us at gps-toolbox@ngs.noaa.gov. To comment on any of the source code discussed here, or to download source code, visit our website at . This column is edited by Stephen Hilla, National Geodetic Survey, NOAA, Silver Spring, Maryland, and Mike Craymer, Geodetic Survey Division, Natural Resources Canada, Ottawa, Ontario, Canada.