Barbara Romanowicz

Bio: Barbara Romanowicz is an academic researcher from University of California, Berkeley. The author has contributed to research in topics: Mantle (geology) & Seismic tomography. The author has an hindex of 67, co-authored 284 publications receiving 14950 citations. Previous affiliations of Barbara Romanowicz include Collège de France & Centre national de la recherche scientifique.

Broad plumes rooted at the base of the Earth's mantle beneath major hotspots

Abstract: Plumes of hot upwelling rock rooted in the deep mantle have been proposed as a possible origin of hotspot volcanoes, but this idea is the subject of vigorous debate. On the basis of geodynamic computations, plumes of purely thermal origin should comprise thin tails, only several hundred kilometres wide, and be difficult to detect using standard seismic tomography techniques. Here we describe the use of a whole-mantle seismic imaging technique--combining accurate wavefield computations with information contained in whole seismic waveforms--that reveals the presence of broad (not thin), quasi-vertical conduits beneath many prominent hotspots. These conduits extend from the core-mantle boundary to about 1,000 kilometres below Earth's surface, where some are deflected horizontally, as though entrained into more vigorous upper-mantle circulation. At the base of the mantle, these conduits are rooted in patches of greatly reduced shear velocity that, in the case of Hawaii, Iceland and Samoa, correspond to the locations of known large ultralow-velocity zones. This correspondence clearly establishes a continuous connection between such zones and mantle plumes. We also show that the imaged conduits are robustly broader than classical thermal plume tails, suggesting that they are long-lived, and may have a thermochemical origin. Their vertical orientation suggests very sluggish background circulation below depths of 1,000 kilometres. Our results should provide constraints on studies of viscosity layering of Earth's mantle and guide further research into thermochemical convection.

The three‐dimensional shear velocity structure of the mantle from the inversion of body, surface and higher‐mode waveforms

Abstract: SUMMARY We present a 3-D model of shear heterogeneity (Vsh) in the whole mantle, derived from the inversion of hand-picked body, surface and higher-mode waveforms. The forward and inverse problems are formulated using the non-linear asymptotic coupling theory, the zeroth-order asymptotic expansion of a Born seismogram computed by normal mode summation and including coupling across mode branches. We invert iteratively for mantle heterogeneity and for centroid moment tensors until convergence. Our model, SAW24B16, is parametrized laterally in spherical harmonics up to degree 24 and radially in 16 cubic b-splines with knots spaced to reflect the data sampling with depth. The power spectrum of the model is dominated by heterogeneity in the boundary layers at low degrees: degrees 5 and 6 near the surface and degree 2 near the core–mantle boundary. The rest of the lower mantle is dominated by degree 3. We find significant heterogeneity up to spherical harmonic degree 24 in the transition zone and the uppermost part of the lower mantle. Our model displays three slow domains extending continuously between uppermost and lowermost mantle, one originating under Africa and two in the Pacific. The results of our modelling also suggest an interaction near the surface between the anomalies in the Pacific and the network of mid-ocean ridges in the East Pacific, and between the African anomaly and the low-velocity zone in the Red Sea, the East African Rift, the Mid-Indian Ridge and the Mid-Atlantic Ridge. These anomalies appear to ascend into the upper mantle without altering their shape near the 670 km discontinuity. Fast anomalies accumulate near the 670 km discontinuity under the Southern Kurile, the Japan trench and the Izu arc, with little or no penetration into the lower mantle, in agreement with results from a variety of regional studies. Cross-sections through subduction zones such as the Marianas, Tonga–Kermadec, Java, the palaeosubduction zone beneath North America and the Peru–Chile trench show continuous fast structures penetrating into the lower mantle, with the depth of penetration varying from region to region. The distribution of both fast and slow anomalies near 670 km thus suggests that the discontinuity does not act as a strong barrier to upgoing and downgoing flow between the upper and the lower mantle.

Lithospheric layering in the North American craton

