Extension of forward modeling phase-screen code to for AVO analysis in isotropic and anisotropic media

The computationally efficient phase-screen forward modeling technique is extended to allow investigation of non-normal ray paths. The code is developed to accommodate all diffracted and converted phases up to critical angle, building on a geometrical construction method. The new approach relies upon pre-scanning the model space to assess the complexity of each screen. The propagating wavefields are then divided as a function of horizontal wavenumber, and each subset is transformed to the spatial domain separately, carrying with it angular information. This allows both accurate 3-D phase corrections and Zoeppritz reflection and transmission coefficients to be applied. The phase-screen code is further developed to handle simple anisotropic media. During phasescreen modeling, propagation is undertaken in the wavenumber domain where exact expressions for anisotropic phase velocities are available. Travel-times and amplitude effects from a range of anisotropic shales are computed and compared with previous literature

Structure and dynamics of the Ecuador Fracture Zone, Panama Basin

In this study, multiple geophysical data types are used to investigate the structure and dynamics of the Ecuador Fracture Zone – a complex multi-stranded strike-slip fault system located in the Panama Basin. Gravity modelling reveals a 25-30 km-wide region of ∼3 km-thick, low-density crust beneath this system, and an anomalously low-density region in the uppermost mantle. Along both edges, the transition to the ‘normal’ structure and thickness oceanic crust formed at both the Ecuador and Costa Rica Rifts is abrupt. Within the Ecuador Fracture Zone itself, normal faults bound the median ridges. These faults traverse the entire thickness of accumulated sediment and offset the seabed, while sediment layer geometries document multiple phases of relative uplift, with the most recent phase still ongoing. Active extensional faulting, with an approximately spreading ridge-parallel strike, is also observed in 6-7 Ma Costa Rica Rift crust. The median ridges and the transverse ridge at the eastern edge of the Ecuador Fracture Zone also have contrasting crustal density structures. Both median ridges have a lower density crust than between the intervening valleys, while the transverse ridge crust has an equivalent thickness and density structure to that formed at the Costa Rica Rift. The active median valley basement-cutting normal faults allow seawater ingress and alternation of the crustal footwall, and also flow to mantle depth where, based on gravity modelling, 30-50% serpentinization of mantle peridotite occurs. The resulting serpentinite-driven buoyancy acts as the primary control on the observed median ridge relative vertical tectonism. In contrast, the relative uplift of the transverse ridge results from lithospheric flexure in response to a change in spreading direction between the Ecuador and Costa Rica Rifts. Contrary to the widely accepted assumption that fracture zones are tectonically inactive systems, the Ecuador Fracture Zone provides evidence of extension, serpentinization due to ongoing hydrothermal circulation, and relative uplift

Redefining relative biological effectiveness in the context of the EQDX formalism: implications for alpha-particle emitter therapy.

Alpha-particle radiopharmaceutical therapy (αRPT) is currently enjoying increasing attention as a viable alternative to chemotherapy for targeting of disseminated micrometastatic disease. In theory, αRPT can be personalized through pre-therapeutic imaging and dosimetry. However, in practice, given the particularities of α-particle emissions, a dosimetric methodology that accurately predicts the thresholds for organ toxicity has not been reported. This is in part due to the fact that the biological effects caused by α-particle radiation differ markedly from the effects caused by traditional external beam (photon or electron) radiation or β-particle emitting radiopharmaceuticals. The concept of relative biological effectiveness (RBE) is used to quantify the ratio of absorbed doses required to achieve a given biological response with alpha particles versus a reference radiation (typically a beta emitter or external beam radiation). However, as conventionally defined, the RBE varies as a function of absorbed dose and therefore a single RBE value is limited in its utility because it cannot be used to predict response over a wide range of absorbed doses. Therefore, efforts are underway to standardize bioeffect modeling for different fractionation schemes and dose rates for both nuclear medicine and external beam radiotherapy. Given the preponderant use of external beams of radiation compared to nuclear medicine in cancer therapy, the more clinically relevant quantity, the 2 Gy equieffective dose, EQD2(α/β), has recently been proposed by the ICRU. In concert with EQD2(α/β), we introduce a new, redefined RBE quantity, named RBE2(α/β), as the ratio of the two linear coefficients that characterize the α particle absorbed dose-response curve and the low-LET megavoltage photon 2 Gy fraction equieffective dose-response curve. The theoretical framework for the proposed new formalism is presented along with its application to experimental data obtained from irradiation of a breast cancer cell line. Radiobiological parameters are obtained using the linear quadratic model to fit cell survival data for MDA-MB-231 human breast cancer cells that were irradiated with either α particles or a single fraction of low-LET (137)Cs γ rays. From these, the linear coefficient for both the biologically effective dose (BED) and the EQD2(α/β) response lines were derived for fractionated irradiation. The standard RBE calculation, using the traditional single fraction reference radiation, gave RBE values that ranged from 2.4 for a surviving fraction of 0.82-6.0 for a surviving fraction of 0.02, while the dose-independent RBE2(4.6) value was 4.5 for all surviving fraction values. Furthermore, bioeffect modeling with RBE2(α/β) and EQD2(α/β) demonstrated the capacity to predict the surviving fraction of cells irradiated with acute and fractionated low-LET radiation, α particles and chronic exponentially decreasing dose rates of low-LET radiation. RBE2(α/β) is independent of absorbed dose for α-particle emitters and it provides a more logical framework for data reporting and conversion to equieffective dose than the conventional dose-dependent definition of RBE. Moreover, it provides a much needed foundation for the ongoing development of an α-particle dosimetry paradigm and will facilitate the use of tolerance dose data available from external beam radiation therapy, thereby helping to develop αRPT as a single modality as well as for combination therapies

