The slab geometry laser. I - Theory
NASA Technical Reports Server (NTRS)
Eggleston, J. M.; Kane, T. J.; Kuhn, K.; Byer, R. L.; Unternahrer, J.
1984-01-01
Slab geometry solid-state lasers offer significant performance improvements over conventional rod-geometry lasers. A detailed theoretical description of the thermal, stress, and beam-propagation characteristics of a slab laser is presented. The analysis includes consideration of the effects of the zig-zag optical path, which eliminates thermal and stress focusing and reduces residual birefringence.
Convection in Slab and Spheroidal Geometries
NASA Technical Reports Server (NTRS)
Porter, David H.; Woodward, Paul R.; Jacobs, Michael L.
2000-01-01
Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.
Numerical quadrature for slab geometry transport algorithms
Hennart, J.P.; Valle, E. del
1995-12-31
In recent papers, a generalized nodal finite element formalism has been presented for virtually all known linear finite difference approximations to the discrete ordinates equations in slab geometry. For a particular angular directions {mu}, the neutron flux {Phi} is approximated by a piecewise function Oh, which over each space interval can be polynomial or quasipolynomial. Here we shall restrict ourselves to the polynomial case. Over each space interval, {Phi} is a polynomial of degree k, interpolating parameters given by in the continuous and discontinuous cases, respectively. The angular flux at the left and right ends and the k`th Legendre moment of {Phi} over the cell considered are represented as.
NASA Astrophysics Data System (ADS)
Wu, J. E.; Suppe, J.; Renqi, L.; Lin, C.; Kanda, R. V.
2013-12-01
The past locations, shapes and polarity of subduction trenches provide first-order constraints for plate tectonic reconstructions. Analogue and numerical models of subduction zones suggest that relative subducting (Vs) and overriding (Vor) plate velocities may strongly influence final subducted slab geometries. Here we have mapped the 3D geometries of subducted slabs in the upper and lower mantle of Asia from global seismic tomography. We have incorporated these slabs into plate tectonic models, which allows us to infer the subducting and overriding plate velocities. We describe two distinct slab geometry styles, ';flat slabs' and ';slab curtains', and show their implications for paleo-trench positions and subduction geometries in plate tectonic reconstructions. When compared to analogue and numerical models, the mapped slab styles show similarities to modeled slabs that occupy very different locations within Vs:Vor parameter space. ';Flat slabs' include large swaths of sub-horizontal slabs in the lower mantle that underlie the well-known northward paths of India and Australia from Eastern Gondwana, viewed in a moving hotspot reference. At India the flat slabs account for a significant proportion of the predicted lost Ceno-Tethys Ocean since ~100 Ma, whereas at Australia they record the existence of a major 8000km by 2500-3000km ocean that existed at ~43 Ma between East Asia, the Pacific and Australia. Plate reconstructions incorporating the slab constraints imply these flat slab geometries were generated when continent overran oceanic lithosphere to produce rapid trench retreat, or in other words, when subducting and overriding velocities were equal (i.e. Vs ~ Vor). ';Slab curtains' include subvertical Pacific slabs near the Izu-Bonin and Marianas trenches that extend from the surface down to 1500 km in the lower mantle and are 400 to 500 km thick. Reconstructed slab lengths were assessed from tomographic volumes calculated at serial cross-sections. The ';slab
Coarse-mesh diffusion synthetic acceleration in slab geometry
Kim, K.S.; Palmer, T.S.
2000-07-01
It has long been known that the success of a diffusion synthetic acceleration (DSA) scheme is very sensitive to the consistency between the discretization of the transport and diffusion acceleration equations. Acceleration schemes involving inconsistent discretizations have been successful, but no prescription is available that determines a priori an allowable degree of inconsistency. It is notable, however, that all current DSA schemes involve diffusion equations discretized on the spatial mesh used to solve the transport equations. Often the solution of a large number of low-order equations is an expensive part of the transport simulation. This motivates the desire to find stable and rapidly convergent acceleration schemes that are discretized on a mesh that is coarse relative to the transport mesh. The authors present here results showing that the low-order diffusion equation can be solved on a mesh coarser (by a factor of 2) than that used for the slab geometry transport equation. Their results show that coarse-mesh DSA is unconditionally stable and is as rapidly convergent as a DSA method discretized on the transport mesh. They have used Adams and Martin's modified four-step acceleration method (M4S) applied to the linear discontinuous (LD) finite element transport equations in slab geometry. To evaluate their procedure, they have performed a Fourier analysis to calculate theoretical spectral radii. They compare this analysis with convergence behavior observed in an implementation code for several model problems.
Slab2 - Providing updated subduction zone geometries and modeling tools to the community
NASA Astrophysics Data System (ADS)
Hayes, G. P.; Hearne, M. G.; Portner, D. E.; Borjas, C.; Moore, G.; Flamme, H.
2015-12-01
The U.S. Geological Survey database of global subduction zone geometries (Slab1.0) combines a variety of geophysical data sets (earthquake hypocenters, moment tensors, active source seismic survey images of the shallow subduction zone, bathymetry, trench locations, and sediment thickness information) to image the shape of subducting slabs in three dimensions, at approximately 85% of the world's convergent margins. The database is used extensively for a variety of purposes, from earthquake source imaging, to magnetotelluric modeling. Gaps in Slab1.0 exist where input data are sparse and/or where slabs are geometrically complex (and difficult to image with an automated approach). Slab1.0 also does not include information on the uncertainty in the modeled geometrical parameters, or the input data used to image them, and provides no means to reproduce the models it described. Currently underway, Slab2 will update and replace Slab1.0 by: (1) extending modeled slab geometries to all global subduction zones; (2) incorporating regional data sets that may describe slab geometry in finer detail than do previously used teleseismic data; (3) providing information on the uncertainties in each modeled slab surface; (4) modifying our modeling approach to a fully-three dimensional data interpolation, rather than following the 2-D to 3-D steps of Slab1.0; (5) migrating the slab modeling code base to a more universally distributable language, Python; and (6) providing the code base and input data we use to create our models, such that the community can both reproduce the slab geometries, and add their own data sets to ours to further improve upon those models in the future. In this presentation we describe our vision for Slab2, and the first results of this modeling process.
NASA Astrophysics Data System (ADS)
Preston, L.; Creager, K.; Crosson, R.; Symons, N.
2002-12-01
The 1998 Wet SHIPS experiment, an onshore-offshore active source experiment within NW Washington State and SW British Columbia, provided travel times for over 1200 secondary arrivals that we interpret as wide-angle reflections from the subducting Juan de Fuca slab. The reflection points are distributed over a broad area under the Olympic Peninsula, WA. Using a nonlinear, iterative inverse procedure incorporating full 3-D ray tracing, we combine these reflected arrivals with 100,000 first arrivals obtained from SHIPS and other active source experiments and from local earthquakes to simultaneously determine smooth 3-D velocity structure, smooth slab geometry and hypocentral positions. The inversion converges stably after many iterations to the same solution over the range of starting models we have tested. The reflector is generally parallel to and within the upper half of the intraslab events, which form a 6-km thick dipping zone. The vast majority of these earthquakes occur where seismic velocities exceed 7.5 km/s. We are especially interested in testing the hypothesis that dehydration embrittlement, associated with the basalt to eclogite phase transformation, is responsible for the intraslab earthquakes. This hypothesis predicts that the majority of the intraslab events should occur within the subducted oceanic crust above a Moho reflector. An unambiguous interpretation requires knowledge of the reflector depth, relative to the earthquakes, with relative errors less than a few km. We present results from a variety of empirical and numerical error and resolution tests including jack-knife and checkerboard analyses. We analyze the interrelationships and trade-offs among the reflector position, intraslab earthquakes, and the fine-scale deep velocity structure to determine error and resolution bounds in this crucial area of our model. In general, the active source experiments provide dense crossing rays with known source locations to give excellent resolution and
62-dB-gain multiple-pass slab geometry Nd:YAG amplifier
NASA Technical Reports Server (NTRS)
Kane, T. J.; Kozlovsky, W. J.; Byer, R. L.
1986-01-01
A flash-lamp-pumped slab Nd:YAG laser amplifier with 62 dB of gain is demonstrated. The amplifier remains linear for energy extraction up to 12 mJ, and the beam spatial quality remains good at repetition rates up to 10 Hz. The slab geometry provides for gain averaging in the zigzag dimension, thus eliminating the undesirable effects of nonuniform pumping. This averaging and the necessity for only a single compact gain element are valuable attributes that may allow for eventual laser-diode array pumping of the slab Nd:YAG amplifier.
Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity
McCrory, Patricia A.; Blair, J. Luke; Waldhause, Felix; Oppenheimer, David H.
2012-01-01
A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.
Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity
NASA Astrophysics Data System (ADS)
McCrory, Patricia A.; Blair, J. Luke; Waldhauser, Felix; Oppenheimer, David H.
2012-09-01
A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.
Seismicity and shallow slab geometry in the central Vanuatu subduction zone
NASA Astrophysics Data System (ADS)
Baillard, Christian; Crawford, Wayne C.; Ballu, Valérie; Régnier, Marc; Pelletier, Bernard; Garaebiti, Esline
2015-08-01
The Vanuatu arc in the southwest Pacific Ocean is one of the world's most seismically active regions, with almost 39 magnitude 7+ earthquakes in the past 43 years. Convergence rates are around 90-120 mm/yr along most of the arc, but drop to 25-43 mm/yr in the central section, probably due to the subduction of the d'Entrecasteaux ridge. We characterize the slab geometry and tectonic state in this central section by analyzing data from a 10 month deployment of 30 seismometers over this section. We located more than 30,000 events (all less than magnitude 5.5), constructed an improved 1-D velocity model, calculated focal mechanisms and cluster geometries, and determined the 3-D geometry of the interplate seismogenic zone. The seismogenic zone has a shallow bulge in front of the d'Entrecasteaux ridge, which could be explained by the ridge's buoyancy contributing to the uplift of the fore-arc islands. The seismogenic zone extends to ~45 km depth, significantly below the 26-27 km depth of the fore-arc Moho, indicating that the upper mantle wedge is not significantly serpentinized, which is consistent with the relatively high thermal parameter of the subducting plate. The maximum width of the seismogenic zone is 80 km, indicating an upper earthquake magnitude limit of Mw 7.85 ± 0.4, assuming standard rupture zone aspect ratios. The data also reveal a double seismic zone, 20 to 30 km below the seismogenic zone, which is presumably caused by flexure of the downgoing plate.
Geometry and deformation of the subducting slab beneath the Altiplano
NASA Astrophysics Data System (ADS)
Kumar, A.; Wagner, L.; Beck, S. L.; Zandt, G.
2011-12-01
One of the largest continental plateaus on Earth - the Altiplano - Puna in the central Andes presents challenges to the scientific community due to the complex geodynamic processes that are responsible for its formation and evolution. Current theories predict two different modes of evolution for Altiplano-Puna plateau: (1) Slow-and-steady rise of the plateau that isostatically compensated horizontal crustal shortening or (2) rapid rise of the plateau due to large scale crustal shortening that induced episodic removal of dense lower lithosphere. Previous geophysical studies along the central and southern portion of the Altiplano suggest that the thickness of the crust varies across the Andean plateau from 70 km below the highest topography in the Eastern and Western Cordilleras to 59-64 km in the central Altiplano. Velocity profiling across the plateau indicates the presence of high velocity layer of mafic lower crust below the Western Cordillera while the central Altiplano and Eastern Cordillera have no high velocity lower crust. We present preliminary earthquake relocations and focal mechanisms using data from the currently deployed network of 50 broadband seismic stations that is part of the NSF-Continental Dynamics-funded project "CAUGHT" (Central Andean Uplift and the Geodynamics of High Topography). This deployment ranges from ~13oS to 18oS across the northern flank of the Andean plateau in Peru and Bolivia. Our relocated events and focal mechanisms provide a much better constrained look at the subducted plate and deforming over-riding crust than have previously been possible using the global network alone. In addition to illuminating processes associated with the formation of the northern Altiplano, the relocated slab events and new focal mechanisms provide new constraints on the nature of the transition from flat-slab subduction beneath the northernmost Altiplano to normal slab dip subduction beneath the rest of the Altiplano to the south. References 1
Casimir effects for classical and quantum liquids in slab geometry: A brief review
Biswas, Shyamal
2015-05-15
We analytically explore Casimir effects for confinement of classical and quantum fluctuations in slab (film) geometry (i) for classical (critical) fluctuations over {sup 4}He liquid around the λ point, and (ii) for quantum (phonon) fluctuations of Bogoliubov excitations over an interacting Bose-Einstein condensate. We also briefly review Casimir effects for confinement of quantum vacuum fluctuations confined to two plates of different geometries.
NASA Astrophysics Data System (ADS)
Gibert, G.; Gerbault, M.; Hassani, R.; Tric, E.
2012-05-01
The aim of this study is to quantify the relationship between the kinematics of subduction, deformation in the overriding plate and the evolution of slab geometry. A 2-D finite element numerical code is used, and a first objective consists in benchmarking previously published analogue models. Far-field plate velocities are applied, and once the subducting plate reaches the 660 km discontinuity, modelled as a rigid base, we obtain two different forms or styles of subduction that depend on the overriding plate velocity vop: if vop > 0, the slab lies forwards on the 660 km discontinuity (style 1), and if vop≤ 0, the slab lies backwards on the discontinuity (style 2). We also obtain a cyclic pattern with the slab folding on itself repeatedly when vsp > 0 and 2vop+vsp > 0 (where vsp is the subducting plate velocity). These conditions result from the analysis of several simulations in which the subduction velocities and plate viscosities are varied. When the slab periodically folds on the 660 km discontinuity, periods of shallow slab dip and compression in the overriding plate are followed by periods of slab steepening and relative extension in the overriding plate. Folding periodicity is controlled by the slab viscosity and subduction velocity. When a low-viscosity zone is included in the overriding plate, the trench motion is effectively decoupled from the overriding plate velocity, therefore allowing it to be directly controlled by the deep dynamics of the slab. For the cyclic style 2 corresponding to forward folding of the slab, the low-viscosity region in the overriding plate increases the stress amplitudes oscillations, the trench motion and the folding periodicity with time. Therefore the strength of the entire overriding plate is shown to directly control the dynamics of subduction. Using the Nazca and South American plate velocities we produce models of cyclic folding with a period of ca. 22 Ma and a minimal dip angle of ca. 10°. Episodic folding of the slab
Error analysis of the quadratic nodal expansion method in slab geometry
Penland, R.C.; Turinsky, P.J.; Azmy, Y.Y.
1994-10-01
As part of an effort to develop an adaptive mesh refinement strategy for use in state-of-the-art nodal diffusion codes, the authors derive error bounds on the solution variables of the quadratic Nodal Expansion Method (NEM) in slab geometry. Closure of the system is obtained through flux discontinuity relationships and boundary conditions. In order to verify the analysis presented, the authors compare the quadratic NEM to the analytic solution of a test problem. The test problem for this investigation is a one-dimensional slab [0,20cm] with L{sup 2} = 6.495cm{sup 2} and D = 0.1429cm. The slab has a unit neutron source distributed uniformly throughout and zero flux boundary conditions. The analytic solution to this problem is used to compute the node-average fluxes over a variety of meshes, and these are used to compute the NEM maximum error on each mesh.
Geometry and State of Stress of the Slab Beneath the North Central Andes
NASA Astrophysics Data System (ADS)
Kumar, A.; Beck, S. L.; Wagner, L. S.; Zandt, G.; Long, M. D.
2012-12-01
The central Andean plateau of southern Peru and Bolivia is one of the largest topographic features on Earth. It has strongly influenced the local and regional climate since the early Miocene by affecting the regional dynamics that control circulation and precipitation. The surface and subsurface processes responsible for the plateau formation and evolution are still unclear. There are two end member models proposed for this uplift: (1) Slow and steady rise since the late Eocene (~40 Ma) with maximum upper crustal shortening between 30 and 10 Ma or (2) rapid surface uplift of ~2.5 km in the late Miocene between 10.3 and 6.7 Ma. The rapid uplift theory argues for the wholesale removal of a thick portion of the lower eclogitic crust and upper mantle lithosphere. A slow and steady uplift of the Andes would suggest a continuous removal of the lower lithosphere or piecemeal delamination, proportional to the rate of shortening. We present earthquake locations and focal mechanisms using data from two ongoing temporary arrays: the network of 50 broadband seismic stations that was part of the NSF-Continental Dynamics-funded project "CAUGHT" (Central Andean Uplift and the Geodynamics of High Topography) and the 40 station NSF- Geophysics funded "PULSE" array (PerU Lithosphere and Slab Experiment). Our new earthquake locations provide an improved insight about the geometry of subducting Nazca slab and also put an upper bound on the thickness of overriding lithosphere. Obvious clustering of intermediate depth earthquakes suggests strong and localized release of tectonic stress in the slab at ~15.5oS. The seismic section drawn from the precisely located slab events provide a better idea about the lateral variations of the slab geometry and the geometry of asthenoshperic corner flow to help understand its geodynamic effect on the lithospheric delamination or ablative subduction process. . Focal mechanisms of the slab events are helpful in understanding the stress state of the
NASA Astrophysics Data System (ADS)
Wada, I.; Wang, K.; He, J.
2013-12-01
In this study, we revisit the effects of along-strike variation in slab geometry and oblique subduction on subduction zone thermal structures. Along-strike variations in slab dip cause changes in the descending rate of the slab and generate trench-parallel pressure gradients that drive trench-parallel mantle flow (e.g., Kneller and van Keken, 2007). Oblique subduction also drives trench-parallel mantle flow. In this study, we use a finite element code PGCtherm3D and examine a range of generic subduction geometries and parameters to investigate the effects of the above two factors. This exercise is part of foundational work towards developing detailed 3-D thermal models for NE Japan, Nankai, and Cascadia to better constrain their 3-D thermal structures and to understand the role of temperature in controlling metamorphic, seismogenic, and volcanic processes. The 3-D geometry of the subducting slabs in the forearc and arc regions are well delineated at these three subduction zones. Further, relatively large compilations of surface heat flow data at these subduction zones make them excellent candidates for this study. At NE Japan, a megathrust earthquake occurred on March 11, 2011; at Nankai and Cascadia, there has been a great effort to constrain the scale of the next subduction thrust earthquake for the purpose of disaster prevention. Temperature influences the slip behavior of subduction faults by (1) affecting the rheology of the interface material and (2) controlling dehydration reactions, which can lead to elevated pore fluid pressure. Beyond the depths of subduction thrust earthquakes, the thermal structure is affected strongly by the pattern of mantle wedge flow. This flow is driven by viscous coupling between the subducting slab and the overriding mantle, and it brings in hot flowing mantle into the wedge. The trench-ward (up-dip) extent of the slab-mantle coupling is thus a key factor that controls the thermal structure. Slab-mantle decoupling at shallow
Complex geometry of the subducted Pacific slab inferred from receiver function
NASA Astrophysics Data System (ADS)
Zhang, Ruiqing; Wu, Qingju; Zhang, Guangcheng
2014-05-01
In recent years, slab tear has received considerable attention and been reported in many arc-arc junctures in Pacific plate subdution zones. From 2009 to 2011, we deployed two portable experiments equipped with CMG-3ESPC seismometers and the recorders of REFTEK-130B in NE China. The two linear seismic arrays were designed nearly parallel, and each of them containing about 60 seismic stations extended about 1200 km from west to east spanning all surface geological terrains of NE China. The south one was firstly set up and continually operated over two year, while the north deployment worked only about one year. By using the teleseismic data collected by these two arrays, we calculate the P receiver functions to map topographic variation of the upper mantle discontinuities. Our sampled region is located where the juncture between the subducting Kuril and Japan slabs reaches the 660-km discontinuity. Distinct variation of the 660-km discontinuity is mapped beneath the regions. A deeper-than-normal 660 km discontinuity is observed locally in the southeastern part of our sampled region. The depression of the 660 km discontinuity may be resulted from an oceanic lithospheric slab deflected in the mantle transition zone, in good agreement with the result of earlier tomographic and other seismic studies in this region. The northeastern portion of our sampled region, however, does not show clearly the deflection of the slab. The variation of the tomography of the 660-km discontinuity in our sampled regions may indicate a complex geometry of the subducted Pacific slab.
NASA Technical Reports Server (NTRS)
Stamnes, K.; Lie-Svendsen, O.; Rees, M. H.
1991-01-01
The linear Boltzmann equation can be cast in a form mathematically identical to the radiation-transport equation. A multigroup procedure is used to reduce the energy (or velocity) dependence of the transport equation to a series of one-speed problems. Each of these one-speed problems is equivalent to the monochromatic radiative-transfer problem, and existing software is used to solve this problem in slab geometry. The numerical code conserves particles in elastic collisions. Generic examples are provided to illustrate the applicability of this approach. Although this formalism can, in principle, be applied to a variety of test particle or linearized gas dynamics problems, it is particularly well-suited to study the thermalization of suprathermal particles interacting with a background medium when the thermal motion of the background cannot be ignored. Extensions of the formalism to include external forces and spherical geometry are also feasible.
Mathews, K.; Sjoden, G.; Minor, B. )
1994-09-01
The exponential characteristic spatial quadrature for discrete ordinates neutral particle transport in slab geometry is derived and compared with current methods. It is similar to the linear characteristic (or, in slab geometry, the linear nodal) quadrature but differs by assuming an exponential distribution of the scattering source within each cell, S(x) = a exp(bx), whose parameters are root-solved to match the known (from the previous iteration) average and first moment of the source over the cell. Like the linear adaptive method, the exponential characteristic method is positive and nonlinear but more accurate and more readily extended to other cell shapes. The nonlinearity has not interfered with convergence. The authors introduce the exponential moment functions,'' a generalization of the functions used by Walters in the linear nodal method, and use them to avoid numerical ill-conditioning. The method exhibits O([Delta]x[sup 4]) truncation error on fine enough meshes; the error is insensitive to mesh size for coarse meshes. In a shielding problem, it is accurate to 10% using 16-mfp-thick cells; conventional methods err by 8 to 15 orders of magnitude. The exponential characteristic method is computationally more costly per cell than current methods but can be accurate with very thick cells, leading to increased computational efficiency on appropriate problems.
Slab Geometry Control on Mantle Flow Regime: A case study from Central South America Subduction Zone
NASA Astrophysics Data System (ADS)
Biryol, C. B.; Beck, S. L.; Zandt, G.; Wagner, L. S.
2013-12-01
The subduction of oceanic lithosphere along convergent plate margins plays an important role in the dynamics of the upper mantle beneath convergent margins and major orogenic belts. Many studies of mantle dynamics show that the flow pattern of the mantle varies greatly between different subduction zones as well as within the same subduction zone. The factors that control such variations are poorly understood and need to be investigated further in order to develop a better understanding of various subduction zone processes such as the deformation of mantle beneath convergent plate margins and transport of melts and volatiles in the mantle wedge above subducting slabs. Earlier studies of mantle flow inferred from seismic anisotropy via shear-wave splitting analysis indicated that the dynamics and deformation of subducting and overriding plates as well as the slab geometry have important roles on mantle flow regime. In an effort to test the significance of these factors in constraining the mantle dynamics along the central South America subduction zone, we carried out a shear-wave splitting analysis. Our study area covers southern Peru and northwestern Bolivia encompassing the northernmost Altiplano plateau where subduction of the Nazca plate begins to gradually flatten towards the north. The major part of the data for our analysis comes from the CAUGHT temporary seismic deployment (2010 - 2012) with 49 three-component broadband seismometers. In our study we used SKS, SKKS and PKS arrivals from over 80 teleseismic earthquakes, located between the distance-range of 60 to 140 degrees. We determined polarization direction and delay-time of shear-wave arrivals that are polarized into fast and slow components and split in time. The resultant fast polarization directions indicate the direction of mantle flow beneath the study area and the delay-times show the strength and depth extend of the associated seismic anisotropy. The results of our analysis revealed a
Simulations of Coulomb systems with slab geometry using an efficient 3D Ewald summation method
NASA Astrophysics Data System (ADS)
dos Santos, Alexandre P.; Girotto, Matheus; Levin, Yan
2016-04-01
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.
Spatial Treatment of the Slab-geometry Discrete Ordinates Equations Using Artificial Neural Networks
Brantley, P S
2001-03-23
An artificial neural network (ANN) method is developed for treating the spatial variable of the one-group slab-geometry discrete ordinates (S{sub N}) equations in a homogeneous medium with linearly anisotropic scattering. This ANN method takes advantage of the function approximation capability of multilayer ANNs. The discrete ordinates angular flux is approximated by a multilayer ANN with a single input representing the spatial variable x and N outputs representing the angular flux in each of the discrete ordinates angular directions. A global objective function is formulated which measures how accurately the output of the ANN approximates the solution of the discrete ordinates equations and boundary conditions at specified spatial points. Minimization of this objective function determines the appropriate values for the parameters of the ANN. Numerical results are presented demonstrating the accuracy of the method for both fixed source and incident angular flux problems.
NASA Astrophysics Data System (ADS)
Wada, Ikuko; He, Jiangheng; Hasegawa, Akira; Nakajima, Junichi
2015-09-01
We develop a 3-D thermal model for the Northeast Japan subduction margin, using a realistic slab geometry for the subducting Pacific plate, and investigate the effects of oblique subduction and 3-D slab geometry on the mantle wedge flow pattern and the thermal structure. In the Tohoku region, the mantle wedge flow pattern is nearly two-dimensional resulting in a thermal structure similar to those obtained by a 2-D model, owing to the simple slab geometry and subduction nearly perpendicular to the margin. However, in Hokkaido, oblique subduction leads to 3-D mantle wedge flow with northerly inflow and west-northwestward outflow and also results in lower temperatures in the shallow part of the mantle wedge than in Tohoku due to lower sinking rate of the slab. Between Hokkaido and Tohoku, the slab has a hinge-like shape due to a relatively sharp change in the dip direction. In this hinge zone, northerly mantle inflow from Hokkaido and westerly mantle inflow from Tohoku converge, discouraging inflow from northwest and resulting in a cooler mantle wedge. The model-predicted mantle wedge flow patterns are consistent with observed seismic anisotropy and may explain the orientations of volcanic cross-chains. The predicted 3-D thermal structure correlates well with the along-arc variations in the location of the frontal arc volcanoes and help to provide new insights into the surface heat flow pattern and the down-dip extent of interplate earthquakes.
Community interactive webtool to retrieve Greenland glacier data for 1-D geometry
NASA Astrophysics Data System (ADS)
Perrette, Mahé
2015-04-01
Marine-terminating, outlet glaciers are challenging to include in conventional Greenland-wide ice sheet models because of the large variation in scale between model grid size (typically 10 km) and outlet glacier width (typically 1-5km), making it a subgrid scale feature. A possible approach to tackle this problem is to use one-dimensional flowline models for the individual glaciers (e.g. Nick et al., 2013, Nature; Enderlin et al 2013a,b, The Cryosphere). Here we present a python- and javascript- based webtool to prepare data required to feed in or validate a flowline model. It is designed primarily to outline the glacier geometry and returns relevant data averaged over cross-sections. The tool currently allows to: visualize 2-D ice sheet data (zoom/pan), quickly switch between datasets (e.g. ice thickness, bedrock elevation, surface velocity) interpolated / transformed on a common grid. draw flowlines from user-input seeds on the map, calculated from a vector field of surface velocity, as an helpful guide for point 3 interactively draw glacier outline (side and middle lines) on top of the data mesh the outlined glacier domain in the horizontal plane extract relevant data into a 1-D longitudinal profile download the result as a netCDF file The project is hosted on github to encourage collaboration, under the open-source MIT Licence. The server-side is written in python (open-source) using the web-framework flask, and the client-side (javascript) makes use of the d3 library for interactive figures. For now it only works locally in a web browser (start server: "python runserver.py"). Data need to be downloaded separately from the original sources. See the README file in the project for information how to use it. Github projects: https://github.com/perrette/webglacier1d (main) https://github.com/perrette/dimarray (dependency)
The dynamical Cooperative Lamb Shift in a system of two-level atoms in a slab-geometry
NASA Astrophysics Data System (ADS)
Friedberg, Richard; Manassah, Jamal T.
2009-09-01
Using the eigenmode analysis, we compute the Cooperative Lamb Shift ( CLS) as a function of time from the vector photon model for a system of two-level atoms in a slab-geometry for forward and backward emission in two initial states of weak excitation and complete inversion.
Westerly, David C.; Mo, Xiaohu; Tomé, Wolfgang A.; Mackie, Thomas R.; DeLuca, Paul M.
2013-01-01
Purpose: Pencil beam algorithms are commonly used for proton therapy dose calculations. Szymanowski and Oelfke [“Two-dimensional pencil beam scaling: An improved proton dose algorithm for heterogeneous media,” Phys. Med. Biol. 47, 3313–3330 (2002)10.1088/0031-9155/47/18/304] developed a two-dimensional (2D) scaling algorithm which accurately models the radial pencil beam width as a function of depth in heterogeneous slab geometries using a scaled expression for the radial kernel width in water as a function of depth and kinetic energy. However, an assumption made in the derivation of the technique limits its range of validity to cases where the input expression for the radial kernel width in water is derived from a local scattering power model. The goal of this work is to derive a generalized form of 2D pencil beam scaling that is independent of the scattering power model and appropriate for use with any expression for the radial kernel width in water as a function of depth. Methods: Using Fermi-Eyges transport theory, the authors derive an expression for the radial pencil beam width in heterogeneous slab geometries which is independent of the proton scattering power and related quantities. The authors then perform test calculations in homogeneous and heterogeneous slab phantoms using both the original 2D scaling model and the new model with expressions for the radial kernel width in water computed from both local and nonlocal scattering power models, as well as a nonlocal parameterization of Molière scattering theory. In addition to kernel width calculations, dose calculations are also performed for a narrow Gaussian proton beam. Results: Pencil beam width calculations indicate that both 2D scaling formalisms perform well when the radial kernel width in water is derived from a local scattering power model. Computing the radial kernel width from a nonlocal scattering model results in the local 2D scaling formula under-predicting the pencil beam width by as
NASA Astrophysics Data System (ADS)
Lopes Cardozo, David; Holdsworth, Peter C. W.
2016-04-01
The magnetization probability density in d = 2 and 3 dimensional Ising models in slab geometry of volume L\\paralleld-1× {{L}\\bot} is computed through Monte-Carlo simulation at the critical temperature and zero magnetic field. The finite-size scaling of this distribution and its dependence on the system aspect-ratio ρ =\\frac{{{L}\\bot}}{{{L}\\parallel}} and boundary conditions are discussed. In the limiting case ρ \\to 0 of a macroscopically large slab ({{L}\\parallel}\\gg {{L}\\bot} ) the distribution is found to scale as a Gaussian function for all tested system sizes and boundary conditions.
High energy efficient solid state laser sources. [slab geometry laser sources
NASA Technical Reports Server (NTRS)
Byer, R. L.
1983-01-01
Slab glass performance studies demonstate 18 J of output at 2 Hz with 2.3% wall plug efficiency. The goal is to achieve 10 J per pulse at 10 Hz and 3% wall plug efficiency during the next annual period. The slab concept was extended to Nd:YAG and to Nd:GGG. To date over 80 W of CW output power at 2% efficiency was generated in slab Nd:YAG. A multiplexed slab Nd:YAG pre-amplifier was invented and a Nd:YAG oscillator was demonstrated with 100kHz linewidth for eventual use in wind velocity measurements.
Quality of Slab Track Construction - Track Alignment Design and Track Geometry
NASA Astrophysics Data System (ADS)
Šestáková, Janka
2015-05-01
The slab track superstructure design (without ballast) is a perspective construction especially for building tunnels and bridges in the modernized sections of railway tracks in Slovakia. Monitoring of the structure described in this article is focused on the transition areas between standard structure with ballast and slab track construction.
Net Rotation of the Lithosphere Induced by Slabs, Plate Geometry, and Keels
NASA Astrophysics Data System (ADS)
Gérault, M.; Becker, T. W.; Kaus, B. J.; Faccenna, C.; Moresi, L. N.
2010-12-01
Most absolute plate motion reconstructions that infer crustal velocities with respect to a deep mantle reference frame (e.g. hotspots) contain a net rotation of the lithosphere around an Euler pole located in the Southwest Indian Ocean. This net rotation, however, differs in the predicted amplitudes by a factor of ~4 between the different reference frames. Global numerical models in three dimensions (3-D) have been able to reproduce the correct direction of net rotation, but only with amplitudes toward the low end of those predicted by plate motion reconstructions in hotspots reference frames. Several recent studies have suggested that geologic features of the Pacific Ocean domain may have a dominant influence on trench and global plate motions. Here, we use a 2-D cylindrical unstructured finite element model to study the influence of subducting slabs, oceanic ridge position, and continental keels on the net rotation of the lithosphere and trench migration. The role of viscosity is also examined, in particular the effect of the asthenosphere and slab viscosity. The fastest predicted net rotation is as large as the values predicted by plate motion reconstructions in the Pacific hotspots reference frames (e.g. HS-3). Such rapid net rotations are induced primarily by asymmetric slab dip angles, in particular a flat slab beneath the Andes and a steep slab in the western Pacific. Asymmetric ridge position can also promote the net rotation of the lithosphere, though to a smaller extent. The effect of the surrounding continental keels is less clear, unlike what has been found in earlier global studies. Keels can either promote or decrease the net rotation of the lithosphere depending on their width and the distance between them and the slab hinge. Several models yield slab advance in the western Pacific, and simultaneously, slab retreat in the eastern Pacific. While these results have been obtained using models with free-slip boundary conditions, a few studies have
The Geometry of the Subducting Slabs Beneath the PRVI Microplate Based on 3D Tomography
NASA Astrophysics Data System (ADS)
Xu, X.; Keller, G. R.; Holland, A. A.; Keranen, K. M.; Li, H.
2011-12-01
The Puerto Rico and the Virgin Islands (PRVI) microplate is located between two subduction zones, with the Puerto Rico trench to the north and the Muertos trough to the south. The Puerto Rico trench is caused by southward-directed subduction of the North American Plate, and the Muertos trough is the northern boundary of the Caribbean Plate. There is no active volcanism on Puerto Rico; however, earthquake depths and seismic tomography imply that the slab of Caribbean plate continues northward beneath Puerto Rico. Puerto Rico overlies these two slabs with extension to both the west (Mona Passage) and southeast (Anacapa Passage). The cause of the extension is unknown, but GPS measurements show that Puerto Rico is experiencing anti-clockwise rotation, and the extension associated with the Anacapa Passage may be produced by this rotation. To the west, it is debated whether the Mona Passage is a boundary between two micro-plates or simple a local rift basin. To address the sources of the extension and the cause of the rotation, we are investigating if the deep structures can be the dynamic source for the observed kinematic movements. We collected data on earthquakes occurring between 2009-2011 in the PRVI region and relocated them using the SEISAN code provided by the Institute of Solid Earth Physics, University of Bergen. The FMTOMO code from Australian National University was used for 3D tomography from P and S wave arrival times. By comparing the relocated epicenters and the 3D tomography results, the subducting slabs were identified. When integrated with the results of previous studies, the geometric model of the slabs is a critical key to understanding the evolution of the PRVI microplate in the past and the future.
NASA Astrophysics Data System (ADS)
Kincaid, C. R.; MacDougall, J. G.; Druken, K. A.; Fischer, K. M.
2010-12-01
Understanding patterns in plate scale mantle flow in subduction zones is key to models of thermal structure, dehydration reactions, volatile distributions and magma generation and transport in convergent margins. Different patterns of flow in the mantle wedge can generate distinct signatures in seismological observables. Observed shear wave fast polarization directions in several subduction zones are inconsistent with predictions of simple 2-D wedge corner flow. Geochemical signatures in a number of subduction zones also indicate 3-D flow and entrainment patterns in the wedge. We report on a series of laboratory experiments on subduction driven flow to characterize spatial and temporal variability in 3-D patterns in flow and shear-induced finite strain. Cases focus on how rollback subduction, along-strike dip changes in subducting plates and evolving gaps or tears in subduction zones control temporal-spatial patterns in 3-D wedge flow. Models utilize a glucose working fluid with a temperature dependent viscosity to represent the upper 2000 km of the mantle. Subducting lithosphere is modeled with two rubber-reinforced continuous belts. Belts pass around trench and upper/lower mantle rollers. The deeper rollers can move laterally to allow for time varying dip angle. Each belt has independent speed control and dip adjustment, allowing for along-strike changes in convergence rate and the evolution of slab gaps. Rollback is modeled using a translation system to produce either uniform and asymmetric lateral trench motion. Neutral density finite strain markers are distributed throughout the fluid and used as proxies for tracking the evolution of anisotropy through space and time in the evolving flow fields. Particle image velocimetry methods are also used to track time varying 3-D velocity fields for directly calculating anisotropy patterns. Results show that complex plate motions (rollback, steepening) and morphologies (gaps) in convergent margins produce flows with
NASA Astrophysics Data System (ADS)
Pavlis, G. L.; Wang, Y.
2015-12-01
A significant number of P and S wave tomography models have been produced in the past decade using various subsets of data from the Earthscope USArray and different inversion algorithms. We focus here on published tomography results that span large portions of the final footprint of the USArray. We use 3D visualization techniques to search for common features in different tomography models. We also compare tomography results to features seen in our current generation wavefield images. Recent innovations of our plane wave migration method have yielded what is arguably the highest resolution image ever produced of the mantle in the vicinity of the transition zone. The new results reveal a rich collection of coherent, dipping structures seen throughout the upper mantle and transition zone. These dipping interfaces are judged significant according to a coherence metric. We treat these surfaces as strain markers to assess proposed models for geometry of the 3D geometry of the Farallon Slab under North America. We find the following geologic interpretations are well supported by independent results: 1. The old Farallon under eastern North America and below the base of transition zone is universally seen as a high velocity anomaly. 2. All results support a simple, 3D kinematic model of the updip limit of the Farallon slab window that follows a track from Cape Mendocino, across Nevada, and northern Arizona and New Mexico. 3. All models show a strong low-velocity mantle under the southwestern U.S. 4. A low-velocity features is universally seen related to the Yellowstone-Snake River system. Shorter wavelength features observed in different tomography models are inconsistent showing that the theme of this session is very important to understand what features are in current results are real. Isopach maps of the thickness of the transition show a systematic difference in transition zone thickness in the western and eastern US. The transition zone thickens in the eastern US in
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Févotte, François; Lathuilière, Bruno; Plagne, Laurent
2014-06-01
The past few years have been marked by a noticeable increase in the interest in 3D whole-core heterogeneous deterministic neutron transport solvers for reference calculations. Due to the extremely large problem sizes tackled by such solvers, they need to use adapted numerical methods and need to be efficiently implemented to take advantage of the full computing power of modern systems. As for numerical methods, one possible approach consists in iterating over resolutions of 2D and 1D MOC problems by taking advantage of prismatic geometries. The MICADO solver, developed at EDF R&D, is a parallel implementation of such a method in distributed and shared memory systems. However it is currently unable to use SIMD vectorization to leverage the full computing power of modern CPUs. In this paper, we describe our first effort to support vectorization in MICADO, typically targeting Intel© SSE CPUs. Both the 2D and 1D algorithms are vectorized, allowing for high expected speedups for the whole spatial solver. We present benchmark computations, which show nearly optimal speedups for our vectorized implementation on the TAKEDA case.
NASA Astrophysics Data System (ADS)
Kuraptsev, A. S.; Sokolov, I. M.; Fofanov, Ya. A.
2016-03-01
We analyze resonant light scattering from a dense and disordered ensemble of motionless point-like scatters with uniform (on average) spatial distribution of the density. The average interatomic distance is considered comparable with the resonant wavelength and the mean free path of photon. The inhomogeneity of dipole-dipole interaction near the surface is discussed. Angular distribution of the light scattered from a medium in a slab geometry is calculated. The total reflected light power P and the reflectivity R depending on the optical thickness of a medium bz are analyzed. It is shown that in the case of small optical thickness (bz ≪ 1) P(bz) ∝ bz2 and R(bz) ∝ bz. With further increasing of the optical thickness we observe oscillations of the dependencies P(bz) and R(bz). The attenuation coefficient of these oscillations is discussed.
NASA Astrophysics Data System (ADS)
Albrecht, G. F.; Eggleston, J. M.; Petr, R. A.
1986-01-01
After a brief overview of slab geometry work published to date, zig-zag optical path slab laser development toward a 100-Hz slab YAG system ('Centurion') and a 10-Hz dual slab glass system ('Gemini') is described. Some of the major diagnostic work performed and in the process is also described to illustrate the major mechanisms responsible for beam distortion observed in these types of slab lasers.
Gyrokinetic Vlasov-Poisson simulation in slab geometry using the conservative IDO scheme
NASA Astrophysics Data System (ADS)
Imadera, Kenji; Kishimoto, Yasuaki; Li, Jiquan; Saito, Daisuke; Utsumi, Takayuki
2008-11-01
We have introduced the IDO-CF (Conservative Form of Interpolated Differential Operator) scheme [1], which is one of the multi-moment schemes and has been applied to various CFD problems, in solving a Vlasov-Poisson system. The IDO scheme is found to be efficient in capturing a sharp domain interface like shock propagation, and in introducing dissipations like particle collision and also external source/sink terms. Furthermore, the IDO-CF scheme has exact mass conservation properties, so that we can apply it to the problems that need long time scale simulations. We first apply the scheme in studying the nonlinear Landau damping and two-stream instability. We have investigated the conservation property of the total mass, energy and entropy, and found that the IDO-CF scheme allows stable simulation over many bounce periods keeping higher accuracy than other multi-moment schemes. We have also developed a gyrokinetic full-f Vlasov code with the IDO-CF scheme in studying the slab ITG driven turbulence. [1] Y.Imai et al., J. Comput. Phys. 227, 2263(2008).
Stability of a two-volume MRxMHD model in slab geometry
NASA Astrophysics Data System (ADS)
Tuen, Li Huey; Hole, Matthew; Dewar, Robert; Dennis, Graham
2015-11-01
Ideal MHD models are inadequate to describe various physical attributes of toroidal plasmas with no continuous symmetry, such as magnetic islands and stochastic regions. A new MHD formulation, MRxMHD, uses only a finite number of ideal-MHD flux surface interfaces, with relaxed plasma regions in between, thus allowing a stepped-pressure profile, magnetic islands, and stochastic regions. In toroidally asymmetric plasma, the existence of interfaces in MRxMHD is contingent on the field pitch, or rotational transform, of flux surfaces being irrational; a KAM argument shows that some good equilibrium flux surfaces continue to exist for small perturbations to an integrable system (foliated by flux surfaces), provided that the rotational transforms on both sides of each interface are sufficiently irrational. Building upon the MRxMHD stability model by Hole, we study the effects of irrationality of the rotational transform at interfaces in MRxMHD on plasma stability. Investigating the plasma stability of a two-volume MRxMHD periodic-slab, we find that the 2D system stability conditions are dependent on the interface and resonance layer magnetic field pitch at minimised energy states. While tearing instabilities exist at low order rational resonances, investigating instability of high-order rationals requires study of pressure-driven instabilities.
Li, Daijun; Ma, Zhe; Haas, Rüdiger; Schell, Alexander; Simon, Janosch; Diart, Robert; Shi, Peng; Hu, Peixin; Loosen, Peter; Du, Keming
2007-05-15
We demonstrate a diode-pumped electro-optical Q-switched slab laser with a high optical efficiency, high pulse energy, and short pulse width with two Nd:YLF crystals inside one resonator. The single compact slab resonator can generate a 1D top-hat beam at both the far field and the near field. With a slab-geometry-design lithium triborate (LBO) crystal, efficient critical phase-matching second-harmonic generation for a 1D top-hat beam with multiple transverse modes is achieved. PMID:17440558
Sjoden, G.E.
1992-03-01
A new discrete ordinates spatial quadrature scheme is presented for solving neutral particle transport problems. This new scheme, called the exponential characteristic method, is developed here in slab geometry with isotropic scattering. This method uses a characteristic integration of the Boltzmann transport equation with an exponential function as the assumed from of the source distribution, continuous across each spatial cell. The exponential source function is constructed to globally conserve zeroth and first spatial source moments and is non-negative. Characteristic integration ensures non-negative fluxes and flux moments. Numerical testing indicates that convergence of the exponential characteristic scheme is fourth order in the limit of vanishingly thin cells. Highly accurate solutions to optically thick problems can result using this scheme with very coarse meshes. Comparing accuracy and computational cost with existing spatial quadrature schemes (diamond difference, linear discontinuous, linear characteristic, linear adaptive, etc.), the exponential characteristic scheme typically performed best. This scheme is expected to be expandable to two dimensions in a straight forward manner. Due to the high accuracies achievable using coarse meshes, this scheme may allow researchers to obtain solutions to transport problems once thought too large or too difficult to be adequately solved conventional computer systems.
Menezes, W. A.; Filho, H. A.; Barros, R. C.
2013-07-01
A generalization of the spectral Green's function (SGF) method is developed for multigroup, fixed-source, slab-geometry discrete ordinates (S{sub N}) problems with anisotropic scattering. The offered SGF method with the one-node block inversion (NBI) iterative scheme converges numerical solutions that are completely free from spatial truncation errors for multigroup slab-geometry S{sub N} problems with scattering anisotropy of order L, provided L < N. As a coarse-mesh numerical method, the SGF method generates numerical solutions that generally do not give detailed information on the problem solution profile, as the grid points can be located considerably away from each other. Therefore, presented here is a technique for the spatial reconstruction of the coarse-mesh solution generated by the multigroup SGF method. Numerical results are given to illustrate the method's accuracy. (authors)
Schnack, D. D.; Cheng, J.; Parker, S. E.; Barnes, D. C.
2013-06-15
We perform linear stability studies of the ion temperature gradient (ITG) instability in unsheared slab geometry using kinetic and extended magnetohydrodynamics (MHD) models, in the regime k{sub ∥}/k{sub ⊥}≪1. The ITG is a parallel (to B) sound wave that may be destabilized by finite ion Larmor radius (FLR) effects in the presence of a gradient in the equilibrium ion temperature. The ITG is stable in both ideal and resistive MHD; for a given temperature scale length L{sub Ti0}, instability requires that either k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} be sufficiently large. Kinetic models capture FLR effects to all orders in either parameter. In the extended MHD model, these effects are captured only to lowest order by means of the Braginskii ion gyro-viscous stress tensor and the ion diamagnetic heat flux. We present the linear electrostatic dispersion relations for the ITG for both kinetic Vlasov and extended MHD (two-fluid) models in the local approximation. In the low frequency fluid regime, these reduce to the same cubic equation for the complex eigenvalue ω=ω{sub r}+iγ. An explicit solution is derived for the growth rate and real frequency in this regime. These are found to depend on a single non-dimensional parameter. We also compute the eigenvalues and the eigenfunctions with the extended MHD code NIMROD, and a hybrid kinetic δf code that assumes six-dimensional Vlasov ions and isothermal fluid electrons, as functions of k{sub ⊥}ρ{sub i} and ρ{sub i}/L{sub Ti0} using a spatially dependent equilibrium. These solutions are compared with each other, and with the predictions of the local kinetic and fluid dispersion relations. Kinetic and fluid calculations agree well at and near the marginal stability point, but diverge as k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} increases. There is good qualitative agreement between the models for the shape of the unstable global eigenfunction for L{sub Ti0}/ρ{sub i}=30 and 20. The results quantify how far
Cisneros, Anselmo T.; Ilas, Dan
2013-01-01
The Advanced High-Temperature Reactor (AHTR) is a 3400 MWth fluoride-salt-cooled high-temperature reactor (FHR) that uses TRISO particle fuel compacted into slabs rather than spherical or cylindrical fuel compacts. Simplified methods are required for parametric design studies such that analyzing the entire feasible design space for an AHTR is tractable. These simplifications include fuel homogenization techniques to increase the speed of neutron transport calculations in depletion analysis and equilibrium depletion analysis methods to analyze systems with multi-batch fuel management schemes. This paper presents three elements of significant novelty. First, the Reactivity-Equivalent Physical Transformation (RPT) methodology usually applied in systems with coated-particle fuel in cylindrical and spherical geometries has been extended to slab geometries. Secondly, based on this newly developed RPT method for slab geometries, a methodology that uses Monte Carlo depletion approaches was further developed to search for the maximum discharge burnup in a multi-batch system by iteratively estimating the beginning of equilibrium cycle (BOEC) composition and sampling different discharge burnups. This Iterative Equilibrium Depletion Search (IEDS) method fully defines an equilibrium fuel cycle (keff, power, flux, and composition evolutions) but is computationally demanding, although feasible on single-processor workstations. Finally, an analytical method, the Non-Linear Reactivity Model, was developed by expanding the linear reactivity model to include an arbitrary number of higher order terms so that single-batch depletion results could be extrapolated to estimate the maximum discharge burnup and BOEC keff in systems with multi-batch fuel management schemes. Results from this method were benchmarked against equilibrium depletion analysis results using the IEDS.
NASA Technical Reports Server (NTRS)
Egan, Michael P.; Leung, Chun Ming; Spagna, George F., Jr.
1988-01-01
The program solves the radiation transport problem in a dusty medium with one-dimensional planar, spherical or cylindrical geometry. It determines self-consistently the effects of multiple scattering, absorption, and re-emission of photons on the temperature of dust grains and the characteristics of the internal radiation field. The program can treat radiation field anisotropy, linear anisotropic scattering, and multi-grain components. The program output consists of the dust-temperature distribution, flux spectrum, surface brightness at each frequency and the observed intensities (involving a convolution with a telescope beam pattern).
Lloyd, R.C.
1986-12-01
A series of critical experiments was completed with mixed plutonium-uranium solutions having Pu/(Pu + U) ratios of approximately 0.5. These experiments were a part of the Criticality Data Development Program between the United States Department of Energy (USDOE), and the Power Reactor and Nuclear Fuel Development Corporation (PNC) of Japan. A complete description of, and data from, the experiments are included in this report. The experiments were performed with mixed plutonium-uranium solutions in cylindrical and slab geometries and included measurements with a water reflector, a concrete reflector, and without an added reflector. The concentration was varied from 112 to 332 g (Pu + U)/liter. The ratio of plutonium to total heavy metal (plutonium plus uranium) was 52% for all experiments.
Lloyd, RC
1988-04-01
A series of critical experiments was completed with mixed plutonium-uranium solutions having Pu/(Pu + U) ratios of approximately 0.4. These experiments were a part of the Criticality Data Development Program between the United States Department of Energy (USDOE), and the Power Reactor and Nuclear Fuel Development Corporation (PNC) of Japan. A complete description of, and data from, the experiments are included in this report. The experiments were performed with mixed plutonium-uranium solutions in cylinqrical and slab geometries and included measurements with a water reflector, a concrete reflector, and without an added reflector. The concentration was varied from 105 to 436 g (Pu + U)/liter. The ratio of plutonium to total heavy metal (plutonium plus uranium) was 0.4 for all experiments.
NASA Astrophysics Data System (ADS)
Maginot, Peter G.; Morel, Jim E.; Ragusa, Jean C.
2012-08-01
We present a new nonlinear spatial finite-element method for the linearized Boltzmann transport equation with Sn angular discretization in 1-D and 2-D Cartesian geometries. This method has two central characteristics. First, it is equivalent to the linear-discontinuous (LD) Galerkin method whenever that method yields a strictly non-negative solution. Second, it always satisfies both the zeroth and first spatial moment equations. Because it yields the LD solution when that solution is non-negative, one might interpret our method as a classical fix-up to the LD scheme. However, fix-up schemes for the LD equations derived in the past have given up solution of the first moment equations when the LD solution is negative in order to satisfy positivity in a simple manner. We present computational results comparing our method in 1-D to the strictly non-negative linear exponential-discontinuous method and to the LD method. We present computational results in 2-D comparing our method to a recently developed LD fix-up scheme and to the LD scheme. It is demonstrated that our method is a valuable alternative to existing methods.
Pohl, B.A.; Keeton, S.C.
1997-09-01
R. C. Lloyd of PNL has completed and published a series of critical experiments with mixed plutonium- uranium nitrate solutions (Reference 1). This series of critical experiments was part of an extensive program jointly sponsored by the U. S. Department of Energy (DOE) and the Power Reactor and Nuclear Fuel Development Corporation (PNC) of Japan and was carried out in the mid-1980`s. The experiments evaluated here (published as Report PNL-6327) were performed with mixed plutonium- uranium nitrate solution in a variable thickness slab tank with two 106.7 cm square sides and a width that could be varied from 7.6 to 22.8 cm. The objective of these experiments was to obtain experimental data to permit the validation of computer codes for criticality calculations and of cross-section data to minimize the uncertainties inherent therein, so that facility safety, efficiency, and reliability could be enhanced. The concentrations of the solution were about 105, 293, and 435 g(Pu+U)/liter with a ratio of plutonium to total heavy metal (plutonium plus uranium) of about 0. 40 for all eight experiments. Four measurements were made with a water reflector, and four with no reflector. Following the publication of the initial PNL reports, considerable effort was devoted to an extensive reevaluation of this series of experiments by a collaboration of researchers from ORNL, PNL, and PNC (Reference 2). Their work resulted in a more accurate description of the ``as built`` hardware configuration and the materials specifications. For the evaluations in this report, the data published in Reference 2 by Smolen et al. is selected to supersede the original PNL report. Eight experiments have been evaluated and seven (063, 064, 071, 072, 074, 075, and 076) provide benchmark criticality data. Experiment 073 could not achieve criticality within vessel height limitations.
NASA Astrophysics Data System (ADS)
Shea, Thomas; Krimer, Daniel; Costa, Fidel; Hammer, Julia
2014-05-01
One of the achievements in recent years in volcanology is the determination of time-scales of magmatic processes via diffusion in minerals and its addition to the petrologists' and volcanologists' toolbox. The method typically requires one-dimensional modeling of randomly cut crystals from two-dimensional thin sections. Here we address the question whether using 1D (traverse) or 2D (surface) datasets exploited from randomly cut 3D crystals introduces a bias or dispersion in the time-scales estimated, and how this error can be improved or eliminated. Computational simulations were performed using a concentration-dependent, finite-difference solution to the diffusion equation in 3D. The starting numerical models involved simple geometries (spheres, parallelepipeds), Mg/Fe zoning patterns (either normal or reverse), and isotropic diffusion coefficients. Subsequent models progressively incorporated more complexity, 3D olivines possessing representative polyhedral morphologies, diffusion anisotropy along the different crystallographic axes, and more intricate core-rim zoning patterns. Sections and profiles used to compare 1, 2 and 3D diffusion models were selected to be (1) parallel to the crystal axes, (2) randomly oriented but passing through the olivine center, or (3) randomly oriented and sectioned. Results show that time-scales estimated on randomly cut traverses (1D) or surfaces (2D) can be widely distributed around the actual durations of 3D diffusion (~0.2 to 10 times the true diffusion time). The magnitude over- or underestimations of duration are a complex combination of the geometry of the crystal, the zoning pattern, the orientation of the cuts with respect to the crystallographic axes, and the degree of diffusion anisotropy. Errors on estimated time-scales retrieved from such models may thus be significant. Drastic reductions in the uncertainty of calculated diffusion times can be obtained by following some simple guidelines during the course of data
Fevotte, F.; Lathuiliere, B.
2013-07-01
The large increase in computing power over the past few years now makes it possible to consider developing 3D full-core heterogeneous deterministic neutron transport solvers for reference calculations. Among all approaches presented in the literature, the method first introduced in [1] seems very promising. It consists in iterating over resolutions of 2D and ID MOC problems by taking advantage of prismatic geometries without introducing approximations of a low order operator such as diffusion. However, before developing a solver with all industrial options at EDF, several points needed to be clarified. In this work, we first prove the convergence of this iterative process, under some assumptions. We then present our high-performance, parallel implementation of this algorithm in the MICADO solver. Benchmarking the solver against the Takeda case shows that the 2D-1D coupling algorithm does not seem to affect the spatial convergence order of the MOC solver. As for performance issues, our study shows that even though the data distribution is suited to the 2D solver part, the efficiency of the ID part is sufficient to ensure a good parallel efficiency of the global algorithm. After this study, the main remaining difficulty implementation-wise is about the memory requirement of a vector used for initialization. An efficient acceleration operator will also need to be developed. (authors)
Benchmarks and models for 1-D radiation transport in stochastic participating media
Miller, D S
2000-08-21
Benchmark calculations for radiation transport coupled to a material temperature equation in a 1-D slab and 1-D spherical geometry binary random media are presented. The mixing statistics are taken to be homogeneous Markov statistics in the 1-D slab but only approximately Markov statistics in the 1-D sphere. The material chunk sizes are described by Poisson distribution functions. The material opacities are first taken to be constant and then allowed to vary as a strong function of material temperature. Benchmark values and variances for time evolution of the ensemble average of material temperature energy density and radiation transmission are computed via a Monte Carlo type method. These benchmarks are used as a basis for comparison with three other approximate methods of solution. One of these approximate methods is simple atomic mix. The second approximate model is an adaptation of what is commonly called the Levermore-Pomraning model and which is referred to here as the standard model. It is shown that recasting the temperature coupling as a type of effective scattering can be useful in formulating the third approximate model, an adaptation of a model due to Su and Pomraning which attempts to account for the effects of scattering in a stochastic context. This last adaptation shows consistent improvement over both the atomic mix and standard models when used in the 1-D slab geometry but shows limited improvement in the 1-D spherical geometry. Benchmark values are also computed for radiation transmission from the 1-D sphere without material heating present. This is to evaluate the performance of the standard model on this geometry--something which has never been done before. All of the various tests demonstrate the importance of stochastic structure on the solution. Also demonstrated are the range of usefulness and limitations of a simple atomic mix formulation.
NASA Astrophysics Data System (ADS)
Vyawahare, Vishwesh A.; Nataraj, P. S. V.
2013-07-01
In this paper, we report the development and analysis of some novel versions and approximations of the fractional-order (FO) point reactor kinetics model for a nuclear reactor with slab geometry. A systematic development of the FO Inhour equation, Inverse FO point reactor kinetics model, and fractional-order versions of the constant delayed neutron rate approximation model and prompt jump approximation model is presented for the first time (for both one delayed group and six delayed groups). These models evolve from the FO point reactor kinetics model, which has been derived from the FO Neutron Telegraph Equation for the neutron transport considering the subdiffusive neutron transport. Various observations and the analysis results are reported and the corresponding justifications are addressed using the subdiffusive framework for the neutron transport. The FO Inhour equation is found out to be a pseudo-polynomial with its degree depending on the order of the fractional derivative in the FO model. The inverse FO point reactor kinetics model is derived and used to find the reactivity variation required to achieve exponential and sinusoidal power variation in the core. The situation of sudden insertion of negative reactivity is analyzed using the FO constant delayed neutron rate approximation. Use of FO model for representing the prompt jump in reactor power is advocated on the basis of subdiffusion. Comparison with the respective integer-order models is carried out for the practical data. Also, it has been shown analytically that integer-order models are a special case of FO models when the order of time-derivative is one. Development of these FO models plays a crucial role in reactor theory and operation as it is the first step towards achieving the FO control-oriented model for a nuclear reactor. The results presented here form an important step in the efforts to establish a step-by-step and systematic theory for the FO modeling of a nuclear reactor.
Development of a hybrid deterministic/stochastic method for 1D nuclear reactor kinetics
NASA Astrophysics Data System (ADS)
Terlizzi, Stefano; Rahnema, Farzad; Zhang, Dingkang; Dulla, Sandra; Ravetto, Piero
2015-12-01
A new method has been implemented for solving the time-dependent neutron transport equation efficiently and accurately. This is accomplished by coupling the hybrid stochastic-deterministic steady-state coarse-mesh radiation transport (COMET) method [1,2] with the new predictor-corrector quasi-static method (PCQM) developed at Politecnico di Torino [3]. In this paper, the coupled method is implemented and tested in 1D slab geometry.
Development of a hybrid deterministic/stochastic method for 1D nuclear reactor kinetics
Terlizzi, Stefano; Dulla, Sandra; Ravetto, Piero; Rahnema, Farzad; Zhang, Dingkang
2015-12-31
A new method has been implemented for solving the time-dependent neutron transport equation efficiently and accurately. This is accomplished by coupling the hybrid stochastic-deterministic steady-state coarse-mesh radiation transport (COMET) method [1,2] with the new predictor-corrector quasi-static method (PCQM) developed at Politecnico di Torino [3]. In this paper, the coupled method is implemented and tested in 1D slab geometry.
Thorkelson, D.J.; Taylor, R.P. )
1989-09-01
The geometry and geologic implications of subducted spreading ridges are topics that have bedeviled earth scientists ever since the recognition of plate tectonics. As a consequence of subduction of the Kula-Farallon and East Pacific rises, slab windows formed and migrated beneath the North American Cordillera. The probable shape and extent of these windows, which represent the asthenosphere-filled gaps between two separating, subducting oceanic plates, are depicted from the Late Cretaceous to the present. Possible effects of the existence and migration of slab windows on the Cordillera at various times include cessation of arc volcanism and replacement by rift or plate-edge volcanism; lithospheric uplift, attenuation, and extension; and increased intensity of compressional tectonism. Eocene extensional tectonism and alkaline magmatism in southern British Columbia and the northwestern United States were facilitated by slab-window development.
Goedecke, Niels; Maul, Christof; Kauczok, Sebastian; Gericke, Karl-Heinz; Chichinin, Alexey I.
2009-08-07
The bimolecular reaction O({sup 1}D)+N{sub 2}O{yields}NO+NO was photoinitiated in the (N{sub 2}O){sub 2} dimer at a wavelength of 193 nm and was investigated by three-dimensional (3D) velocity map imaging. State selective 3D momentum vector distributions were monitored and analyzed. For the first time, kinetic energy resolution and stereodynamic information about the reaction under constrained geometry conditions is available. Directly observable NO products exhibit moderate vibrational excitation and are rotationally and translationally cold. Speed and spatial distributions suggest a pronounced backward scattering of the observed products with respect to the direction of motion of the O({sup 1}D) atom. Forward scattered partner products, which are not directly detectable are also translationally cold, but carry very large internal energy as vibration or rotation. The results confirm and extend previous studies on the complex initiated reaction system. The restricted geometry of the van der Waals complex seems to favor an abstraction reaction of the terminal nitrogen atom by the O({sup 1}D) atom, which is in striking contrast to the behavior observed for the unrestricted gas phase reaction under bulk conditions.
NASA Astrophysics Data System (ADS)
Manassah, Jamal T.
2016-08-01
Using the expansion in the eigenmodes of 1-D Lienard-Wiechert kernel, the temporal and spectral profiles of the radiation emitted by a fully inverted collection of two-level atoms in a sub-wavelength slab geometry are computed. The initial number of amplifying modes determine the specific regime of radiation. In particular, the temporal profile of the field intensity is oscillatory and the spectral profile is non-Lorentzian with two unequal height peaks in a narrow band centered at the slab thickness value at which the real parts of the lowest order odd and even eigenvalues are equal.
NASA Technical Reports Server (NTRS)
Spurrier, Francis R. (Inventor); DeZubay, Egon A. (Inventor); Murray, Alexander P. (Inventor); Vidt, Edward J. (Inventor)
1984-01-01
Slab-shaped high efficiency catalytic reformer configurations particularly useful for generation of fuels to be used in fuel cell based generation systems. A plurality of structures forming a generally rectangular peripheral envelope are spaced about one another to form annular regions, an interior annular region containing a catalytic bed and being regeneratively heated on one side by a hot comubstion gas and on the other side by the gaseous products of the reformation. An integrally mounted combustor is cooled by impingement of incoming oxidant.
Spurrier, Francis R.; DeZubay, Egon A.; Murray, Alexander P.; Vidt, Edward J.
1984-02-07
Slab-shaped high efficiency catalytic reformer configurations particularly useful for generation of fuels to be used in fuel cell based generation systems. A plurality of structures forming a generally rectangular peripheral envelope are spaced about one another to form annular regions, an interior annular region containing a catalytic bed and being regeneratively heated on one side by a hot comubstion gas and on the other side by the gaseous products of the reformation. An integrally mounted combustor is cooled by impingement of incoming oxidant.
Spurrier, Francis R.; DeZubay, Egon A.; Murray, Alexander P.; Vidt, Edward J.
1985-03-12
Slab-shaped high efficiency catalytic reformer configurations particularly useful for generation of fuels to be used in fuel cell based generation systems. A plurality of structures forming a generally rectangular peripheral envelope are spaced about one another to form annular regions, an interior annular region containing a catalytic bed and being regeneratively heated on one side by a hot combustion gas and on the other side by the gaseous products of the reformation. An integrally mounted combustor is cooled by impingement of incoming oxidant.
NASA Astrophysics Data System (ADS)
Santos, Frederico P.; Filho, Hermes Alves; Barros, Ricardo C.
2013-10-01
The scattering source iterative (SI) scheme is traditionally applied to converge fine-mesh numerical solutions to fixed-source discrete ordinates (SN) neutron transport problems. The SI scheme is very simple to implement under a computational viewpoint. However, the SI scheme may show very slow convergence rate, mainly for diffusive media (low absorption) with several mean free paths in extent (low leakage). In this work we describe an acceleration technique based on an improved initial guess for the scattering source distribution within the slab. In other words, we use as initial guess for the fine-mesh scattering source, the coarse-mesh solution of the neutron diffusion equation with special boundary conditions to account for the classical SN prescribed boundary conditions, including vacuum boundary conditions. Therefore, we first implement a spectral nodal method that generates coarse-mesh diffusion solution that is completely free from spatial truncation errors, then we reconstruct this coarse-mesh solution within each spatial cell of the discretization grid, to further yield the initial guess for the fine-mesh scattering source in the first SN transport sweep (forward and backward) across the spatial grid. We consider a number of numerical experiments to illustrate the efficiency of the offered diffusion synthetic acceleration (DSA) technique.
Hauck, Cory D; Alldredge, Graham; Tits, Andre
2012-01-01
We present a numerical algorithm to implement entropy-based (M{sub N}) moment models in the context of a simple, linear kinetic equation for particles moving through a material slab. The closure for these models - as is the case for all entropy-based models - is derived through the solution of constrained, convex optimization problem. The algorithm has two components. The first component is a discretization of the moment equations which preserves the set of realizable moments, thereby ensuring that the optimization problem has a solution (in exact arithmetic). The discretization is a second-order kinetic scheme which uses MUSCL-type limiting in space and a strong-stability-preserving, Runge-Kutta time integrator. The second component of the algorithm is a Newton-based solver for the dual optimization problem, which uses an adaptive quadrature to evaluate integrals in the dual objective and its derivatives. The accuracy of the numerical solution to the dual problem plays a key role in the time step restriction for the kinetic scheme. We study in detail the difficulties in the dual problem that arise near the boundary of realizable moments, where quadrature formulas are less reliable and the Hessian of the dual objection function is highly ill-conditioned. Extensive numerical experiments are performed to illustrate these difficulties. In cases where the dual problem becomes 'too difficult' to solve numerically, we propose a regularization technique to artificially move moments away from the realizable boundary in a way that still preserves local particle concentrations. We present results of numerical simulations for two challenging test problems in order to quantify the characteristics of the optimization solver and to investigate when and how frequently the regularization is needed.
Exponentially-convergent Monte Carlo for the 1-D transport equation
Peterson, J. R.; Morel, J. E.; Ragusa, J. C.
2013-07-01
We define a new exponentially-convergent Monte Carlo method for solving the one-speed 1-D slab-geometry transport equation. This method is based upon the use of a linear discontinuous finite-element trial space in space and direction to represent the transport solution. A space-direction h-adaptive algorithm is employed to restore exponential convergence after stagnation occurs due to inadequate trial-space resolution. This methods uses jumps in the solution at cell interfaces as an error indicator. Computational results are presented demonstrating the efficacy of the new approach. (authors)
Tunable one-dimensional photonic crystal slabs
NASA Astrophysics Data System (ADS)
Beccherelli, R.; Bellini, B.; Zografopoulos, D.; Kriezis, E.
2007-05-01
A 1D photonic crystal slab based on preferential etching of commercially available silicon-on-insulator wafers is presented. Compared to dry etching, anisotropic wet etching is more tolerant to errors as it is self-stopping on crystallographic {111} planes and it produces a more precise geometry with symmetries and homothetic properties, with surface roughness close to 1 nm. The resulting grooves are infiltrated by low viscosity liquid crystal having large positive optical anisotropy. The use of slanted grooves provides advantages: first of all the complete filling of slanted grooves is simplified when compared to vertical walls structures. Furthermore alignment is significantly facilitated. Indeed the liquid crystal molecules tend to align with their long axis along the submicron grooves. Therefore by forcing reorientation out of a rest position, the liquid crystal presents a choice of refractive indices to the propagating optical field. The liquid crystal behavior is simulated by a finite element method, and coupled to a finite difference time domain method. We investigate different photonic crystal configurations. Large tunability of bandgap edge for TE polarization is demonstrated when switching the liquid crystal with an applied voltage. We have also studied the use of the same device geometry as a very compact microfluidic refractometric sensor.
Improving Seismic Constraints on Subduction Zone Geometry
NASA Astrophysics Data System (ADS)
Syracuse, E. M.; Abers, G. A.; Fischer, K. M.; van Keken, P. E.; Kneller, E. A.; Rychert, C. A.
2007-12-01
Accurate slab geometries are necessary for 3D flow modeling, and for understanding the variations in temperature and melting geometry between different subduction zones. Recent studies have shown that the depth to slab beneath arc volcanoes varies by as much as a factor of two between subduction zones, but these results are based on teleseismic earthquake catalogs with potentially large errors. When available, local seismic arrays provide better constraints. The TUCAN array (Tomography Under Costa Rica and Nicaragua) deployed 48 three component broadband PASSCAL instruments for 18 months with station spacing of 10-50 km across the Central America arc. This dataset provides some of the best control anywhere for ground-truth comparison of teleseismic catalogs in steeply dipping subduction zones. Joint inversion of TUCAN arrival times for velocity and hypocenters illuminate a 10-15 km thick Wadati-Benioff zone (WBZ), with absolute hypocenter uncertainties of 1-5 km. Besides providing accurate hypocenters, the tomographic images provide independent constraints on melting and temperature, through the imaging of low Vp (7.5-7.8 km/s) and highly attenuating (40
NASA Astrophysics Data System (ADS)
Vanacore, E. A.; Lopez, A. M.; Huerfano Moreno, V. A.
2015-12-01
The determination of earthquake locations are dependent on the velocity model selected. Consequently, the refinement and updating of the velocity models used at the local and regional network level is a critical component for network efficiency through location accuracy. With the expansion of broadband instruments within the Puerto Rico -Virgin Islands region, updating the velocity model is a current long term goal of the Puerto Rico Seismic Network (PRSN). As a first step to this long term goal, receiver functions of ~20 broadband stations with data between 2010 and 2015 were calculated using iterative time domain deconvolution. The receiver function analysis not only provides insight into the crustal velocity structure but also leads to a better understanding of the region's larger tectonic structure. Preliminary results of the receiver function analysis exhibit evidence of a "slab signal"; the receiver function backazimuth sweeps for some stations particularly on the northern side of the island contain a strong P to S conversion at approximately 7 seconds which likely corresponds to the top of the slab beneath Puerto Rico. This strong slab signal implies that simple 1-D analyses of the data (e.g. H-K stacking) may lead to misleading results. To further understand the crustal structure of PRVI, we employ a 3D common-conversion-point analysis. This analysis yields a Moho beneath the island between 32-42km and a possible southward dipping slab structure between 60-80km depth. Further analysis is needed to determine the 2D or 3D velocity structure of Puerto Rico and the surrounding environs such as waveform modeling. Given the current geometry of the available array, detailed imaging of the slab and mantle wedge beneath Puerto Rico is limited. PRSN is currently seeking to install denser temporary networks to improve local imaging that will help understand the nature of the crust, mantle wedge and slab structure beneath the island as well as the structure's influence
Linearly exact parallel closures for slab geometry
Ji, Jeong-Young; Held, Eric D.; Jhang, Hogun
2013-08-15
Parallel closures are obtained by solving a linearized kinetic equation with a model collision operator using the Fourier transform method. The closures expressed in wave number space are exact for time-dependent linear problems to within the limits of the model collision operator. In the adiabatic, collisionless limit, an inverse Fourier transform is performed to obtain integral (nonlocal) parallel closures in real space; parallel heat flow and viscosity closures for density, temperature, and flow velocity equations replace Braginskii's parallel closure relations, and parallel flow velocity and heat flow closures for density and temperature equations replace Spitzer's parallel transport relations. It is verified that the closures reproduce the exact linear response function of Hammett and Perkins [Phys. Rev. Lett. 64, 3019 (1990)] for Landau damping given a temperature gradient. In contrast to their approximate closures where the vanishing viscosity coefficient numerically gives an exact response, our closures relate the heat flow and nonvanishing viscosity to temperature and flow velocity (gradients)
Linearly exact parallel closures for slab geometry
NASA Astrophysics Data System (ADS)
Ji, Jeong-Young; Held, Eric D.; Jhang, Hogun
2013-08-01
Parallel closures are obtained by solving a linearized kinetic equation with a model collision operator using the Fourier transform method. The closures expressed in wave number space are exact for time-dependent linear problems to within the limits of the model collision operator. In the adiabatic, collisionless limit, an inverse Fourier transform is performed to obtain integral (nonlocal) parallel closures in real space; parallel heat flow and viscosity closures for density, temperature, and flow velocity equations replace Braginskii's parallel closure relations, and parallel flow velocity and heat flow closures for density and temperature equations replace Spitzer's parallel transport relations. It is verified that the closures reproduce the exact linear response function of Hammett and Perkins [Phys. Rev. Lett. 64, 3019 (1990)] for Landau damping given a temperature gradient. In contrast to their approximate closures where the vanishing viscosity coefficient numerically gives an exact response, our closures relate the heat flow and nonvanishing viscosity to temperature and flow velocity (gradients).
Evolution and diversity of subduction zones controlled by slab width.
Schellart, W P; Freeman, J; Stegman, D R; Moresi, L; May, D
2007-03-15
Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges. PMID:17361181
History vs. snapshot: how slab morphology relates to slab age evolution
NASA Astrophysics Data System (ADS)
Garel, Fanny; Goes, Saskia; Davies, Rhodri; Davies, Huw; Lallemand, Serge; Kramer, Stephan; Wilson, Cian
2016-04-01
The age of the subducting plate at the trench ("slab age") spans a wide range, from less than 10 Myr in Central and South America to 150 Myr in the Marianas. The morphology of subducting slab in the upper mantle is also very variable, from slabs stagnating at the top of the lower mantle to slabs penetrating well beyond 1000 km depth. People have looked rather unsucessfully for correlations between slab morphology and subduction parameters, including age at the trench, on the basic assumption that old (thick) plates are likely to generate a large slab pull force that would influence slab dip. Thermo-mechanical models reveal complex feedbacks between temperature, strain rate and rheology, and are able to reproduce the evolution of plate ages as a function of time, subducting plate velocity and trench velocity. In particular, we show how initially young subducting plates can rapidly age at the surface because of a slow sinking velocity. As a consequence, different slab morphologies can exhibit similar ages at the trench provided that subduction history is different. We illustrate how models provide insights into Earth subduction zones for which we have to consider their history (evolution of trench velocity, relative plate ages at time of initiation) in order to unravel their present-day geometry.
Detecting lower mantle slabs beneath Asia and the Aleutians
NASA Astrophysics Data System (ADS)
Schumacher, L.; Thomas, C.
2016-03-01
To investigate the descend of subducted slabs we search for and analyse seismic arrivals that reflected off the surface of the slab. In order to distinguish between such arrivals and other seismic phases, we search for waves that reach a seismic array with a backazimuth deviating from the theoretical backazimuth of the earthquake. Source-receiver combinations are chosen in a way that their great circle paths do not intersect the slab region, hence the direct arrivals can serve as reference. We focus on the North and Northwest Pacific region by using earthquakes from Japan, the Philippines and the Hindu Kush area recorded at North American networks (e.g. USArray, Alaska and Canada). Using seismic array techniques for analysing the data and record information on slowness, backazimuth and travel time of the observed out-of-plane arrivals we use these measurements to trace the wave back through a 1D velocity model to its scattering/reflection location. We find a number of out-of-plane reflections. Assuming only single scattering, most out-of-plane signals have to travel as P-to-P phases and only a few as S-to-P phases, due to the length of the seismograms we processed. The located reflection points present a view of the 3D structures within the mantle. In the upper mantle and the transition zone they correlate well with the edges of fast velocity regions in tomographic images. We also find reflection points in the mid- and lower mantle and their locations generally agree with fast velocities mapped by seismic tomography models suggesting that in the subduction regions we map, slabs enter the lower mantle. To validate our approach, we calculate and process synthetic seismograms for 3D wave field propagation through a model containing a slab-like heterogeneity. We show, that depending on the source-receiver geometry relative to the reflection plane, it is indeed possible to observe and back-trace out-of-plane signals.
Detecting lower-mantle slabs beneath Asia and the Aleutians
NASA Astrophysics Data System (ADS)
Schumacher, L.; Thomas, C.
2016-06-01
To investigate the descend of subducted slabs we search for and analyse seismic arrivals that reflected off the surface of the slab. In order to distinguish between such arrivals and other seismic phases, we search for waves that reach a seismic array with a backazimuth deviating from the theoretical backazimuth of the earthquake. Source-receiver combinations are chosen in a way that their great circle paths do not intersect the slab region, hence the direct arrivals can serve as reference. We focus on the North and Northwest Pacific region by using earthquakes from Japan, the Philippines and the Hindu Kush area recorded at North American networks (e.g. USArray, Alaska and Canada). Using seismic array techniques for analysing the data and record information on slowness, backazimuth and traveltime of the observed out-of-plane arrivals we use these measurements to trace the wave back through a 1-D velocity model to its scattering/reflection location. We find a number of out-of-plane reflections. Assuming only single scattering, most out-of-plane signals have to travel as P-to-P phases and only a few as S-to-P phases, due to the length of the seismograms we processed. The located reflection points present a view of the 3-D structures within the mantle. In the upper mantle and the transition zone they correlate well with the edges of fast velocity regions in tomographic images. We also find reflection points in the mid- and lower mantle and their locations generally agree with fast velocities mapped by seismic tomography models suggesting that in the subduction regions we map, slabs enter the lower mantle. To validate our approach, we calculate and process synthetic seismograms for 3-D wave field propagation through a model containing a slab-like heterogeneity. We show, that depending on the source-receiver geometry relative to the reflection plane, it is indeed possible to observe and back-trace out-of-plane signals.
The gravitational field of an infinite flat slab
NASA Astrophysics Data System (ADS)
Fulling, S. A.; Bouas, J. D.; Carter, H. B.
2015-08-01
We study Einstein's equations with a localized plane-symmetric source, with close attention to gauge freedom/fixing and to listing all physically distinct solutions. In the vacuum regions there are only two qualitatively different solutions, one curved and one flat; in addition, on each of the two sides there is a free parameter describing how the slab is embedded into the vacuum region. Surprisingly, for a generic slab source the solution must be curved on one side and flat on the other. We treat infinitely thin slabs in full detail and indicate how thick slabs can increase the variety of external geometry pairs. Positive energy density seems to force external geometries with curvature singularities at some distance from the slab; we speculate that such singularities occur in regions where the solution cannot be physically relevant anyway.
Along-strike translation of a fossil slab
NASA Astrophysics Data System (ADS)
Pikser, Jacob E.; Forsyth, Donald W.; Hirth, Greg
2012-05-01
The Isabella anomaly is a high seismic velocity anomaly beneath the southern Central Valley of California that has previously been interpreted to represent a lithospheric drip or delaminated Sierra Nevada root. However, recent work suggests that the anomaly is a remnant slab, left over from Cenozoic subduction, attached to the Pacific plate underneath the edge of the North American plate. This hypothesis requires the slab to translate hundreds of kilometers along-strike while remaining intact and attached to the Pacific plate despite drag from the surrounding asthenosphere and overriding lithosphere. We use 3-D finite element models to simulate this scenario, and calculate the viscosity ratio between the slab and the surrounding asthenosphere required for the slab to translate undeformed. The required viscosity ratio increases with increasing downdip extent of the slab, and decreases with increasing slab dip; for geometries approximating that of our proposed slab, the required ratio ranges from 102 to 104. Given the thermal and petrologic history of the slab, mantle flow laws predict viscosity contrasts greater than or equal to these requirements. As such, we conclude that along-strike translation of a remnant slab is feasible, and serves as a possible explanation of the Isabella anomaly. The significance of this finding extends beyond our general understanding of subduction dynamics, in that the presence of such a slab could have implications for the water budget of the San Andreas Fault and its role in aseismic slip.
1-D EQUILIBRIUM DISCRETE DIFFUSION MONTE CARLO
T. EVANS; ET AL
2000-08-01
We present a new hybrid Monte Carlo method for 1-D equilibrium diffusion problems in which the radiation field coexists with matter in local thermodynamic equilibrium. This method, the Equilibrium Discrete Diffusion Monte Carlo (EqDDMC) method, combines Monte Carlo particles with spatially discrete diffusion solutions. We verify the EqDDMC method with computational results from three slab problems. The EqDDMC method represents an incremental step toward applying this hybrid methodology to non-equilibrium diffusion, where it could be simultaneously coupled to Monte Carlo transport.
Deepest hypocentral distributions associated with stagnant slabs and penetrated slabs
NASA Astrophysics Data System (ADS)
Fukao, Y.; Obayashi, M.
2013-12-01
We constructed a new P-wave tomographic model of the mantle, GAP_P4, using more than ten millions of travel time data, including waveform-based differential travel times from ocean bottoms, to all of which the finite frequency kernels were applied in the inversion. Based on this model, we made a systematic survey for subducted slab images around the Circum Pacific. This survey revealed a progressive lateral variation of slab configuration along arc(s), where a subducted slab is in general in one or two of the following four stages: I. slab stagnant above the 660, II. slab penetrating the 660, III. slab trapped in the uppermost lower mantle (660 to ˜1000 km in depth), and IV. slab descending well into the deep lower mantle. The majority of the slab images are either at stage I or III. We interpret I to IV as the successive stages of slab subduction through the transition region with the 660 at the middle. There is a remarkable correlation of the slab configuration with the deepest shock hypocentral distribution. Subhorizontal distributions of deepest shocks are associated with stagnant slabs in the transition zone (slabs at stage I). Their focal depths are limited to shallower than ˜620 km. Steeply dipping deepest shock distributions are associated with penetrating slabs across the 660-km discontinuity or trapped slabs below it (slabs at stages II and III). Their focal depths extend well beyond ˜620 km. There are no cases of association of either a stagnant slab (at stage I) with subvertical distribution of deepest shocks or a trapped slab (at stage II or III) with their subhorizontal distribution. Only steeply dipping slabs appear to penetrate the 660 to be trapped in the uppermost lower mantle. The along-arc variations of stagnant-slab configuration and deepest shock distribution beneath the Bonin arc indicate a process of how the slab begins to penetrate the 660-km discontinuity after the slab stagnation. Those beneath the Java arc and Kermadec arc commonly
Was there a Laramide "flat slab"?
NASA Astrophysics Data System (ADS)
Jones, C. H.
2014-12-01
Slab-continent interactions drive most non-collisional orogenies; this has led us to usually anticipate that temporal changes or spatial variations in orogenic style are related to changes in the slab, most especially in the slab's dip. This is most dramatically evident for orogenies in the foreland, well away from the trench, such as the Laramide orogeny. However, the physical means of connecting slab geometry to crustal deformation remain obscure. Dickinson and Snyder (1978) and Bird (1984) laid out a conceptually elegant means of creating foreland deformation from shear between a slab and overriding continental lithosphere, but such strong shear removed all of the continental lithosphere in the western U.S. when included in a numerical simulation of flat slab subduction (Bird, 1988), a removal in conflict with observations of volcanic rocks and xenoliths in many locations. Relying on an increase in edge normal stresses results, for the Laramide, in requiring the little-deformed Colorado Plateau to either be unusually strong or to have risen rapidly enough and high enough to balance edge stresses with body forces. Early deformation in the Plateau rules out unusual strength, and the accumulation and preservation of Late Cretaceous near-sea level sedimentary rocks makes profound uplift unlikely (though not impossible). Relying on comparisons with the Sierras Pampeanas is also fraught with problems: the Sierras are not separated from the Andean fold-and-thrust belt by several hundred kilometers of little-deformed crust, nor were they buried under kilometers of marine muds as were large parts of the Laramide foreland. We have instead suggested that some unusual interactions of an obliquely subducting plate with a thick Archean continental root might provide a better explanation than a truly flat slab (Jones et al., 2011). From this, and given that several flat-slab segments today are not associated with foreland orogenesis and noting that direct evidence for truly
Mantle Response to a Slab Gap and Three-dimensional Slab Interaction in Central America
NASA Astrophysics Data System (ADS)
Jadamec, M. A.; Fischer, K. M.
2013-12-01
Seismically constrained global slab geometries suggest the Middle America-South American subduction system contains a gap on the order of 500 km separating the east-dipping Cocos and Nazca slabs at depth (Gudmundsson and Sambridge, 1998; Syracuse and Abers, 2006; Hayes et al., 2012). The location of the gap correlates with tectonic features impinging on the Pacific side of the Middle America trench, in particular the incoming young buoyant oceanic lithosphere and oceanic ridges associated with the Galapagos hotspot and Cocos-Nazca spreading center (Protti et al., 1994; Mann et al., 2007; Muller et al., 2008). Moreover, geochemical studies focusing on the arc chemistry in the Central American volcanic front argue for a slab window of some kind in this region (Johnston and Thorkelson, 1997; Abratis and Worner, 2001; Hoernle et al., 2008). We use high-resolution, three-dimensional (3D) geodynamic modeling of the Middle America-South American subduction system to investigate the role of the incoming young oceanic lithosphere and a gap between the Cocos and Nazca slabs in controlling mantle flow velocity and geochemical signatures beneath Central America. The geodynamic models are geographically referenced with the geometry and thermal structure for the overriding and subducting plates based on geological and geophysical observables and constructed with the multi-plate subduction generator code, SlabGenerator (Jadamec and Billen, 2010; Jadamec et al., 2012; Jadamec and Billen, 2012). The viscous flow simulations are solved using the mantle convection finite-element code, CitcomCU (Zhong, 2006), modified by Jadamec and Billen (2010) to take into account the experimentally derived flow law for olivine and allow for variable 3D plate interface geometries and magnitudes of inter-plate coupling. The 3D numerical models indicate the gap between the Cocos and Nazca slabs serves as a conduit for Pacific-Cocos mantle to pass into the Caribbean, with toroidal flow around the
NASA Astrophysics Data System (ADS)
Taylor, Marika
2006-03-01
Two charge BPS horizon free supergravity geometries are important in proposals for understanding black hole microstates. In this paper we construct a new class of geometries in the NS1-P system, corresponding to solitonic strings carrying fermionic as well as bosonic condensates. Such geometries are required to account for the full microscopic entropy of the NS1-P system. We then briefly discuss the properties of the corresponding geometries in the dual D1-D5 system.
ERIC Educational Resources Information Center
Suitor, Cheryl
2012-01-01
In science class, fourth graders investigate the structure of plants and leaves from trees and how the process of photosynthesis turns sunlight into sugar proteins. In this article, the author fuses art and science for a creative and successful clay slab project in her elementary art classroom. (Contains 1 online resource.)
Development of common conversion point stacking of receiver functions for detecting subducted slabs
NASA Astrophysics Data System (ADS)
Abe, Y.; Ohkura, T.; Hirahara, K.; Shibutani, T.
2010-12-01
In subduction zones, the subducting slabs are thought to convey fluid into the mantle wedge to cause arc volcanism (Hasegawa et al., 2008. Iwamori, 2007). Kawakatsu & Watada (2007) examined the Pacific slab subducting beneath northeast Japan with receiver function (RF) analysis, and revealed where the hydrated oceanic crust and the serpentinized mantle wedge exist. In the other subduction zones, it is also essential to examine subducting slabs for better understanding of water transportation and volcanic activities. In this study, we develop a new method to migrate RFs in order to examine subducting slabs with high dip angle (Abe et al., submitted to GJI) and apply this method to examination of the Philippine Sea slab (PHS). The RF technique is one of the useful methods to obtain seismic velocity discontinuities. Ps phases converted at discontinuities in a teleseismic coda can be detected by RF analysis. RFs are usually converted to depth domain assuming a 1-d velocity structure, and the geometry of discontinuities is obtained (e.g. Yamauchi et al., 2003). In subduction zones, however, subducting slabs usually dip, and we should take into account the refraction of seismic waves at dipping interfaces. Therefore, we use the multi-stage fast marching method (FMM, de Kool et al., 2006) to convert RFs into depth domain. We stack transverse RFs, since polarity of them does not change depending on their dip angles and they are better at detecting phases converted at dipping interfaces than radial RFs. We have confirmed that this method works properly with synthetic test. We apply our method to waveform data observed in Kyushu, Japan, where PHS is subducting toward WNW and the Wadati-Benioff zone dips at 30° at depths up to 80 km, and dips at 70° at depths between 80 km and 170 km. We obtain a vertical section, on which RF amplitude is projected, across central part of Kyushu perpendicular to the depth contour of the Wadati-Benioff zone. On the section, positive peaks of
Decay-ratio calculation in the frequency domain with the LAPUR code using 1D-kinetics
Munoz-Cobo, J. L.; Escriva, A.; Garcia, C.; Berna, C.
2012-07-01
This paper deals with the problem of computing the Decay Ratio in the frequency domain codes as the LAPUR code. First, it is explained how to calculate the feedback reactivity in the frequency domain using slab-geometry i.e. 1D kinetics, also we show how to perform the coupling of the 1D kinetics with the thermal-hydraulic part of the LAPUR code in order to obtain the reactivity feedback coefficients for the different channels. In addition, we show how to obtain the reactivity variation in the complex domain by solving the eigenvalue equation in the frequency domain and we compare this result with the reactivity variation obtained in first order perturbation theory using the 1D neutron fluxes of the base case. Because LAPUR works in the linear regime, it is assumed that in general the perturbations are small. There is also a section devoted to the reactivity weighting factors used to couple the reactivity contribution from the different channels to the reactivity of the entire reactor core in point kinetics and 1D kinetics. Finally we analyze the effects of the different approaches on the DR value. (authors)
NASA Astrophysics Data System (ADS)
Scire, A. C.; Zandt, G.; Beck, S. L.; Bishop, B.; Biryol, C. B.; Wagner, L. S.; Long, M. D.; Minaya, E.; Tavera, H.
2014-12-01
The modern central Peruvian Andes are dominated by a laterally extensive region of flat slab subduction. The Peruvian flat slab extends for ~1500 km along the strike of the Andes, correlating with the subduction of the Nazca Ridge in the south and the theorized Inca Plateau in the north. We have used data from the CAUGHT and PULSE experiments for finite frequency teleseismic P- and S-wave tomography to image the Nazca slab in the upper mantle below 95 km depth under central Peru between 10°S and 18°S as well as the surrounding mantle. Since the slab inboard of the subducting Nazca Ridge is mostly aseismic, our results provide important constraints on the geometry of the subducting Nazca slab in this region. Our images of the Nazca slab suggest that steepening of the slab inboard of the subducting Nazca Ridge locally occurs ~100 km further inland than was indicated in previous studies. The region where we have imaged the steepening of the Nazca slab inboard of the Nazca Ridge correlates with the location of the Fitzcarrald Arch, a long wavelength upper plate topographic feature which has been suggested to be a consequence of ridge subduction. When the slab steepens inboard of the flat slab region, it does so at a very steep (~70°) angle. The transition from the Peruvian flat slab to the more normally dipping slab south of 16°S below Bolivia is characterized by an abrupt bending of the slab anomaly in the mantle in response to the shift from flat to normal subduction. The slab anomaly appears to be intact south of the Nazca Ridge with no evidence for tearing of the slab in response to the abrupt change in slab dip. A potential tear in the slab is inferred from an observed offset in the slab anomaly north of the Nazca Ridge extending subparallel to the ridge axis between 130 and 300 km depth. A high amplitude (-5-6%) slow S-wave velocity anomaly is observed below the projection of the Nazca Ridge. This anomaly appears to be laterally confined to the mantle
Along-strike translation of a fossil slab
NASA Astrophysics Data System (ADS)
Eichenbaum-Pikser, J. M.; Forsyth, D. W.; Hirth, G.
2011-12-01
The Isabella anomaly is a high seismic velocity anomaly beneath the southern Central Valley of California. Breaking from previous interpretations of the anomaly as a lithospheric drip (Zandt and Carrigan, 1993) or delaminated Sierra Nevada root (Zandt et al., 2004), Forsyth et al. (2011) propose that it is a remnant slab, left over from Cenozoic subduction, attached to the Monterey microplate and translating along-strike with the Pacific plate underneath the edge of the North American plate. This hypothesis requires the slab to translate hundreds of kilometers along strike while remaining intact and attached to the Pacific plate despite drag from the surrounding asthenosphere and overriding lithosphere. Using COMSOL Multiphysics, we design 3-D finite element fluid flow models to simulate this scenario, and calculate the viscosity ratio required between the slab and the surrounding asthenosphere in order for the slab to translate undeformed. The ratio needed increases with downdip extent of the slab, and decreases with slab dip; for geometries approximating that of our proposed slab, it ranges from 10^2 to 10^4. Given the thermal and hydrological history of the slab, mantle flow laws predict viscosity contrasts greater than or equal to these requirements. As such, we conclude that along-strike translation of a remnant slab is entirely feasible, and serves as possible explanation of the Isabella anomaly. The significance of this finding extends beyond our general understanding of subduction dynamics, in that the presence of such a slab could have interesting implications for the water budget of the San Andreas Fault and its role in aseismic slip.
Lithosphere-Mantle Interactions Associated with Flat-Slab Subduction
NASA Astrophysics Data System (ADS)
Gerault, M.; Becker, T. W.; Husson, L.; Humphreys, E.
2014-12-01
Episodes of flat-slab subduction along the western margin of the Americas may have lead to the formation of intra-continental basins and seas, as well as mountain belts and continental plateaux. Here, we explore some of the consequences of a flat slab morphology, linking dynamic topography and stress patterns in continents to slab and mantle dynamics. Using a 2-D cylindrical code, we develop general models and apply them to the North and South America plates. The results are primarily controlled by the coupling along the slab-continent interface (due to geometry and viscosity), the viscosity of the mantle wedge, and the buoyancy of the subducted lithosphere. All models predict broad subsidence, large deviatoric stresses, and horizontal compression above the tip of the flat slab and the deep slab hinge. In models where the slab lays horizontally for hundreds of kilometers, overriding plate compression focuses on both ends of the flat segment, where normal-dip subduction exerts a direct downward pull. In between, a broad low-stress region gets uplifted proportionally to the amount of coupling between the slab and the continent. Anomalously buoyant seafloor enhances this effect but is not required. The downward bending of the flat slab extremities causes its upper part to undergo extension and the lower part to compress. These results have potential for explaining the existence of relatively undeformed, uplifted regions surrounded by mountain belts, such as in the western U.S. and parts of the Andes. Adequately modeling topography and stress in the unusual setting of southwestern Mexico requires a low-viscosity subduction interface and mantle wedge. Our results are only partially controlled by the buoyancy of the subducting plate, suggesting that the viscosity and the morphology of the slab are important, and that the often-used low resolution and "Stokeslet" models may be missing substantial effects.
Tulip, J.
1988-01-12
A gas slab laser is described comprising: first and second elongated electrodes each including a planar light reflecting surface disposed so as to form a light guide only in a plane perpendicular to the planar surface and to define a gas discharge gap therebetween; a laser gas disposed in the gap; and means for applying a radio frequency current between the first and second electrodes to establish a laser-exciting discharge in the laser gas.
Estimation of dielectric slab permittivity using a flared coaxial line
NASA Astrophysics Data System (ADS)
Shin, Dong H.; Eom, Hyo J.
2003-04-01
Estimation of dielectric slab permittivity is considered by using a flared coaxial line. A problem of reflection from a flared coaxial line that radiates into a dielectric slab with a flange is solved. A flared coaxial line is modeled with multiply stepped coaxial lines with different inner and outer conductors. A set of simultaneous equations for the modal coefficients is constituted based on the boundary conditions. Computations are performed to illustrate the reflection behavior in terms of the coaxial line geometry, frequency, and permittivity of a dielectric slab. Nomograms are developed to estimate the permittivity from the measured reflection coefficients. The utility of a flared coaxial line for the determination of slab permittivity is discussed.
The Role of Subducting Ridges in the Formation of Flat Slabs: Insights from the Peruvian Flat Slab
NASA Astrophysics Data System (ADS)
Knezevic Antonijevic, Sanja; Wagner, Lara; Kumar, Abhash; Beck, Susan; Long, Maureen; Zandt, George; Eakin, Caroline M.
2015-04-01
Flattening of the subducting plate is often used to explain various geological features removed far from the subducting margins, including basement-cored uplifts, the cessation of arc volcanism, ignimbrite flare-ups, and the formation of high plateaus and ore deposits [Humphreys et al., 2003; Gutscher et al., 2000; Rosenbaum et al., 2005, Kay and Mpodozis, 2001]. Today, flat slab subduction is observed in central Chile and Peru, representing the modern analogues to the immense paleo-flat slab that subducted beneath the North American continent during the Laramide orogeny (80-55 Ma) [English et al., 2003]. However, how flat slabs form and what controls their inboard and along-strike extent is still poorly understood. To better understand modern and paleo-flat slabs, we focus on the Peruvian flat slab, where the Nazca plate starts to bend at ~90 km depth and travels horizontally for several hundred kilometers beneath the South American plate. Earlier studies propose a correlation between the flat slab and the subducting Nazca Ridge that has been migrating to the south over the past 11 ~Ma [Hampel et al., 2004, Gutscher et al., 2003]. Combining 3D shear wave velocity structure and Rayleigh wave phase anisotropy between ~10° and 18° S, we find that the flat slab has the greatest inboard extent along the track of the subducting Nazca Ridge. North of the ridge track, where the flat slab was initially formed, the flat slab starts to sag, tear and re-initiate steep slab subduction, allowing inflow of warm asthenosphere. Based on our new constraints on the geometry of the subducted plate, we find that the subduction of buoyant oceanic features with overthickened oceanic crust plays a vital role in the formation of flat slabs. We further develop a model of temporal evolution of the Peruvian flab slab that forms as a result of the combined effects of the subducting ridge, trench retreat, and suction forces. Once the buoyant ridge subducts to ~90 km depth, it will fail to
Efficient 3D/1D self-consistent integral-equation analysis of ICRH antennae
NASA Astrophysics Data System (ADS)
Maggiora, R.; Vecchi, G.; Lancellotti, V.; Kyrytsya, V.
2004-08-01
This work presents a comprehensive account of the theory and implementation of a method for the self-consistent numerical analysis of plasma-facing ion-cyclotron resonance heating (ICRH) antenna arrays. The method is based on the integral-equation formulation of the boundary-value problem, solved via a weighted-residual scheme. The antenna geometry (including Faraday shield bars and a recess box) is fairly general and three-dimensional (3D), and the plasma is in the one-dimensional (1D) 'slab' approximation; finite-Larmor radius effects, as well as plasma density and temperature gradients, are considered. Feeding via the voltages in the access coaxial lines is self-consistently accounted throughout and the impedance or scattering matrix of the antenna array obtained therefrom. The problem is formulated in both the dual space (physical) and spectral (wavenumber) domains, which allows the extraction and simple handling of the terms that slow the convergence in the spectral domain usually employed. This paper includes validation tests of the developed code against measured data, both in vacuo and in the presence of plasma. An example of application to a complex geometry is also given.
Tomographic imaging of the Cascadia subduction zone: Constraints on the Juan de Fuca slab
NASA Astrophysics Data System (ADS)
Chen, Chuanxu; Zhao, Dapeng; Wu, Shiguo
2015-04-01
We used 40,343 P-wave arrival times from 1883 local earthquakes and 105,455 P-wave arrivals from 6361 teleseismic events to study the detailed structure of the Cascadia subduction zone. We conducted tomographic inversions using a starting velocity model which includes the high-velocity subducting Juan de Fuca slab as a priori information. A number of such slab-constrained inversions are conducted by changing the slab thickness and the velocity contrast between the slab and the surrounding mantle. Our optimal 3-D velocity model fits the data much better than that determined by an inversion with a 1-D homogeneous starting model. Our results show that the subducting Juan de Fuca slab has a thickness of 30-50 km and a P-wave velocity of 1-3% higher than that of the surrounding mantle. Beneath the northern and southern parts of the Cascadia, P-wave velocity is lower in the slab and along the slab interface, which may reflect a more hydrated slab and more active slab dehydration there. The lateral velocity variations may indicate different degrees of slab dehydration and forearc mantle serpentinization. The segmentation in episodic tremor and slip (ETS) is also spatially coincident with the velocity heterogeneities, indicating that the ETS occurrence and recurrence interval are controlled by fluid activity in and around the mantle wedge corner.
Slab laser development at MSNW - The Gemini and Centurion systems
NASA Astrophysics Data System (ADS)
Eggleston, J. M.; Albrecht, G. F.
Two, zig-zag-optical-path, slab-geometry, solid-state lasers, referred to as Gemini and Centurion, are described. The Nd:glass laser (Gemini) uses a pump geometry in which the flash lamps are located between two slabs in the same laser head. The dimensions and functions of the glass slabs are studied and the single-sided pumping of the Nd:glass laser is examined. The system is verified using the Nd:YAG laser system (Centurion). The Centurion system uses four flash lamps to pump a single 6 mm x 2 cm x 15.5 cm Nd:YAG slab; the reflector structure of the system is analyzed. The thermal-optical model for the Nd:glass laser and the Trace 3D, a three-dimensional flashlamp-slab coupling code, are evaluated. The oscillation performance and defocusing of a single-pass beam are measured; it is observed that the single-sided pump output is 30 percent more efficient than the standard configuration and no major defocusing effect is detected. The use of the Trace 3D code to design a reflector system for Gemini is discussed.
Double seismic zone and dehydration embrittlement of the subducting slab
NASA Astrophysics Data System (ADS)
Yamasaki, Tadashi; Seno, Tetsuzo
2003-04-01
Dehydration embrittlement of metamorphosed oceanic crust and mantle in the subducting slab may be responsible for the occurrence of intermediate-depth earthquakes. We explore the possibility that this hypothesis can explain the morphology of the double seismic zones observed in northeast Japan, southwest Japan, northeast Taiwan, northern Chile, Cape Mendocino, and eastern Aleutians. We calculate transient temperature structures of slabs based on geologically estimated subduction histories of these regions. We then determine dehydration loci of metamorphosed oceanic crust and serpentinized mantle using experimentally derived phase diagrams. The depth range of the dehydration loci of metamorphosed oceanic crust and serpentine is dependent on slab age. The dehydration loci of serpentine produce a double-layered structure. Because the upper dehydration loci of serpentine are mostly located in the wedge mantle above the slab, we regard the upper plane seismicity representing dehydration embrittlement in the oceanic crust, and we fix the slab geometry so that the upper plane seismicity is just below the upper surface of the slab. We find that the lower plane seismicity is located at the lower dehydration loci of serpentine, which indicates that the morphology of the double seismic zones is consistent with the dehydration embrittlement.
Sub-slab sampling has become an integral part of vapor intrusion investigations. It is now recommended in guidance documents developed by EPA and most states. A method for sub-slab probe installation was devised in 2002, presented at conferences through 2005, and finally docume...
[Thermoluminescence Slab Dosimeter].
Shinsho, Kiyomitsu; Koba, Yusuke; Tamatsu, Satoshi; Sakurai, Noboru; Wakabayashi, Genichiro; Fukuda, Kazusige
2013-01-01
In 1953 F. Daniels et al. used the property of thermoluminescence in dosimetry for the first time. Since then, numerous TLD have been developed. 2D TLD was investigated for the first time in 1972 by P Broadhead. However, due to excessive fading, difficulties with handling and the time required for measurements, development stalled. At the current time, the majority of TLD are used in small scale, localized dosimetry with a wide dynamic range and personal dosimeters for exposure management. Urushiyama et. al. have taken advantage of the commoditization of CCD cameras in recent years--making large area, high resolution imaging easier--to introduce and develop a 2D TLD. It is expected that these developments will give rise to a new generation of applications for 2D TL dosimetry. This paper introduces the "TL Slab Dosimeter" developed jointly by Urushiyama et. al. and our team, its measurement system and several typical usage scenarios. PMID:24893451
NASA Astrophysics Data System (ADS)
MacDougall, Julia G.; Kincaid, Chris; Szwaja, Sara; Fischer, Karen M.
2014-05-01
Observed seismic anisotropy and geochemical anomalies indicate the presence of 3-D flow around and above subducting slabs. To investigate how slab geometry and velocity affect mantle flow, we conducted a set of experiments using a subduction apparatus in a fluid-filled tank. Our models comprise two independently adjustable, continuous belts to represent discrete sections of subducting slabs that kinematically drive flow in the surrounding glucose syrup that represents the upper mantle. We analyse how slab dip (ranging from 30° to 80°), slab dip difference between slab segments (ranging from 20° to 50°), rates of subduction (4-8 cm yr-1) and slab/trench rollback (0-3 cm yr-1) affect mantle flow. Whiskers were used to approximate mineral alignment induced by the flow, as well as to predict directions of seismic anisotropy. We find that dip variations between slab segments generate 3-D flow in the mantle wedge, where the path lines of trenchward moving mantle material above the slab are deflected towards the slab segment with the shallower dip. The degree of path line deflection increases as the difference in slab dip between the segments increases, and, for a fixed dip difference, as slab dip decreases. In cases of slab rollback and large slab dip differences, we observe intrusion of subslab material through the gap and into the wedge. Flow through the gap remains largely horizontal before eventual downward entrainment. Whisker alignment in the wedge flow is largely trench-normal, except near the lateral edges of the slab where toroidal flow dominates. In addition, whisker azimuths located above the slab gap deviate most strongly from trench-normal orientations when slab rollback does not occur. Such flow field complexities are likely sufficient to affect deep melt production and shallow melt delivery. However, none of the experiments produced flow fields that explain the trench-parallel shear wave splitting fast directions observed over broad arc and backarc
Arc Interrupted: The birth, life, and death of the Peruvian flat slab
NASA Astrophysics Data System (ADS)
Wagner, Lara; Knezevic Antonijevic, Sanja; Kumar, Abhash; Beck, Susan; Long, Maureen; Zandt, George; Tavera, Hernando; Minaya, Estela
2016-04-01
The Peruvian flat slab is a unique natural laboratory for investigating the temporal evolution of flat slab subduction and its associated thermal, tectonic, and seismic implications. This is because the flat slab has been hypothesized to have first formed further north (at approximately the latitude of Lima, Peru), but broadened to the south over the past 11 Ma. This means that areas further to the north represent an older, more evolved flat slab setting, whereas the southernmost edge of the modern flat slab reflects conditions experienced by a newly formed flat slab. Here we present findings from a suite of recent temporary broadband seismic deployments that spanned this entire region. Results from intermediate depth earthquake locations, surface wave tomography (ballistic and ambient Rayleigh wave), and Rayleigh wave anisotropy all indicate that the flat slab did indeed first form further to the north and broadened to the south, along with the southward migration of the Nazca ridge. Subsequently, a trench-parallel tear developed in the older portions of the flat slab north of the ridge, resulting in a resumption of normal subduction geometry where once a flat slab had existed. This tear allows for an interchange of mantle material from beneath the slab to the south to above the slab to the north. This mantle flow has significant thermal implications, both beneath the flat slab and in the lower continental crust located above the relatively newly formed tear. Our results provide unique constraints on the thermal and tectonic evolution of this unusual subduction geometry that may help us to understand better subduction zone processes everywhere.
Diffractive properties of imaginary-part photonic crystal slab
2012-01-01
The diffraction spectra of imaginary-part photonic crystal (IPPC) slabs are analyzed using the scattering-matrix method. By investigating the thickness dependence of the diffraction, we find a remarkable red shift of central wavelength of the diffraction spectrum, which obviously distinguishes from the phenomenon of spectral hole. We observe that diffraction efficiency can be enhanced more than 20-fold by optimizing the geometry parameters. These imply that the diffraction spectra of the IPPC slab can be controlled at will and used to guide the design to achieve useful nanoscale devices. PMID:22720871
Coherent thermal conductance of 1-D photonic crystals
NASA Astrophysics Data System (ADS)
Tschikin, Maria; Ben-Abdallah, Philippe; Biehs, Svend-Age
2012-10-01
We present an exact calculation of coherent thermal conductance in 1-D multilayer photonic crystals using the S-matrix method. In particular, we study the thermal conductance in a bilayer structure of Si/vacuum or Al2O3/vacuum slabs by means of the exact radiative heat flux expression. Based on the results obtained for the Al2O3/vacuum structure we show by comparison with previous works that the material losses and (localized) surface modes supported by the inner layers play a fundamental role and cannot be omitted in the definition of thermal conductance. Our results could have significant implications in the conception of efficient thermal barriers.
New Packaging for Amplifier Slabs
Riley, M.; Thorsness, C.; Suratwala, T.; Steele, R.; Rogowski, G.
2015-03-18
The following memo provides a discussion and detailed procedure for a new finished amplifier slab shipping and storage container. The new package is designed to maintain an environment of <5% RH to minimize weathering.
Rheologic Controls on the Dynamic Evolution of Slabs in the Upper Mantle
NASA Astrophysics Data System (ADS)
Billen, M.; Hirth, G.
2004-12-01
Subduction of tectonic plates is characterized by long-lived subduction zones, asymmetric subduction and slab dip angles of 25--80o in the upper mantle. Several mechanisms proposed to explain the variation in observed dip include large-scale mantle flow, trench roll-back, and interaction of the slab with the transition zone. Previous dynamic models of subduction that include only Newtonian viscosity and moderately strong slabs generally fail to predict subduction angles less than 60--90o at shallow depths (100--300 km). We find that the observed characteristics of subduction are reproduced by viscous flow models, in which the rheologic structure is consistent with experimentally determined flow laws for Newtonian and non-Newtonian visco-plastic deformation of olivine. The properties of the models required to match the observed characteristics of slabs are: non-Newtonian viscosity in the mantle producing a weak mantle wedge (1018--1019~Pa s), a stiff slab interior (1025~Pa s) limited by a plastic yield criterion and a weak plate boundary shear zone (1020--1021~Pa s). The shallow slab dip reaches a minimum of 25--30o for high convergence rates and a stiff slab, without trench roll-back or relative motion of the entire lithosphere with respect to the mantle, suggesting these other mechanisms are not the primary controls on slab geometry. The deep slab dip (350--650 km) decreases as the slab penetrates the stiffer (x10), Newtonian viscosity lower mantle, eventually stabilizing the upper mantle slab geometry.
Wu, R.S.; Xie, X.B.
1994-12-31
The theory of spatial distribution of seismic energy density in one dimensional (1D) random media derived in part 1 (Wu, 1993) is tested by numerical experiments using a full wave propagation matrix method. The geometry of numerical experiment mimics the configuration of zero-offset VSP (Vertical Seismic Profiling) along a borehole. A procedure of octave-band frequency averaging is applied to the measured data to reduce fluctuation of spatial energy distribution, so that stable estimations of medium parameters can be achieved without resorting to ensemble averaging. Results from Monte-Carlo numerical experiments for both infinite random media and finite random slabs with or without bottom reflections show good agreement for dark-to-gray (weak to intermediate scattering compared with absorption) media. When scattering is very strong (when backscattering-absorption ratio S{sub b} > 3), results from single realization fluctuate substantially. However, most the practical situations of sedimentary rocks in the crust fall into the validity region of the energy transfer theory.
NASA Astrophysics Data System (ADS)
Burckel, David Bruce
wavelength to period ratio. These optical results indicated a need for a deeper understanding of the confinement/guiding mechanisms in such waveguide structures. A simplification of the full 2-D problem to a more tractable "tilted 1-D" geometry led to the proposal of a new waveguide geometry, Generalized Transverse Bragg Waveguides (GTBW), as well as a new propagation mode characterized by spatial variation in both the transverse direction as well as the direction of propagation. GTBW demonstrate many of the same dispersion tunability traits exhibited in complete bandgap photonic crystal waveguides, under more modest fabrication demands, and moreover provide much insight into photonic crystal waveguide modes of all types. Generalized Transverse Bragg Waveguides are presented in terms of the standard physical properties associated with waveguides, including the dispersion relation, expressions for the spatial field profile, and the concepts of phase and group velocity. In addition, the proposal of at least one obvious application, semiconductor optical amplifiers, is offered.
Cryogenic, high power, near diffraction limited, Yb:YAG slab laser.
Ganija, Miftar; Ottaway, David; Veitch, Peter; Munch, Jesper
2013-03-25
A cryogenic slab laser that is suitable for scaling to high power, while taking full advantage of the improved thermo-optical and thermo-mechanical properties of Yb:YAG at cryogenic temperatures is described. The laser uses a conduction cooled, end pumped, zigzag slab geometry resulting in a near diffraction limited, robust, power scalable design. The design and the initial characterization of the laser up to 200W are presented. PMID:23546080
Evidence for a coherent southeastern edge of the Farallon slab window from USArray tomography models
NASA Astrophysics Data System (ADS)
Panessa, Alexander
Before the San Andreas Fault developed, subduction of the ancient Farallon plate dominated North America's west coast. As the Pacific-Farallon ridge reached the continent and Farallon subduction became restricted to separate regions, an area with no slab formed under southwestern North America. The northern boundary of this "slab window" is well-defined in recent high-resolution tomography models derived from USArray data, but the geometry of the southeastern edge of the slab window has received less attention. I will show that this eastern edge exists as a coherent boundary separating the subducted Cocos plate from the slab window. I adapt stage pole data for the Cocos plate into a single surface slab model designed to track the motion of the plate through time. This model of the Cocos plate geometry is consistent with recently published tomography models derived from combinations of USArray and global seismic data. The tomography results suggest the subducted slab remains intact along the southeastern edge of the Farallon slab window.
Down-dip and Along-strike Stress in Subducting Slabs
NASA Astrophysics Data System (ADS)
Petricca, P.; Carminati, E.; Basili, R.; Doglioni, C.
2014-12-01
The reconstructed stress fields for intermediate depths in slabs (between 100 and 300 km) remain enigmatic. Slab pull should enhance downdip tension, but downdip compression is often observed.Based on hotspot tracks (indicating a relative motion between the lithosphere and the mantle) several researchers proposed a westward drift of the lithosphere with debated velocities (2-13 cm/yr). Available data indicates that dominant down-dip tension occurs along E-directed subductions (e.g., Chile) sustained by the eastward relative mantle flow and down-dip compression along W-directed subductions (e.g., Tonga) encroached by the mantle motion. In this perspective, 3D viscoelastic models (FEM) were performed to investigate the origin of stress in slabs at intermediate depths. Geometries inspired to real subduction zones, were used to analyze the effect of the principal geodynamic forces (slab pull, mantle flow, plate convergence), the correlation between subduction polarity (i.e., how the mantle flow approaches the slab) and geometry, the nature of down-dip and along-strike stresses (poorly investigated in the literature).We conclude that, although the stress in slabs is controlled by several geodynamic processes, down-dip compression is favored by mantle flow opposing the slab, whereas down-dip tension is enhanced by mantle flow sustaining the slab. These predictions are in agreement with available geophysical observations.Substantial along-strike deformation affects intermediate and deeper regions of the slab when mantle flow is enhanced. Down-dip and along-strike axes form couples that tend to rotate approaching the slab sides. The larger the lateral extent of the plate the most this effect is evident. According to these findings, along-strike tension or compression comes out in response to the flexure (forward or backward) of the lithosphere depending on the slab longitudinal curvature. Along-strike tension occurs in case of slabs curved towards the mantle flow direction
A wrinkly phononic crystal slab
NASA Astrophysics Data System (ADS)
Bayat, Alireza; Gordaninejad, Faramarz
2015-03-01
The buckling induced surface instability is employed to propose a tunable phononic crystal slab composed of a stiff thin film bonded on a soft elastomer. Wrinkles formation is used to generate one-dimensional periodic scatterers at the surface of a finitely thick slab. Wrinkles' pattern change and corresponding stress is employed to control wave propagation triggered by a compressive strain. Simulation results show that the periodic wrinkly structure can be used as a transformative phononic crystal which can switch band diagram of the structure in a reversible behavior. Results of this study provide opportunities for the smart design of tunable switch and elastic wave filters at ultrasonic and hypersonic frequency ranges.
Phase conjugated slab laser designator
Chandra, S.; Paul, J.L.
1989-06-06
A laser designator is described comprising a laser pump means; a high power phase conjugated slab laser amplifier formed of GSGG:Cr:Nd as a lasing material on one side of the pump means; a low power rod shaped laser oscillator on the opposite side of the pump means; and a first plurality of optical reflecting and refracting means for directing and shaping a laser beam which surrounds the pump means and passes through the rod and slab; and a telescope means coupled to the beam to direct it to a distant target.
NASA Astrophysics Data System (ADS)
Guillaume, Benjamin; Moroni, Monica; Funiciello, Francesca; Martinod, Joseph; Faccenna, Claudio
2010-12-01
We present dynamically self-consistent mantle-scale laboratory models that have been conducted to improve our understanding of the influence of slab window opening on subduction dynamics, mantle flow and associated dynamic topography over geological time scales. The adopted setup consists of a two-layer linearly viscous system simulating the subduction of a fixed plate of silicone (lithosphere) under negative buoyancy in a viscous layer of glucose syrup (mantle). Our experimental setting is also characterized by a constant-width rectangular window located at the center of a laterally confined slab, modeling the case of the interaction of a trench-parallel spreading ridge with a wide subduction zone. We found that the opening of a slab window does not produce consistent changes of the geometry and the kinematics of the slab. On the contrary, slab-induced mantle circulation, quantified both in the vertical and horizontal sections using the Feature Tracking image analysis technique, is strongly modified. In particular, rollback subduction and the opening of the slab window generate a complex mantle circulation pattern characterized by the presence of poloidal and toroidal components, with the importance of each evolving according to kinematic stages. Mantle coming from the oceanic domain floods through the slab window, indenting the supra-slab mantle zone and producing its deformation without any mixing between mantle portions. The opening of the slab window and the upwelling of sub-slab mantle produce a regional-scale non-isostatic topographic uplift of the overriding plate that would correspond to values ranging between ca. 1 and 5 km in nature. Assuming that our modeling results can be representative of the natural behavior of subduction zones, we compared them to the tectonics and volcanism of the Patagonian subduction zone. We found that the anomalous backarc volcanism that has been developing since the middle Miocene could result from the lateral flow of sub-slab
Seismicity and state of stress in the central and southern Peruvian flat slab
NASA Astrophysics Data System (ADS)
Kumar, Abhash; Wagner, Lara S.; Beck, Susan L.; Long, Maureen D.; Zandt, George; Young, Bissett; Tavera, Hernando; Minaya, Estella
2016-05-01
We have determined the Wadati-Benioff Zone seismicity and state of stress of the subducting Nazca slab beneath central and southern Peru using data from three recently deployed local seismic networks. Our relocated hypocenters are consistent with a flat slab geometry that is shallowest near the Nazca Ridge, and changes from steep to normal without tearing to the south. These locations also indicate numerous abrupt along-strike changes in seismicity, most notably an absence of seismicity along the projected location of subducting Nazca Ridge. This stands in stark contrast to the very high seismicity observed along the Juan Fernandez ridge beneath central Chile where, a similar flat slab geometry is observed. We interpret this as indicative of an absence of water in the mantle beneath the overthickened crust of the Nazca Ridge. This may provide important new constraints on the conditions required to produce intermediate depth seismicity. Our focal mechanisms and stress tensor inversions indicate dominantly down-dip extension, consistent with slab pull, with minor variations that are likely due to the variable slab geometry and stress from adjacent regions. We observe significantly greater variability in the P-axis orientations and maximum compressive stress directions. The along strike change in the orientation of maximum compressive stress is likely related to slab bending and unbending south of the Nazca Ridge.
How mantle slabs drive plate tectonics.
Conrad, Clinton P; Lithgow-Bertelloni, Carolina
2002-10-01
The gravitational pull of subducted slabs is thought to drive the motions of Earth's tectonic plates, but the coupling between slabs and plates is not well established. If a slab is mechanically attached to a subducting plate, it can exert a direct pull on the plate. Alternatively, a detached slab may drive a plate by exciting flow in the mantle that exerts a shear traction on the base of the plate. From the geologic history of subduction, we estimated the relative importance of "pull" versus "suction" for the present-day plates. Observed plate motions are best predicted if slabs in the upper mantle are attached to plates and generate slab pull forces that account for about half of the total driving force on plates. Slabs in the lower mantle are supported by viscous mantle forces and drive plates through slab suction. PMID:12364804
NASA Astrophysics Data System (ADS)
Scire, Alissa
The Nazca-South America convergent margin is marked by the presence of the Andean mountain belt, which stretches along the 8000-km long western margin of the South American plate. The subduction zone is characterized by significant along-strike changes in both upper plate structure and slab geometry that make it an ideal region to study the relationship between the subducting slab, the surrounding mantle, and the overriding plate. This dissertation summarizes the results of three finite frequency teleseismic tomography studies of the central Nazca-South America subduction zone which improve our understanding of how along-strike variations in the Andean mountain belt and the subducting Nazca plate interact with each other and with the surrounding mantle. This is accomplished by first focusing on two smaller adjacent regions of the central Andes to explore upper mantle variations and then by using a combined dataset, which covers a larger region, to image the deeply subducted Nazca slab to investigate the fate of the slab. The first study focuses on the central Andean upper mantle under the Altiplano-Puna Plateau where normally dipping subduction of the Nazca plate is occurring (18° to 28°S). The shallow mantle under the Eastern Cordillera is generally fast, consistent with either underthrusting of the Brazilian cratonic lithosphere from the east or a localized "curtain" of delaminating material. Additional evidence for delamination is seen in the form of high amplitude low velocities under the Puna Plateau, consistent with proposed asthenospheric influx following lithospheric removal. In the second study, we explore the transition between normal and flat subduction along the north edge of the Altiplano Plateau (8° to 21°S). We find that the Peruvian flat slab extends further inland along the projection of the Nazca Ridge than was previously proposed and that when re-steepening of the slab occurs, the slab dips very steeply (˜70°) down through the mantle
Vibration characteristics of floating slab track
NASA Astrophysics Data System (ADS)
Kuo, Chen-Ming; Huang, Cheng-Hao; Chen, Yi-Yi
2008-11-01
Coupled equilibrium equations of suspended wheels and floating slab track system were solved with the fourth-order Runge-Kutta method to obtain the deflections, vibration velocities, and wheel-rail contact forces. The program was validated through several aspects. Cases with various vehicle speed, slab mass, and stiffness of slab bearing were analyzed to reveal the effects of slab bearing on track responses. The correlation between wheel-rail resonance and train speed was also discussed. It was found that rail deflections increase significantly as train speed increases. Although large slab mass may lower tuning frequency, it could also result in higher wheel-rail contact force and rail deflections. The floating slab track is effective in isolating loading above 10 Hz, which might present in some railway sections with irregularities. Adopting floating slab track for vibration control for environment along the railway may cause concerns about ride quality and track damages.
NASA Astrophysics Data System (ADS)
Dougherty, Sara L.; Clayton, Robert W.
2014-04-01
The morphology of the transition from flat to normal subduction in eastern central Mexico is explored using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. Observations of a sharp transition in slab dip near the abrupt end of the Trans-Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of a thin ultra-slow velocity layer (USL) imaged atop the Cocos slab in recent studies along the Meso America Subduction Experiment array is examined here using additional data. We find an end to this USL which is coincident with the western boundary of a zone of decreased seismicity and the end of the TMVB near the sharp transition in slab dip. Waveform modeling of the 2-D structure in this region using a finite difference algorithm provides constraints on the velocity and geometry of the slab's seismic structure and confirms the location of the USL. Analysis of intraslab seismicity patterns reveals clustering, sudden increase in depth, variable focal mechanism orientations and faulting types, and alignment of source mechanisms along the sharp transition in slab dip. The seismicity and structural evidence suggests a possible tear in the South Cocos slab. This potential tear, together with the tear along the Orozco Fracture Zone to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments.
Preface: Deep Slab and Mantle Dynamics
NASA Astrophysics Data System (ADS)
Suetsugu, Daisuke; Bina, Craig R.; Inoue, Toru; Wiens, Douglas A.
2010-11-01
We are pleased to publish this special issue of the journal Physics of the Earth and Planetary Interiors entitled "Deep Slab and Mantle Dynamics". This issue is an outgrowth of the international symposium "Deep Slab and Mantle Dynamics", which was held on February 25-27, 2009, in Kyoto, Japan. This symposium was organized by the "Stagnant Slab Project" (SSP) research group to present the results of the 5-year project and to facilitate intensive discussion with well-known international researchers in related fields. The SSP and the symposium were supported by a Grant-in-Aid for Scientific Research (16075101) from the Ministry of Education, Culture, Sports, Science and Technology of the Japanese Government. In the symposium, key issues discussed by participants included: transportation of water into the deep mantle and its role in slab-related dynamics; observational and experimental constraints on deep slab properties and the slab environment; modeling of slab stagnation to constrain its mechanisms in comparison with observational and experimental data; observational, experimental and modeling constraints on the fate of stagnant slabs; eventual accumulation of stagnant slabs on the core-mantle boundary and its geodynamic implications. This special issue is a collection of papers presented in the symposium and other papers related to the subject of the symposium. The collected papers provide an overview of the wide range of multidisciplinary studies of mantle dynamics, particularly in the context of subduction, stagnation, and the fate of deep slabs.
Numerical Models of Subduction and Slab Detachment: Placing a Lower Bound on the Strength of Slabs
NASA Astrophysics Data System (ADS)
Billen, M. I.; Andrews, E.
2007-12-01
Subduction provides the main driving force for the motion of tectonic plates at the Earth's surface through slab- pull and sinking-induced flow in the surrounding mantle. The ability of the slab to directly transmit slab-pull forces to the tectonic plate at the surface depends on the minimum strength and rheology (e.g., viscous, plastic) of the slab. Previous models have shown that observations including the state of stress in slabs, dynamic topography and the geoid above slabs, the evolution of slab and the kinematic history of subduction can be well-matched by a variety of models with either low viscosity (i.e., 100-1,000 times more viscous than the surrounding mantle) or high viscosity slabs (i.e., more than 10,000 times more viscous than the surrounding mantle). However, in many of the models in which a good match to observations is found for low viscosity slabs, the maximum slab viscosity is imposed as a cut-off value, which forces the entire slab to have a more or less uniform viscosity independent of strain-rate or stress magnitude, rather than a plastic yielding-type rheology. We present numerical models demonstrating that when the non-Newtonian viscosity of the upper mantle and plastic yielding behavior of slabs are taken into account, the minimum yield strength that allows for continuous subduction is approximately 300- 500 MPa, which leads to high viscosity slabs with some localized lower viscosity regions. A yield stress of 10-100 MPa is required to form uniformly low viscosity slabs, but these slabs detach from the subducting plate, due to localized weakening, when the slab reaches a length of 200-300 km, even when subduction is facilitated by a low viscosity shear zone and kinematically-imposed surface velocities. In contrast, detachment of higher strength slabs in fully-dynamic models only occurs when the shear zone is removed and prevents further subduction.
Slab melting versus slab dehydration in subduction-zone magmatism.
Mibe, Kenji; Kawamoto, Tatsuhiko; Matsukage, Kyoko N; Fei, Yingwei; Ono, Shigeaki
2011-05-17
The second critical endpoint in the basalt-H(2)O system was directly determined by a high-pressure and high-temperature X-ray radiography technique. We found that the second critical endpoint occurs at around 3.4 GPa and 770 °C (corresponding to a depth of approximately 100 km in a subducting slab), which is much shallower than the previously estimated conditions. Our results indicate that the melting temperature of the subducting oceanic crust can no longer be defined beyond this critical condition and that the fluid released from subducting oceanic crust at depths greater than 100 km under volcanic arcs is supercritical fluid rather than aqueous fluid and/or hydrous melts. The position of the second critical endpoint explains why there is a limitation to the slab depth at which adakitic magmas are produced, as well as the origin of across-arc geochemical variations of trace elements in volcanic rocks in subduction zones. PMID:21536910
Slab melting versus slab dehydration in subduction-zone magmatism
Mibe, Kenji; Kawamoto, Tatsuhiko; Matsukage, Kyoko N.; Fei, Yingwei; Ono, Shigeaki
2011-01-01
The second critical endpoint in the basalt-H2O system was directly determined by a high-pressure and high-temperature X-ray radiography technique. We found that the second critical endpoint occurs at around 3.4 GPa and 770 °C (corresponding to a depth of approximately 100 km in a subducting slab), which is much shallower than the previously estimated conditions. Our results indicate that the melting temperature of the subducting oceanic crust can no longer be defined beyond this critical condition and that the fluid released from subducting oceanic crust at depths greater than 100 km under volcanic arcs is supercritical fluid rather than aqueous fluid and/or hydrous melts. The position of the second critical endpoint explains why there is a limitation to the slab depth at which adakitic magmas are produced, as well as the origin of across-arc geochemical variations of trace elements in volcanic rocks in subduction zones. PMID:21536910
Thermal runaway in microwave heated isothermal slabs, cylinders, and spheres
NASA Astrophysics Data System (ADS)
Vriezinga, C. A.
1998-01-01
The absorption of electromagnetic energy within a microwave heated isothermal slab, cylinder, and sphere is analyzed and compared to each other. It is shown that the absorbed heat oscillates as a function of temperature, regardless of the geometry of the irradiated object. It is possible to formulate this behavior in a simple mathematical equation, which proves that the oscillation is basically caused by resonance of the electromagnetic waves within the object. This oscillation, combined with the heat loss, is found to be responsible for thermal runaway phenomenon in isothermal objects. Based on such an observation, a general rule to prevent thermal runaway has been developed.
Three-dimensional Numerical Models of Mantle Flow Through the Cocos-Nazca Slab Gap
NASA Astrophysics Data System (ADS)
Jadamec, M.; Fischer, K. M.
2013-05-01
Global slab geometry models suggest a 350 km to 1000 km spacing between the southern extent of the Cocos slab and the northern extent of the Nazca slab (Gudmundsson and Sambridge, 1998; Syracuse and Abers, 2006; Hayes et al., 2012). The apparent gap between the east-dipping Cocos and Nazca slabs at depth correlates to several tectonic features on the Pacific side of Central and northern South America that may limit subduction, namely the (a) Panama Fault zone, (b) incoming young lithosphere associated with the Cocos-Nazca spreading center, and (c) the Cocos, Coiba, Malpelo, and Carnegie ridges associated with the Galapogos hotspot and Cocos-Nazca spreading center (Protti et al., 1994; Johnston and Thorkelson, 1997; Gutscher et al., 1999; Abratis and Worner, 2001; Sdrolias and Muller, 2006; Mann et al., 2007; Gazel et al., 2011). In addition, on the Caribbean side of Central and northern South America, seismic data suggest that part of the Caribbean plate is subducting and dipping in a direction opposite to the Cocos and Nazca slabs (van der Hilst and Mann, 1994; Camacho et al., 2010). We construct high-resolution three-dimensional numerical models of the Cocos-Nazca subduction system to test the effects of a slab gap and variable overriding plate thickness on surface plate motion and mantle flow. The 3D tectonic configuration is generated with SlabGenerator (Jadamec and Billen, 2010, 2012) and the mantle convection code CitcomCU is used to solve for the viscous flow (Moresi and Solomatov, 1995; Zhong, 2006). The negative thermal buoyancy of the slabs drive the flow. No driving velocities are applied to the plates or any of the slabs in the model. The detailed geometries of the Cocos and Nazca slabs are constructed from seismicity and seismic tomography (Protti et al., 1994; Colombo et al., 1997; Gudmundsson and Sambridge, 1998; Rogers et al., 2002; Husen et al., 2003; Syracuse and Abers, 2006; Syracuse et al., 2008; Dzierma et al., 2011). Seismic tomography
High Average Power Nd:YAG Slab Laser
NASA Astrophysics Data System (ADS)
Kasai, Takeshi; Sindo, Yoshihiko; Haga, Keiji
1989-07-01
A slab geometry Nd:YAG laser with a zigzag optical path is described. The dimensions of the Nd:YAG slab are 5.6 x 18.4 x 153.9 mm, and Nei' ion concentration is 1.1 at.%. Two krypton flashlamps, one located on each side of the YAG slab, are used for pumping. The conditions for normal pulsed operation were as follows: the repetition rate was from 5 to 27 pps, and the pulse durations were 4 and 9.9 ms. With the above conditions, a maximum average output power of 500 W was obtained with an efficiency of 2 %, the slope efficiency being 2.4 %. The beam divergence was estimated to be 10x25 mrad. The stability of the laser output power was about +/-1.5 %. Another oscillator that includes intra-cavity cylindrical lenses, was also designed. Using this resonator configuration reduced the beam divergence to about 7.6 x8.2 mrad. The preliminary laser processing experiment was attemped using this laser oscillator.
Dynamic uplift during slab flattening
NASA Astrophysics Data System (ADS)
Dávila, Federico M.; Lithgow-Bertelloni, Carolina
2015-09-01
Subduction exerts a strong control on surface topography and is the main cause of large vertical motions in continents, including past events of large-scale marine flooding and tilting. The mechanism is dynamic deflection: the sinking of dense subducted lithosphere gives rise to stresses that directly pull down the surface. Here we show that subduction does not always lead to downward deflections of the Earth's surface. Subduction of young lithosphere at shallow angles (flat subduction) leaves it neutrally or even positively buoyant with respect to underlying mantle because the lithosphere is relatively warm compared with older lithosphere, and because the thickened and hence drier oceanic crust does not undergo the transformation of basalt to denser eclogite. Accounting for neutrally buoyant flat segments along with large variations in slab morphology in the South American subduction zone explains along-strike and temporal changes in dynamic topography observed in the geologic record since the beginning of the Cenozoic. Our results show that the transition from normal subduction to slab flattening generates dynamic uplift, preventing back-arc marine flooding.
Lee, Myoung-Jae; Jung, Young-Dae
2015-02-15
The nonthermal and geometric effects on the propagation of the surface dust acoustic waves are investigated in a Lorentzian dusty plasma slab. The symmetric and anti-symmetric dispersion modes of the dust acoustic waves are obtained by the plasma dielectric function with the spectral reflection conditions the slab geometry. The variation of the nonthermal and geometric effects on the symmetric and the anti-symmetric modes of the surface plasma waves is also discussed.
Opening and closing slab windows in congested subduction zones
NASA Astrophysics Data System (ADS)
Moresi, Louis
2013-04-01
Subduction zones often try to swallow buoyant material which is embedded in the oceanic lithosphere: plume material or hotspot residues, oceanic plateaux, and fragments of continental material. This often results in the formation of a slab window and it has been shown (Mason et al, 2010; Betts et al, 2012) that this window strongly influences the subsequent evolution of the slab and the advance/retreat rate of the trench. The buoyant material typically pushes the trench into a local state of advance, and the creation of the slab window allows the rest of the trench to retreat as the mantle behind the slab flows in through the window. This situation is inherently unstable: if the buoyancy anomaly is finite in size, then the retreating trench will soon move behind the anomaly and juxtapose negatively buoyant oceanic lithosphere with active subduction. This creates the potential to close the slab window and, in doing so, transfer the buoyant material to the over-riding plate. Models show that this closure of the window initially occurs through a lateral rollback process followed by a catastrophic re-initiation of subduction behind the colliding buoyant anomaly. This rollback leaves a characteristic, tightly rolled remnant in the mantle and significant rotation in the over-riding plate and the newly-docked block. The over-riding plate is thrown into extension perpendicular to the original orientation of the trench. This same situation applies at the late-stages of a closing ocean due to the passive margin geometry and the presence of debris collected from the closing ocean floor and it seems likely that these models can also be applied to the complicated geometry of subduction in such environments. Mason, W. G.; Moresi, L.; Betts, P. G. & Miller, M. S. Three-dimensional numerical models of the influence of a buoyant oceanic plateau on subduction zones Tectonophysics, 2010, 483, 71-79 P. Betts, W. Mason, L. Moresi, The influence of mantle plumes on subduction zone
Subducting Slabs: Jellyfishes in the Earth's Mantle
NASA Astrophysics Data System (ADS)
Loiselet, C.; Braun, J.; Husson, L.; Le Carlier de Veslud, C.; Thieulot, C.; Yamato, P.; Grujic, D.
2010-12-01
The constantly improving resolution of geophysical data, seismic tomography and seismicity in particular, shows that the lithosphere does not subduct as a slab of uniform thickness but is rather thinned in the upper mantle and thickened around the transition zone between the upper and lower mantle. This observation has traditionally been interpreted as evidence for the buckling and piling of slabs at the boundary between the upper and lower mantle, where a strong contrast in viscosity may exist and cause resistance to the penetration of slabs into the lower mantle. The distribution and character of seismicity reveal, however, that slabs undergo vertical extension in the upper mantle and compression near the transition zone. In this paper, we demonstrate that during the subduction process, the shape of low viscosity slabs (1 to 100 times more viscous than the surrounding mantle) evolves toward an inverted plume shape that we coin jellyfish. Results of a 3D numerical model show that the leading tip of slabs deform toward a rounded head skirted by lateral tentacles that emerge from the sides of the jellyfish head. The head is linked to the body of the subducting slab by a thin tail. A complete parametric study reveals that subducting slabs may achieve a variety of shapes, in good agreement with the diversity of natural slab shapes evidenced by seismic tomography. Our work also suggests that the slab to mantle viscosity ratio in the Earth is most likely to be lower than 100. However, the sensitivity of slab shapes to upper and lower mantle viscosities and densities, which remain poorly constrained by independent evidence, precludes any systematic deciphering of the observations.
Subducting slabs: Jellyfishes in the Earth's mantle
NASA Astrophysics Data System (ADS)
Loiselet, Christelle; Braun, Jean; Husson, Laurent; Le Carlier de Veslud, Christian; Thieulot, Cedric; Yamato, Philippe; Grujic, Djordje
2010-08-01
The constantly improving resolution of geophysical data, seismic tomography and seismicity in particular, shows that the lithosphere does not subduct as a slab of uniform thickness but is rather thinned in the upper mantle and thickened around the transition zone between the upper and lower mantle. This observation has traditionally been interpreted as evidence for the buckling and piling of slabs at the boundary between the upper and lower mantle, where a strong contrast in viscosity may exist and cause resistance to the penetration of slabs into the lower mantle. The distribution and character of seismicity reveal, however, that slabs undergo vertical extension in the upper mantle and compression near the transition zone. In this paper, we demonstrate that during the subduction process, the shape of low viscosity slabs (1 to 100 times more viscous than the surrounding mantle) evolves toward an inverted plume shape that we coin jellyfish. Results of a 3D numerical model show that the leading tip of slabs deform toward a rounded head skirted by lateral tentacles that emerge from the sides of the jellyfish head. The head is linked to the body of the subducting slab by a thin tail. A complete parametric study reveals that subducting slabs may achieve a variety of shapes, in good agreement with the diversity of natural slab shapes evidenced by seismic tomography. Our work also suggests that the slab to mantle viscosity ratio in the Earth is most likely to be lower than 100. However, the sensitivity of slab shapes to upper and lower mantle viscosities and densities, which remain poorly constrained by independent evidence, precludes any systematic deciphering of the observations.
Geomorphic Response to Flat Slab Subduction along the Eastern Foothills of the Colombian Andes
NASA Astrophysics Data System (ADS)
Veloza, G.; Taylor, M. H.; Gosse, J. C.; Mora, A.; Becker, T. W.
2013-12-01
It is thought that in northwest South America flat slab subduction plays a key role in the recent development of the eastern Colombian Andes. Here we show that the geomorphic response to flat slab subduction is presently occurring >500 km inboard of the subduction zone plate boundary. The Llanos basin located along the eastern edge of the Colombian Andes is experiencing active uplift along the seismically active Cusiana, Yopal, Paz de Ariporo and Tame thrust faults, which we refer to as the Llanos Foothills thrust system (LFTS). The LFTS is comprised of east-directed thrust faults that are listric in geometry with shallowly west-dipping decollements. Locally, actively growing north-south plunging folds are cored by blind thrust faults, and are being incised by antecedent east-flowing streams. Using a combination of field-based observations on the geometry of faulted and folded fluvial terraces, and geochronology from terrestrial cosmogenic nuclides, we show that the fluvial terraces have been uplifted, and locally, incised >200 meters at incision rates exceeding 3 mm/yr. The field observations in combination with earthquakes and geodynamic simulations can be reconciled by flat slab subduction, but it is presently unknown whether the flat slab has a Caribbean or Nazca plate affinity. Different geodynamic scenarios can be tested to understand how the leading edge of the flat slab interacts with the South American craton, and how that interaction controls upper crustal deformation.
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2015-12-01
Observations of seismicity and seismic tomography provide constraints on the geometry of slabs within mantle, while compression/tension axis derived from moment tensor solutions provide constraints on the internal deformation of slabs. However, since these observations provide only a somewhat blurred or incomplete snapshot of the slab in time, it is difficult to directly relate these observations to the evolution of the slab geometry and the forces acting on and within the slab. In contrast, plate tectonic reconstructions provide time-dependent constraints on the surface motion of plates and the trench at subduction zones, which are related to the dynamical evolution of the slab. We use 2D geodynamical simulations of subduction to explore the relationship between dynamical process within the deforming slab and the observations of surface plate motion and the state-of-stress in slabs. Specifically we utilize models that include the extended Boussinesq approximation (shear heating and latent heat terms in the energy equation), a layered lithosphere with pyrolite, harzburgite and basalt/eclogite, compositionally-dependent phase transitions, and a composite rheology with yielding. The models employ a weak crustal layer that decouples the overriding and subducting plates and allows for dynamically determined trench motion. Here we show that, 1) multiple phase transitions increase slab folding, 2) ridge push significantly increases trench retreat, and 3) strength of the weak crustal layer influences slab detachment. Compared to past studies a more realistic treatment of the phase transitions makes trench retreat more difficult to generate: a weaker plate may encourage slab retreat but detaches once the slab tip crosses into the transition zone due to the rapid increase in slab density. As suggested by previous studies, slab folding within the transition zone changes the direction of forces on the slab and causes periodic changes from trench retreat to trench advance. We
Vertical slab sinking and westward subduction offshore of Mesozoic North America
NASA Astrophysics Data System (ADS)
Sigloch, Karin; Mihalynuk, Mitchell G.
2013-04-01
Subducted slabs in the mantle, as imaged by seismic tomography, preserve a record of ancient subduction zones. Ongoing debate concerns how direct this link is. How long ago did each parcel of slab subduct, and where was the trench located relative to the imaged slab position? Resolving these questions will benefit paleogeographic reconstructions, and restrict the range of plausible rheologies for mantle convection simulations. We investigate one of the largest and best-constrained Mesozoic slab complexes, the "Farallon" in the transition zone and lower mantle beneath North America. We quantitatively integrate observations from whole-mantle P-wave tomography, global plate reconstructions, and land geological evidence from the North American Cordillera. These three data sets permit us to test the simplest conceivable hypothesis for linking slabs to paleo-trenches: that each parcel of slab sank only vertically shortly after entering the trench That is, we test whether within the limits of tomographic resolution, all slab material lies directly below the location where it subducted beneath its corresponding arc. Crucially and in contrast to previous studies, we do not accept or impose an Andean-style west coast trench (Farallon-beneath-continent subduction) since Jurassic times, as this scenario is inconsistent with many geological observations. Slab geometry alone suggests that trenches started out as intra-oceanic because tomography images massive, linear slab "walls" in the lower mantle, extending almost vertically from about 800 km to 2000+ km depth. Such steep geometries would be expected from slabs sinking vertically beneath trenches that were quasi-stationary over many tens of millions of years. Intra-oceanic trenches west of Mesozoic North America could have been stationary, whereas a coastal Farallon trench could not, because the continent moved westward continuously as the Atlantic opened. Overlap of North American west-coast positions, as reconstructed in a
Synthetic Hydrogen Spectra of Oscillating Prominence Slabs Immersed in the Solar Corona
NASA Astrophysics Data System (ADS)
Zapiór, M.; Oliver, R.; Ballester, J. L.; Heinzel, P.
2016-08-01
We study the behavior of Hα and Hβ spectral lines and their spectral indicators in an oscillating solar prominence slab surrounded by the solar corona, using an MHD model combined with a 1D radiative transfer code taken in the line of sight perpendicular to the slab. We calculate the time variation of the Doppler shift, half-width, and maximum intensity of the Hα and Hβ spectral lines for different modes of oscillation. We find a non-sinusoidal time dependence of some spectral parameters with time. Because Hα and Hβ spectral indicators have different behavior for different modes, caused by differing optical depths of formation and different plasma parameter variations in time and along the slab, they may be used for prominence seismology, especially to derive the internal velocity field in prominences.
Looking east inside of the 44" slab mill building at ...
Looking east inside of the 44" slab mill building at the red hot slabs being conveyed to the hot beds. - U.S. Steel Edgar Thomson Works, 44" Slab Mill, Along Monongahela River, Braddock, Allegheny County, PA
Ionin, Andrei A; Kozlov, A Yu; Seleznev, L V; Sinitsyn, D V
2009-03-31
A compact capacitive transverse RF-discharge-pumped slab CO laser with cryogenically cooled electrodes, which operates both in the cw and repetitively pulsed regimes, is fabricated. The laser operation is studied in the free running multifrequency regime at the vibrational - rotational transitions of the fundamental (V + 1 {yields} V) vibrational bands of the CO molecule in the spectral region from 5.1 to 5.4 {mu}m. Optimal operation conditions (gas mixture composition and pressure, RF pump parameters) are determined. It is shown that only gas mixtures with a high content of oxygen (up to 20% with respect to the concentration of CO molecules) can be used as an active medium of this laser. It is demonstrated that repetitively pulsed pumping is more efficient compared to cw pumping. In this case, quasi-cw lasing regime can be obtained. The maximum average output power of {approx}12 W was obtained for this laser operating on fundamental bands and its efficiency achieved {approx}14 %. The frequency-selective operation regime of the slab RF-discharge-pumped CO laser was realised at {approx} 100 laser lines in the spectral region from 5.0 to 6.5 {mu}m with the average output power of up to several tens of milliwatts in each line. Lasing at the transitions of the first vibrational overtone (V + 2 {yields} V) of the CO molecule is obtained in the spectral region from 2.5 to 3.9 {mu}m. The average output power of the overtone laser achieved 0.3 W. All the results were obtained without the forced gas mixture exchange in the discharge chamber. Under fixed experimental conditions, repetitively pulsed lasing (with fluctuations of the output characteristics no more than {+-}10 %) was stable for more than an hour. (lasers)
Slab stagnation and detachment under northeast China
NASA Astrophysics Data System (ADS)
Honda, Satoru
2016-03-01
Results of tomography models around the Japanese Islands show the existence of a gap between the horizontally lying (stagnant) slab extending under northeastern China and the fast seismic velocity anomaly in the lower mantle. A simple conversion from the fast velocity anomaly to the low-temperature anomaly shows a similar feature. This feature appears to be inconsistent with the results of numerical simulations on the interaction between the slab and phase transitions with temperature-dependent viscosity. Such numerical models predict a continuous slab throughout the mantle. I extend previous analyses of the tomography model and model calculations to infer the origins of the gap beneath northeastern China. Results of numerical simulations that take the geologic history of the subduction zone into account suggest two possible origins for the gap: (1) the opening of the Japan Sea led to a breaking off of the otherwise continuous subducting slab, or (2) the western edge of the stagnant slab is the previous subducted ridge, which was the plate boundary between the extinct Izanagi and the Pacific plates. Origin (2) suggesting the present horizontally lying slab has accumulated since the ridge subduction, is preferable for explaining the present length of the horizontally lying slab in the upper mantle. Numerical models of origin (1) predict a stagnant slab in the upper mantle that is too short, and a narrow or non-existent gap. Preferred models require rather stronger flow resistance of the 660-km phase change than expected from current estimates of the phase transition property. Future detailed estimates of the amount of the subducted Izanagi plate and the present stagnant slab would be useful to constrain models. A systematic along-arc variation of the slab morphology from the northeast Japan to Kurile arcs is also recognized, and its understanding may constrain the 3D mantle flow there.
On the Viability of Slab Melting
NASA Astrophysics Data System (ADS)
Van Hunen, J.; Bouilhol, P.; Magni, V.; Maunder, B. L.
2014-12-01
Melting subducted mafic crust is commonly assumed to be the main process leading to silicic melts with an adakitic signature, which may form Archaean granitoids and generate early continental crust. Alternatively, melting of the overriding lower mafic crust and near-Moho depth fractional crystallisation of mantle melts can form differentiated magmas with an adakitic signature. Previous work shows how only very young slabs melt through dehydration melting, or depict melting of dry eclogites via water addition from deeper slab dehydration. Alternatively, underplated subducted material via delamination and diapirism may be important in the generation of felsic continental crust. We quantify subduction dehydration and melting reactions in a warm subduction system using a thermo-mechanical subduction model with a thermodynamic database. We find that even young (hot) slabs dehydrate before reaching their solidus, which suppresses any slab dehydration melting and creates significant amounts of mantle wedge melting irrespective of slab age. Significant slab crust melting is only achieved in young slabs via water present melting if metamorphic fluids from the subducted mantle flux through the dry eclogites. These slab melts, however, interfere with massive mantle wedge melting and unlikely to participate in the overriding plate felsic magmatism, unlike the shallower, primitive mantle wedge melts. We also explore the conditions for delaminating the mafic subducted crust. For a wide range of ages, the uppermost part of the subducted slab might delaminate to form compositionally buoyant plumes that rise through the mantle wedge. Thick crust on young slabs (as perhaps representative for a hotter, early Earth) may delaminate entirely and reside in the mantle wedge. Under such conditions, this ponded crust might melts subsequently, forming "adakitic" felsic melts contributing to a significant amount of the overriding plate crustal volumes.
Characterization of Yb:YAG active slab media based on a layered structure with different doping
NASA Astrophysics Data System (ADS)
Lapucci, A.; Ciofini, M.; Esposito, L.; Ferrara, P.; Gizzi, L. A.; Hostaša, J.; Labate, L.; Pirri, A.; Toci, G.; Vannini, M.
2013-05-01
Slabs with non-uniform doping distribution are studied with the aim of reducing thermal deformations in high-energy high-average-power Yb:YAG slab systems. We present a numerical analysis based on Finite Element Mesh (FEM) methods suitable to model non-uniform devices. The thermal variation of the refractive index, the end-faces deformations and the photo-elastic effect have been calculated in order to estimate the total thermal-lens effect. The stress distributions are also obtained. Some results of this numerical approach are compared to experimental thermal lens measurements in a simple geometry for both uniform and structured samples, in order to validate the numerical procedures. Finally we compare numerical simulations for different single- or double-sided pumping and cooling geometries. They show that structured slabs can reduce thermal gradients with respect to uniformly doped means with comparable absorption and geometry. This means a reduction of thermal lens effect and thus an increase of maximum allowed pump power loading. Previous literature reports some work made with structured slabs where higher doping was located in layers with lower pump radiation levels, in order to get a more uniform absorption. Interestingly our modeling indicates that reduced thermal effects are instead obtained when a higher doping is located close to the cooled surfaces.
Exact image theory for the problem of dielectric/magnetic slab
NASA Technical Reports Server (NTRS)
Lindell, I. V.
1987-01-01
Exact image method, recently introduced for the exact solution of electromagnetic field problems involving homogeneous half spaces and microstrip-like geometries, is developed for the problem of homogeneous slab of dielectric and/or magnetic material in free space. Expressions for image sources, creating the exact reflected and transmitted fields, are given and their numerical evaluation is demonstrated. Nonradiating modes, guided by the slab and responsible for the loss of convergence of the image functions, are considered and extracted. The theory allows, for example, an analysis of finite ground planes in microstrip antenna structures.
NASA Astrophysics Data System (ADS)
Breitsprecher, Katrin; Thorkelson, Derek J.
2009-01-01
The Patagonian slab window is a subsurface tectonic feature resulting from subduction of the Nazca-Antarctic spreading-ridge system (Chile Rise) beneath southern South America. The geometry of the slab window had not been rigorously defined, in part because of the complex nature of the history of ridge subduction in the southeast Pacific region, which includes four interrelated spreading-ridge systems since 20 Ma: first, the Nazca-Phoenix ridge beneath South America, then simultaneous subduction of the Nazca-Antarctic and the northern Phoenix-Antarctic spreading-ridge systems beneath South America, and the southern Phoenix-Antarctic spreading-ridge system beneath Antarctica. Spreading-ridge paleo-geographies and rotation poles for all relevant plate pairs (Nazca, Phoenix, Antarctic, South America) are available from 20 Ma onward, and form the mathematical basis of our kinematic reconstruction of the geometry of the Patagonia and Antarctic slab windows through Neogene time. At approximately 18 Ma, the Nazca-Phoenix-Antarctic oceanic (ridge-ridge-ridge) triple junction enters the South American trench; we recognize this condition as an unstable quadruple junction. Heat flow at this junction and for some distance beneath the forearc would be considerably higher than is generally recognized in cases of ridge subduction. From 16 Ma onward, the geometry of the Patagonia slab window developed from the subduction of the trailing arms of the former oceanic triple junction. The majority of the slab window's areal extent and geometry is controlled by the highly oblique (near-parallel) subduction angle of the Nazca-Antarctic ridge system, and by the high contrast in relative convergence rates between these two plates relative to South America. The very slow convergence rate of the Antarctic slab is manifested by the shallow levels achieved by the slab edge (< 45 km); thus no point on the Antarctic slab is sufficiently deep to generate "normal" mantle-derived arc-type magmas
Cretaceous Arctic magmatism: Slab vs. plume? Or slab and plume?
NASA Astrophysics Data System (ADS)
Gottlieb, E. S.; Miller, E. L.; Andronikov, A. V.; Brumley, K.; Mayer, L. A.; Mukasa, S. B.
2010-12-01
Tectonic models for the Cretaceous paleogeographic evolution of the Arctic Ocean and its adjacent landmasses propose that rifting in the Amerasia Basin (AB) began in Jura-Cretaceous time, accompanied by the development of the High Arctic Large Igneous Province (HALIP). During the same timespan, deformation and slab-related magmatism, followed by intra-arc rifting, took place along the Pacific side of what was to become the Arctic Ocean. A compilation and comparison of the ages, characteristics and space-time variation of circum-Arctic magmatism allows for a better understanding of the role of Pacific margin versus Arctic-Atlantic plate tectonics and the role of plume-related magmatism in the origin of the Arctic Ocean. In Jura-Cretaceous time, an arc built upon older terranes overthrust the Arctic continental margins of North America and Eurasia, shedding debris into foreland basins in the Brooks Range, Alaska, across Chukotka, Russia, to the Lena Delta and New Siberian Islands region of the Russian Arctic. These syn-tectonic sediments have some common sources (e.g., ~250-300 Ma magmatic rocks) as determined by U-Pb detrital zircon geochronology. They are as young as Valanginian-Berriasian (~136 Ma, Gradstein et al., 2004) and place a lower limit on the age of formation of the AB. Subsequent intrusions of granitoid plutons, inferred to be ultimately slab-retreat related, form a belt along the far eastern Russian Arctic continental margin onto Seward Peninsula and have yielded a continuous succession of zircon U-Pb ages from ~137-95 Ma (n=28) and a younger suite ~91-82 Ma (n=16). All plutons dated were intruded in an extensional tectonic setting based on their relations to wall-rock deformation. Regional distribution of ages shows a southward migration of the locus of magmatism during Cretaceous time. Basaltic lavas as old as 130 Ma and as young as 80 Ma (40Ar/39Ar)) erupted across the Canadian Arctic Islands, Svalbard and Franz Josef Land and are associated with
NASA Astrophysics Data System (ADS)
Eakin, C. M.; Long, M. D.
2013-05-01
Flat or shallow subduction is a relatively widespread global occurrence but the dynamics remain poorly understood. In particular, the interaction between flat-slabs and the surrounding mantle flow has yet to be studied in detail. In this study measurements of seismic anisotropy are utilized to investigate mantle flow beneath the Peruvian flat-slab segment, the largest present-day region of flat-subduction. A detailed shear wave splitting analysis is conducted at a long-running seismic station (NNA) located near Lima, Peru. Measurements of apparent splitting parameters (fast direction φ and delay time δt) for SKS, ScS and local S phases from 86 events were obtained. Well defined frequency dependence and back-azimuthal variability is observed, indicating the likely presence of complex anisotropy. Forward modeling the observations with two or three layers of anisotropy reveals a likely layer with a trench-parallel fast axis overlying a layer(s) with a more trench-normal fast geometry. In order to further constrain the anisotropic geometry, source-side splitting from events originating within the slab measured at distant stations have been analyzed. Beneath the flat-slab segment, trench-normal fast splitting directions in the sub-slab mantle were found and likely trench-parallel anisotropy within the slab itself. This sub-slab pattern contradicts observations from elsewhere in South America for 'normal' (i.e. more steeply dipping) slab conditions. It is similar, however, to inferences from other shallowly dipping subduction zones around the world. While there is an apparent link between slab dip and the surrounding mantle flow, at least beneath Peru, the precise nature of the relationship requires further investigation.
NASA Technical Reports Server (NTRS)
1945-01-01
Vought F4U-1D Corsair: In February and March of 1945 this Corsair was examined in the NACA's 30 x 60 Full Scale Tunnel at Langley Field. The F4U-1D has rockets mounted on its wings for this test. After installation and during testing, the wings would be lowered to their flight position.
Sausage oscillations of coronal plasma slabs
NASA Astrophysics Data System (ADS)
Hornsey, C.; Nakariakov, V. M.; Fludra, A.
2014-07-01
Context. Sausage oscillations are observed in plasma non-uniformities of the solar corona as axisymmetric perturbations of the non-uniformity. Often, these non-uniformities can be modelled as field-aligned slabs of the density enhancement. Aims: We perform parametric studies of sausage oscillations of plasma slabs, aiming to determine the dependence of the oscillation period on its parameters, and the onset of leaky and trapped regimes of the oscillations. Methods: Slabs with smooth transverse profiles of the density of a zero-beta plasma are perturbed by an impulsive localised perturbation of the sausage symmetry. In particular, the slab can contain an infinitely thin current sheet in its centre. The initial value problem is then solved numerically. The numerical results are subject to spectral analysis. The results are compared with analytical solutions for a slab with a step-function profile and also with sausage oscillations of a plasma cylinder. Results: We established that sausage oscillations in slabs generally have the same properties as in plasma cylinders. In the trapped regime, the sausage oscillation period increases with the increase in the longitudinal wavelength. In the leaky regime, the dependence of the period on the wavelength experiences saturation, and the period becomes independent of the wavelength in the long-wavelength limit. In the leaky regime the period is always longer than in the trapped regime. The sausage oscillation period in a slab is always longer than in a cylinder with the same transverse profile. In slabs with steeper transverse profiles, sausage oscillations have longer periods. The leaky regime occurs at shorter wavelengths in slabs with smoother profiles.
Detecting slab structure beneath the Mediterranean
NASA Astrophysics Data System (ADS)
Miller, Meghan S.; Sun, Daoyuan; Piana Agostinetti, Nicola
2013-04-01
The presence of subducted slabs in the Mediterranean has been well documented with seismic tomography, however, these images, which are produced by smoothed, damped inversions, underestimate the sharpness of the structures. The position and extent of the slabs and the presence possible tears or gaps in the subducted lithosphere are still debated, yet the shape and location these structures are important for kinematic reconstructions and evolution of the entire subduction zone system. Extensive distribution of broadband seismic instrumentation in the Mediterranean (Italian National Seismic Network in Italy and the NSF-PICASSO project in Spain and Morocco) has allowed us to use alternative methodologies to detect the position of the slabs and slab tears beneath the Central and Western Mediterranean. Using S receiver functions we are able to identify S-to-p conversions from the bottom of the subducted slab and a lack of these signals where there are gaps or tears in the slab. We also analyze broadband waveforms for changes in P wave coda from deep (> 300 km depth) local earthquakes. The waveform records for stations in southern Italy and around the Betic-Rif show large amplitude, high frequency (f > 5 Hz) late arrivals with long coda after relatively low-frequency onset. High frequency arrivals are the strongest from events whose raypaths travel within the slab to the stations where they are recorded allowing for mapping of where the subducted material is located within the upper mantle. These two methods, along with inferring the slab position from fast P-wave velocity perturbations in tomography and intermediate depth seismicity, provide additional geophysical evidence to aid in interpretation of the complex, segmented slab structure beneath the Mediterranean.
Andean flat subduction maintained by slab tunneling
NASA Astrophysics Data System (ADS)
Schepers, Gerben; van Hinsbergen, Douwe; Kosters, Martha; Boschman, Lydian; McQuarrie, Nadine; Spakman, Wim
2016-04-01
In two segments below the Andean mountain belt, the Nazca Plate is currently subducting sub-horizontally below South America over a distance of 200-300 km before the plate bends into the mantle. Such flat slab segments have pronounced effects on orogenesis and magmatism and are widely believed to be caused by the downgoing plate resisting subduction due to its local positive buoyancy. In contrast, here we show that flat slabs primarily result from a local resistance against rollback rather than against subduction. From a kinematic reconstruction of the Andean fold-thrust belt we determine up to ~390 km of shortening since ~50 Ma. During this time the South American Plate moved ~1400 km westward relative to the mantle, thus forcing ~1000 km of trench retreat. Importantly, since the 11-12 Ma onset of flat slab formation, ~1000 km of Nazca Plate subduction occurred, much more than the flat slab lengths, which leads to our main finding that the flat slabs, while being initiated by arrival of buoyant material at the trench, are primarily maintained by locally impeded rollback. We suggest that dynamic support of flat subduction comes from the formation of slab tunnels below segments with the most buoyant material. These tunnels trap mantle material until tearing of the tunnel wall provides an escape route. Fast subduction of this tear is followed by a continuous slab and the process can recur during ongoing rollback of the 7000 km wide Nazca slab at segments with the most buoyant subducting material, explaining the regional and transient character of flat slabs. Our study highlights the importance of studying subduction dynamics in absolute plate motion context.
Numerical Modeling of Flat Slab Formation in Central Chile
NASA Astrophysics Data System (ADS)
Manea, M.; Perez-Gussinye, M.; Manea, V.; Fernandez, M.
2009-12-01
CitcomS Version 2.0.1 available from the Computational Infrastructure for Geodynamics (CIG) (http://geodynamics.org), we integrate the models from 20 Ma to the present day. Our preliminary results show that when an oceanic plate subducts beneath thin continental lithosphere, the angle of subduction is steep. However, when the oceanic plate subducts beneath or in the vicinity of a thick and highly viscous continent, shallower slab geometry is developed. Continuing trenchward motion of the thick continental lithosphere results in flat subduction.
Evolution of attached and detached slabs and their associated mantle dynamics
NASA Technical Reports Server (NTRS)
Hsui, Albert T.
1992-01-01
possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.
Dual-mode characteristics of the Buneman instability in a bounded slab plasma
Hong, Woo-Pyo; Jung, Young-Dae
2015-06-15
The dual-mode characteristics of the Buneman instability are investigated in a slab plasma, including the geometric effects. The dual symmetric and anti-symmetric dispersion modes of the Buneman instability are obtained by the plasma dielectric function with the spectral reflection conditions for the slab geometry. The result shows that the magnitudes of the growth rate for the symmetric mode are always greater than those for the anti-symmetric mode. It is also found that the geometric effect suppresses the position of the maximum growth rate for the Buneman instability in bounded slab plasmas since the maximum conditions for the symmetric and anti-symmetric modes of the Buneman instability are given by 0.60
Non-uniform black strings and the critical dimension in the 1/D expansion
NASA Astrophysics Data System (ADS)
Suzuki, Ryotaku; Tanabe, Kentaro
2015-10-01
Non-uniform black strings (NUBS) are studied by the large D effective theory approach. By solving the near-horizon geometry in the 1 /D expansion, we obtain the effective equation for the deformed horizon up to the next-to-next-to-leading order (NNLO) in 1 /D. We also solve the far-zone geometry by the Newtonian approximation. Matching the near and far zones, the thermodynamic variables are computed in the 1 /D expansion. As the result, the large D analysis gives a critical dimension D * ≃ 13 .5 at which the translation-symmetry-breaking phase transition changes between first and second order. This value of D * agrees perfectly, within the precision of the 1 /D expansion, with the result previously obtained by E. Sorkin through the numerical resolution. We also compare our NNLO results for the thermodynamics of NUBS to earlier numerical calculations, and find good agreement within the expected precision.
NASA Astrophysics Data System (ADS)
Sigloch, K.; Mihalynuk, M. G.
2014-12-01
How rapidly slabs sink, which trajectories they follow, and how they deform in the process, presents an inferential challenge to geophysics. Mantle rheologies remain highly uncertain, and seismic tomography can merely offer present-day snapshots of a process defined by temporal evolution. Thus observational constraints on slab sinking have tended to remain non-unique. Subduction zones are complex litho-consumers whose time-variant activity can be reconstructed from geological observations on paleo-arcs, but the association of arcs to their subducted, tomographically imaged lithosphere is iffy. Except for young slabs that can be reliably linked with coeval paleo-arc activity a priori, deeper geological time information cannot be exploited with certainty. As long as slab geometries remain "undated", few constraints on slab sinking behavior and hence mantle rheology can be extracted. Sigloch & Mihalynuk (2013) demonstrated a quantitative method to tighten constraints on slab sinking in the lower mantle by investigating the least ambiguous slab geometries observed. Extremely massive and almost vertical slab walls should have been deposited by vertical sinking beneath (intra-oceanic) trenches that remained stationary for a long time (~100 m.y.). We showed how this hypothesis of vertical sinking can be tested quantitatively and successfully, making only minimal assumptions on mantle rheology, and with proper error propagation for all observations (tomography, plate reconstructions, geology). Here the discussion of sinking trajectories and rates is extended to more challenging geometries. Dipping slabs in the lower mantle, and laterally extensive "stagnant slabs" in the transition zone can also be rendered dateable and trackable by (re-)investigation of their paleo-trenches. We discuss examples and link to recent geodynamic modeling of viscous sheet sinking. Reference: Sigloch K & Mihalynuk MG (2013), Intra-oceanic subduction shaped the assembly of Cordilleran North
Investigation of Subducting Slab beneath Northeastern Taiwan from Sp Converted Phases
NASA Astrophysics Data System (ADS)
Wu, W.
2004-12-01
The converted waves from upper plate boundary of the subducting slab denoted as Sp phase recorded by the Central Weather Bureau Seismographic Network (CWBSN) for events in northern Taiwan provide us an advantageous opportunity to construct the geometry of upper boundary of descending Philippine Sea Plate (PSP). The events were chosen for the region with latitude of 24.4 degree N - 25.2 degree N, longitude of 121.5 degree E - 122.2 degree E, depth from 50 -300 km, and time period of January 1991 to Mary 2003 for magnitude greater than 4.0. We apply the polarization filter analysis, particle motion diagram and theoretical travel time to the digital 3-compoent short period seismograms to identify the conceivable Sp converted phases. The initial velocity model with subducting slab was used for theoretical travel time calculation based on the recent tomographic results. The model is constructed firstly by considering a fixed velocity contrasts among the discontinuities. The dip angle and latitude location of slab are allowed to change to obtain the best travel time to the data. Several tests had been made for the inversion. Considering variance reduction of the all observed data to the constructed models by grid searching technique, variance reduction is not in satisfactory. Further studies by grouping the events in north-south and northwestern-southeastern trend were examined. Although the geometry of subducting slab might be difficult to model due to disadvantage in tectonic setting of northeastern Taiwan, which is located to northern subducting slab, the most possible image of the slab from the clear identified converted phases is sketched.
Photocatalytic, highly hydrophilic porcelain stoneware slabs
NASA Astrophysics Data System (ADS)
Raimondo, M.; Guarini, G.; Zanelli, C.; Marani, F.; Fossa, L.; Dondi, M.
2011-10-01
Photocatalytic, highly hydrophilic industrial porcelain stoneware large slabs were realized by deposition of nanostructured TiO2 coatings. Different surface finishing and experimental conditions were considered in order to assess the industrial feasibility. Photocatalytic and wetting behaviour of functionalized slabs mainly depends on surface phase composition in terms of anatase/rutile ratio, this involving - as a key issue - the deposition of TiO2 on industrially sintered products with an additional annealing step to strengthen coatings' performances and durability.
ERIC Educational Resources Information Center
Desseyn, H. O.; And Others
1985-01-01
Compares linear-nonlinear and planar-nonplanar geometry through the valence-shell electron pairs repulsion (V.S.E.P.R.), Mulliken-Walsh, and electrostatic force theories. Indicates that although the V.S.E.P.R. theory has more advantages for elementary courses, an explanation of the best features of the different theories offers students a better…
Development of 1D Liner Compression Code for IDL
NASA Astrophysics Data System (ADS)
Shimazu, Akihisa; Slough, John; Pancotti, Anthony
2015-11-01
A 1D liner compression code is developed to model liner implosion dynamics in the Inductively Driven Liner Experiment (IDL) where FRC plasmoid is compressed via inductively-driven metal liners. The driver circuit, magnetic field, joule heating, and liner dynamics calculations are performed at each time step in sequence to couple these effects in the code. To obtain more realistic magnetic field results for a given drive coil geometry, 2D and 3D effects are incorporated into the 1D field calculation through use of correction factor table lookup approach. Commercial low-frequency electromagnetic fields solver, ANSYS Maxwell 3D, is used to solve the magnetic field profile for static liner condition at various liner radius in order to derive correction factors for the 1D field calculation in the code. The liner dynamics results from the code is verified to be in good agreement with the results from commercial explicit dynamics solver, ANSYS Explicit Dynamics, and previous liner experiment. The developed code is used to optimize the capacitor bank and driver coil design for better energy transfer and coupling. FRC gain calculations are also performed using the liner compression data from the code for the conceptual design of the reactor sized system for fusion energy gains.
Three-dimensional necking during viscous slab detachment
NASA Astrophysics Data System (ADS)
Tscharner, M.; Schmalholz, S. M.; Duretz, T.
2014-06-01
We study the three-dimensional (3-D) deformation during detachment of a lithospheric slab with simple numerical models using the finite element method. An initially vertical layer of power law viscous fluid mimics the slab and is surrounded by a linear or power law viscous fluid representing asthenospheric mantle. We quantify the impact of slab size and shape (symmetric/asymmetric) on slab detachment and identify two processes that control the lateral (i.e., along-trench) slab deformation: (1) the horizontal deflection of the lateral, vertical slab sides (> 100 km with velocities up to 16 mm/yr) and (2) the propagation of localized thinning (necking) inside the slab (with velocities >9 cm/yr). The lateral propagation velocity is approximately constant during slab detachment. Larger slabs (here wider than approximately 300 km) detach with rates similar to those predicted by 2-D models, whereas smaller slabs detach slower. Implications for geodynamic processes and interpretations of seismic tomography are discussed.
NASA Astrophysics Data System (ADS)
Faccenda, M.; Burlini, L.; Gerya, T.; Mancktelow, N.
2012-12-01
In this presentation we summarize the results of a project started in 2007 from a brilliant intuition of Luigi Burlini that suggested an additional anisotropy source for the interpretation of seismic anisotropy patterns observed at subduction zones. Such an anisotropic body located in the upper part of the slab would result from the hydration of the oceanic plate at the trench outer-rise. The natural continuation of the project was to understand the mechanical processes behind slab hydration and the fluid flow patterns established during slab dehydration. In both cases, we found that tectonic pressure gradients due to the bending and unbending of the subducting oceanic plate are fundamental in driving fluid flow. This last part of the project led to the other two chapters of the final trilogy about the long route of water in the slab. This trilogy is here described in detail and a chronologically ordered series of events presented below. The first episode is related to slab hydration occurring during bending at the trench-outer rise. Here, fluids are driven downward along active normal faults by bending-related, sub-hydrostatic pressure gradients. Water can percolate down to 15-20 km below the seafloor, triggering hydrothermal reactions and the formation of hydrous minerals. This results in an elongated pattern of mostly trench-dipping hydrated faults with a strike parallel to the trench and whose orientation below the forearc becomes sub-vertical. The second episode is related to the geophysical implications of a hydrated slab below the forearc. Indeed, both the subvertical layering of closely spaced hydrated and dry levels (SPO) and the syn-deformational, fault-parallel alignment of highly anisotropic minerals (CPO) such as serpentine and talc may contribute to the SKS splitting patterns observed in the forearc. We suggest that the upper part of the slab may have a strong seismic anisotropy that can be approximated by a transverse isotropic body with a sub
NASA Astrophysics Data System (ADS)
Cembranos, J. A. R.; Dobado, A.; Maroto, A. L.
Extra-dimensional theories contain additional degrees of freedom related to the geometry of the extra space which can be interpreted as new particles. Such theories allow to reformulate most of the fundamental problems of physics from a completely different point of view. In this essay, we concentrate on the brane fluctuations which are present in brane-worlds, and how such oscillations of the own space-time geometry along curved extra dimensions can help to resolve the Universe missing mass problem. The energy scales involved in these models are low compared to the Planck scale, and this means that some of the brane fluctuations distinctive signals could be detected in future colliders and in direct or indirect dark matter searches.
Contribution of Elasticity in Slab Bending
NASA Astrophysics Data System (ADS)
Fourel, L.; Goes, S. D. B.; Morra, G.
2014-12-01
Previous studies have shown that plate rheology exerts a dominant control on the shape and velocity of subducting plates. Here, we perform a systematic investigation of the, often disregarded, role of elasticity in slab bending at the trench, using simple, yet fully dynamic, set of 2.5D models where an elastic, visco-elastic or visco-elasto-plastic plate subducts freely into a purely viscous mantle. We derive a scaling relationship between the bending radius of visco-elastic slabs and the Deborah number, De, which is the ratio of Maxwell time over deformation time. We show that De controls the ratio of elastically stored energy over viscously dissipated energy and find that at De exceeding 10-2, it requires substantially less energy to bend a visco-elastic slab to the same shape as a purely viscous slab with the same viscosity (90% less for De=0.1). Elastically stored energy at higher De facilitates slab unbending and hence favours retreating modes of subduction, while trench advance only occurs for some cases with De<10-2. We use our scaling relation to estimate apparent Deborah numbers, Deapp, from a global compilation of subduction-zone parameters. Values range from 10-3 to >1, where most zones have low Deapp<10-2, but a few young plates have Deapp>0.1. Slabs with Deapp ≤ 10-2 either have very low viscosities, ≤10 times mantle viscosity, or they may be yielding, in which case our apparent Deborah number may underestimate actual De by up to an order of magnitude. If a significant portion of the low Deapp slabs yield, then elastically stored energy may actually be important in quite a large number of subduction zones. Interestingly, increasing Deapp correlates with increasing proportion of larger seismic events (b-value) in both instrumental and historic catalogues, indicating that increased contribution of elasticity may facilitate rupture in larger, less frequent earthquakes.
Cracking behavior of structural slab bridge decks
NASA Astrophysics Data System (ADS)
Baah, Prince
Bridge deck cracking is a common problem throughout the United States, and it affects the durability and service life of concrete bridges. Several departments of transportation (DOTs) in the United States prefer using continuous three-span solid structural slab bridges without stringers over typical four-lane highways. Recent inspections of such bridges in Ohio revealed cracks as wide as 0.125 in. These measured crack widths are more than ten times the maximum limit recommended in ACI 224R-01 for bridge decks exposed to de-icing salts. Measurements using digital image correlation revealed that the cracks widened under truck loading, and in some cases, the cracks did not fully close after unloading. This dissertation includes details of an experimental investigation of the cracking behavior of structural concrete. Prism tests revealed that the concrete with epoxy-coated bars (ECB) develops the first crack at smaller loads, and develops larger crack widths compared to the corresponding specimens with uncoated (black) bars. Slab tests revealed that the slabs with longitudinal ECB developed first crack at smaller loads, exhibited wider cracks and a larger number of cracks, and failed at smaller ultimate loads compared to the corresponding test slabs with black bars. To develop a preventive measure, slabs with basalt and polypropylene fiber reinforced concrete were also included in the test program. These test slabs exhibited higher cracking loads, smaller crack widths, and higher ultimate loads at failure compared to the corresponding slab specimens without fibers. Merely satisfying the reinforcement spacing requirements given in AASHTO or ACI 318-11 is not adequate to limit cracking below the ACI 224R-01 recommended maximum limit, even though all the relevant design requirements are otherwise met. Addition of fiber to concrete without changing any steel reinforcing details is expected to reduce the severity and extent of cracking in reinforced concrete bridge decks.
Accidents due to falls from roof slabs.
Rudelli, Bruno Alves; Silva, Marcelo Valerio Alabarce da; Akkari, Miguel; Santili, Claudio
2013-01-01
CONTEXT AND OBJECTIVE Falls from the roof slabs of houses are accidents of high potential severity that occur in large Brazilian cities and often affect children and adolescents. The aims of this study were to characterize the factors that predispose towards this type of fall involving children and adolescents, quantify the severity of associated lesions and suggest preventive measures. DESIGN AND SETTING Descriptive observational prospective longitudinal study in two hospitals in the metropolitan region of São Paulo. METHODS Data were collected from 29 cases of falls from roof slabs involving children and adolescents between October 2008 and October 2009. RESULTS Cases involving males were more prevalent, accounting for 84%. The predominant age group was schoolchildren (7 to 12 years old; 44%). Leisure activities were most frequently being practiced on the roof slab at the time of the fall (86%), and flying a kite was the most prevalent game (37.9%). In 72% of the cases, the children were unaccompanied by an adult responsible for them. Severe conditions such as multiple trauma and traumatic brain injuries resulted from 79% of the accidents. CONCLUSION Falls from roof slabs are accidents of high potential severity, and preventive measures aimed towards informing parents and guardians about the dangers and risk factors associated with this type of accident are needed, along with physical protective measures, such as low walls around the slab and gates with locks to restrict free access to these places. PMID:23903263
NASA Astrophysics Data System (ADS)
Tarlow, S.; Tan, E.; Billen, M. I.
2015-12-01
At the Ryukyu subduction zone, seismic anisotropy observations suggest that there may be strong trench-parallel flow within the mantle wedge driven by complex 3D slab geometry. However, previous simulations have either failed to account for 3D flow or used the infinite strain axis (ISA) approximation for LPO, which is known to be inaccurate in complex flow fields. Additionally, both the slab depth and shape of the Ryukyu slab are contentious. Development of strong trench-parallel flow requires low viscosity to decouple the mantle wedge from entrainment by the sinking slab. Therefore, understanding the relationship between seismic anisotropy and the accompanying flow field will better constrain the material and dynamic properties of the mantle near subduction zones. In this study, we integrate a kinematic model for calculation of LPO (D-Rex) into a buoyancy-driven, instantaneous 3D flow simulation (ASPECT), using composite non-Newtonian rheology to investigate the dependence of LPO on slab geometry and depth at the Ryukyu Trench. To incorporate the 3D flow effects, the trench and slab extends from the southern tip of Japan to the western edge of Taiwan and the model region is approximately 1/4 of a spherical shell extending from the surface to the core-mantle boundary. In the southern-most region we vary the slab depth and shape to test for the effects of the uncertainties in the observations. We also investigate the effect of adding locally hydrated regions above the slab that affect both the mantle rheology and development of LPO through the consequent changes in mantle flow and dominate (weakest) slip system. We characterize how changes in the simulation conditions affect the LPO within the mantle wedge, subducting slab and sub-slab mantle and relate these to surface observations of seismic anisotropy.
NASA Astrophysics Data System (ADS)
Bernal-Olaya, R.; Mann, P.; Vargas, C. A.; Koulakov, I.
2013-12-01
We define the length and geometry of eastward and southeastward-subducting slabs beneath northwestern South America in Colombia using ~100,000 earthquake events recorded by the Colombian National Seismic Network from 1993 to 2012. Methods include: hypocenter relocation, compilation of focal mechanisms, and P and S wave tomographic calculations performed using LOTOS and Seisan. The margins of Colombia include four distinct subduction zones based on slab dip: 1) in northern Colombia, 12-16-km-thick oceanic crust subducts at a modern GPS rate of 20 mm/yr in a direction of 110 degrees at a shallow angle of 8 degrees; as a result of its low dip, Pliocene-Pleistocene volcanic rocks are present 400 km from the frontal thrust; magmatic arc migration to the east records 800 km of subduction since 58 Ma ago (Paleocene) with shallow subduction of the Caribbean oceanic plateau starting ~24-33 Ma (Miocene); at depths of 90-150 km, the slab exhibits a negative velocity anomaly we associate with pervasive fracturing; 2) in the central Colombia-Panama area, we define an area of 30-km-thick crust of the Panama arc colliding/subducting at a modern 30/mm in a direction of 95 degrees; the length of this slab shows subduction/collision initiated after 20 Ma (Middle Miocene); we call this feature the Panama indenter since it has produced a V-shaped indentation of the Colombian margin and responsible for widespread crustal deformation and topographic uplift in Colombia; an incipient subduction area is forming near the Panama border with intermediate earthquakes at an eastward dip of 70 degrees to depths of ~150 km; this zone is not visible on tomographic images; 3) a 250-km-wide zone of Miocene oceanic crust of the Nazca plate flanking the Panama indenter subducts at a rate of 25 mm/yr in a direction of 55 degrees and at a normal dip of 40 degrees; the length of this slab suggests subduction began at ~5 Ma; 4) the Caldas tear defines a major dip change to the south where a 35 degrees
1D ferrimagnetism in homometallic chains
NASA Astrophysics Data System (ADS)
Coronado, E.; Gómez-García, C. J.; Borrás-Almenar, J. J.
1990-05-01
The magnetic properties of the cobalt zigzag chain Co(bpy)(NCS)2 (bpy=2,2'-bipyridine) are discussed on the basis of an Ising-chain model that takes into account alternating Landé factors. It is emphasized, for the first time, that a homometallic chain containing only one type of site can give rise to a 1D ferrimagneticlike behavior.
DESIGN PACKAGE 1D SYSTEM SAFETY ANALYSIS
L.R. Eisler
1995-02-02
The purpose of this analysis is to systematically identify and evaluate hazards related to the Yucca Mountain Project Exploratory Studies Facility (ESF) Design Package 1D, Surface Facilities, (for a list of design items included in the package 1D system safety analysis see section 3). This process is an integral part of the systems engineering process; whereby safety is considered during planning, design, testing, and construction. A largely qualitative approach was used since a radiological System Safety analysis is not required. The risk assessment in this analysis characterizes the accident scenarios associated with the Design Package 1D structures/systems/components in terms of relative risk and includes recommendations for mitigating all identified risks. The priority for recommending and implementing mitigation control features is: (1) Incorporate measures to reduce risks and hazards into the structure/system/component (S/S/C) design, (2) add safety devices and capabilities to the designs that reduce risk, (3) provide devices that detect and warn personnel of hazardous conditions, and (4) develop procedures and conduct training to increase worker awareness of potential hazards, on methods to reduce exposure to hazards, and on the actions required to avoid accidents or correct hazardous conditions. The scope of this analysis is limited to the Design Package 1D structures/systems/components (S/S/Cs) during normal operations excluding hazards occurring during maintenance and ''off normal'' operations.
Fluorescence losses from Yb:YAG slab lasers.
Chen, Ying; Rapaport, Alexandra; Chung, Te-yuan; Chen, Bin; Bass, Michael
2003-12-20
We report on the distribution of fluorescence that can be emitted through the surfaces of a ytterbium-doped yttrium aluminum garnet (Yb:YAG) slab-shaped high-power solid-state laser. Slab shapes considered include parallel or antiparallel Brewster endfaced slabs and rectangular parallelepiped slabs. We treat cases in which all the faces of these slabs are in air, or with water or another coating on the largest faces. The fraction of the fluorescence emitted through each face, its distribution over that face, and the directions in which it travels are shown to be important to the design of high-power slab lasers. PMID:14717292
NASA Astrophysics Data System (ADS)
Horton, B. K.; Fuentes, F.; McKenzie, N. R.; Constenius, K. N.; Alvarado, P. M.
2014-12-01
Debate persists over the effects of flat-slab subduction on the kinematics of overriding plate deformation and the evolution of retroarc sedimentary basins. In western Argentina, major spatial and temporal variations in the geometry of the subducting Nazca slab since ~15 Ma provide opportunities to evaluate the late Cenozoic response of the Andean fold-thrust belt and foreland basin to subhorizontal subduction. Preliminary results from several structural and sedimentary transects spanning the frontal thrust belt and foreland basin system between 31°S and 35°S reveal Oligocene-middle Miocene hinterland exhumation during normal-slab subduction followed thereafter by progressive slab shallowing with initial rapid cratonward propagation of ramp-flat thrust structures (prior to basement-involved foreland uplifts) and accompanying wholesale exhumation and recycling of the early Andean foreland basin (rather than regional dynamic subsidence). Detrital zircon U-Pb geochronologic data prove instrumental for revealing shifts in thrust-belt exhumation, defining depositional ages within the foreland basin, and constraining the timing of activity along frontal thrust structures. In both the San Juan (31-32°S) and Malargüe (34-35°S) segments of the fold-thrust belt, geochronological results for volcaniclastic sandstones and syndeformational growth strata are consistent with a major eastward advance in shortening at 12-9 Ma. This episode of rapid thrust propagation precedes the reported timing of Sierras Pampeanas basement-involved foreland uplifts and encompasses modern regions of both normal- and flat-slab subduction, suggesting that processes other than slab dip (such as inherited crustal architecture, critical wedge dynamics, and arc magmatism) are additional regulators of thrust-belt kinematics and foreland basin evolution.
Fatigue of concrete beams and slabs
NASA Astrophysics Data System (ADS)
Roesler, Jeffrey Raphael
Traditionally, simply supported concrete beam (SSB) fatigue results have been used to characterize the fatigue resistance of fully supported concrete slabs (FSS). SSB concrete fatigue tests have been assumed to be equivalent to the fatigue resistance of concrete slabs in the field. The effect specimen size, boundary conditions, and loading configurations have on the fatigue of concrete beams and slabs have not been considered in the design of concrete pavements against fatigue. A laboratory study was undertaken to determine if the fatigue behavior of FSS and SSB were similar. A fully supported beam (FSB) was also tested under repeated loading, since it represented an intermediate specimen size and testing configuration between SSB and FSS. The best way to present fatigue results for all specimens was a stress ratio (S) to number of cycles to failure (N) curve (S-N curve). SSB fatigue behavior was similar to results obtained from the literature. FSB had similar fatigue behavior to SSB. The fatigue curve derived from repeated loading of FSS was 30 percent higher than the SSB fatigue curve. This suggested for a given number of cycles to failure, FSS could take a 30 percent higher bending stress as compared to SSB and FSB. The concrete modulus of rupture from a FSS test configuration was 30 percent greater than the concrete modulus of rupture from a SSB test setup. If the concrete modulus of rupture from a FSS test configuration was used in the slab's stress ratio, the slab's fatigue curve was the same as the SSB and FSB. This meant concrete behaved the same under fatigue loading, irrespective of specimen size and test configuration, as long as the correct concrete modulus of rupture was used in the stress ratio. Strain gages, attached to all specimens tested, indicated cracking in concrete occurred in a narrow band. Regions of high plastic strain were found in the plane of cracking, while adjacent areas experienced decreases in strain levels with cracking. Strain
NASA Astrophysics Data System (ADS)
Bina, Craig; Cizkova, Hana
2014-05-01
Subducting slabs may exhibit buckling instabilities and consequent folding behavior in the mantle transition zone for various combinations of dynamical parameters, accompanied by temporal variations in dip angle, plate velocity, and trench retreat. Parameters governing such behavior include both viscous forces (slab and mantle rheology) and buoyancy forces (slab thermal structure and mineral phase relations). 2D numerical experiments show that many parameter sets lead to slab deflection at the base of the transition zone, typically accompanied by quasi-periodic oscillations (consistent with previous scaling analyses) in largely anticorrelated plate and rollback velocities, resulting in undulating stagnant slabs as buckle folds accumulate subhorizontally atop the lower mantle. Slab interactions with mantle phase transitions are important components of this process (Bina and Kawakatsu, 2010; Čížková and Bina, 2013). For terrestrial parameter sets, trench retreat is found to be nearly ubiquitous, and trench advance is quite rare - due to both rheological structure and ridge-push effects (Čížková and Bina, 2013). Recent analyses of global plate motions indicate that significant trench advance is also rare on Earth, being largely restricted to the Izu-Bonin arc (Matthews et al., 2013). Consequently, we explore the conditions necessary for terrestrial trench advance through dynamical models involving the unusual geometry associated with the Philippine Sea region. Detailed images of buckled stagnant slabs are difficult to resolve due to smoothing effects inherent in seismic tomography, but velocity structures computed for compositionally layered slabs, using laboratory data on relevant mineral assemblages, can be spatially low-pass filtered for comparison with tomographic images of corresponding resolution. When applied to P-wave velocity anomalies from stagnant slab material beneath northeast China, model slabs which undulate due to compound buckling fit
Quadratic Finite Element Method for 1D Deterministic Transport
Tolar, Jr., D R; Ferguson, J M
2004-01-06
In the discrete ordinates, or SN, numerical solution of the transport equation, both the spatial ({und r}) and angular ({und {Omega}}) dependences on the angular flux {psi}{und r},{und {Omega}}are modeled discretely. While significant effort has been devoted toward improving the spatial discretization of the angular flux, we focus on improving the angular discretization of {psi}{und r},{und {Omega}}. Specifically, we employ a Petrov-Galerkin quadratic finite element approximation for the differencing of the angular variable ({mu}) in developing the one-dimensional (1D) spherical geometry S{sub N} equations. We develop an algorithm that shows faster convergence with angular resolution than conventional S{sub N} algorithms.
Effective theory of black holes in the 1/D expansion
NASA Astrophysics Data System (ADS)
Emparan, Roberto; Shiromizu, Tetsuya; Suzuki, Ryotaku; Tanabe, Kentaro; Tanaka, Takahiro
2015-06-01
The gravitational field of a black hole is strongly localized near its horizon when the number of dimensions D is very large. In this limit, we can effectively replace the black hole with a surface in a background geometry (e.g. Minkowski or Anti-deSitter space). The Einstein equations determine the effective equations that this `black hole surface' (or membrane) must satisfy. We obtain them up to next-to-leading order in 1/ D for static black holes of the Einstein-(A)dS theory. To leading order, and also to next order in Minkowski backgrounds, the equations of the effective theory are the same as soap-film equations, possibly up to a redshift factor. In particular, the Schwarzschild black hole is recovered as a spherical soap bubble. Less trivially, we find solutions for `black droplets', i.e. black holes localized at the boundary of AdS, and for non-uniform black strings.
NASA Astrophysics Data System (ADS)
Jones, C. H.; Reeg, H.; Zandt, G.; Gilbert, H. J.; Owens, T. J.; Stachnik, J. C.
2013-12-01
The high-wavespeed Isabella anomaly in the upper mantle at the southwestern edge of the Sierra Nevada has been interpreted as a convective (Rayleigh-Taylor) instability (or drip), the remains of a fragment of the Farallon plate (a slab) or the product of delamination of lithosphere from the east or south. P-wave tomography using 29,186 picks from portable deployments from 1988, 1997, and the SNEP deployment of 2005-7 and surrounding TA and permanent broadband stations was run from a variety of starting models. Some models started from a 1-D earth model, some from the Moschetti et al. (JGR 2010) 3-D S-wave model, some from the Gilbert et al. (Geosphere, 2012) 3-D SV-wave model. S-wave models were converted to P-wavespeeds using the regression suggested by Brocher (BSSA 2005). In some cases the upper levels of the 3-D models were fixed and only wavespeeds below ~55 km were allowed to change. Because of the relatively poor vertical resolution of the teleseismic body wave tomography and the bias towards a minimal model variance, the resulting images in the upper lithosphere vary considerably between models, producing results resembling slabs, drips, and delaminations for the Isabella anomaly. Thus the shape of the shallow part of the anomaly does not reliably determine its origin. Deeper (100-250 km) parts of the anomaly are consistent between different inversions with a ≥4% fast body dipping 60-70° east. Vertical integrals of wavespeed anomaly are relatively insensitive to the models; use of such an integral over the body from 95 to 245 km depth yields an equivalent volume of 7 × 1 x106 km3 at a mean anomaly of +1%. We expect from geological considerations that a volume on average 5% fast of 0.9-1.6 x 106 km3 was removed from under the southern Sierra, equivalent to a volume of 4.4-8.4 x 106 km3 at 1% fast, in close agreement with the equivalent volume of the Isabella anomaly. We prefer some kind of 3-D convective removal for the Isabella anomaly, noting that the
Validation of 3D/1D Analysis of ICRF Antennas
NASA Astrophysics Data System (ADS)
Milanesio, D.; Lancellotti, V.; Kyrytsya, V.; Maggiora, R.; Vecchi, G.; Parisot, A.; Wukitch, S. J.
2004-11-01
An innovative tool has been realized for the 3D/1D simulation of Ion Cyclotron Radio Frequency (ICRF), i.e. accounting for antennas in a realistic 3D geometry and with an accurate 1D plasma model. The approach to the problem is based on an integral-equation formulation for the self-consistent evaluation of the current distribution on the conductors. The environment has been subdivided in two coupled region: the plasma region and the vacuum region. The two problems are linked by means of electromagnetic current distribution on the aperture between the two regions. The plasma enters the formalism via a surface impedance matrix for this reason any plasma model can be used. The source term directly models the TEM mode of the coax feeding the antenna and the current in the coax is determined self-consistently, giving the input impedance/admittance of the antenna itself. The suite, called TOPICA, has been used in the design of various ICRF antennas and also for the performance prediction of the ALCATOR C-MOD D and E antenna. An extensive set of comparisons between measured and simulated antenna parameters during ALCATOR C-MOD operation will be presented.
NASA Astrophysics Data System (ADS)
Wu, J. E.; Suppe, J.; Kanda, R. V.
2012-12-01
Subducted slabs were mapped in the mantle under East Asia using MITP08 global seismic tomography (Li et al., 2008), Benioff zone seismicities and published local tomography. 3D gridded slab surfaces were constructed by manually picking and correlating the midpoint of fast seismic anomalies along variable cross-section orientations. The mapped slabs were structurally 'unfolded' and restored to the spherical Earth surface to assess their pre-subduction geometries. Gplates software was used to constrain plate tectonic reconstructions using the unfolded slabs. The unfolded SE Asia upper mantle slabs reveal a 'picture puzzle' fit along their edges that suggests a larger NE Indo-Australian ocean once existed that included the Philippine Sea, Molucca Sea and Celebes Sea. Deeper lower mantle detached slabs indicate an early to mid-Cenozoic 'East Asia Sea' between east Sundaland and the Pacific that stretched from the Ryukyu Islands north of present-day Taiwan southward to Sulawesi. The unfolded slab constraints produced gap and overlap incompatibilities when used in published plate tectonic reconstructions. Here a plate tectonic reconstruction incorporating the unfolded slab constraints is proposed that has the Philippine Sea, Molucca Sea and Celebes Sea clustered at the northern margin of Australia during the early Cenozoic. At the mid-Cenozoic these plates moved NNE with 'Australia-like' plate motions and overrode the 'East Asia Sea'. Plate motions were accommodated by N-S transforms at the eastern margin of Sundaland. Between 25 to 15 Ma the Philippine Sea, Molucca Sea and Celebes Sea plates were fragmented from the greater Indo-Australian ocean. The Philippine Sea was captured by the Pacific plate and now has Pacific-like westward motions.
NASA Astrophysics Data System (ADS)
Eakin, Caroline M.; Long, Maureen D.
2013-09-01
Flat or shallow subduction is a relatively widespread global occurrence, but the dynamics remain poorly understood. In particular, the interaction between flat slabs and the surrounding mantle flow has yet to be studied in detail. Here we present measurements of seismic anisotropy to investigate mantle flow beneath the Peruvian flat-slab segment, the largest present-day region of flat subduction. We conduct a detailed shear wave splitting analysis at a long-running seismic station (NNA) located near Lima, Peru. We present measurements of apparent splitting parameters (fast direction φ and delay time δt) for SKS, ScS, and local S phases from 80 events. We observe well-defined frequency dependence and backazimuthal variability, indicating the likely presence of complex anisotropy. Forward modeling the observations with two or three layers of anisotropy reveals a likely layer with a trench-normal fast direction underlying a layer with a more trench-oblique (to trench-subparallel) fast direction. In order to further constrain the anisotropic geometry, we analyzed the source-side splitting from events originating within the slab measured at distant stations. Beneath the flat-slab segment, we found trench-normal fast splitting directions in the subslab mantle, while within the dipping portion of the slab further to the east, likely trench-subparallel anisotropy within the slab itself. This subslab pattern contradicts observations from elsewhere in South America for "normal" (i.e., more steeply dipping) slab conditions. It is similar, however, to inferences from other shallowly dipping subduction zones around the world. While there is an apparent link between slab dip and the surrounding mantle flow, at least beneath Peru, the precise nature of the relationship remains to be clarified.
NASA Astrophysics Data System (ADS)
Jadamec, M. A.; Fischer, K. M.
2014-12-01
We present a series of three-dimensional (3D), high-resolution, end-member tectonic configurations of the Central American plate system and use these to solve for the 3D viscous mantle flow and surface plate motions. The 3D geodynamic models test the relative control of the viscosity structure (Newtonian versus Composite), subducting plate geometry (continuous slab versus Cocos-Nazca slab gap), and overriding plate thickness (uniform versus laterally variable) on the predicted motion of the Cocos and Nazca plates and the slab-induced 3D flow field in the upper mantle. Models using the composite viscosity formulation result in increased surface plate motions, which better fit the observed motion of the Cocos and Nazca plates. This is particularly significant because these 3D regional models contain the entire Cocos plate, suggesting the importance of the non-linear rheology in models that aim to predict surface plate motions. Faster flow velocities occur in models using the composite viscosity due to the decreased resistance to subduction and reduced viscous support of the slab as the mantle surrounding the slab undergoes non-linear weakening. A zone of partial decoupling between the uppermost mantle and lithosphere, thus, naturally develops due to the composite viscosity formulation. Models that include a gap between the Cocos and Nazca slabs better fit the mantle flow pathways interpreted from the geochemical signatures, as material is brought from beneath the Cocos plate around the slab edge and northward into the mantle wedge beneath Central America. The mantle-lithosphere decoupling is enhanced in models with the slab gap, wherein the mantle flow field contains both counter-clockwise toroidal flow around the Cocos slab edge and clockwise toroidal flow around the northern Nazca slab edge, both of which are non-parallel to surface motions. The models also demonstrate that overriding plate thickness places a control on both the predicted surface motion and
Preparation of 1D nanostructures using biomolecules
NASA Astrophysics Data System (ADS)
Pruneanu, Stela; Olenic, Liliana; Barbu Tudoran, Lucian; Kacso, Irina; Farha Al-Said, Said A.; Hassanien, Reda; Houlton, Andrew; Horrocks, Benjamin R.
2009-08-01
In this paper we have shown that one-dimensional (1D) particle arrays can be obtained using biomolecules, like DNA or amino-acids. Nano-arrays of silver and gold were prepared in a single-step synthesis, by exploiting the binding abilities of λ-DNA and L-Arginine. The morphology and optical properties of these nanostructures were investigated using AFM, TEM and UV-Vis absorption spectroscopy.
A Simple Vertical Slab Gel Electrophoresis Apparatus.
ERIC Educational Resources Information Center
Carter, J. B.; And Others
1983-01-01
Describes an inexpensive, easily constructed, and safe vertical slab gel kit used routinely for sodium dodecyl sulphate-polyacrylamide gel electrophoresis research and student experiments. Five kits are run from a single transformer. Because toxic solutions are used, students are given plastic gloves and closely supervised during laboratory…
Oceanic slab melting and mantle metasomatism.
Scaillet, B; Prouteau, G
2001-01-01
Modern plate tectonic brings down oceanic crust along subduction zones where it either dehydrates or melts. Those hydrous fluids or melts migrate into the overlying mantle wedge trigerring its melting which produces arc magmas and thus additional continental crust. Nowadays, melting seems to be restricted to cases of young (< 50 Ma) subducted plates. Slab melts are silicic and strongly sodic (trondhjemitic). They are produced at low temperatures (< 1000 degrees C) and under water excess conditions. Their interaction with mantle peridotite produces hydrous metasomatic phases such as amphibole and phlogopite that can be more or less sodium rich. Upon interaction the slab melt becomes less silicic (dacitic to andesitic), and Mg, Ni and Cr richer. Virtually all exposed slab melts display geochemical evidence of ingestion of mantle material. Modern slab melts are thus unlike Archean Trondhjemite-Tonalite-Granodiorite rocks (TTG), which suggests that both types of magmas were generated via different petrogenetic pathways which may imply an Archean tectonic model of crust production different from that of the present-day, subduction-related, one. PMID:11838241
Transition on the Geometry of the Cocos Plate in Central-Southern Mexico.
NASA Astrophysics Data System (ADS)
Rodríguez-Domínguez, M. Á.; Perez-Campos, X.; Valencia-Cabrera, D.; Clayton, R. W.; Cordoba-Montiel, F.; Valdes-Gonzales, C. M.; Brudzinski, M. R.; Cabral-Cano, E.; Arciniega-Ceballos, A.
2014-12-01
The tectonic setting, produced by the interaction between the Cocos and North American plates, follows complex geometries along the Pacific coast. Previous studies in central Mexico showed that the slab dips nearly horizontally before steeply subducting into the continental mantle; in contrast, in southern Mexico, the slabs dips under the continental plate at a constant angle. Receiver functions from four seismic networks: GECO (Geometry of Cocos), SSN (Servicio Sismológico Nacional), OxNet (Oaxaca Network) and UV (Universidad Veracruzana) are used to study the crustal structure underneath the stations, and image the subducting Cocos plate in order to define the geometry and the transition angle in central-southern Mexico.
Centrosome Positioning in 1D Cell Migration
NASA Astrophysics Data System (ADS)
Adlerz, Katrina; Aranda-Espinoza, Helim
During cell migration, the positioning of the centrosome and nucleus define a cell's polarity. For a cell migrating on a two-dimensional substrate the centrosome is positioned in front of the nucleus. Under one-dimensional confinement, however, the centrosome is positioned behind the nucleus in 60% of cells. It is known that the centrosome is positioned by CDC42 and dynein for cells moving on a 2D substrate in a wound-healing assay. It is currently unknown, however, if this is also true for cells moving under 1D confinement, where the centrosome position is often reversed. Therefore, centrosome positioning was studied in cells migrating under 1D confinement, which mimics cells migrating through 3D matrices. 3 to 5 μm fibronectin lines were stamped onto a glass substrate and cells with fluorescently labeled nuclei and centrosomes migrated on the lines. Our results show that when a cell changes directions the centrosome position is maintained. That is, when the centrosome is between the nucleus and the cell's trailing edge and the cell changes direction, the centrosome will be translocated across the nucleus to the back of the cell again. A dynein inhibitor did have an influence on centrosome positioning in 1D migration and change of directions.
Coalescence phenomena in 1D silver nanostructures
NASA Astrophysics Data System (ADS)
Gutiérrez-Wing, C.; Pérez-Alvarez, M.; Mondragón-Galicia, G.; Arenas-Alatorre, J.; Gutiérrez-Wing, M. T.; Henk, M. C.; Negulescu, I. I.; Rusch, K. A.
2009-07-01
Different coalescence processes on 1D silver nanostructures synthesized by a PVP assisted reaction in ethylene glycol at 160 °C were studied experimentally and theoretically. Analysis by TEM and HRTEM shows different defects found on the body of these materials, suggesting that they were induced by previous coalescence processes in the synthesis stage. TEM observations showed that irradiation with the electron beam eliminates the boundaries formed near the edges of the structures, suggesting that this process can be carried out by the application of other means of energy (i.e. thermal). These results were also confirmed by theoretical calculations by Monte Carlo simulations using a Sutton-Chen potential. A theoretical study by molecular dynamics simulation of the different coalescence processes on 1D silver nanostructures is presented, showing a surface energy driven sequence followed to form the final coalesced structure. Calculations were made at 1000-1300 K, which is near the melting temperature of silver (1234 K). Based on these results, it is proposed that 1D nanostructures can grow through a secondary mechanism based on coalescence, without losing their dimensionality.
NASA Astrophysics Data System (ADS)
Margirier, A.; Robert, X.; Laurence, A.; Gautheron, C.; Bernet, M.; Simon-Labric, T.; Hall, S. R.
2015-12-01
Processes driving surface uplift in the Andes are still debated and the role of subduction processes as slab flattening on surface uplift and relief building in the Andes is not well understood. Some of the highest Andean summits, the Cordillera Blanca (6768 m) and the Cordillera Negra (5187 m), are located above a present flat subduction zone (3-15°S), in northern Peru. In this area, both the geometry and timing of the flattening of the slab are well constrained (Gutscher et al., 1999; Rosenbaum et al., 2005). This region is thus a perfect target to explore the effect of slab flattening on the Andean topography and uplift. We obtained new apatite (U-Th)/He and fission-track ages from three vertical profiles located in the Cordillera Blanca and the Cordillera Negra. Time-temperature paths obtained from inverse modeling of the thermochronological data indicates a Middle Miocene cooling for both Cordillera Negra profiles. We interpret it as regional exhumation in the Cordillera Occidental starting in Middle Miocene, synchronous with the onset of the subduction of the Nazca ridge (Rosenbaum et al., 2005). We propose that the Nazca ridge subduction at 15 Ma and onset of slab flattening in northern Peru drove regional positive dynamic topography and thus enhanced exhumation in the Cordillera Occidental. This study provides new evidence of the impact subduction processes and associated dynamic topography on paleogeography and surface uplift in the Andes.
Effect of reentrant cone geometry on energy transport in intense laser-plasma interactions
Lancaster, K. L.; Sherlock, M.; Heathcote, R.; Green, J. S.; Norreys, P. A.; Gregory, C. D.; Hakel, P.; Akli, K. U.; Hey, D. S.; Stephens, R. B.; Beg, F. N.; Chen, S. N.; Wei, M. S.; Yabuuchi, T.; Freeman, R. R.; Highbarger, K.; Van Woerkom, L.; Weber, R. L.; Habara, H.; Key, M. H.
2009-10-15
The energy transport in cone-guided low-Z targets has been studied for laser intensities on target of 2.5x10{sup 20} W cm{sup -2}. Extreme ultraviolet (XUV) imaging and transverse optical shadowgraphy of the rear surfaces of slab and cone-slab targets show that the cone geometry strongly influences the observed transport patterns. The XUV intensity showed an average spot size of 65{+-}10 {mu}m for slab targets. The cone slabs showed a reduced spot size of 44{+-}10 {mu}m. The shadowgraphy for the aforementioned shots demonstrate the same behavior. The transverse size of the expansion pattern was 357{+-}32 {mu}m for the slabs and reduced to 210{+-}30 {mu}m. A transport model was constructed which showed that the change in transport pattern is due to suppression of refluxing electrons in the material surrounding the cone.
INTERIOR VIEW, LOOKING NORTHEAST, WITH OVENS AND SLAB BEING PROCESSED ...
INTERIOR VIEW, LOOKING NORTHEAST, WITH OVENS AND SLAB BEING PROCESSED THROUGH PRESS/STECKLER (RIGHT). HEATED SLABS MAKE SEVERAL PASSES THROUGH THE STECKLER WITH THICKNESS OF THE SLAB DECREASED UNTIL REQUIRED GAGE IS REACHED. - Central Iron Foundry, Hot Strip Mill Building, 1700 Holt Road, Holt, Tuscaloosa County, AL
Necessity of the Ridge for the Flat Slab Subduction: Insights from the Peruvian Flat Slab
NASA Astrophysics Data System (ADS)
Knezevic Antonijevic, S.; Wagner, L. S.; Beck, S. L.; Long, M. D.; Zandt, G.; Tavera, H.
2014-12-01
Flattening of the subducting plate has been linked to the formation of various geological features, including basement-cored uplifts, the cessation of arc volcanism, ignimbrite flare-ups, and the formation of high plateaus and ore deposits [Humphreys et al., 2003; Gutscher et al., 2000; Rosenbaum et al., 2005]. However, the mechanism responsible for the slab flattening is still poorly understood. Here we focus on the Peruvian flat slab, where the Nazca plate starts to bend at ~80 km depth and travels horizontally for several hundred kilometers, at which point steep subduction resumes. Based on a 1500 km long volcanic gap and intermediate depth seismicity patterns, the Peruvian flat slab appears to have the greatest along-strike extent and, therefore, has been suggested as a modern analogue to the putative flat slab during the Laramide orogeny in the western United States (~80-55 Ma). Combining 3D shear wave velocity structure and Rayleigh wave phase anisotropy between ~10° and 18° S, we find that the subducting Nazca plate is not uniformly flat along the entire region, but fails to the north of the subducting Nazca Ridge. Our results show that, in combination with trench retreat, rapid overriding plate motion, and/or presence of a thick cratonic root, the subduction of buoyant overthickened oceanic crust, such as the Nazca Ridge, is necessary for the formation and sustainability of flat slabs. This finding has important implications for the formation of flat slabs both past and present.
Magnetotelluric imaging of the subducting slab in Cascadia with constraints from seismology
NASA Astrophysics Data System (ADS)
Yang, B.; Egbert, G. D.; Kelbert, A.; Humphreys, E.
2015-12-01
We present results from three-dimensional (3D) inversion of long-period magnetotelluric (MT) data from Cascadia, using seismological constraints on plate geometry and back-arc structure, to refine 3D images of electrical resistivity across this subduction zone. For this study we employed the impedances and vertical transfer functions from 144 sites from the EarthScope Transportable Array, along with data from previous higher density MT profiles from Cascadia (EMSLAB, CAFE-MT etc.). Morphological parameters for the subducting Juan de Fuca and Gorda plates (e.g. upper boundary and thickness) were extracted from McCrory et al (2012) and Schmandt and Humphreys (2010) seismological models and used to define a resistive subducting slab structure in 3D. This was then either used as a prior model, or fixed (both resistivity and geometry) during the MT inversion. By imposing constraints on the geometry of the slab (which is otherwise imaged as an amorphous broad resistive zone) we improve recovery and resolution of subduction related conductivity features. The constrained inversions also allowed us to test sensitivity of the MT data to variants on slab geometry, such as the proposed slab "tear" near the Oregon-Washington border suggested by some seismic tomography models, and to explore consistency of the MT data with seismic models, which suggest segmentation of back-arc upwelling. Three zones of substantially reduced resistivity were found, all exhibiting significant along-strike variability. In the forearc, an N-S stripe of high conductivity (10 ohm-m or less) was found just above the plate interface, near the tip of the mantle wedge. This conductive feature is spatially coincident with mapped locations of episodic tremor and slip, and likely represents aqueous fluids associated with slab dehydration. To the east, a second, clearly separated, N-S elongate zone of similarly high conductivity occurs in the mid-lower crust and upper mantle beneath the modern arc, again
Overriding plate thickness control on subducting slab curvature
NASA Astrophysics Data System (ADS)
Holt, A.; Buffett, B. A.; Becker, T. W.
2014-12-01
The curvature of subducting lithosphere controls deformation due to bending at the trench, which results in a force that dissipates gravitational potential energy and may affect seismic coupling. We use 2-D, thermo-mechanical subduction models to explore the dependence of the radius of curvature on the thickness of the subducting and overriding plates for models with both viscous and effectively plastic lithospheric rheologies. Such a plastic rheology has been shown to reproduce the bending stresses/moment computed using a kinematic strain rate description and a laboratory derived composite rheology. Laboratory and numerical models show that the bending geometry of subducting slabs with a viscous rheology is strongly dependent on slab thickness; thicker plates have a larger radius of curvature. However, the curvature of subducting plates on Earth, illuminated by the distribution of earthquake hypocenters, shows little to no dependence on the plate thickness or age. Such an observation is instead compatible with plates that have a plastic rheology. Indeed, our numerical models show that the radius of curvature of viscous plates has a stronger dependence on subducting plate thickness than in equivalent plastic models. In viscous plates, the bending moment produces a torque, which balances the torque exerted by buoyancy. However, for the plastic plate case the bending moment saturates at a maximum value and so cannot balance the gravitational torque. The saturation of bending moment means that, (a) the radius of curvature of the bending region is not constrained by this torque balance, and, (b) other forces are required to balance the gravitational torque. We explore the role that the overriding plate could play in controlling the subducting plate curvature in plastic plate models where the bending stresses have saturated. For such plates, we find that increasing the thickness of the overriding plate causes the radius of curvature to increase. The same correlation is
Geometry Transition in the Cocos Plate, from Flat-Steep to Constant Dip: Smooth or Abrupt?
NASA Astrophysics Data System (ADS)
Perez-Campos, X.; Clayton, R. W.; Brudzinski, M. R.; Valdés-González, C. M.; Cabral-Cano, E.; Arciniega-Ceballos, A.; Córdoba-Montiel, F.
2013-05-01
Subduction of the Cocos Plate beneath North America has a variable and complex behavior along the Middle-American Trench. Initially, its geometry was delineated from regional seismicity. In the last 10 years, seismic experiments have illuminated some details in the geometry. They have reported, from NW to SE an abrupt dip transition, from 50 to 26°, as the result of a tear that splits Cocos North from Cocos South; then there is a smooth transition to a horizontal geometry under central Mexico. Further southeast, under the Isthmus of Tehuantepec, the Cocos plate shows a constant ~26° subduction dip. This last transition has been assumed to be smooth from the sparse seismicity in the region. A first glimpse of the slab geometry under Oaxaca, shows the slab continues to be flat at least until 97.5°W longitude, where the slab suddenly changes to a ~55° dip to the northeast. This occurs at a distance of ~75 km from the Pico de Orizaba volcano, which is a similar distance as the active Popocatepetl volcano from the place where the slab dives into the mantle along the Meso-American Subduction Experiment line, in central Mexico. East of this region, receiver function images show an abrupt change in the geometry and length of the slab.
A 1-D dusty plasma photonic crystal
Mitu, M. L.; Ticoş, C. M.; Toader, D.; Banu, N.; Scurtu, A.
2013-09-21
It is demonstrated numerically that a 1-D plasma crystal made of micron size cylindrical dust particles can, in principle, work as a photonic crystal for terahertz waves. The dust rods are parallel to each other and arranged in a linear string forming a periodic structure of dielectric-plasma regions. The dispersion equation is found by solving the waves equation with the boundary conditions at the dust-plasma interface and taking into account the dielectric permittivity of the dust material and plasma. The wavelength of the electromagnetic waves is in the range of a few hundred microns, close to the interparticle separation distance. The band gaps of the 1-D plasma crystal are numerically found for different types of dust materials, separation distances between the dust rods and rod diameters. The distance between levitated dust rods forming a string in rf plasma is shown experimentally to vary over a relatively wide range, from 650 μm to about 1350 μm, depending on the rf power fed into the discharge.
Slab photonic crystals with dimer colloid bases
Riley, Erin K.; Liddell Watson, Chekesha M.
2014-06-14
The photonic band gap properties for centered rectangular monolayers of asymmetric dimers are reported. Colloids in suspension have been organized into the phase under confinement. The theoretical model is inspired by the range of asymmetric dimers synthesized via seeded emulsion polymerization and explores, in particular, the band structures as a function of degree of lobe symmetry and degree of lobe fusion. These parameters are varied incrementally from spheres to lobe-tangent dimers over morphologies yielding physically realizable particles. The work addresses the relative scarcity of theoretical studies on photonic crystal slabs with vertical variation that is consistent with colloidal self-assembly. Odd, even and polarization independent gaps in the guided modes are determined for direct slab structures. A wide range of lobe symmetry and degree of lobe fusion combinations having Brillouin zones with moderate to high isotropy support gaps between odd mode band indices 3-4 and even mode band indices 1-2 and 2-3.
Subduction zone earthquakes and stress in slabs
NASA Technical Reports Server (NTRS)
Vassiliou, M. S.; Hager, B. H.
1988-01-01
Simple viscous fluid models of subducting slabs are used to explain observations of the distribution of earthquakes as a function of depth and the orientation of stress axes of deep (greater than 300 km) and intermediate (70-300 km) earthquakes. Results suggest the following features in the distribution of earthquakes with depth: (1) an exponential decrease from shallow depths down to 250 to 300 km, (2) a minimum near 250 to 300 km, and (3) a deep peak below 300 km. Many shallow subducting slabs show only the first characteristic, while deeper extending regions tend to show all three features, with the deep peak varying in position and intensity. These data, combined with the results on the stress orientations of various-depth earthquakes, are consistent with the existence of a barrier of some sort at 670-km depth and a uniform viscosity mantle above this barrier.
Effects of Thermodynamic Properties on Slab Evolution
NASA Astrophysics Data System (ADS)
Wada, I.; King, S. D.; Caddick, M. J.; Ross, N.
2012-12-01
We perform a series of numerical experiments to investigate the effects of thermodynamic properties on the geometrical evolution of subducting slabs. We calculate density (ρ), thermal expansivity (α), and heat capacity (cp) of mantle mineral assemblages of a harzburgite composition over a range of pressure and temperature conditions applicable to the Earth's mantle, using the thermodynamic database of Stixrude and Lithgow-Bertelloni [2011] and the thermodynamic calculation code Perple_X [Connolly, 2009]. Following Nakagawa et al. [2009], we assume that thermal diffusivity (κ) follows a power-law relationship with density (κ=κ0(ρ/ρ0)3, where κ0 and ρ0 are reference diffusivity and density, respectively). The calculations show that ρ, α, and κ change significantly along mantle geotherms; the ranges of their values are 3300-5100 km/m3, 1.5-3.5 10-5/K, and 3-17 W/m K, respectively. The change in cp is small (< 5%). We incorporate the pressure and temperature (PT) dependence of these thermodynamic properties into a 2-D finite element code with compressible convection formulations under the truncated anelastic liquid approximation [Lee and King, 2009] and develop a dynamic subduction model with kinematic boundary conditions. In the model, we use a composite mantle rheology that accounts for both diffusion and dislocation creep with flow law parameterization of wet olivine [Hirth and Kohlstedt, 2003]. Following Billen and Hirth [2007] and Lee and King [2011], we adjust the flow law parameter values for the lower mantle to test the effects of viscosity contrast between the upper and lower mantle on slab evolution. We use a reference model with a constant ρ, κ α, and cp, which is equivalent to using the incompressible extended Bousisnesq approximation. Preliminary results show that incorporating PT-dependent ρ enhances the vigor of the buoyancy driven flow compared to the reference model. Further, lithostatic pressure at a given depth is higher than in the
The effect of metastable pyroxene on the slab dynamics
NASA Astrophysics Data System (ADS)
Agrusta, Roberto; Hunen, Jeroen; Goes, Saskia
2014-12-01
Seismic studies show that some subducting slabs penetrate straight into the lower mantle, whereas others seem to flatten near the base of the mantle transition zone. Slab stagnation is often attributed to an increase in viscosity and phase transformations in the olivine system. However, recent mineral physics studies showed that due to extremely low transformational diffusion rates, low-density metastable pyroxene may persist into the transition zone in cool slabs. Here we use a dynamically fully self-consistent subduction model to investigate the influence of metastable pyroxene on the dynamics of subducting oceanic lithosphere. Our results show that metastable pyroxene affects slab buoyancy at least as much as olivine metastability. However, unlike metastable olivine, which can inhibit slab penetration in the lower mantle only for cold, old, and fast slabs, metastable pyroxene is likely to also affect sinking of relatively young and slow slabs.
A numerical study of the interaction between the mantle wedge, subducting slab, and overriding plate
NASA Astrophysics Data System (ADS)
Eberle, Michael A.; Grasset, Olivier; Sotin, Christophe
2002-12-01
We have formulated a numerical model with strongly temperature-dependent viscosity to calculate thermal structure and flow-field in subduction zones. One important particularity of the model is that the overriding plate is not fixed over its whole thickness in order to allow material exchange between the wedge and the upper lithosphere. Numerical problems due to very high-viscosity contrasts are avoided by coupling a finite difference method and a finite element method for solving the energy conservation equation and the Stockes equation, respectively. In this model, a temperature decrease from 1400 to 1300 °C increases the viscosity by an order of magnitude. We study the temperature structure and the velocity field of the subducting slab and mantle wedge. Surface heat flow, velocity anomalies, and geometry of the partial melting zone are also calculated. To study the effect that boundary conditions play on the interaction between the mantle wedge, overriding plate and subducting plate, we examine models with both fixed and free-slip conditions applied to the overriding plate. When the overriding plate is allowed to move laterally (free-slip), the subducting slab is thick, and both the temperature field and the convective motions in the mantle wedge are similar to those observed when using constant viscosity numerical models or analytical corner flow models. If the surface of the overriding plate is fixed, the subducting slab is thin and the mantle wedge impinges upon the overriding plate forming a high-temperature nose between the overriding plate and subducting lithosphere. Furthermore, viscous decoupling occurs implicitly at shallow depth between the slab and the wedge because hot material from the wedge is entrained close to the trench. In that case, the subducting slab tectonically erodes the lower lithosphere of the overriding plate leading to high-temperatures, low seismic velocities, high attenuation and high heat flow beneath volcanic arc, in agreement
The 3-D dynamics of slab break-off and implications for continental collision zones
NASA Astrophysics Data System (ADS)
van Hunen, Jeroen; Allen, Mark
2010-05-01
Some of the world best studied mountain ranges are a result of continental collision, such as the Himalayas, Zagros mountains, and the Alps. Continental collision forms the last stage of the closure of an oceanic basin, and leads to the slow-down or complete cessation of the subduction process. Previously subducted slab material will experience a period of thermal warming (Gerya et al., 2004) and/or a larger tensile stress, and will eventually weaken, yield and sink into the mantle. This process has potentially important implications for the thermal and stress regime of the overlying convergence zone, and has been held responsible for various phenomena such as late-stage magmatism (Davies and von Blanckenburg, 1995) and surface uplift or depression (van der Meulen et al., 1998, Buiter et al., 2002). Even though the collision process itself is relatively short-lived compared to the preceding oceanic subduction, its remnants are often preserved, and probably provide a valuable window into the plate tectonic process during the Proterozoic and perhaps the Archaean (e.g. Calvert et al., 1995). The three-dimensional nature of this break-off process has previously been discussed with conceptual models. E.g. slab break-off has been suggested to propagate laterally through an advancing tear (Wortel and Spakman, 2000). In this study we present 3D numerical results of the evolution of slab break-off. We focus on the development and evolution of a laterally migrating slab tear, and present results on the sensitivity of this process to the geometry of the closing oceanic basin, the tensile stresses in and the rheological properties of the slab, and the thermal state of the surrounding mantle. By comparing our numerical results to previously published analogue results (Regard et al., 2004) and various tomographic, structural, and magmatic observations of well-studied subduction collision systems, we are able to extract valuable insights in to the dynamics and strength of
The Effect of Subducting Slabs in Global Shear Wave Tomography
NASA Astrophysics Data System (ADS)
Lu, Chang; Grand, Stephen P.
2016-03-01
Subducting slabs create strong short wavelength seismic anomalies in the upper mantle where much of Earth's seismicity is located. As such, they have the potential to bias longer wavelength seismic tomography models. To evaluate the effect of subducting slabs in global tomography, we performed a series of inversions using a global synthetic shear wave travel time dataset for a theoretical slab model based on predicted thermal anomalies within slabs. The spectral element method was applied to predict the travel time anomalies produced by the 3D slab model for paths corresponding to our current data used in actual tomography models. Inversion tests have been conducted first using the raw travel time anomalies to check how well the slabs can be imaged in global tomography without the effect of earthquake mislocation. Our results indicate that most of the slabs can be identified in the inversion result but with smoothed and reduced amplitude. The recovery of the total mass anomaly in slab regions is about 88%. We then performed another inversion test to investigate the effect of mislocation caused by subducting slabs. We found that source mislocation largely removes slab signal and significantly degrades the imaging of subducting slabs - potentially reducing the recovery of mass anomalies in slab regions to only 41%. We tested two source relocation procedures - an iterative relocation inversion and joint relocation inversion. Both methods partially recover the true source locations and improve the inversion results, but the joint inversion method worked significantly better than the iterative method. In all of our inversion tests, the amplitude of artifact structures in the lower mantle caused by the incorrect imaging of slabs (up to ˜0.5% S velocity anomalies) are comparable to some large scale lower mantle heterogeneities seen in global tomography studies. Based on our inversion tests, we suggest including a-priori subducting slabs in the starting models in global
The effect of subducting slabs in global shear wave tomography
NASA Astrophysics Data System (ADS)
Lu, Chang; Grand, Stephen P.
2016-05-01
Subducting slabs create strong short wavelength seismic anomalies in the upper mantle where much of Earth's seismicity is located. As such, they have the potential to bias longer wavelength seismic tomography models. To evaluate the effect of subducting slabs in global tomography, we performed a series of inversions using a global synthetic shear wave traveltime data set for a theoretical slab model based on predicted thermal anomalies within slabs. The spectral element method was applied to predict the traveltime anomalies produced by the 3-D slab model for paths corresponding to our current data used in actual tomography models. Inversion tests have been conducted first using the raw traveltime anomalies to check how well the slabs can be imaged in global tomography without the effect of earthquake mislocation. Our results indicate that most of the slabs can be identified in the inversion result but with smoothed and reduced amplitude. The recovery of the total mass anomaly in slab regions is about 88 per cent. We then performed another inversion test to investigate the effect of mislocation caused by subducting slabs. We found that source mislocation largely removes slab signal and significantly degrades the imaging of subducting slabs-potentially reducing the recovery of mass anomalies in slab regions to only 41 per cent. We tested two source relocation procedures-an iterative relocation inversion and joint relocation inversion. Both methods partially recover the true source locations and improve the inversion results, but the joint inversion method worked significantly better than the iterative method. In all of our inversion tests, the amplitudes of artefact structures in the lower mantle caused by the incorrect imaging of slabs (up to ˜0.5 per cent S velocity anomalies) are comparable to some large-scale lower-mantle heterogeneities seen in global tomography studies. Based on our inversion tests, we suggest including a-priori subducting slabs in the
Fast Waves in Smooth Coronal Slab
NASA Astrophysics Data System (ADS)
Lopin, I.; Nagorny, I.
2015-03-01
This work investigates the effect of transverse density structuring in coronal slab-like waveguides on the properties of fast waves. We generalized previous results obtained for the exponential and Epstein profiles to the case of an arbitrary transverse density distribution. The criteria are given to determine the possible (trapped or leaky) wave regime, depending on the type of density profile function. In particular, there are plasma slabs with transverse density structuring that support pure trapped fast waves for all wavelengths. Their phase speed is nearly equal to the external Alfvén speed for the typical parameters of coronal loops. Our findings are obtained on the basis of Kneser’s oscillation theorem. To confirm the results, we analytically solved the wave equation evaluated at the cutoff point and the original wave equation for particular cases of transverse density distribution. We also used the WKB method and obtained approximate solutions of the wave equation at the cutoff point for an arbitrary transverse density profile. The analytic results were supplemented by numerical solutions of the obtained dispersion relations. The observed high-quality quasi-periodic pulsations of flaring loops are interpreted in terms of the trapped fundamental fast-sausage mode in a slab-like coronal waveguide.
Implementing slab solar water heating system
NASA Astrophysics Data System (ADS)
Raveendran, S. K.; Shen, C. Q.
2015-08-01
Water heating contributes a significant part of energy consumption in typical household. One of the most employed technologies today that helps in reducing the energy consumption of water heating would be conventional solar water heating system. However, this system is expensive and less affordable by most family. The main objective of this project is to design and implement an alternative type of solar water heating system that utilize only passive solar energy which is known as slab solar water heating system. Slab solar water heating system is a system that heat up cold water using the solar radiance from the sun. The unique part of this system is that it does not require any form of electricity in order to operate. Solar radiance is converted into heat energy through convection method and cold water will be heated up by using conduction method [1]. The design of this system is governed by the criteria of low implementation cost and energy saving. Selection of material in the construction of a slab solar water heating system is important as it will directly affect the efficiency and performance of the system. A prototype has been built to realize the idea and it had been proven that this system was able to provide sufficient hot water supply for typical household usage at any given time.
1D-VAR Retrieval Using Superchannels
NASA Technical Reports Server (NTRS)
Liu, Xu; Zhou, Daniel; Larar, Allen; Smith, William L.; Schluessel, Peter; Mango, Stephen; SaintGermain, Karen
2008-01-01
Since modern ultra-spectral remote sensors have thousands of channels, it is difficult to include all of them in a 1D-var retrieval system. We will describe a physical inversion algorithm, which includes all available channels for the atmospheric temperature, moisture, cloud, and surface parameter retrievals. Both the forward model and the inversion algorithm compress the channel radiances into super channels. These super channels are obtained by projecting the radiance spectra onto a set of pre-calculated eigenvectors. The forward model provides both super channel properties and jacobian in EOF space directly. For ultra-spectral sensors such as Infrared Atmospheric Sounding Interferometer (IASI) and the NPOESS Airborne Sounder Testbed Interferometer (NAST), a compression ratio of more than 80 can be achieved, leading to a significant reduction in computations involved in an inversion process. Results will be shown applying the algorithm to real IASI and NAST data.
A two- and three-dimensional numerical comparison study of slab detachment
NASA Astrophysics Data System (ADS)
Thieulot, Cedric; Buiter, Susanne; Brune, Sascha; Davies, Rhodri; Duretz, Thibault; Gerbault, Muriel; Glerum, Anne; Quinteros, Javier; Schmalholz, Stefan; Spakman, Wim
2015-04-01
Subduction is likely to be the most studied phenomenon in Numerical Geodynamics. Over the past 20 years, hundreds of publications have focused on its various aspects (influence of the rheology and thermal state of the plates, slab-mantle coupling, roll-back, mantle wedge evolution, buoyancy changes due to phase change, ...) and results were obtained with a variety of codes. Slab detachment has recently received some attention but remains a field worth exploring due to its profound influence on dynamic topography, mantle flow and subsequent stress state of the plates, and is believed to have occured in the Zagros, Carpathians and beneath eastern Anatolia, to name only a few regions. Following the work of Schmalholz (2011), we propose a two- and three-dimensional numerical benchmark of slab detachment. The geometry is simple: a power-law T-shaped plate including an already subducted slab overlies the mantle whose viscosity is either linear or power-law. Boundary conditions are free-slip on the top and the bottom of the domain, and no-slip on the sides. When the system evolves in time, the slab stretches out vertically and shows buoyancy-driven necking, until it finally detaches. The benchmark is subdivided into several sub-experiments with gradually increase in complexity (free surface, coupling of the rheology with temperature, ...). An array of objective measurements is recorded throughout the simulation such as the width of the necked slab over time and the exact time of detachment. The experiments will be run in two-dimensions and repeated in three-dimensional, the latter case being designed so as to allow both poloidal and toroidal flow. We show results obtained with a multitude of Finite Element and Finite Difference codes, using either compositional fields, level sets or tracers to track the compositions. A good agreement is found for most of the measurements in the two-dimensional case, and preliminary three-dimensional measurements will be shown. Schmalholz
Slab flattening driving regional uplift in the Cordilleras Blanca and Negra, Western Andes
NASA Astrophysics Data System (ADS)
Margirier, Audrey; Audin, Laurence; Robert, Xavier; Bernet, Matthias; Gautheron, Cécile
2015-04-01
The Andean range topographic evolution is known to have had a strong impact on regional climate by building an orographic barrier that preserved its western flank from the south Atlantic moisture. Even if largely invoked, the impact of subduction processes on the uplift and relief building is not yet well understood in the Andes. The northern Peru is characterized by a present day flat subduction zone (3-15°S), where both the geometry and temporal evolution of the flat-slab are well constrained. The subduction of two buoyant anomalies, the Nazca ridge and the lost Inca plateau controlled the slab flattening. The highest Peruvian peaks in the Cordillera Blanca (6768 m), and the Cordillera Negra (5187 m) are located just above the flat-slab segment. Both ranges trend parallel to the subduction zone and are separated by the NW-SE Rio Santa valley. The Cordillera Blanca batholith emplaced at 8-5 Ma and renders of an abnormal magmatic activity over a planar subduction. This area is a perfect target to explore the impact of slab flattening on the topography and uplift in the Occidental Cordillera of the Andes. We present new AHe and AFT data from three vertical profiles located in both the Cordilleras Blanca and Negra. We compare time-temperature paths obtained from inverse modeling of the thermochronological data with the timing of the slab flattening, the arrival of the Nazca ridge and magmatism. Our thermochronological data evidences a regional exhumation in the Occidental Cordillera from ~10 Ma. We propose that the Nazca ridge subduction below the Occidental Cordillera (11 Ma) and slab flattening (8 Ma) drive the Occidental Cordillera uplift and thus exhumation. We evidence the important contribution of the magmatism in the Cordillera Blanca exhumation and high relief building in the Occidental Cordillera. Our new thermochronological data highlight the control of both the subduction processes and magmatism on the paleogeography and uplift in the Andes. Finally, the
NASA Astrophysics Data System (ADS)
Bina, C. R.; Čížková, H.
2015-12-01
Recent work indicates that subducting slabs may exhibit buckling instabilities and consequent folding behavior in the mantle transition zone for various dynamical parameters, accompanied by temporal variations in dip angle, plate velocity, and trench retreat. Governing parameters include both viscous (rheological) and buoyancy (thermo-petrological) forces. 2D numerical experiments show that many parameter sets lead to slab deflection at the base of the transition zone, typically accompanied by quasi-periodic oscillations in largely anticorrelated plate and rollback velocities, resulting in undulating stagnant slabs as buckle folds accumulate subhorizontally atop the lower mantle. Slab petrology, of mantle phase transitions and hydrated crust, is a dominant factor in this process (Čížková and Bina, 2013). For terrestrial parameter sets, trench retreat is found to be nearly ubiquitous and trench advance quite rare, largely due to rheological and ridge-push effects. Recently updated analyses of global plate motions indicate that significant trench advance is also rare on Earth, being largely restricted to the Izu-Bonin arc (Matthews et al., 2013). Thus, we explore conditions necessary for terrestrial trench advance through dynamical models involving the unusual geometry of the Philippine Sea region. Our 2D modeling of such geometries, in which distal subduction of the overriding plate overprints an opposed slab-pull force on the usual ridge-push at the trench, yields persistent trench advance interrupted by episodes of back-arc extension, demonstrating that trench advance can occur for terrestrial rheologies in such special geometries (Čížková and Bina, 2015).
Viscous Dissipation and Criticality of Subducting Slabs
NASA Astrophysics Data System (ADS)
Riedel, Mike; Karato, Shun; Yuen, Dave
2016-04-01
Rheology of subducting lithosphere appears to be complicated. In the shallow part, deformation is largely accomodated by brittle failure, whereas at greater depth, at higher confining pressures, ductile creep is expected to control slab strength. The amount of viscous dissipation ΔQ during subduction at greater depth, as constrained by experimental rock mechanics, can be estimated on the basis of a simple bending moment equation [1,2] 2ɛ˙0(z) ∫ +h/2 2 M (z) = h ṡ ‑h/2 4μ(y,z)y dy , (1) for a complex multi-phase rheology in the mantle transition zone, including the effects of a metastable phase transition as well as the pressure, temperature, grain-size and stress dependency of the relevant creep mechanisms; μ is here the effective viscosity and ɛ˙0(z) is a (reference) strain rate. Numerical analysis shows that the maximum bending moment, Mcrit, that can be sustained by a slab is of the order of 1019 Nm per m according to Mcrit˜=σp ∗h2/4, where σp is the Peierl's stress limit of slab materials and h is the slab thickness. Near Mcrit, the amount of viscous dissipation grows strongly as a consequence of a lattice instability of mantle minerals (dislocation glide in olivine), suggesting that thermo-mechanical instabilities become prone to occur at places where a critical shear-heating rate is exceeded, see figure. This implies that the lithosphere behaves in such cases like a perfectly plastic solid [3]. Recently available detailed data related to deep seismicity [4,5] seems to provide support to our conclusion. It shows, e.g., that thermal shear instabilities, and not transformational faulting, is likely the dominating mechanism for deep-focus earthquakes at the bottom of the transition zone, in accordance with this suggested "deep criticality" model. These new findings are therefore briefly outlined and possible implications are discussed. References [1] Riedel, M. R., Karato, S., Yuen, D. A. Criticality of Subducting Slabs. University of Minnesota
Prediction of car cabin environment by means of 1D and 3D cabin model
NASA Astrophysics Data System (ADS)
Fišer, J.; Pokorný, J.; Jícha, M.
2012-04-01
Thermal comfort and also reduction of energy requirements of air-conditioning system in vehicle cabins are currently very intensively investigated and up-to-date issues. The article deals with two approaches of modelling of car cabin environment; the first model was created in simulation language Modelica (typical 1D approach without cabin geometry) and the second one was created in specialized software Theseus-FE (3D approach with cabin geometry). Performance and capabilities of this tools are demonstrated on the example of the car cabin and the results from simulations are compared with the results from the real car cabin climate chamber measurements.
Asymmetric photonic resonances in GaN slab waveguide for mid infrared selective filters.
Vangala, S R; Avrutsky, I; Keiffer, P; Nader, N; Walker, D; Cleary, J W; Hendrickson, J R
2014-10-01
We demonstrate a spectrally selective reflector that exploits asymmetric photonic resonances of a 1D photonic crystal. The proposed spectrally selective reflector has a very simple structure - essentially just a single high-index slab of GaN, properly perforated, and supported by a transparent sapphire substrate. With the proper 1D array design, nearly 100% reflection is achieved with a narrow spectral width between 10 cm⁻¹ - 18 cm⁻¹, while the background reflection remains low across the entire mid-IR range. The reflection peak can be tuned over a large wavelength span based on physical parameters. Resonant transmission dips in the experimentally measured spectra corroborate the device theory and simulation, exhibiting the narrowband low-background mid-IR reflection as predicted. PMID:25322049
Dynamic buckling of subducting slabs reconciles geological and geophysical observations
NASA Astrophysics Data System (ADS)
Lee, Changyeol; King, Scott D.
2011-12-01
Ever since the early days of the development of plate tectonic theory, subduction zones have been engrained in geological thinking as the place where steady, linear slabs descend into the mantle at a constant, uniform dip angle beneath volcanic arcs. However, growing evidence from geological and geophysical observations as well as analog and numerical modeling indicates that subducting slabs buckle in a time-dependent manner, in contrast to the steady-state, linear cartoons that dominate the literature. To evaluate the implication of time-dependent slab buckling of geological events, we conduct a series of 2-D numerical dynamic/kinematic subduction experiments by varying the viscosity increase across the 660 km discontinuity and the strength of the subducting slab. Our results show that slab buckling is a universal figure in all the experiments when rate of the trench migration ( Vtrench) is relatively slow ( Vtrench| < 2 cm/a) and viscosity increases across the 660 km discontinuity are greater than a factor of 30. Slab buckling is expressed as alternate shallowing and steepening dip of the subducting slab (from ~ 40 to ~ 100°) which is correlated with increasing and decreasing convergent rate of the incoming oceanic plate toward the trench. Further, the slab buckling in our experiments is consistent with the previously developed scaling laws for slab buckling; using reasonable parameters from subducted slabs the buckling amplitude and period are ~ 400 km and ~ 25 Myr, respectively. The slab buckling behavior in our experiments explains a variety of puzzling geological and geophysical observations. First, the period of slab buckling is consistent with short periodic variations (~ 25 Myr) in the motions of the oceanic plates that are anchored by subduction zones. Second, the scattered distributions of slab dips (from ~ 20 to ~ 90°) in the upper mantle are snapshots of time-dependent slab dip. Third, the current compressional and extensional stress environments in
Buoyancy, bending, and seismic visibility in deep slab stagnation
NASA Astrophysics Data System (ADS)
Bina, Craig R.; Kawakatsu, Hitoshi; Suetsugu, D.; Bina, C.; Inoue, T.; Wiens, D.; Jellinek, M.
2010-11-01
The petrological consequences of deep subhorizontal deflection ("stagnation") of subducting slabs should affect both apparent seismic velocity structures and slab morphology. We construct kinematic thermal models of stagnant slabs and perform thermodynamic modeling of the consequent perturbation of high-pressure phase transitions in mantle minerals, focusing upon Japan as our study area. We calculate associated thermo-petrological buoyancy forces and bending moments which (along with other factors such as viscosity variations and rollback dynamics) may contribute to slab deformation. We consider effects of variations in depth of stagnation, post-stagnation dip angle, phase transition sharpness, transition triplication due to multiple intersection of geotherms with phase boundaries, and potential persistence of metastable phases due to kinetic hindrance. We also estimate seismic velocity anomalies, as might be imaged by seismic tomography, and corresponding seismic velocity gradients, as might be imaged by receiver-function analysis. We find that buoyant bending moment gradients of petrological origin at the base of the transition zone may contribute to slab stagnation. Such buoyancy forces vary with the depth at which stagnation occurs, so that slabs may seek an equilibrium slab stagnation depth. Metastable phase bending moment gradients further enhance slab stagnation, but they thermally decay after ∱/4600•700 km of horizontal travel, potentially allowing stagnant slabs to descend into the lower mantle. Stagnant slabs superimpose zones of negative velocity gradient onto a depressed 660-km seismic discontinuity, affecting the seismological visibility of such features. Seismologically resolvable details should depend upon both stagnation depth and the nature of the imaging technique (travel-time tomography vs. boundary-interaction phases). While seismic tomography appears to yield images of stagnant slabs, discontinuity topography beneath Japan resolved by
NASA Astrophysics Data System (ADS)
Nishikawa, T.; Ide, S.
2014-12-01
There are clear variations in maximum earthquake magnitude among Earth's subduction zones. These variations have been studied extensively and attributed to differences in tectonic properties in subduction zones, such as relative plate velocity and subducting plate age [Ruff and Kanamori, 1980]. In addition to maximum earthquake magnitude, the seismicity of medium to large earthquakes also differs among subduction zones, such as the b-value (i.e., the slope of the earthquake size distribution) and the frequency of seismic events. However, the casual relationship between the seismicity of medium to large earthquakes and subduction zone tectonics has been unclear. Here we divide Earth's subduction zones into over 100 study regions following Ide [2013] and estimate b-values and the background seismicity rate—the frequency of seismic events excluding aftershocks—for subduction zones worldwide using the maximum likelihood method [Utsu, 1965; Aki, 1965] and the epidemic type aftershock sequence (ETAS) model [Ogata, 1988]. We demonstrate that the b-value varies as a function of subducting plate age and trench depth, and that the background seismicity rate is related to the degree of slab bending at the trench. Large earthquakes tend to occur relatively frequently (lower b-values) in shallower subduction zones with younger slabs, and more earthquakes occur in subduction zones with deeper trench and steeper dip angle. These results suggest that slab buoyancy, which depends on subducting plate age, controls the earthquake size distribution, and that intra-slab faults due to slab bending, which increase with the steepness of the slab dip angle, have influence on the frequency of seismic events, because they produce heterogeneity in plate coupling and efficiently inject fluid to elevate pore fluid pressure on the plate interface. This study reveals tectonic factors that control earthquake size distribution and seismicity rate, and these relationships between seismicity and
NASA Astrophysics Data System (ADS)
Scire, Alissa; Zandt, George; Beck, Susan; Long, Maureen; Wagner, Lara; Minaya, Estela; Tavera, Hernando
2016-01-01
Two arrays of broad-band seismic stations were deployed in the north central Andes between 8° and 21°S, the CAUGHT array over the normally subducting slab in northwestern Bolivia and southern Peru, and the PULSE array over the southern part of the Peruvian flat slab where the Nazca Ridge is subducting under South America. We apply finite frequency teleseismic P- and S-wave tomography to data from these arrays to investigate the subducting Nazca plate and the surrounding mantle in this region where the subduction angle changes from flat north of 14°S to normally dipping in the south. We present new constraints on the location and geometry of the Nazca slab under southern Peru and northwestern Bolivia from 95 to 660 km depth. Our tomographic images show that the Peruvian flat slab extends further inland than previously proposed along the projection of the Nazca Ridge. Once the slab re-steepens inboard of the flat slab region, the Nazca slab dips very steeply (˜70°) from about 150 km depth to 410 km depth. Below this the slab thickens and deforms in the mantle transition zone. We tentatively propose a ridge-parallel slab tear along the north edge of the Nazca Ridge between 130 and 350 km depth based on the offset between the slab anomaly north of the ridge and the location of the re-steepened Nazca slab inboard of the flat slab region, although additional work is needed to confirm the existence of this feature. The subslab mantle directly below the inboard projection of the Nazca Ridge is characterized by a prominent low-velocity anomaly. South of the Peruvian flat slab, fast anomalies are imaged in an area confined to the Eastern Cordillera and bounded to the east by well-resolved low-velocity anomalies. These low-velocity anomalies at depths greater than 100 km suggest that thick mantle lithosphere associated with underthrusting of cratonic crust from the east is not present. In northwestern Bolivia a vertically elongated fast anomaly under the Subandean Zone
NASA Astrophysics Data System (ADS)
Cotte, N.; Perfettini, H.; Rousset, B.; Kostoglodov, V.; Lasserre, C.; Walpersdorf, A.; Cabral, E.
2014-12-01
Since the first cGPS have been installed in 1997 in the Guerrero-Oaxaca region of the Central American subduction zone in Mexico, several slow slip events (SSE) have been observed there. Further studies showed that the SSE characteristics are disparate in Guerrero and Oaxaca in terms of duration, amplitude and depth. In particular, recent investigations show that the slip during SSE intrudes the seismogenic zone in Guerrero, but barely in Oaxaca. The SSE lateral and depth extent has important implications on the seismic cycle. It affects the redistribution of stress (increasing or decreasing it) within the seismogenic zone, thus the timing and magnitude of future earthquakes. It is thus important to constrain a degree of coupling on the subduction interface and how the coupling evolves in space and time. This can be achieved by refining coupling models based on the observations and exploring a sensitivity of the models to variations of their main parameters.To better understand the effect of the SSE on the seismic cycle along the Mexican subduction zone, we analyze cGPS data and invert the time series, using a method based on a principal component analysis: PCAIM (Principal Component Analysis-based Inversion Method). For both Guerrero and Oaxaca areas, we first discuss the geometry of the plates' interface. Then we investigate in particular the effect of possible lateral and down-dip variations of the subduction interface on coupling models. Subsequent modeling using the PCAIM of the sequential SSEs that occurred in the past decade permitted to compare our results with the inferences of previously published works. We present the first analysis of the latest large SSE that started in January 2014 in Guerrero. We show that the PCAIM allows to describe a complex slip history in space and time of different SSEs and to provide a comparative analysis of behavior of the individual SSEs. We then discuss the potential interplay between SSE in Guerrero and Oaxaca, and the
Thermal management in inertial fusion energy slab amplifiers
NASA Astrophysics Data System (ADS)
Sutton, Steven B.; Albrecht, George F.
1995-12-01
As the technology associated with the development of solid-state drivers for inertial fusion energy (IFE) has evolved, increased emphasis has been placed on the development of an efficient approach for managing the waste heat generated in the laser media. This paper addresses the technical issues associated with the gas cooling of large aperture slabs, where the laser beam propagates through the cooling fluid. It is shown that the major consequence of proper thermal management is the introduction of simple wedge, or beam steering, into the system. Achieving proper thermal management requires careful consideration of the geometry, cooling fluid characteristics, cooling flow characteristics, as well as the thermal/mechanical/optical characteristics of the laser media. Particularly important are the effects of cooling rate variation and turbulent scattering on the system optical performance. Helium is shown to have an overwhelming advantage with respect to turbulent scattering losses. To mitigate cooling rate variations, we introduce the concept of flow conditioning. Finally, optical path length variations across the aperture are calculated. A comparison of two laser materials (S-FAP and YAG) shows the benefit of a nearly a-thermal material on optical variations in the system.
Tunable surface plasmon-polaritons in a gyroelectric slab sandwiched between two graphene layers
NASA Astrophysics Data System (ADS)
Xu, Guoding; Cao, Ming; Liu, Chang; Sun, Jian; Pan, Tao
2016-05-01
We study numerically the properties of surface plasmon-polaritons (SPPs) in a gyroelectric slab sandwiched between two graphene layers, where the external static magnetic field is applied in the Voigt geometry. It is shown that the dispersion characteristics and propagation lenghts of the SPPs for both the optical and the acoustic branches can be tuned flexibly by the external magnetic field and graphene's chemical potential, and that the nonreciprocal properties of the SPPs caused by the external magnetic field are rather obvious. The results provide a method for adjusting and improving the dispersion and propagation properties of the SPPs, which might be helpful for the design of the related plasmonic devices.
Photonic crystal slab quantum cascade detector
Reininger, Peter Schwarz, Benedikt; Harrer, Andreas; Zederbauer, Tobias; Detz, Hermann; Maxwell Andrews, Aaron; Gansch, Roman; Schrenk, Werner; Strasser, Gottfried
2013-12-09
In this Letter, we demonstrate the design, fabrication, and characterization of a photonic crystal slab quantum cascade detector (PCS-QCD). By employing a specifically designed resonant cavity, the performance of the photodetector is improved in three distinct ways. The PCS makes the QCD sensitive to surface normal incident light. It resonantly enhances the photon lifetime inside the active zone, thus increasing the photocurrent significantly. And, the construction form of the device inherently decreases the noise. Finally, we compare the characteristics of the PCS-QCD to a PCS - quantum well infrared photodetector and outline the advantages for certain fields of applications.
Automated inspection of hot steel slabs
Martin, Ronald J.
1985-01-01
The disclosure relates to a real time digital image enhancement system for performing the image enhancement segmentation processing required for a real time automated system for detecting and classifying surface imperfections in hot steel slabs. The system provides for simultaneous execution of edge detection processing and intensity threshold processing in parallel on the same image data produced by a sensor device such as a scanning camera. The results of each process are utilized to validate the results of the other process and a resulting image is generated that contains only corresponding segmentation that is produced by both processes.
Automated inspection of hot steel slabs
Martin, R.J.
1985-12-24
The disclosure relates to a real time digital image enhancement system for performing the image enhancement segmentation processing required for a real time automated system for detecting and classifying surface imperfections in hot steel slabs. The system provides for simultaneous execution of edge detection processing and intensity threshold processing in parallel on the same image data produced by a sensor device such as a scanning camera. The results of each process are utilized to validate the results of the other process and a resulting image is generated that contains only corresponding segmentation that is produced by both processes. 5 figs.
Hybrid slab-microchannel gel electrophoresis system
Balch, Joseph W.; Carrano, Anthony V.; Davidson, James C.; Koo, Jackson C.
1998-01-01
A hybrid slab-microchannel gel electrophoresis system. The hybrid system permits the fabrication of isolated microchannels for biomolecule separations without imposing the constraint of a totally sealed system. The hybrid system is reusable and ultimately much simpler and less costly to manufacture than a closed channel plate system. The hybrid system incorporates a microslab portion of the separation medium above the microchannels, thus at least substantially reducing the possibility of non-uniform field distribution and breakdown due to uncontrollable leakage. A microslab of the sieving matrix is built into the system by using plastic spacer materials and is used to uniformly couple the top plate with the bottom microchannel plate.
Slab coupled optical fiber sensor calibration
NASA Astrophysics Data System (ADS)
Whitaker, B.; Noren, J.; Chadderdon, S.; Wang, W.; Forber, R.; Selfridge, R.; Schultz, S.
2013-02-01
This paper presents a method for calibrating slab coupled optical fiber sensors (SCOS). An automated system is presented for selecting the optimal laser wavelength for use in SCOS interrogation. The wavelength calibration technique uses a computer sound card for both the creation of the applied electric field and the signal detection. The method used to determine the ratio between the measured SCOS signal and the applied electric field is also described along with a demonstration of the calibrated SCOS involving measuring the dielectric breakdown of air.
Diode-pumped high-efficiency high-brightness Q-switched ND:YAG slab laser.
Armandillo, E; Norrie, C; Cosentino, A; Laporta, P; Wazen, P; Maine, P
1997-08-01
A high-efficiency diode-pumped Q-switched Nd:YAG oscillator designed for future spaceborne applications has been demonstrated and characterized. The laser is based on a side-pumped slab geometry and uses an unstable resonator with a radially variable-reflectivity output coupler. The laser provides an output pulse energy of 100 mJ at a 100-Hz repetition rate, with a near-diffraction-limited beam and an overall electrical optical efficiency exceeding 6%. PMID:18185784
NASA Astrophysics Data System (ADS)
Bishop, B.; Beck, S. L.; Zandt, G.; Scire, A. C.; Wagner, L. S.; Long, M. D.; Tavera, H.
2014-12-01
In central Peru the combination of an unusually shallow Wadati-Benioff zone and lack of arc volcanism are indicators of flat slab subduction and are associated with both the ongoing subduction of the Nazca Ridge and the prior subduction of the Inca Plateau. Data from the PULSE experiment has allowed us to better constrain the geometry of the southern half of the Peruvian flat slab through analysis of teleseismic receiver functions, Pn and Sn phases from regional intermediate (>100 km depth) and deep (>500 km depth) earthquakes at the margins of the flat slab region, and teleseismic tomgraphy. We observe a low velocity anomaly below the subducted portion of the Nazca Ridge in the teleseismic S-wave tomography. Utilizing both Pn and Sn phases from regional intermediate and deep earthquakes at the margins of the flat slab, we have found significant travel time delays for propagation paths passing through this anomaly, confirming the presence of this low velocity anomaly under the flat slab. This anomaly likely contributes to the buoyancy of this segment of the flat slab, increasing the coupling with the upper plate. Both the teleseismic tomographic and our receiver function results indicate that the southern segment of the Peruvian flat slab extends locally more than 100 km further inboard than previous estimates. As the shallow portion of the slab inboard of the subducting Nazca Ridge is largely aseismic, these new results help to better constrain the geometry of the Peruvian flat slab as it re-subducts back into the mantle. Between 10°S and 16°S the subducted oceanic crust along the inboard projection of the Nazca Ridge lies at a depth of 60 km to 70 km while subducted crust immediately north and south of the ridge projection lies at depths of 80 km to 90 km suggesting the slab is sinking north and south of the ridge. The unusually shallow depth of the slab along the ridge's projection may indicate that the subducted Nazca Plate is coupled to the South American
Observation of the Early Transition from Slab to mixed Slab-Toroidal ETG Turbulence
NASA Astrophysics Data System (ADS)
Balbaky, Abed; Sokolov, Vladimir; Sen, Amiya K.
2014-10-01
Parametric studies of the transition between the slab branch of electron temperature gradient (ETG) mode and the mixed slab-toroidal branch of the ETG mode in CLM are reported. CLM was operated in a mirror machine configuration with a cell length of 50--100 cm, and a mirror ratio of 1--2. For typical CLM parameters and a mode localized at r = 2 cm this provides a range for inverse radius of curvature Rc - 1 between 0 and .006 cm-1. Under normal conditions theory predicts transition between slab and toroidal modes would occur when the parameter k| |Rc / 2k⊥ ρ ~ 1. Recent experiments have obtained an experimental scaling of mode amplitude and frequency as a function of Rc - 1. They indicate that even for much more modest levels of k| |Rc / 2k⊥ ρ ~ . 1 , there are substantial increases in saturated mode, up to 5 times larger than the pure slab mode. Changes in real frequency in the mode are generally small, on the order of <5%. This research was supported by the Department of Electrical Engineering of Columbia University.
Thermally induced birefringence in Nd:YAG slab lasers
Ostermeyer, Martin; Mudge, Damien; Veitch, Peter J.; Munch, Jesper
2006-07-20
We study thermally induced birefringence in crystalline Nd:YAG zigzag slab lasers and the associated depolarization losses. The optimum crystallographic orientation of the zigzag slab within the Nd:YAG boule and photoelastic effects in crystalline Nd:YAG slabs are briefly discussed. The depolarization is evaluated using the temperature and stress distributions, calculated using a finite element model, for realistically pumped and cooled slabs of finite dimensions. Jones matrices are then used to calculate the depolarization of the zigzag laser mode. We compare the predictions with measurements of depolarization, and suggest useful criteria for the design of the gain media for such lasers.
MULTIPLE SETS OF TWIN SLABS ON THE RUN OUT. THE ...
MULTIPLE SETS OF TWIN SLABS ON THE RUN OUT. THE RUN OUT INCLUDES THE TRAVELING TORCH WHICH CUTS SLABS TO DESIRED LENGTH, AN IDENTIFICATION SYSTEM TO INDICATE HEAT NUMBER AND TRACE IDENTITY OF EVERY SLAB, AND A DEBURRING DEVICE TO SMOOTH SLABS. AT LEFT OF ROLLS IS THE DUMMY BAR. DUMMY BAR IS INSERTED UP THROUGH CONTAINMENT SECTION INTO MOLD PRIOR TO START OF CAST. WHEN STEEL IS INTRODUCED INTO MOLD IT CONNECTS WITH BAR AS CAST BEGINS, AT RUN OUT DUMMY BAR DISCONNECTS AND IS STORED - U.S. Steel, Fairfield Works, Continuous Caster, Fairfield, Jefferson County, AL
MULTIPLE SETS OF TWIN SLABS ON THE RUN OUT. THE ...
MULTIPLE SETS OF TWIN SLABS ON THE RUN OUT. THE RUN OUT INCLUDES THE TRAVELING TORCH WHICH CUTS SLABS TO DESIRED LENGTH, AN IDENTIFICATION SYSTEM TO INDICATE HEAT NUMBER AND TRACE IDENTITY OF EVERY SLAB, AND A DEBURRING DEVICE TO SMOOTH SLABS. AT LEFT OF ROLLS IS THE DUMMY BAR. DUMMY BAR IS INSERTED UP THROUGH CONTAINMENT SECTION INTO MOLD PRIOR TO START OF CAST. WHEN STEEL IS INTRODUCED INTO MOLD IT CONNECTS WITH BAR AS CAST BEGINS, AT RUN OUT DUMMY BAR DISCONNECTS AND IS STORED. - U.S. Steel, Fairfield Works, Continuous Caster, Fairfield, Jefferson County, AL
Characteristics of the surface plasma wave in a self-gravitating magnetized dusty plasma slab
Lee, Myoung-Jae; Jung, Young-Dae
2015-11-15
The dispersion properties of surface dust ion-acoustic waves in a self-gravitating magnetized dusty plasma slab are investigated. The dispersion relation is derived by using the low-frequency magnetized dusty dielectric function and the surface wave dispersion integral for the slab geometry. We find that the self-gravitating effect suppresses the frequency of surface dust ion-acoustic wave for the symmetric mode in the long wavelength regime, whereas it hardly changes the frequency for the anti-symmetric mode. As the slab thickness and the wave number increase, the surface wave frequency slowly decreases for the symmetric mode but increases significantly for the anti-symmetric mode. The influence of external magnetic field is also investigated in the case of symmetric mode. We find that the strength of the magnetic field enhances the frequency of the symmetric-mode of the surface plasma wave. The increase of magnetic field reduces the self-gravitational effect and thus the self-gravitating collapse may be suppressed and the stability of dusty objects in space is enhanced.
Stripe to slab confinement for the linearization of macromolecules in nanochannels.
Benková, Zuzana; Námer, Pavol; Cifra, Peter
2015-03-21
We investigated the recently suggested advantageous analysis of chain linearization experiments with macromolecules confined in a stripe-like channel (Huang and Battacharya, EPL, 2014, 106, 18004) using Monte Carlo simulations. The enhanced chain extension in a stripe, which is due to the significant excluded volume interactions between the monomers in two dimensions, weakens considerably on transition to an experimentally feasible slit-like channel. Based on the chain extension-confinement strength dependence and the structure factor behavior for a chain in a stripe, we infer the excluded volume regime (de Gennes regime) typical for two-dimensional systems. On widening of the stripe in a direction perpendicular to the stripe plane, i.e. on the transition to the slab geometry, the advantageous chain extension decreases and a Gaussian regime is observed for not very long semiflexible chains. The evidence for pseudo-ideality in confined chains is based on four indicators: the extension curves, variation of the extension with the persistence length P, estimated limits for the regimes in the investigated systems, and the structure factor behavior. The slab behavior can be observed when the two-dimensional stripe (originally of a one-monomer thickness) reaches a reduced thickness D larger than approximately D/P ≈ 0.2 in the third dimension. This maximum height of a slab at which the advantage of a stripe is retained is very low and has implications for DNA linearization experiments. PMID:25651783
Polarization effects in active Fresnel rhomb zig-zag slab amplifier
Bikmatov, R.G.; Chernyak, V.M.; Ignat`ev, L.P.; Kuznetsov, V.G.; Pergament, M.I.; Smirnov, R.V.; Sokolov, V.I.; Hunt, J.T.; Manes, K.
1997-01-27
The concept to use a slab as active element, working in zig-zag geometry, and also as Fresnel rhomb, seems to be rather attractive. However, in this case different depolarization effects in active element arc of crucial importance. We have carried out the estimations of depolarization effects arising both due to mechanical loading of an active element at its fastening and due to thermooptical distortions. To check up these rigid requirements to depolarization (0.1 % - 0.01 %) careful measurements of depolarization effects and their sources are being carried out. Mechanical loading gives one of the main contributions in depolarization at fastening of active element. Using model experiments with glass Fresnel rhomb under mechanical loading we have measured depolarization effects. It is proposed to use additional glass plate to compensate beam depolarization in zig-zag slab. The received results allow to expect successful use of the slab amplifier as a Fresnel rhomb providing rather high quality of optical material of active clement.
Slab pull and indentation tectonics: insights from 3D laboratory experiments
NASA Astrophysics Data System (ADS)
Regard, Vincent; Faccenna, Claudio; Martinod, Joseph; Bellier, Olivier
2005-03-01
We investigate, using 3D laboratory experiments, how the dynamics of indentation process are affected by the evolution at depth of the oceanic and continental subductions. Lithospheric plates are modelled by sand-silicone plates floating on glucose syrup, and the density contrast between oceanic and continental lithospheric plates and asthenosphere is reproduced. Analogue experiments model the convergence between two lithospheric plates, a small continent indenting a large continental plate. We show that the surface deformation in front of the indenter and above the oceanic subduction zone depends on the behaviour of the slab below the collision zone. Slab break-off following the subduction of the small continent favours the indentation process, because it results in an increasing compression in front of the indenter, and extension above the neighbouring oceanic subduction, both of them being responsible for the appearance of the indenter-like geometry of the plate boundary. When the slab does not deform significantly at depth, in contrast, the closure of the oceanic domain in front of the indenter is followed by a longer period of continental subduction, during which the tectonic regime within the wide continent remains quite homogeneous. Comparing the presented analogue experiments with the subductions both part of the Arabian indenter within Eurasia, our results suggest that the different tectonic regime on both sides of the Arabia indenter may partly result from the probable occurrence of a detachment at depth under eastern Anatolia.
Plate deformation at depth under northern California: Slab gap or stretched slab?
ten Brink, U.S.; Shimizu, N.; Molzer, P.C.
1999-01-01
Plate kinematic interpretations for northern California predict a gap in the underlying subducted slab caused by the northward migration of the Pacific-North America-Juan de Fuca triple junction. However, large-scale decompression melting and asthenospheric upwelling to the base of the overlying plate within the postulated gap are not supported by geophysical and geochemical observations. We suggest a model for the interaction between the three plates which is compatible with the observations. In this 'slab stretch' model the Juan de Fuca plate under coastal northern California deforms by stretching and thinning to fill the geometrical gap formed in the wake of the northward migrating Mendocino triple junction. The stretching is in response to boundary forces acting on the plate. The thinning results in an elevated geothermal gradient, which may be roughly equivalent to a 4 Ma oceanic lithosphere, still much cooler than that inferred by the slab gap model. We show that reequilibration of this geothermal gradient under 20-30 km thick overlying plate can explain the minor Neogene volcanic activity, its chemical composition, and the heat flow. In contrast to northern California, geochemical and geophysical consequences of a 'true' slab gap can be observed in the California Inner Continental Borderland offshore Los Angeles, where local asthenospheric upwelling probably took place during the Miocene as a result of horizontal extension and rotation of the overlying plate. The elevated heat flow in central California can be explained by thermal reequilibration of the stalled Monterey microplate under the Coast Ranges, rather than by a slab gap or viscous shear heating in the mantle.
Radiative transfer model for contaminated rough slabs.
Andrieu, François; Douté, Sylvain; Schmidt, Frédéric; Schmitt, Bernard
2015-11-01
We present a semi-analytical model to simulate the bidirectional reflectance distribution function (BRDF) of a rough slab layer containing impurities. This model has been optimized for fast computation in order to analyze massive hyperspectral data by a Bayesian approach. We designed it for planetary surface ice studies but it could be used for other purposes. It estimates the bidirectional reflectance of a rough slab of material containing inclusions, overlaying an optically thick media (semi-infinite media or stratified media, for instance granular material). The inclusions are assumed to be close to spherical and constituted of any type of material other than the ice matrix. It can be any other type of ice, mineral, or even bubbles defined by their optical constants. We assume a low roughness and we consider the geometrical optics conditions. This model is thus applicable for inclusions larger than the considered wavelength. The scattering on the inclusions is assumed to be isotropic. This model has a fast computation implementation and thus is suitable for high-resolution hyperspectral data analysis. PMID:26560577
The Equivalent Slab Thickness of Mars' Ionosphere
NASA Astrophysics Data System (ADS)
Mendillo, M.; Lawler, G.; Narvaez, C.; Kofman, W.; Mouginot, J.; Morgan, D.; Gurnett, D.
2014-04-01
The integral with height of an electron density profile, called the ionospheric total electron content (TEC), is dominated by plasma near the height of maximum electron density (Nmax). The ratio τ = TEC/Nmax has a unit of distance, representing the thickness of a slab of plasma of uniform density (Nmax) with the same TEC. At Earth, the parameter τ has been found to vary far less than either Nmax or TEC, and thus models of τ can be used to generate values of TEC or Nmax when only one is observed. For an ionospheric layer dominated by photo-chemical processes, τ has also been related to the scale height of the neutral gas (H = kT/mg) that is ionized by sunlight. The MARSIS radio science package on the Mars Express satellite has produced large independent data sets of TEC and Nmax. We have used them to form slab thickness patterns versus phase of the solar cycle and solar zenith angle. The overall sample average for daytime (SZA < 90o) conditions is < τ >day = 71 ± 21 km, and for SZA between 90o - 100o, < τ >terminator = 55 ± 25 km. We will report on the possible use of τ patterns to infer characteristics of the martian ionosphere and thermosphere.
1D-1D Coulomb drag in a 6 Million Mobility Bi-layer Heterostructure
NASA Astrophysics Data System (ADS)
Bilodeau, Simon; Laroche, Dominique; Xia, Jian-Sheng; Lilly, Mike; Reno, John; Pfeiffer, Loren; West, Ken; Gervais, Guillaume
We report Coulomb drag measurements in vertically-coupled quantum wires. The wires are fabricated in GaAs/AlGaAs bilayer heterostructures grown from two different MBE chambers: one at Sandia National Laboratories (1.2M mobility), and the other at Princeton University (6M mobility). The previously observed positive and negative drag signals are seen in both types of devices, demonstrating the robustness of the result. However, attempts to determine the temperature dependence of the drag signal in the 1D regime proved challenging in the higher mobility heterostructure (Princeton), in part because of difficulties in aligning the wires within the same transverse subband configuration. Nevertheless, this work, performed at the Microkelvin laboratory of the University of Florida, is an important proof-of-concept for future investigations of the temperature dependence of the 1D-1D drag signal down to a few mK. Such an experiment could confirm the Luttinger charge density wave interlocking predicted to occur in the wires. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under Contract DE-AC04-94AL8500.
Subduction Zone Science - Examples of Seismic Images of the Central Andes and Subducting Nazca Slab
NASA Astrophysics Data System (ADS)
Beck, S. L.; Zandt, G.; Scire, A. C.; Ward, K. M.; Portner, D. E.; Bishop, B.; Ryan, J. C.; Wagner, L. S.; Long, M. D.
2015-12-01
Subduction has shaped large regions of the Earth and constitute over 55,000 km of convergent plate margin today. The subducting slabs descend from the surface into the lower mantle and impacts earthquake occurrence, surface uplift, arc volcanism and mantle convection as well as many other processes. The subduction of the Nazca plate beneath the South America plate is one example and constitutes the largest present day ocean-continent convergent margin system and has built the Andes, one of the largest actively growing mountain ranges on Earth. This active margin is characterized by along-strike variations in arc magmatism, upper crustal shortening, crustal thickness, and slab geometry that make it an ideal region to study the relationship between the subducting slab, the mantle wedge, and the overriding plate. After 20 years of portable seismic deployments in the Central Andes seismologists have combined data sets and used multiple techniques to generate seismic images spanning ~3000 km of the South American subduction zone to ~800 km depth with unprecedented resolution. For example, using teleseismic P- waves we have imaged the Nazca slab penetrating through the mantle transition zone (MTZ) and into the uppermost lower mantle. Our tomographic images show that there is significant along-strike variation in the morphology of the Nazca slab in the upper mantle, MTZ, and the lower mantle, including possible tears, folding, and internal deformation. Receiver function studies and surface wave tomography have revealed major changes in lithospheric properties in the Andes. Improved seismic images allow us to more completely evaluate tectonic processes in the formation and uplift of the Andes including: (1) overthickened continental crust driven by crustal shortening, (2) changes in slab dip and coupling with the overlying plate (3) localized lithospheric foundering, and (4) large-scale mantle and crustal melting leading to magmatic addition and/or crustal flow. Although
Okal, E.A.; Kirby, S.H.
1998-01-01
Previous work has suggested that many of the deep earthquakes beneath the Fiji Basin occur in slab material that has been detached and foundered to the bottom of the transition zone or has been laid down by trench migration in a similar recumbent position. Since nowhere else in the Earth do so many earthquakes occur in slabs stagnated in the transition zone, these earthquakes merit closer study. Accordingly, we have assembled from historical and modern data a comprehensive catalogue of the relocated hypocenters and focal mechanisms of well-located deep events in the geographic area between the bottoms of the main Vanuatu and Tonga Wadati-Benioff zones. Two regions of deep seismogenesis are recognized there: (i) 163 deep shocks have occurred north of 15??S in the Vityaz Group from 1949 to 1996. These seismological observations and the absence of other features characteristic of active subduction suggest that the Vityaz group represents deep failure in a detached slab that has foundered to a horizontal orientation near the bottom of the transition zone. (ii) Another group of nearly 50 'outboard' deep shocks occur between about 450 and 660 km depth, west of the complexly buckled and offset western edge of the Tonga Wadati-Benioff zone. Their geometry is in the form of two or possibly three small-circle arcs that roughly parallel the inferred motion of Tonga trench migration. Earthquakes in the southernmost of these arcs occur in a recumbent high-seismic-wavespeed slab anomaly that connects both to the main inclined Tonga anomaly to the east and a lower mantle anomaly to the west [Van der Hilst, R., 1995. Complex morphology of subducted lithosphere in the mantle beneath the Tonga trench. Nature, Vol. 374, pp. 154-157.]. Both groups show complexity in their focal mechanisms. The major question raised by these observations is the cause of this apparent temporary arrest in the descent of the Tonga slab into the lower mantle. We approach these questions by considering the
Van, Anh T; Aksoy, Murat; Holdsworth, Samantha J; Kopeinigg, Daniel; Vos, Sjoerd B; Bammer, Roland
2014-01-01
Purpose To propose a method for mitigating slab boundary artifacts in 3D multislab diffusion imaging with no or minimal increases in scan time. Methods The multislab acquisition was treated as parallel imaging acquisition where the slab profiles acted as the traditional receiver sensitivity profiles. All the slabs were then reconstructed simultaneously along the slab direction using Cartesian-based sensitivity encoding (SENSE) reconstruction. The slab profile estimation was performed using either a Bloch simulation or a calibration scan. Results Both phantom and in vivo results showed negligible slab boundary artifacts after reconstruction using the proposed method. The performance of the proposed method is comparable to the state-of-the-art slab combination method without the scan time penalty that depends on the number of acquired volumes. The obtained g-factor map of the SENSE reconstruction problem showed a maximum g-factor of 1.7 in the region of interest. Conclusion We proposed a novel method for mitigating slab boundary artifacts in 3D diffusion imaging by treating the multislab acquisition as a parallel imaging acquisition and reconstructing all slabs simultaneously using Cartesian SENSE. Unlike existing methods, the scan time increase, if any, does not scale with the number of image volumes acquired. PMID:24691843
ERIC Educational Resources Information Center
Cukier, Mimi; Asdourian, Tony; Thakker, Anand
2012-01-01
Geometry provides a natural window into what it is like to do mathematics. In the world of geometry, playful experimentation is often more fruitful than following a procedure, and logic plus a few axioms can open new worlds. Nonetheless, teaching a geometry course in a way that combines both rigor and play can be difficult. Many geometry courses…
ERIC Educational Resources Information Center
Kuntz, Gilles
The first section of this paper on World Wide Web applications related to dynamic geometry addresses dynamic geometry and teaching, including the relationship between dynamic geometry and direct manipulation, key features of dynamic geometry environments, the importance of direct engagement of the learner using construction software for…
Understanding 1D Electrostatic Dust Levitation
NASA Astrophysics Data System (ADS)
Hartzell, C. M.; Scheeres, D. J.
2011-12-01
Electrostatically-dominated dust motion has been hypothesized since the Lunar Horizon Glow was observed by the Surveyor spacecraft. The hypothesized occurence of this phenomenon was naturally extended to asteroids due to their small gravities. Additionally, it has been suggested that the dust ponds observed on Eros by the NEAR mission may be created by electrostatically-dominated dust transport. Previous attempts to numerically model dust motion on the Moon and Eros have been stymied by poorly understood dust launching mechanisms. As a result, the initial velocity and charge of dust particles used in numerical simulations may or may not have any relevance to the actual conditions occurring in situ. It has been seen that properly tuned initial states (velocity and charge) result in dust particles levitating above the surface in both 1D and 2D simulations. Levitation is of interest to planetary scientists since it provides a way to quickly redistribute the surface dust particles over a body. However, there is currently no method to predict whether or not a certain initial state will result in levitation. We have developed a method to provide constraints on the initial states that result in levitation as a function of dust particle size and central body gravity. Additionally, our method can be applied to several models of the plasma sheath. Thus, we limit the guesswork involved in determining which initial conditions result in levitation. We provide a more detailed understanding of levitation phenomena couched in terms of the commonly recognized spring-mass system. This method of understanding dust motion removes the dependency on the launching mechanism, which remains fraught with controversy. Once a feasible dust launching mechanism is identified (be it micrometeoroid bombardment or electrostatic lofting), our method will allow the community to quickly ascertain if dust levitation will occur in situ or if it is simply a numerical artifact. In addition to
NASA Astrophysics Data System (ADS)
Tian, You; Zhao, Dapeng
2012-06-01
We used 190,947 high-quality P-wave arrival times from 8421 local earthquakes and 1,098,022 precise travel-time residuals from 6470 teleseismic events recorded by the EarthScope/USArray transportable array to determine a detailed three-dimensional P-wave velocity model of the crust and mantle down to 1000 km depth under the western United States (US). Our tomography revealed strong heterogeneities in the crust and upper mantle under the western US. Prominent high-velocity anomalies are imaged beneath Idaho Batholith, central Colorado Plateau, Cascadian subduction zone, stable North American Craton, Transverse Ranges, and Southern Sierra Nevada. Prominent low-velocity anomalies are imaged at depths of 0-200 km beneath Snake River Plain, which may represent a small-scale convection beneath the western US. The low-velocity structure deviates variably from a narrow vertical plume conduit extending down to ˜1000 km depth, suggesting that the Yellowstone hotspot may have a lower-mantle origin. The Juan de Fuca slab is imaged as a dipping high-velocity anomaly under the western US. The slab geometry and its subducted depth vary in the north-south direction. In the southern parts the slab may have subducted down to >600 km depth. A "slab hole" is revealed beneath Oregon, which shows up as a low-velocity anomaly at depths of ˜100 to 300 km. The formation of the slab hole may be related to the Newberry magmatism. The removal of flat subducted Farallon slab may have triggered the vigorous magmatism in the Basin and Range and southern part of Rocky Mountains and also resulted in the uplift of the Colorado Plateau and Rocky Mountains.
Reconciling surface plate motions with rapid three-dimensional mantle flow around a slab edge.
Jadamec, Margarete A; Billen, Magali I
2010-05-20
The direction of tectonic plate motion at the Earth's surface and the flow field of the mantle inferred from seismic anisotropy are well correlated globally, suggesting large-scale coupling between the mantle and the surface plates. The fit is typically poor at subduction zones, however, where regional observations of seismic anisotropy suggest that the direction of mantle flow is not parallel to and may be several times faster than plate motions. Here we present three-dimensional numerical models of buoyancy-driven deformation with realistic slab geometry for the Alaska subduction-transform system and use them to determine the origin of this regional decoupling of flow. We find that near a subduction zone edge, mantle flow velocities can have magnitudes of more than ten times the surface plate motions, whereas surface plate velocities are consistent with plate motions and the complex mantle flow field is consistent with observations from seismic anisotropy. The seismic anisotropy observations constrain the shape of the eastern slab edge and require non-Newtonian mantle rheology. The incorporation of the non-Newtonian viscosity results in mantle viscosities of 10(17) to 10(18) Pa s in regions of high strain rate (10(-12) s(-1)), and this low viscosity enables the mantle flow field to decouple partially from the motion of the surface plates. These results imply local rapid transport of geochemical signatures through subduction zones and that the internal deformation of slabs decreases the slab-pull force available to drive subducting plates. PMID:20485433
Imaging the slab beneath central Chile using the Spectral Elements Method and adjoint techniques
NASA Astrophysics Data System (ADS)
Mercerat, E. D.; Nolet, G.; Marot, M.; Deshayes, P.; Monfret, T.
2010-12-01
This work focuses on imaging the subducting slab beneath Central Chile using novel inversion techniques based on the adjoint method and accurate wave propagation simulations using the Spectral Elements Method. The study area comprises the flat slab portion of the Nazca plate between 29 S and 34 S subducting beneath South America. We will use a database of regional seismicity consisting of both crustal and deep slab earthquakes with magnitude 3 < Mw < 6 recorded by different temporary and permanent seismological networks. Our main goal is to determine both the kinematics and the geometry of the subducting slab in order to help the geodynamical interpretation of such particular active margin. The Spectral Elements Method (SPECFEM3D code) is used to generate the synthetic seismograms and it will be applied for the iterative minimization based on adjoint techniques. The numerical mesh is 600 km x 600 km in horizontal coordinates and 220 km depth. As a first step, we are faced to well-known issues concerning mesh generation (resolution, quality, absorbing boundary conditions). In particular, we must evaluate the influence of free surface topography, as well as the MOHO and other geological interfaces in the synthetic seismograms. The initial velocity model from a previous travel-time tomography study, is linearly interpolated to the Gauss-Lobatto-Legendre grid. The comparison between the first forward simulations (up to 4 seconds minimum period) validate the initial velocity model of the study area, although many features not reproduced by the initial model have already been identified. Next step will concentrate in the comparison between finite-frequency kernels calculated by travel-time methods with ones based on adjoint methods, in order to highlight advantages and disadvantages in terms of resolution, accuracy, but also computational cost.
DESIGN OF A CONCRETE SLAB FOR STORAGE OF SNF AND HLW CASKS
J. Bisset
2005-02-14
This calculation documents the design of the Spent Nuclear Fuel (SNF) and High-Level Waste (HLW) Cask storage slab for the Aging Area. The design is based on the weights of casks that may be stored on the slab, the weights of vehicles that may be used to move the casks, and the layout shown on the sketch for a 1000 Metric Ton of Heavy Metal (MTHM) storage pad on Attachment 2, Sht.1 of the calculation 170-C0C-C000-00100-000-00A (BSC 2004a). The analytical model used herein is based on the storage area for 8 vertical casks. To simplify the model, the storage area of the horizontal concrete modules and their related shield walls is not included. The heavy weights of the vertical storage casks and the tensile forces due to pullout at the anchorages will produce design moments and shear forces that will envelope those that would occur in the storage area of the horizontal modules. The design loadings will also include snow and live loads. In addition, the design will also reflect pertinent geotechnical data. This calculation will document the preliminary thickness and general reinforcing steel requirements for the slab. This calculation also documents the initial design of the cask anchorage. Other slab details are not developed in this calculation. They will be developed during the final design process. The calculation also does not include the evaluation of the effects of cask drop loads. These will be evaluated in this or another calculation when the exact cask geometry is known.
XCHEM-1D: A Heat Transfer/Chemical Kinetics Computer Program for multilayered reactive materials
Gross, R.J.; Baer, M.R.; Hobbs, M.L.
1993-10-01
An eXplosive CHEMical kinetics code, XCHEM, has been developed to solve the reactive diffusion equations associated with thermal ignition of energetic materials. This method-of-lines code uses stiff numerical methods and adaptive meshing to resolve relevant combustion physics. Solution accuracy is maintained between multilayered materials consisting of blends of reactive components and/or inert materials. Phase change and variable properties are included in one-dimensional slab, cylindrical and spherical geometries. Temperature-dependent thermal properties have been incorporated and the modification of thermal conductivities to include decomposition effects are estimated using solid/gas volume fractions determined by species fractions. Gas transport properties, including high pressure corrections, have also been included. Time varying temperature, heat flux, convective and thermal radiation boundary conditions, and layer to layer contact resistances have also been implemented.
3D/1D Analysis of ICRF Antennas
NASA Astrophysics Data System (ADS)
Maggiora, Riccardo; Lancellotti, Vito; Vecchi, Giuseppe
2003-10-01
An innovative tool has been realized for the 3D/1D simulation of Ion Cyclotron Radio Frequency (ICRF), i.e. accounting for antennas in a realistic 3D geometry and with an accurate 1D plasma model. The approach to the problem is based on an integral-equation formulation for the self-consistent evaluation of the current distribution on the conductors. The environment has been subdivided in two coupled region: the plasma region and the vacuum region. The two problems are linked by means of a magnetic current (electric field) distribution on the aperture between the two regions. In the vacuum region all the calculations are executed in the spatial domain while in the plasma region an extraction in the spectral domain of some integrals is employed that permits to significantly reduce the integration support and to obtain a high numerical efficiency leading to the practical possibility of using a large number of sub-domain (rectangular or triangular) basis functions on each solid conductor of the system. The plasma enters the formalism of the plasma region via a surface impedance matrix; for this reason any plasma model can be used; at present the FELICE code has been adopted, that affords density and temperature profiles, and FLR effects. The source term directly models the TEM mode of the coax feeding the antenna and the current in the coax is determined self-consistently, giving the input impedance/admittance of the antenna itself. Calculation of field distributions (both magnetic and electric), useful for sheath considerations, is included. This tool has been implemented in a suite, called TOPICA, that is modular and applicable to ICRF antenna structures of arbitrary shape. This new simulation tool can assist during the detailed design phase and for this reason can be considered a "Virtual Prototyping Laboratory" (VPL). The TOPICA suite has been tested against assessed codes and against measurements and data of mock-ups and existing antennas. The VPL is being used in
Integrated multicolor detector utilizing 1D photonic bandgap filter with wedge-shaped defect
NASA Astrophysics Data System (ADS)
Jaksic, Zoran S.; Petrovic, Radomir; Randjelovic, Danijela; Dankovic, Tatjana; Djuric, Zoran G.; Ehrfeld, Wolfgang; Schmidt, Andreas; Hecker, Karl H.
1999-03-01
We propose a single-chip multicolor photodetector for micrometers range based on a linear IR semiconductor detector array with an integrated 1D photonic bandgap (PBG) filter. A wedge- shaped defect slab is introduced into the filer instead of one of the layers. The bandgap of the photonic crystal coincides with the spectral sensitivity range of the photodetector array, while the built-in defect gives a transmission peak within the same range. The defect thickness varies along the array length and thus shifts the transmission peak wavelength. The optimized photonic bandgap filter including defect is designed using the transfer matrix method. The peak frequency is tuned by choosing the geometrical parameters of the wedge-shaped defect. In our experiments, thin alternating Si and SiO2 films are sputtered onto the array surface, thus forming a 1D PBG structure. The defect is fabricated by gradually changing the middle Si layer thickness over the width of the array. Its wedge-forming is performed by micromachining or, alternatively, by in-situ oblique deposition within the sputtering system and, possibly, subsequent chemomechanical polishing. The characteristics of the finished PBG structure are measured using an IR spectrophotometer. An increase of the number of PBG layers improves the confinement of transmission peaks and thus decreases the crosstalk between the array elements. Although our multicolor detector is designed for the (3-5) micrometers atmospheric window, it can be straightforward redesigned for any other optical range.
Seismicity and Geometry Properties of the Hellenic Subduction Zone
NASA Astrophysics Data System (ADS)
Papadimitriou, E. E.; Karakostas, V. G.; Vallianatos, F.; Makropoulos, C.; Drakatos, G.
2014-12-01
Recent seismicity and fault plane solutions of earthquakes that occurred along the Hellenic Arc-Trench system are engaged for approximating the geometry of the subducted oceanic plate. Seismicity and focal mechanisms confirm the gentle subduction (~15o-20o) of the oceanic crust reaching a depth of 20 km at a distance of 100 km from the trench. The slab is then bending at larger angles, and in particular at ~45o up to the depth of 80 km and at ~65o up to the depth of 180 km, when seismicity ceased. This geometry of the slab is shown in a bunch of cross sections normal to the convergence strike, up to ~25o (east Crete Island). To the east the sparse inslab seismicity reveals an almost vertical dipping of the lower part (from 80 km downdip) of the descending slab. The slab interface that accommodates hazardous earthquakes is clearly nonplanar with the main seismic moment release taking place on its up-dip side. The fore-arc, upper plate seismicity, is remarkably low in comparison with both subduction and back arc seismicity, and confined inside a seismogenic layer having a width not exceeding the 20km. Offshore seismicity is spatially variable forming distinctive streaks thus revealing that parts of the oceanic crust are probably slipped aseismically. This observation along with the fact that coupling in the Hellenic arc is only about a tenth of the plate motion, imply the presence of areas of lower and higher coupling across the subduction interface. Areas of high coupling imply areas of the slab interface subjected to high normal forces and correlate with earthquake asperities. Although asperity distributions vary substantially through time, identification of such characteristics in the seismogenesis can have a significant impact in the seismic hazard assessment. This research has been co-funded by the European Union (European Social Fund) and Greek national resources under the framework of the "THALES Program: SEISMO FEAR HELLARC" project.
42. PRESSING A SLAB OF CLAY ONTO A MOSAIC MOLD. ...
42. PRESSING A SLAB OF CLAY ONTO A MOSAIC MOLD. THE MOLD, WHICH HAS A RAISED DESIGN, LEAVES AND OUTLINE IN THE SLAB, THE PIECES THUS DEFINED, ARE THEN CUT APART TO BE FIRED SEPARATELY AND REASSEMBLED. - Moravian Pottery & Tile Works, Southwest side of State Route 313 (Swamp Road), Northwest of East Court Street, Doylestown, Bucks County, PA
Scattering of electromagnetic waves from a turbulent plasma slab.
NASA Technical Reports Server (NTRS)
Liu, C. H.
1972-01-01
Scattering of electromagnetic waves from a turbulent plasma slab is studied. Part of the effects of the multiple scattering is taken into account. The reflection coefficient is found to be increased and its variation with respect to the slab thickness is smoothed out by the random scattering.
Learning Geometry through Dynamic Geometry Software
ERIC Educational Resources Information Center
Forsythe, Sue
2007-01-01
In this article, the author investigates effective teaching and learning of geometrical concepts using dynamic geometry software (DGS). Based from her students' reactions to her project, the author found that her students' understanding of the concepts was better than if they had learned geometry through paper-based tasks. However, mixing computer…
Flat slab deformation caused by interplate suction force
NASA Astrophysics Data System (ADS)
Ma, Yiran; Clayton, Robert W.
2015-09-01
We image the structure at the southern end of the Peruvian flat subduction zone, using receiver function and surface wave methods. The Nazca slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes the dipping subduction. The flat slab closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the initial half of the flat slab, and the velocity resumes to normal values before the slab steepens again, indicating the resumption of dehydration and ecologitization. Two prominent midcrust structures are revealed in the 70 km thick crust under the Central Andes: molten rocks beneath the Western Cordillera and the underthrusting Brazilian Shield beneath the Eastern Cordillera.
Nonlinear Shear Instabilities in an Infinite Slab
NASA Astrophysics Data System (ADS)
Nepveu, M.
1982-08-01
The dynamical evolution of an infinite slab moving in denser and noisy (turbulent) surroundings is investigated with a 2D hydrodynamic code. The applicability of the results to astrophysical jets is discussed. Inviscid beams show internal shocks at regular intervals of a few beamwidths. Kinks are not obvious. In viscous beams shocks are less outspoken, but turbulence is triggered with maximum scales of a few beamwidths. These viscous beams broaden. Linear stability analysis may hold up to a few e-folding times, although the seed disturbance field is not infinitesimal. The computations suggest that viscous astrophysical beams may become blurred quite suddenly and may give rise to sudden change in radiation patterns (NGC 1265).
Hybrid slab-microchannel gel electrophoresis system
Balch, J.W.; Carrano, A.V.; Davidson, J.C.; Koo, J.C.
1998-05-05
A hybrid slab-microchannel gel electrophoresis system is described. The hybrid system permits the fabrication of isolated microchannels for biomolecule separations without imposing the constraint of a totally sealed system. The hybrid system is reusable and ultimately much simpler and less costly to manufacture than a closed channel plate system. The hybrid system incorporates a microslab portion of the separation medium above the microchannels, thus at least substantially reducing the possibility of non-uniform field distribution and breakdown due to uncontrollable leakage. A microslab of the sieving matrix is built into the system by using plastic spacer materials and is used to uniformly couple the top plate with the bottom microchannel plate. 4 figs.