Abstract: How cratons-extremely stable continental areas of the Earth's crust-formed and remained largely unchanged for more than 2,500 million years is much debated. Recent studies of seismic-wave receiver function data have detected a structural boundary under continental cratons at depths too shallow to be consistent with the lithosphere-asthenosphere boundary, as inferred from seismic tomography and other geophysical studies. Here we show that changes in the direction of azimuthal anisotropy with depth reveal the presence of two distinct lithospheric layers throughout the stable part of the North American continent. The top layer is thick ( approximately 150 km) under the Archaean core and tapers out on the surrounding Palaeozoic borders. Its thickness variations follow those of a highly depleted layer inferred from thermo-barometric analysis of xenoliths. The lithosphere-asthenosphere boundary is relatively flat (ranging from 180 to 240 km in depth), in agreement with the presence of a thermal conductive root that subsequently formed around the depleted chemical layer. Our findings tie together seismological, geochemical and geodynamical studies of the cratonic lithosphere in North America. They also suggest that the horizon detected in receiver function studies probably corresponds to the sharp mid-lithospheric boundary rather than to the more gradual lithosphere-asthenosphere boundary.

Global anisotropy and the thickness of continents.

Abstract: For decades there has been a vigorous debate about the depth extent of continental roots. The analysis of heat-flow, mantle-xenolith and electrical-conductivity data all indicate that the coherent, conductive part of continental roots (the 'tectosphere') is at most 200-250 km thick. Some global seismic tomographic models agree with this estimate, but others suggest that a much thicker zone of high velocities lies beneath continental shields, reaching a depth of at least 400 km. Here we show that this disagreement can be reconciled by taking into account seismic anisotropy. We show that significant radial anisotropy, with horizontally polarized shear waves travelling faster than those that are vertically polarized, is present under most cratons in the depth range 250-400 km--similar to that found under ocean basins at shallower depths of 80-250 km. We propose that, in both cases, the anisotropy is related to shear in a low-viscosity asthenospheric channel, located at different depths under continents and oceans. The seismically defined 'tectosphere' is then at most 200-250 km thick under old continents. The 'Lehmann discontinuity', observed mostly under continents at about 200-250 km, and the 'Gutenberg discontinuity', observed under oceans at depths of about 60-80 km, may both be associated with the bottom of the lithosphere, marking a transition to flow-induced asthenospheric anisotropy.

Global mantle shear velocity model developed using nonlinear asymptotic coupling theory

Abstract: We present a three-dimensional shear velocity model of the whole mantle developed using S H waveform data. The model is expressed horizontally in terms of spherical harmonics up to degree 12, and vertically in terms of Legendre polynomials up to degrees 5 and 7 in the upper and lower mantle, respectively. What distinguishes this model from other tomographic models published to date is (1) the theoretical normal mode-based wave propagation approach, where we include across branch mode coupling terms in order to model the body wave sensitivity to structure along the path more accurately; (2) the wave-packet weighting scheme which allows to balance contributions from high-amplitude and low-amplitude phases, increasing the resolution in some parts of the mantle. We also relax the constraints on the Moho depth, which is allowed to vary in the inversion, thus absorbing some uncertainties in crustal structure. The resulting model is generally in good agreement with other recent global mantle S velocity models and with some regional models. The rms profile with depth has more power than other models in the upper mantle/lower mantle transition region and the zone of increased power and low degree structure near the base of the mantle is confined to the last 500 km in depth. This model provides a particularly good fit to the non-hydrostatic geoid through harmonic degree 12 (79% variance reduction), as well as good fits to observed splitting functions of S velocity sensitive mantle modes, indicating that both large-scale and small-scale features are really well constrained.

Phd by thesis

Abstract: Deposits of clastic carbonate-dominated (calciclastic) sedimentary slope systems in the rock record have been identified mostly as linearly-consistent carbonate apron deposits, even though most ancient clastic carbonate slope deposits fit the submarine fan systems better. Calciclastic submarine fans are consequently rarely described and are poorly understood. Subsequently, very little is known especially in mud-dominated calciclastic submarine fan systems. Presented in this study are a sedimentological core and petrographic characterisation of samples from eleven boreholes from the Lower Carboniferous of Bowland Basin (Northwest England) that reveals a >250 m thick calciturbidite complex deposited in a calciclastic submarine fan setting. Seven facies are recognised from core and thin section characterisation and are grouped into three carbonate turbidite sequences. They include: 1) Calciturbidites, comprising mostly of highto low-density, wavy-laminated bioclast-rich facies; 2) low-density densite mudstones which are characterised by planar laminated and unlaminated muddominated facies; and 3) Calcidebrites which are muddy or hyper-concentrated debrisflow deposits occurring as poorly-sorted, chaotic, mud-supported floatstones. These