Crustal structure of the French Guiana margin, West Equatorial Atlantic

Geophysical data from the Amazon Cone Experiment are used to determine the structure and evolution of the French Guiana and Northeast Brazil continental margin, and to better understand the origin and development of along-margin segmentation. A 427-km-long combined multichannel reflection and wide-angle refraction seismic profile acquired across the southern French Guiana margin is interpreted, where plate reconstructions suggest a rift-type setting. The resulting model shows a crustal structure in which 35–37-km-thick pre-rift continental crust is thinned by a factor of 6.4 over a distance of ∼70 km associated with continental break-up and the initiation and establishment of seafloor spreading. The ocean–continent boundary is a transition zone up to 45 km in width, in which the two-layered oceanic-type crustal structure develops. Although relatively thin at 3.5–5.0 km, such thin oceanic crust appears characteristic of the margin as a whole. There is no evidence of rift-related magmatism, either as seaward-dipping sequences in the reflection data or as a high velocity region in the lower crust in the P-wave velocity model, and as a such the margin is identified as non-volcanic in type. However, there is also no evidence of the rotated fault block and graben structures characteristic of rifted margins. Consequently, the thin oceanic crust, the rapidity of continental crustal thinning and the absence of characteristic rift-related structures leads to the conclusion that the southern French Guiana margin has instead developed in an oblique rift setting, in which transform motion also played a significant role in the evolution of the resulting crustal structure and along-margin segmentation in structural style

Evolution of Heat Flow, Hydrothermal Circulation and Permeability on the Young Southern Flank of the Costa Rica Rift

We analyze 67 new conductive heat flow measurements on the southern flank of the Costa Rica Rift (CRR). Heat flow measurements cover five sites ranging in oceanic crustal age between approximately 1.6 and 5.7 Ma, and are co-located with a high-resolution multi-channel seismic line that extends from slightly north of the first heat flow site (1.6 Ma) to beyond ODP Hole 504B in 6.9 Ma crust. For the five heat flow sites, the mean observed conductive heat flow is ≈ 85 mWm−2. This value is approximately 30 per cent of the mean lithospheric heat flux expected from a half-space conductive cooling model, indicating that hydrothermal processes account for about 70 per cent of the heat loss. The advective heat loss fraction varies from site to site and is explained by a combination of outcrop to outcrop circulation through exposed basement outcrops and discharge through faults. Super-critical convection in Layer 2A extrusives occurs between 1.6 and 3.5 Ma, and flow through a thinly-sedimented basement high occurs at 4.6 Ma. Advective heat loss diminishes rapidly between ≈ 4.5 and ≈ 5.7 Ma, which contrasts with plate cooling reference models that predict a significant deficit in conductive heat flow up to ages ≈ 65 ± 10 Ma. At ≈ 5.7 Ma the CRR topography is buried under sediment with an average thickness ≈ 150 m, and hydrothermal circulation in the basement becomes sub-critical or perhaps marginally critical. The absence of significant advective heat loss at ≈ 5.7 Ma at the CRR is thus a function of both burial of basement exposure under the sediment load and a reduction in basement permeability that possibly occurs as result of mineral precipitation and original permeability at the time of formation. Permeability is a non-monotonic function of age along the southern flank of the CRR, in general agreement with seismic velocity tomography interpretations that reflect variations in the degree of ridge-axis magma supply and tectonic extension. Hydrothermal circulation in the young oceanic crust at southern flank of CRR is affected by the interplay and complex interconnectedness of variations in permeability, sediment thickness, topographical structure, and tectonic and magmatic activities with age

