Magma-Driven Multiple Dike Propagation and Fracture Toughness of Crustal Rocks

Dike swarms consisting of tens to thousands of subparallel dikes are commonly observed at Earth's surface, raising the possibility of simultaneous propagation of two or more dikes at various stages of a swarm's development. The behavior of multiple propagating dikes differs from that of a single dike owing to the interacting stress fields associated with each dike. We analyze an array of parallel, periodically spaced dikes that grow simultaneously from an overpressured source into a semi-infinite, linear elastic host rock. To simplify the analysis, we assume steady state (constant velocity) magma flow and dike propagation. We use a perturbation method to analyze the coupled, nonlinear problem of multiple dike propagation and magma transport. The stress intensity factor at the dike tips and the opening displacements of the dike surfaces are calculated. The numerical results show that dike spacing has a profound effect on the behavior of dike propagation. The stress intensity factors at the tips of parallel dikes decrease with a decrease in dike spacing and are significantly smaller than that for a single dike with the same length. The reduced stress intensity factor indicates that, compared to a single dike, propagation of parallel dikes is more likely to be arrested under otherwise the same conditions. It also implies that fracture toughness of the host rock in a high confining pressure environment may not be as high as inferred from the propagation of a single dike. Our numerical results suggest fracture toughness values on the order of 100 MPa root m. The opening displacements for parallel dikes are smaller than that for a single dike, which results in higher magma pressure gradients in parallel dikes and lower flux of magma transport.

Surface "Waves" on Byrd Glacier, Antarctica

Byrd Glacier has one of the largest ice catchment areas in Antarctica, delivers more ice to the Ross Ice Shelf than any other ice stream, and is the fastest of these ice streams. A force balance, combined with a mass balance, demonstrates that stream flow in Byrd Glacier is transitional from sheet flow in East Antarctica to shelf flow in the Ross Ice Shelf. The longitudinal pulling stress, calculated along an ice flowband from the force balance, is linked to variations of ice thickness, to the ratio of the basal water pressure to the ice overburden pressure where Byrd Glacier is grounded, and is reduced by an ice-shelf buttressing stress where Byrd Glacier is floating. Longitudinal tension peaks at pressure-ratio maxima in grounded ice and close to minima in the ratio of the pulling stress to the buttressing stress in floating ice. The longitudinal spacing of these tension peaks is rather uniform and, for grounded ice, the peaks occur at maxima in surface slope that have no clear relation to the bed slope. This implies that the maxima in surface slope constitute a "wave train" that is related to regular variations in ice-bed coupling, not primarily to bed topography. It is unclear whether these surface "waves" are "standing waves" or are migrating either upslope or downslope, possibly causing the grounding line to either retreat or advance. Deciding which is the case will require obtaining bed topography in the map plane, a new map of surface topography, and more sophisticated modeling that includes ice flow linked to subglacial hydrology in the map plane.

Twentieth Century Delta δC Variability in Surface Water Dissolved Inorganic Carbon Recorded by Coralline Algae in the Northern North Pacific Ocean and the Bering Sea

The oxygen isotopic composition and Mg/Ca ratios in the skeletons of long-lived coralline algae record ambient seawater temperature over time. Similarly, the carbon isotopic composition in the skeletons record delta(13)C values of ambient seawater dissolved inorganic carbon. Here, we measured delta(13)C in the coralline alga Clathromorphum nereostratum to test the feasibility of reconstructing the intrusion of anthropogenic CO(2) into the northern North Pacific Ocean and Bering Sea. The delta(13)C was measured in the high Mgcalcite skeleton of three C. nereostratum specimens from two islands 500 km apart in the Aleutian archipelago. In the records spanning 1887 to 2003, the average decadal rate of decline in delta(13)C values increased from 0.03% yr(-1) in the 1960s to 0.095% yr(-1) in the 1990s, which was higher than expected due to solely the delta(13)C-Suess effect. Deeper water in this region exhibits higher concentrations of CO(2) and low delta(13)C values. Transport of deeper water into surface water (i.e., upwelling) increases when the Aleutian Low is intensified. We hypothesized that the acceleration of the delta(13)C decline may result from increased upwelling from the 1960s to 1990s, which in turn was driven by increased intensity of the Aleutian Low. Detrended delta(13)C records also varied on 4-7 year and bidecadal timescales supporting an atmospheric teleconnection of tropical climate patterns to the northern North Pacific Ocean and Bering Sea manifested as changes in upwelling.