New Empirical Relationships among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement

Abstract: Source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude ( M ), surface rupture length, subsurface rupture length, downdip rupture width, rupture area, and maximum and average displacement per event. The resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment, moment magnitude, subsurface rupture length, downdip rupture width, and average surface displacement. Each source parameter is classified as reliable or unreliable, based on our evaluation of the accuracy of individual values. Only the reliable source parameters are used in the final analyses. In comparing source parameters, we note the following trends: (1) Generally, the length of rupture at the surface is equal to 75% of the subsurface rupture length; however, the ratio of surface rupture length to subsurface rupture length increases with magnitude; (2) the average surface displacement per event is about one-half the maximum surface displacement per event; and (3) the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement. Thus, for most earthquakes in this data base, slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface. Log-linear regressions between earthquake magnitude and surface rupture length, subsurface rupture length, and rupture area are especially well correlated, showing standard deviations of 0.25 to 0.35 magnitude units. Most relationships are not statistically different (at a 95% significance level) as a function of the style of faulting: thus, we consider the regressions for all slip types to be appropriate for most applications. Regressions between magnitude and displacement, magnitude and rupture width, and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area. The large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data. Separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions. Regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults, and to estimate the expected surface displacement along a fault for a given size earthquake.

Geologic Evolution of the Himalayan-Tibetan Orogen

Abstract: A review of the geologic history of the Himalayan-Tibetan orogen suggests that at least 1400 km of north-south shortening has been absorbed by the orogen since the onset of the Indo-Asian collision at about 70 Ma. Significant crustal shortening, which leads to eventual construction of the Cenozoic Tibetan plateau, began more or less synchronously in the Eocene (50–40 Ma) in the Tethyan Himalaya in the south, and in the Kunlun Shan and the Qilian Shan some 1000–1400 km in the north. The Paleozoic and Mesozoic tectonic histories in the Himalayan-Tibetan orogen exerted a strong control over the Cenozoic strain history and strain distribution. The presence of widespread Triassic flysch complex in the Songpan-Ganzi-Hoh Xil and the Qiangtang terranes can be spatially correlated with Cenozoic volcanism and thrusting in central Tibet. The marked difference in seismic properties of the crust and the upper mantle between southern and central Tibet is a manifestation of both Mesozoic and Cenozoic tectonics. The form...

The Mechanics of Earthquakes and Faulting

Abstract: This essential reference for graduate students and researchers provides a unified treatment of earthquakes and faulting as two aspects of brittle tectonics at different timescales. The intimate connection between the two is manifested in their scaling laws and populations, which evolve from fracture growth and interactions between fractures. The connection between faults and the seismicity generated is governed by the rate and state dependent friction laws - producing distinctive seismic styles of faulting and a gamut of earthquake phenomena including aftershocks, afterslip, earthquake triggering, and slow slip events. The third edition of this classic treatise presents a wealth of new topics and new observations. These include slow earthquake phenomena; friction of phyllosilicates, and at high sliding velocities; fault structures; relative roles of strong and seismogenic versus weak and creeping faults; dynamic triggering of earthquakes; oceanic earthquakes; megathrust earthquakes in subduction zones; deep earthquakes; and new observations of earthquake precursory phenomena.

Current plate motions

Abstract: A global plate motion model, named NUVEL-1, which describes current plate motions between 12 rigid plates is described, with special attention given to the method, data, and assumptions used Tectonic implications of the patterns that emerged from the results are discussed It is shown that wide plate boundary zones can form not only within the continental lithosphere but also within the oceanic lithosphere; eg, between the Indian and Australian plates and between the North American and South American plates Results of the model also suggest small but significant diffuse deformation of the oceanic lithosphere, which may be confined to small awkwardly shaped salients of major plates