Stochastic Heterogeneity Mapping around a Mediterranean salt lens

We present the first application of Stochastic Heterogeneity Mapping based on the band-limited von Kármán function to a seismic reflection stack of a Mediterranean water eddy (meddy), a large salt lens of Mediterranean water. This process extracts two stochastic parameters directly from the reflectivity field of the seismic data: the Hurst number, which ranges from 0 to 1, and the correlation length (scale length). Lower Hurst numbers represent a richer range of high wavenumbers and correspond to a broader range of heterogeneity in reflection events. The Hurst number estimate for the top of the meddy (0.39) compares well with recent theoretical work, which required values between 0.25 and 0.5 to model internal wave surfaces in open ocean conditions based on simulating a Garrett-Munk spectrum (GM76) slope of −2. The scale lengths obtained do not fit as well to seismic reflection events as those used in other studies to model internal waves. We suggest two explanations for this discrepancy: (1) due to the fact that the stochastic parameters are derived from the reflectivity field rather than the impedance field the estimated scale lengths may be underestimated, as has been reported; and (2) because the meddy seismic image is a two-dimensional slice of a complex and dynamic three-dimensional object, the derived scale lengths are biased to the direction of flow. Nonetheless, varying stochastic parameters, which correspond to different spectral slopes in the Garrett-Munk spectrum (horizontal wavenumber spectrum), can provide an estimate of different internal wave scales from seismic data alone. We hence introduce Stochastic Heterogeneity Mapping as a novel tool in physical oceanography

Influence of enhanced melt supply on upper crustal structure at a mid-ocean ridge discontinuity: A three-dimensional seismic tomographic study of 9°N East Pacific Rise

We present a three-dimensional upper crustal model of the 9°03′N overlapping spreading center (OSC) on the East Pacific Rise that assists in understanding the relationship between melt sills and upper crustal structure at a ridge discontinuity with enhanced melt supply at crustal levels. Our P wave velocity model obtained from tomographic inversion of ∼70,000 crustal first arrival travel times suggests that the geometry of extrusive emplacement are significantly different beneath the overlapping spreading limbs. Extrusive volcanic rocks above the western melt sill are inferred to be thin (∼250 m). More extensive accumulation of extrusives is inferred to the west than to the east of the western melt sill. The extrusive layer inferred above the eastern melt sill thickens from ∼350 (at the neovolcanic axis) to 550 m (to the west of the melt sill). Volcanic construction is likely to be significant in the formation of ridge crest morphology at the OSC, particularly at the tip of the eastern limb. On the basis of our interpretation of the velocity model, we propose that enhanced magma supply at crustal levels at the OSC may provide an effective mechanism for the migration of ridge discontinuities. This “dynamic magma supply model” may explain the commonly observed nonsteady migration pattern of ridge discontinuities by attributing this to the temporal fluctuations in melt availability to the overlapping spreading limbs

Abrupt episode of mid-Cretaceous ocean acidification triggered by massive volcanism

Large igneous province volcanic activity during the mid-Cretaceous approximately 94.5 million years ago triggered a global-scale episode of reduced marine oxygen levels known as Oceanic Anoxic Event 2. It has been hypothesized that this geologically rapid degassing of volcanic carbon dioxide altered seawater carbonate chemistry, affecting marine ecosystems, geochemical cycles, and sedimentation. Here, we report on two sites drilled by the International Ocean Discovery Program offshore of southwest Australia that exhibit clear evidence for suppressed pelagic carbonate sedimentation in the form of a stratigraphic interval barren of carbonate, recording ocean acidification during the event. We then use the osmium isotopic composition of bulk sediments to directly link this protracted ~600- kiloyear shoaling of the marine calcite compensation depth to the onset of volcanic activity. This decrease in marine pH was prolonged by biogeochemical feedbacks in highly productive regions that elevated heterotrophic respiration of carbon dioxide to the water column. A compilation of mid- Cretaceous marine stratigraphic records reveals a contemporaneous decrease of sedimentary carbonate content at continental slope sites globally. Thus, we contend that changes in marine carbonate chemistry are a primary ecological stress and important consequence of rapid emission of carbon dioxide during many large igneous province eruptions in the geologic past

Are the magnetic fields of millisecond pulsars ~ 10^8 G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are $\sim 10^{8}$G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The $\sim 10^{8}$ G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have $< 10^{11}$ G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit $10^{11}$ G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict $\sim 10^{8}$ G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted $\sim 10^{8}$ G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields $\gg 10^{8}$ G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields $> 10^{8}$ G with their magnetic dipole emission being balanced by low level accretion from their companions.Comment: 16 pages, accept for publication in Astrophysics and Space Scienc