Accelerating DynEarthSol3D on tightly coupled CPU-GPU heterogeneous processors
NASA Astrophysics Data System (ADS)
Ta, Tuan; Choo, Kyoshin; Tan, Eh; Jang, Byunghyun; Choi, Eunseo
2015-06-01
DynEarthSol3D (Dynamic Earth Solver in Three Dimensions) is a flexible, open-source finite element solver that models the momentum balance and the heat transfer of elasto-visco-plastic material in the Lagrangian form using unstructured meshes. It provides a platform for the study of the long-term deformation of earth's lithosphere and various problems in civil and geotechnical engineering. However, the continuous computation and update of a very large mesh poses an intolerably high computational burden to developers and users in practice. For example, simulating a small input mesh containing around 3000 elements in 20 million time steps would take more than 10 days on a high-end desktop CPU. In this paper, we explore tightly coupled CPU-GPU heterogeneous processors to address the computing concern by leveraging their new features and developing hardware-architecture-aware optimizations. Our proposed key optimization techniques are three-fold: memory access pattern improvement, data transfer elimination and kernel launch overhead minimization. Experimental results show that our proposed implementation on a tightly coupled heterogeneous processor outperforms all other alternatives including traditional discrete GPU, quad-core CPU using OpenMP, and serial implementations by 67%, 50%, and 154% respectively even though the embedded GPU in the heterogeneous processor has significantly less number of cores than high-end discrete GPU.
Centroid Moment Tensor Inversion in a 3D heterogeneous Earth: Application to the Australasian region
NASA Astrophysics Data System (ADS)
Hejrani, B.; Tkalcic, H.; Fichtner, A.
2015-12-01
Australia is surrounded by active complex tectonic belts causing significant seismicity. The recent expansion of permanent seismic networks in the Australasian region provides great opportunity to study Earth structure and a great variety of physical mechanisms responsible for earthquakes.On one hand, a better understanding of the Australasian lithosphere, which is now available through tomographic images from full waveform modelling (Fichtner et al. 2010), provides a powerful tool to scrutinize the determination of earthquake source parameters. Even at relatively long periods (40-200s), the 3D effects of regional structure were shown to significantly alter the global centroid moment tensor solutions (Hingee et al. 2012). Thus, we can now explore other types of uncertainties and test the accuracy of global centroid moment tensor (GCMT) solution for the earthquakes in the Australasian region while checking for the systematic inconsistencies in the solutions. This has a significant bearing on tectonic interpretations. For example, azimuth and plunge of fault planes can be investigated in search for systematic biases.On the other hand, the time has come to take a full advantage of the 3D Earth structural model and embrace ongoing advances in computational power and storage. We develop a semi-automated procedure to calculate the Centroid Moment Tensors in a 3D heterogeneous Earth. We utilize the reciprocity theorem to create Green's functions for point sources covering seismogenic zones of Australasia. We focus on improving the capacity of the method to fully complement the existing monitoring tools at Geosciences Australia. Furthermore, we investigate the effects of detailed velocity structure on Centroid location and double-couple percentages. Moreover Azimuth and Plunge of focal mechanisms in GCMT (Global CMT), were investigated in search for any systematic bias.References: Fichtner, A., Kennett, B.L.N., Igel, H., Bunge, H.-P., 2010. Full waveform tomography for
NASA Astrophysics Data System (ADS)
Duru, K.; Dunham, E. M.; Bydlon, S. A.; Radhakrishnan, H.
2014-12-01
Dynamic propagation of shear ruptures on a frictional interface is a useful idealization of a natural earthquake.The conditions relating slip rate and fault shear strength are often expressed as nonlinear friction laws.The corresponding initial boundary value problems are both numerically and computationally challenging.In addition, seismic waves generated by earthquake ruptures must be propagated, far away from fault zones, to seismic stations and remote areas.Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods.We present a numerical method for:a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration;b) dynamic propagation of earthquake ruptures along rough faults; c) accurate propagation of seismic waves in heterogeneous media with free surface topography.We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts finite differences in space. The finite difference stencils are 6th order accurate in the interior and 3rd order accurate close to the boundaries. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme. We have performed extensive numerical experiments using a slip-weakening friction law on non-planar faults, including recent SCEC benchmark problems. We also show simulations on fractal faults revealing the complexity of rupture dynamics on rough faults. We are presently extending our method to rate-and-state friction laws and off-fault plasticity.
NASA Astrophysics Data System (ADS)
Okaya, D. A.; Van Avendonk, H. J.
2013-12-01
Recent anisotropy studies at scales ranging from crust to full mantle have recognized the importance of 3D anisotropy geometry and heterogeneity as well as variability in anisotropic symmetry and orientation (tilt) of the Earth. The strong relationship between seismic anisotropy and geodynamic processes highlights the need to construct realistic Earth models that can explain observations of anisotropy in modern seismic data sets. For example, ray paths through a mantle slab window or a mountain belt may show that the crust or mantle exhibits low-order anisotropy due to a history of deformation and the development of tectonic fabrics. Observed traveltimes might not be fit with simple Transverse Isotropy (TI), so realistic calculations require an Earth model that accurately describes the wave speeds of compressional and shear waves. We have developed an anisotropic traveltime solver that allows for full 3D heterogeneity of anisotropy tensors, degrees of symmetry, and arbitrary orientation. This traveltime solver is based on the robust shortest path method (SPM) and a ray-bending algorithm that were previously applied to isotropic media (e.g., Van Avendonk et al., 2001). Instead of using an isotropic description of the seismic wave velocity, we define the full elastic tensor at each location in the model. The directional seismic velocity can subsequently be extracted using solutions of the Christoffel equations. For computational efficiency, we calculate all directional seismic velocities at each model node before the start of ray tracing. As we calculate a new ray segment, this information is quickly retrieved. We use these directional velocity maps to separately describe the propagation of compressional (P) and shear (S) body waves in anisotropic media and to subsequently calculate their traveltimes. Patterns within the velocity maps represent tensor symmetries and tilts, allowing for the construction of discretized large-scale 3D LPO flow fields or fabric
Energy Science and Technology Software Center (ESTSC)
2013-10-01
Earth3D is a computer code designed to allow fast calculation of seismic rays and travel times through a 3D model of the Earth. LLNL is using this for earthquake location and global tomography efforts and such codes are of great interest to the Earth Science community.
NASA Astrophysics Data System (ADS)
Stocchi, P.; van der Wal, W.; Vermeersen, B. L.; van de Wal, R.; Wu, P. P.
2010-12-01
Constraining the rheology of the upper mantle beneath the Scandinavian shield by means of relative sea-level indicators and instrumental observations demands a realistic spatio-temporal discretization of the late Pleistocene ice coverage over Eurasia. The waxing and waning of ice sheets is testified by the occurrence of surface marks that come in form of datable depositional and erosive features. We use these georeferenced geological evidences to infer the time-dependent ice sheet margins. By adopting both plastic and viscosplastic rheologies, we fill the boundaries with ice mass resting at equilibrium. Thus, in our approach the ice model is not biased by assumptions on the Earth’s rheology. The Fennoscandian ice cap contributes up to ~14 of the ~38 m of equivalent sea level stored over Eurasia, while 82 m are distributed between simplified parabolic ice caps for North America, Greenland and Antarctica in order to conserve the 120 m of global eustatic sea level rise. We apply the equal-volume plastic and viscoplastic versions of our ice model to a 3d spherical Earth model with rheology based on laboratory derived flow laws for mantle rocks. Fixed deformation parameters for diffusion and dislocation creep are used throughout the mantle, except for the top 250 km in which deformation parameters are computed from lateral varying temperature derived from seismology. We vary grain sizes and water content in the top layers to study their impact on glacial isostatic adjustment (GIA) observables. Predictions are compared with GIA observations in Scandinavia to show whether the ice model and laterally varying Earth rheology can explain GIA observations as good as previous models that were based on Newtonian rheology in the Earth.
Pollitz, F.F.
2002-01-01
I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ???0.01-0.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f ??????0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone.
NASA Astrophysics Data System (ADS)
Spada, G.; Antonioli, A.; Cianetti, S.; Giunchi, C.
2006-05-01
The response of the Earth to the melting of the Late Pleistocene ice sheets is commonly studied by spherically layered models, based on well-established analytical methods. In parallel, a few models have been recently proposed to circumvent the limitations imposed by spherical symmetry, and to reproduce the actual structure of the lithosphere and of the upper mantle. Their main outcome is that laterally varying rheological structures may significantly affect various geophysical quantities related to glacial isostatic adjustment (GIA), and particularly post-glacial relative sea-level (RSL) variations and 3-D crustal velocities in formerly ice-covered regions. In this paper, we contribute to the ongoing debate about the role of lithospheric and mantle heterogeneities by new 3-D spherical Newtonian finite elements models and we directly compare their outcomes with publicly available global RSL data. This differs from previous investigations, in that have mainly focused on extensive sensitivity analyses or have considered a limited number of RSL observations from formerly glaciated regions and their periphery. In our study the lithospheric thickness mimics the global structure of the cratons based on geological evidence, and the upper mantle includes a low-viscosity zone beneath the oceanic lithosphere. We use two distinct global surface loads, based upon the ICE1 and ICE3G deglaciation chronologies, respectively. Our main finding is that using all of the available RSL observations in the last 6000 years it is not possible to discern between homogeneous and heterogeneous GIA models. This result, which holds for both ICE1 and ICE3G, suggests that the cumulative effects of laterally varying structures on the synthetic RSL curves cancel out globally, yielding signals that do not significantly differ from those based on the 1-D models. We have also considered specific subsets of the global RSL database, sharing similar geographical settings and distances from the main
NASA Technical Reports Server (NTRS)
Kulikov, anton I.; Doronila, Paul R.; Nguyen, Viet T.; Jackson, Randal K.; Greene, William M.; Hussey, Kevin J.; Garcia, Christopher M.; Lopez, Christian A.
2013-01-01
Eyes on the Earth 3D software gives scientists, and the general public, a realtime, 3D interactive means of accurately viewing the real-time locations, speed, and values of recently collected data from several of NASA's Earth Observing Satellites using a standard Web browser (climate.nasa.gov/eyes). Anyone with Web access can use this software to see where the NASA fleet of these satellites is now, or where they will be up to a year in the future. The software also displays several Earth Science Data sets that have been collected on a daily basis. This application uses a third-party, 3D, realtime, interactive game engine called Unity 3D to visualize the satellites and is accessible from a Web browser.
NASA Technical Reports Server (NTRS)
Hogan, Patrick; Maxwell, Christopher; Kim, Randolph; Gaskins, Tom
2007-01-01
World Wind allows users to zoom from satellite altitude down to any place on Earth, leveraging high-resolution LandSat imagery and SRTM (Shuttle Radar Topography Mission) elevation data to experience Earth in visually rich 3D. In addition to Earth, World Wind can also visualize other planets, and there are already comprehensive data sets for Mars and the Earth's moon, which are as easily accessible as those of Earth. There have been more than 20 million downloads to date, and the software is being used heavily by the Department of Defense due to the code s ability to be extended and the evolution of the code courtesy of NASA and the user community. Primary features include the dynamic access to public domain imagery and its ease of use. All one needs to control World Wind is a two-button mouse. Additional guides and features can be accessed through a simplified menu. A JAVA version will be available soon. Navigation is automated with single clicks of a mouse, or by typing in any location to automatically zoom in to see it. The World Wind install package contains the necessary requirements such as the .NET runtime and managed DirectX library. World Wind can display combinations of data from a variety of sources, including Blue Marble, LandSat 7, SRTM, NASA Scientific Visualization Studio, GLOBE, and much more. A thorough list of features, the user manual, a key chart, and screen shots are available at http://worldwind.arc.nasa.gov.
Visualization of 3D Geological Models on Google Earth
NASA Astrophysics Data System (ADS)
Choi, Y.; Um, J.; Park, M.
2013-05-01
Google Earth combines satellite imagery, aerial photography, thematic maps and various data sets to make a three-dimensional (3D) interactive image of the world. Currently, Google Earth is a popular visualization tool in a variety of fields and plays an increasingly important role not only for private users in daily life, but also for scientists, practitioners, policymakers and stakeholders in research and application. In this study, a method to visualize 3D geological models on Google Earth is presented. COLLAborative Design Activity (COLLADA, an open standard XML schema for establishing interactive 3D applications) was used to represent different 3D geological models such as borehole, fence section, surface-based 3D volume and 3D grid by triangle meshes (a set of triangles connected by their common edges or corners). In addition, we designed Keyhole Markup Language (KML, the XML-based scripting language of Google Earth) codes to import the COLLADA files into the 3D render window of Google Earth. The method was applied to the Grosmont formation in Alberta, Canada. The application showed that the combination of COLLADA and KML enables Google Earth to effectively visualize 3D geological structures and properties.; Visualization of the (a) boreholes, (b) fence sections, (c) 3D volume model and (d) 3D grid model of Grossmont formation on Google Earth
A 3-D perspective of dynamic behaviour of heterogeneous solids
NASA Astrophysics Data System (ADS)
Lu, Yong; Zhou, Rongxin
2015-09-01
The dynamic behaviour of concrete-like materials under high strain rates has been a subject of continuous scrutiny over the years. A prevailing explanation attributes much of the dynamic increase of strength, especially under compression, to the macroscopic inertia confinement. Studies conducted by the authors' group using meso-scale computational models suggest that the heterogeneity of the material composition, in particular the involvement of the aggregates, also plays a sensible part in the process of damage evolution and the increase of the bulk strength under high strain rates, and a detailed investigation into this effect would benefit if a realistic representation of the heterogeneity in 3D can be achieved. This paper presents some recent progress in the development of a 3-D meso-scale computational model incorporating randomly-shaped 3-D aggregate particles, including the general validation of the model, and application in the simulation of the dynamic response of concrete under high strain rate compression.
Quantifying Interparticle Forces and Heterogeneity in 3D Granular Materials.
Hurley, R C; Hall, S A; Andrade, J E; Wright, J
2016-08-26
Interparticle forces in granular materials are intimately linked to mechanical properties and are known to self-organize into heterogeneous structures, or force chains, under external load. Despite progress in understanding the statistics and spatial distribution of interparticle forces in recent decades, a systematic method for measuring forces in opaque, three-dimensional (3D), frictional, stiff granular media has yet to emerge. In this Letter, we present results from an experiment that combines 3D x-ray diffraction, x-ray tomography, and a numerical force inference technique to quantify interparticle forces and their heterogeneity in an assembly of quartz grains undergoing a one-dimensional compression cycle. Forces exhibit an exponential decay above the mean and partition into strong and weak networks. We find a surprising inverse relationship between macroscopic load and the heterogeneity of interparticle forces, despite the clear emergence of two force chains that span the system. PMID:27610890
Upscaling small heterogeneities for seismic wave propagation in 3D complex media
NASA Astrophysics Data System (ADS)
Cupillard, P.; Capdeville, Y.
2012-04-01
Seismic waves propagating in the Earth are affected by different sizes of heterogeneities. When modelling these waves (using numerical methods such as the SEM), taking into account heterogeneities that are much smaller than the minimum wavelength is a challenge because meshing small heterogeneities often requires important efforts and leads to high numerical costs. In this work, we present a technique which allows to upscale the small heterogeneities that can lie in an elastic medium. This technique yields a smooth effective medium and effective equations. We describe its implementation in the 3D case and we show relevant examples.
Towards Forward Modeling of 3D Heterogeneity in D" region
NASA Astrophysics Data System (ADS)
To, A.; Capdeville, Y.; Romanowicz, B.
2002-12-01
The presence of strong lateral heterogeneity in D" is now well documented. While tomographic modeling provides constraints on the large scale patterns, strong variations on shorter scales are best addressed by forward modeling. Appropriate tools are needed for forward modeling that will handle strong 3D heterogeneity, at relatively short periods and including diffracted waves. We use a coupled mode/SEM (Spectral Element Method) to compute synthetic seismograms in 3D models of the D" layer down to 1/12s. This coupled method (Capdeville, 2001) affords faster computations than SEM in cases where heterogeneity can be restricted to a specific layer. We compare them with observed waveforms for several events in the Western Pacific. Observed and synthetic travel time trends are very consistent, although in most cases the observed residuals are significantly larger. Waveform amplitudes are less consistent. In order to understand the origin of the amplitude difference, we test the effect of 3D heterogeneity on Sdiff phase. In particular, the results show opposite trends in the amplitude of Sdiff due to heterogeneity located near the CMB or well above it. This provides constraints on the location of the causative velocity heterogeneity. Because the forward modeling approach requires many iterations, the coupled mode/SEM approach is still computationally intensive. It is more efficient to use a less accurate traditional approach to first get closer to a final model, and only then use coupled mode/SEM to refine the model. Ray theory is the most expedient way to calculate travel times. However, it is an infinite frequency approximation and not appropriate to handle diffracting waves. We show that ray theory predicts larger travel time anomaly for Sdiff phase than the one obtained by coupled mode/SEM. Although it is based on a weak heterogeneity assumption, Non-linear Asymptotic Coupling Theory(NACT) (Li and Romanowicz, 1995) helps to overcome this difficulty. It can handle
Numerical Results of Earth's Core Accumulation 3-D Modelling
NASA Astrophysics Data System (ADS)
Khachay, Yurie; Anfilogov, Vsevolod
2013-04-01
For a long time as a most convenient had been the model of mega impact in which the early forming of the Earth's core and mantle had been the consequence of formed protoplanet collision with the body of Mercurial mass. But all dynamical models of the Earth's accumulation and the estimations after the Pb-Pb system, lead to the conclusion that the duration of the planet accumulation was about 1 milliard years. But isotopic results after the W-Hf system testify about a very early (5-10) million years, dividing of the geochemical reservoirs of the core and mantle. In [1,3] it is shown, that the account of energy dissipating by the decay of short living radioactive elements and first of all Al,it is sufficient for heating even small bodies with dimensions about (50-100) km up to the iron melting temperature and can be realized a principal new differentiation mechanism. The inner parts of the melted preplanets can join and they are mainly of iron content, but the cold silicate fragments return to the supply zone. Only after the increasing of the gravitational radius, the growing area of the future core can save also the silicate envelope fragments. All existing dynamical accumulation models are constructed by using a spherical-symmetrical model. Hence for understanding the further planet evolution it is significant to trace the origin and evolution of heterogeneities, which occur on the planet accumulation stage. In that paper we are modeling distributions of temperature, pressure, velocity of matter flowing in a block of 3D- spherical body with a growing radius. The boundary problem is solved by the finite-difference method for the system of equations, which include equations which describe the process of accumulation, the Safronov equation, the equation of impulse balance, equation Navier-Stocks, equation for above litho static pressure and heat conductivity in velocity-pressure variables using the Businesque approach. The numerical algorithm of the problem solution in
3D Game Content Distributed Adaptation in Heterogeneous Environments
NASA Astrophysics Data System (ADS)
Morán, Francisco; Preda, Marius; Lafruit, Gauthier; Villegas, Paulo; Berretty, Robert-Paul
2007-12-01
Most current multiplayer 3D games can only be played on a single dedicated platform (a particular computer, console, or cell phone), requiring specifically designed content and communication over a predefined network. Below we show how, by using signal processing techniques such as multiresolution representation and scalable coding for all the components of a 3D graphics object (geometry, texture, and animation), we enable online dynamic content adaptation, and thus delivery of the same content over heterogeneous networks to terminals with very different profiles, and its rendering on them. We present quantitative results demonstrating how the best displayed quality versus computational complexity versus bandwidth tradeoffs have been achieved, given the distributed resources available over the end-to-end content delivery chain. Additionally, we use state-of-the-art, standardised content representation and compression formats (MPEG-4 AFX, JPEG 2000, XML), enabling deployment over existing infrastructure, while keeping hooks to well-established practices in the game industry.
3D modeling of carbonates petro-acoustic heterogeneities
NASA Astrophysics Data System (ADS)
Baden, Dawin; Guglielmi, Yves; Saracco, Ginette; Marié, Lionel; Viseur, Sophie
2015-04-01
Characterizing carbonate reservoirs heterogeneity is a challenging issue for Oil & Gas Industry, CO2 sequestration and all kinds of fluid manipulations in natural reservoirs, due to the significant impact of heterogeneities on fluid flow and storage within the reservoir. Although large scale (> meter) heterogeneities such as layers petrophysical contrasts are well addressed by computing facies-based models, low scale (< meter) heterogeneities are often poorly constrained because of the complexity in predicting their spatial arrangement. In this study, we conducted petro-acoustic measurements on cores of different size and diameter (Ø = 1", 1.5" and 5") in order to evaluate anisotropy or heterogeneity in carbonates at different laboratory scales. Different types of heterogeneities which generally occur in carbonate reservoir units (e.g. petrographic, diagenetic, and tectonic related) were sampled. Dry / wet samples were investigated with different ultrasonic apparatus and using different sensors allowing acoustic characterization through a bandwidth varying from 50 to 500 kHz. Comprehensive measurements realized on each samples allowed statistical analyses of petro-acoustic properties such as attenuation, shear and longitudinal wave velocity. The cores properties (geological and acoustic facies) were modeled in 3D using photogrammetry and GOCAD geo-modeler. This method successfully allowed detecting and imaging in three dimensions differential diagenesis effects characterized by the occurrence of decimeter-scale diagenetic horizons in samples assumed to be homogeneous and/or different diagenetic sequences between shells filling and the packing matrix. We then discuss how small interfaces such as cracks, stylolithes and laminations which are also imaged may have guided these differential effects, considering that understanding the processes may be taken as an analogue to actual fluid drainage complexity in deep carbonate reservoir.
NASA Astrophysics Data System (ADS)
Litvinova, T.; Petrova, A.
2009-04-01
The study of magnetic anomaly field structure of the Barents Sea water area along seismic and extended profiles intersecting known fields is carried out. Geomagnetic and density sections down to 40 km depth are constructed. This allowed the estimation of heterogeneities of the Barents Sea water area deep structure. The analysis of geomagnetic and density sections along extended profiles showed the confinedness of oil-and-gas bearing provinces to deep permeable zones characterized by reduced magnetic and density features. Based on the analysis of permeable zones, regional diagnostic features similar to those obtained earlier in oil-and-gas bearing provinces in other regions, for example, in Timan-Pechora, Volga-Urals and Siberian, as well as in the Northern and Norwegian seas water areas, are revealed. The analysis of magnetic and gravity fields over the region area allowed the delineation of weakened zones as intersection areas of weakly magnetic areals with reduced density. Within the Barents Sea water area, permeable areas with lenticular-laminated structure of the upper and lower Earth's crust containing weakly magnetic areals with reduced rock density within the depth range of 8-12 and 15-20 km are revealed. Such ratio of magnetic and density heterogeneities in the Earth's crust is characteristic for zones with proved oil-and-gas content in the European part of the Atlantic Ocean water area. North Kildin field on 1 AR profile is confined to a trough with thick weakly magnetic stratum discontinuously traced to a depth of 6-10 km. At a depth of approximately 15 km, a lens of weakly magnetic and porous formations is observed. Ludlov field in the North Barents trough is confined to a zone of weakly magnetic rocks with reduced density traced to a depth of 8-9 km. Deeper, at Н=15 km, a lenticular areal of weakly magnetic formations with reduced density is observed. The profile transecting the Stockman field shows that it is located in the central part of a permeable
EarthServer - 3D Visualization on the Web
NASA Astrophysics Data System (ADS)
Wagner, Sebastian; Herzig, Pasquale; Bockholt, Ulrich; Jung, Yvonne; Behr, Johannes
2013-04-01
EarthServer (www.earthserver.eu), funded by the European Commission under its Seventh Framework Program, is a project to enable the management, access and exploration of massive, multi-dimensional datasets using Open GeoSpatial Consortium (OGC) query and processing language standards like WCS 2.0 and WCPS. To this end, a server/client architecture designed to handle Petabyte/Exabyte volumes of multi-dimensional data is being developed and deployed. As an important part of the EarthServer project, six Lighthouse Applications, major scientific data exploitation initiatives, are being established to make cross-domain, Earth Sciences related data repositories available in an open and unified manner, as service endpoints based on solutions and infrastructure developed within the project. Clients technology developed and deployed in EarthServer ranges from mobile and web clients to immersive virtual reality systems, all designed to interact with a physically and logically distributed server infrastructure using exclusively OGC standards. In this contribution, we would like to present our work on a web-based 3D visualization and interaction client for Earth Sciences data using only technology found in standard web browsers without requiring the user to install plugins or addons. Additionally, we are able to run the earth data visualization client on a wide range of different platforms with very different soft- and hardware requirements such as smart phones (e.g. iOS, Android), different desktop systems etc. High-quality, hardware-accelerated visualization of 3D and 4D content in standard web browsers can be realized now and we believe it will become more and more common to use this fast, lightweight and ubiquitous platform to provide insights into big datasets without requiring the user to set up a specialized client first. With that in mind, we will also point out some of the limitations we encountered using current web technologies. Underlying the EarthServer web client
Equivalent Body Force Finite Elements Method and 3-D Earth Model Applied In 2004 Sumatra Earthquake
NASA Astrophysics Data System (ADS)
Qu, W.; Cheng, H.; Shi, Y.
2015-12-01
The 26 December 2004 Sumatra-Andaman earthquake with moment magnitude (Mw) of 9.1 to 9.3 is the first great earthquake recorded by digital broadband, high-dynamic-range seismometers and global positioning system (GPS) equipment, which recorded many high-quality geophysical data sets. The spherical curvature is not negligible in far field especially for large event and the real Earth is laterally inhomogeneity and the analytical results still are difficult to explain the geodetic measurements. We use equivalent body force finite elements method Zhang et al. (2015) and mesh the whole earth, to compute global co-seismic displacements using four fault slip models of the 2004 Sumatra earthquake provided by different authors. Comparisons of calculated co-seismic displacements and GPS show that the confidences are well in near field for four models, and the confidences are according to different models. In the whole four models, the Chlieh model (Chlieh et al., 2007) is the best as this slip model not only accord well with near field data but also far field data. And then we use the best slip model, Chlieh model to explore influence of three dimensional lateral earth structure on both layered spherically symmetric (PREM) and real 3-D heterogeneous earth model (Crust 1.0 model and GyPSuM). Results show that the effects of 3-D heterogeneous earth model are not negligible and decrease concomitantly with increasing distance from the epicenter. The relative effects of 3-D crust model are 23% and 40% for horizontal and vertical displacements, respectively. The effects of the 3-D mantle model are much smaller than that of 3-D crust model but with wider impacting area.
CO2 leakage risk in 3D heterogeneous formations
NASA Astrophysics Data System (ADS)
Hou, Z.; Murray, C. J.; Rockhold, M. L.
2011-12-01
In this study we use a stochastic sensitivity analysis framework to evaluate the impact of 3D spatial heterogeneity in permeability on CO2 leakage risk. The leakage is defined as the total mass of CO2 moving into the overburden through the caprock-overburden interface, in both gaseous and liquid (dissolved) phases. The entropy-based framework has the ability to quantify the uncertainty associated with the input parameters/factors in the form of prior pdfs (probability density functions). Effective sampling of the prior pdfs enables us to explore the parameter space and systematically evaluate the individual and combined effects of the factors/parameters of interest on CO2 leakage risk. The parameters that are considered in the study include: mean, variance, and horizontal to vertical spatial anisotropy ratio for caprock permeability, and those same parameters for reservoir permeability. Given the sampled spatial variogram parameters, multiple realizations of permeability fields were generated using GSLIB subroutines. For each permeability field, a numerical simulator STOMP (water-salt-CO2-energy operational mode) is used to simulate the CO2 migration within the reservoir and caprock up to 50 years after injection. Due to intensive computational demand, a scalable version simulator, eSTOMP, is run on the Jaguar supercomputer. We then perform statistical analyses and summarize the relationships between the parameters of interest (mean/variance/anisotropy ratio of caprock/reservoir permeability) and CO2 leakage ratio. We will also present the effects of those parameters on CO2 plume radius and reservoir injectivity.
3D Orbit Visualization for Earth-Observing Missions
NASA Technical Reports Server (NTRS)
Jacob, Joseph C.; Plesea, Lucian; Chafin, Brian G.; Weiss, Barry H.
2011-01-01
This software visualizes orbit paths for the Orbiting Carbon Observatory (OCO), but was designed to be general and applicable to any Earth-observing mission. The software uses the Google Earth user interface to provide a visual mechanism to explore spacecraft orbit paths, ground footprint locations, and local cloud cover conditions. In addition, a drill-down capability allows for users to point and click on a particular observation frame to pop up ancillary information such as data product filenames and directory paths, latitude, longitude, time stamp, column-average dry air mole fraction of carbon dioxide, and solar zenith angle. This software can be integrated with the ground data system for any Earth-observing mission to automatically generate daily orbit path data products in Google Earth KML format. These KML data products can be directly loaded into the Google Earth application for interactive 3D visualization of the orbit paths for each mission day. Each time the application runs, the daily orbit paths are encapsulated in a KML file for each mission day since the last time the application ran. Alternatively, the daily KML for a specified mission day may be generated. The application automatically extracts the spacecraft position and ground footprint geometry as a function of time from a daily Level 1B data product created and archived by the mission s ground data system software. In addition, ancillary data, such as the column-averaged dry air mole fraction of carbon dioxide and solar zenith angle, are automatically extracted from a Level 2 mission data product. Zoom, pan, and rotate capability are provided through the standard Google Earth interface. Cloud cover is indicated with an image layer from the MODIS (Moderate Resolution Imaging Spectroradiometer) aboard the Aqua satellite, which is automatically retrieved from JPL s OnEarth Web service.
Real time 3D and heterogeneous data fusion
Little, C.Q.; Small, D.E.
1998-03-01
This project visualizes characterization data in a 3D setting, in real time. Real time in this sense means collecting the data and presenting it before it delays the user, and processing faster than the acquisition systems so no bottlenecks occur. The goals have been to build a volumetric viewer to display 3D data, demonstrate projecting other data, such as images, onto the 3D data, and display both the 3D and projected images as fast as the data became available. The authors have examined several ways to display 3D surface data. The most effective was generating polygonal surface meshes. They have created surface maps form a continuous stream of 3D range data, fused image data onto the geometry, and displayed the data with a standard 3D rendering package. In parallel with this, they have developed a method to project real-time images onto the surface created. A key component is mapping the data on the correct surfaces, which requires a-priori positional information along with accurate calibration of the camera and lens system.
2D and 3D heterogeneous photonic integrated circuits
NASA Astrophysics Data System (ADS)
Yoo, S. J. Ben
2014-03-01
Exponential increases in the amount of data that need to be sensed, communicated, and processed are continuing to drive the complexity of our computing, networking, and sensing systems. High degrees of integration is essential in scalable, practical, and cost-effective microsystems. In electronics, high-density 2D integration has naturally evolved towards 3D integration by stacking of memory and processor chips with through-silicon-vias. In photonics, too, we anticipate highdegrees of 3D integration of photonic components to become a prevailing method in realizing future microsystems for information and communication technologies. However, compared to electronics, photonic 3D integration face a number of challenges. This paper will review two methods of 3D photonic integration --- fs laser inscription and layer stacking, and discuss applications and future prospects.
Glacial isostatic adjustment on 3-D Earth models: a finite-volume formulation
NASA Astrophysics Data System (ADS)
Latychev, Konstantin; Mitrovica, Jerry X.; Tromp, Jeroen; Tamisiea, Mark E.; Komatitsch, Dimitri; Christara, Christina C.
2005-05-01
We describe and present results from a finite-volume (FV) parallel computer code for forward modelling the Maxwell viscoelastic response of a 3-D, self-gravitating, elastically compressible Earth to an arbitrary surface load. We implement a conservative, control volume discretization of the governing equations using a tetrahedral grid in Cartesian geometry and a low-order, linear interpolation. The basic starting grid honours all major radial discontinuities in the Preliminary Reference Earth Model (PREM), and the models are permitted arbitrary spatial variations in viscosity and elastic parameters. These variations may be either continuous or discontinuous at a set of grid nodes forming a 3-D surface within the (regional or global) modelling domain. In the second part of the paper, we adopt the FV methodology and a spherically symmetric Earth model to generate a suite of predictions sampling a broad class of glacial isostatic adjustment (GIA) data types (3-D crustal motions, long-wavelength gravity anomalies). These calculations, based on either a simple disc load history or a global Late Pleistocene ice load reconstruction (ICE-3G), are benchmarked against predictions generated using the traditional normal-mode approach to GIA. The detailed comparison provides a guide for future analyses (e.g. what grid resolution is required to obtain a specific accuracy?) and it indicates that discrepancies in predictions of 3-D crustal velocities less than 0.1 mm yr-1 are generally obtainable for global grids with ~3 × 106 nodes; however, grids of higher resolution are required to predict large-amplitude (>1 cm yr-1) radial velocities in zones of peak post-glacial uplift (e.g. James bay) to the same level of absolute accuracy. We conclude the paper with a first application of the new formulation to a 3-D problem. Specifically, we consider the impact of mantle viscosity heterogeneity on predictions of present-day 3-D crustal motions in North America. In these tests, the
Joint earthquake source inversions using seismo-geodesy and 3-D earth models
NASA Astrophysics Data System (ADS)
Weston, J.; Ferreira, A. M. G.; Funning, G. J.
2014-08-01
A joint earthquake source inversion technique is presented that uses InSAR and long-period teleseismic data, and, for the first time, takes 3-D Earth structure into account when modelling seismic surface and body waves. Ten average source parameters (Moment, latitude, longitude, depth, strike, dip, rake, length, width and slip) are estimated; hence, the technique is potentially useful for rapid source inversions of moderate magnitude earthquakes using multiple data sets. Unwrapped interferograms and long-period seismic data are jointly inverted for the location, fault geometry and seismic moment, using a hybrid downhill Powell-Monte Carlo algorithm. While the InSAR data are modelled assuming a rectangular dislocation in a homogeneous half-space, seismic data are modelled using the spectral element method for a 3-D earth model. The effect of noise and lateral heterogeneity on the inversions is investigated by carrying out realistic synthetic tests for various earthquakes with different faulting mechanisms and magnitude (Mw 6.0-6.6). Synthetic tests highlight the improvement in the constraint of fault geometry (strike, dip and rake) and moment when InSAR and seismic data are combined. Tests comparing the effect of using a 1-D or 3-D earth model show that long-period surface waves are more sensitive than long-period body waves to the change in earth model. Incorrect source parameters, particularly incorrect fault dip angles, can compensate for systematic errors in the assumed Earth structure, leading to an acceptable data fit despite large discrepancies in source parameters. Three real earthquakes are also investigated: Eureka Valley, California (1993 May 17, Mw 6.0), Aiquile, Bolivia (1998 February 22, Mw 6.6) and Zarand, Iran (2005 May 22, Mw 6.5). These events are located in different tectonic environments and show large discrepancies between InSAR and seismically determined source models. Despite the 40-50 km discrepancies in location between previous geodetic and
Vorticity and upscaled dispersion in 3D heterogeneous porous media
NASA Astrophysics Data System (ADS)
Di Dato, Mariaines; Chiogna, Gabriele; de Barros, Felipe; Bellin, Alberto; Fiori, Aldo
2015-04-01
Modeling flow in porous media is relevant for many environmental, energy and industrial applications. From an environmental perspective, the relevance of porous media flow becomes evident in subsurface hydrology. In general, flow in natural porous media is creeping, yet the large variability in the hydraulic conductivity values encountered in natural aquifers leads to highly heterogeneous flow fields. This natural variability in the conductivity field will affect both dilution rates of chemical species and reactive mixing. A physical consequence of this heterogeneity is also the presence of a various localized kinematical features such as straining, shearing and vorticity in aquifers, which will influence the shape of solute clouds and its fate and transport. This work aims in fundamentally characterizing the vorticity field in spatially heterogeneous flow fields as a function of their statistical properties in order to analyze the impact on transport processes. In our study, three-dimensional porous formations are constructed with an ensemble of N independent, non-overlapping spheroidal inclusions submerged into an homogeneous matrix, of conductivity K0. The inclusions are randomly located in a domain of volume W and are fully characterized by the geometry of spheroid (oblate or prolate), their conductivity K (random and drawn from a given probability density function fκ), the centroid location ¯x, the axes ratio e, the orientation of the rotational axis (α1,α2) and the volume w. Under the assumption of diluted medium, the flow problem is solved analitically by means of only two parameters: the conductivity contrast κ = K/K0 and the volume fraction n = Nw/W . Through the variation of these parameters of the problem, it is possible to approximate the structure of natural heterogeneous porous media. Using a random distribution of the orientation of the inclusions, we create media defined by the same global anisotropy f = Iz/Ix but different micro
Jacob's Interpretation Method Revisited: Accounting for 3-D Spatial Heterogeneity
NASA Astrophysics Data System (ADS)
Sanchez-Vila, X.; Riva, M.; Guadagnini, A.; Carrera, J.
2005-12-01
Traditional approaches to hydraulic test interpretation provide typically individual aquifer parameters, such as hydraulic conductivity (K) and storativity (S) values. The values obtained somehow incorporate some averaging values of aquifer heterogeneity, while the averaging functions are a direct consequence of the method of analysis employed. In recent years most work, casted in a stochastic framework, focused on the relationship between pumping rate and ensemble mean or variance of drawdown, thus having to pre-specify the parameters characterizing the underlying random spatial function. On the contrary, we contend that additional highly relevant information about heterogeneity can be obtained by looking to the spatial distribution of drawdown in individual realizations of the heterogeneous K field, without the need for invoking ergodic arguments. We present an analysis of the spatial distribution of time-dependent drawdown in a tridimensional aquifer produced by constant rate pumping in a fully penetrating well. The aquifer is considered of infinite extension in the x, y directions, and we assume no-flow boundaries in the aquifer top and bottom. The observation point is a fully penetrating piezometer. We consider an unknown spatial distribution of K(x,y,z), and using a perturbation expansion up to second order, we look at the late-time behavior of drawdown at any given observation vertical line. We conclude that: (1) at any given observation line the late-time behavior of drawdown would display a straight line in a drawdown versus log time plot, thus allowing the use of Jacob's method for test interpretation; (2) the slope of the straight line is the same for each observation line, thus providing a global average of K(x,y,z) through the aquifer; (3) the intercept point of the line in the same plot depends on location and is related to connectivity issues between the pumping and observation locations; (4) the intercept value is a weighted function of the local
Reactive Transport from Path3D: A Stream Tube Approach for Heterogeneous Aquifers
NASA Astrophysics Data System (ADS)
LI, L.
2001-05-01
Path3D (Zheng, 1991) is a popular computer program run in series with MODFLOW. Remediation engineers and hydrogeologist use it to track contaminant paths and to estimate solute travel time at heterogeneous sites. In order to predict fate and transport of multiple species at heterogeneous sites, numerical modeling packages, such as MT3D (Zheng,1990) or RT3D (Clement 1997) are often employed. These packages also are based on pre-processing with MODFLOW. However, for complex cases with aquifer heterogeneity, MT3D and RT3D often require very long computer run times. This paper addresses a new, stream-tube, approach that is both highly efficient and accurate to predict multi-species reactive transport at heterogeneous sites with steady flow. Our application of the stream tube approach is different from other stream tube approaches that apply the advection-dispersion-reaction equation in each stream tube (such as Ginn, 2000, Yabusaki, 1998, Charbeneau, 2000). In this work, the authors make use of properties of a linear system,working with decoupled reaction and sorption processes and mixing processes described by residence time distributions (RTDs). RTDs are abstracted from Path3D particle-tracking results and additional temporal and spatial dispersion (not caused by aquifer heterogeneity) is ignored. Reactions, including first order reactions and linear, reversible sorption, are applied through analytical transfer functions (called kinetic response functions). Convolution can then be applied to determine contaminant concentrations at monitoring points, using the RTDs determined from Path3D, kinetic transfer functions (expressed analytically), and expected trends of the source concentration. We are currently testing the approach and noting significant computational advantages for problems in three-dimensions, with first order reaction pathways and different retardation factors. We will demonstrate the method with several examples and compare the performance with MT3D and
The Earth's Seasons in 3-D--Part 1.
ERIC Educational Resources Information Center
Eckroth, Charles A.
1993-01-01
This article describes a teaching aid made from four colored foam balls mounted on a stiff wire circle used to teach about the changing seasons and earth temperature fluctuations. The spheres represent the Earth at the solstice and equinox positions. (MVL)
Heterogeneously Assembled Metamaterials and Metadevices via 3D Modular Transfer Printing
Lee, Seungwoo; Kang, Byungsoo; Keum, Hohyun; Ahmed, Numair; Rogers, John A.; Ferreira, Placid M.; Kim, Seok; Min, Bumki
2016-01-01
Metamaterials have made the exotic control of the flow of electromagnetic waves possible, which is difficult to achieve with natural materials. In recent years, the emergence of functional metadevices has shown immense potential for the practical realization of highly efficient photonic devices. However, complex and heterogeneous architectures that enable diverse functionalities of metamaterials and metadevices have been challenging to realize because of the limited manufacturing capabilities of conventional fabrication methods. Here, we show that three-dimensional (3D) modular transfer printing can be used to construct diverse metamaterials in complex 3D architectures on universal substrates, which is attractive for achieving on-demand photonic properties. Few repetitive processing steps and rapid constructions are additional advantages of 3D modular transfer printing. Thus, this method provides a fascinating route to generate flexible and stretchable 2D/3D metamaterials and metadevices with heterogeneous material components, complex device architectures, and diverse functionalities. PMID:27283594
Heterogeneously Assembled Metamaterials and Metadevices via 3D Modular Transfer Printing.
Lee, Seungwoo; Kang, Byungsoo; Keum, Hohyun; Ahmed, Numair; Rogers, John A; Ferreira, Placid M; Kim, Seok; Min, Bumki
2016-01-01
Metamaterials have made the exotic control of the flow of electromagnetic waves possible, which is difficult to achieve with natural materials. In recent years, the emergence of functional metadevices has shown immense potential for the practical realization of highly efficient photonic devices. However, complex and heterogeneous architectures that enable diverse functionalities of metamaterials and metadevices have been challenging to realize because of the limited manufacturing capabilities of conventional fabrication methods. Here, we show that three-dimensional (3D) modular transfer printing can be used to construct diverse metamaterials in complex 3D architectures on universal substrates, which is attractive for achieving on-demand photonic properties. Few repetitive processing steps and rapid constructions are additional advantages of 3D modular transfer printing. Thus, this method provides a fascinating route to generate flexible and stretchable 2D/3D metamaterials and metadevices with heterogeneous material components, complex device architectures, and diverse functionalities. PMID:27283594
Heterogeneously Assembled Metamaterials and Metadevices via 3D Modular Transfer Printing
NASA Astrophysics Data System (ADS)
Lee, Seungwoo; Kang, Byungsoo; Keum, Hohyun; Ahmed, Numair; Rogers, John A.; Ferreira, Placid M.; Kim, Seok; Min, Bumki
2016-06-01
Metamaterials have made the exotic control of the flow of electromagnetic waves possible, which is difficult to achieve with natural materials. In recent years, the emergence of functional metadevices has shown immense potential for the practical realization of highly efficient photonic devices. However, complex and heterogeneous architectures that enable diverse functionalities of metamaterials and metadevices have been challenging to realize because of the limited manufacturing capabilities of conventional fabrication methods. Here, we show that three-dimensional (3D) modular transfer printing can be used to construct diverse metamaterials in complex 3D architectures on universal substrates, which is attractive for achieving on-demand photonic properties. Few repetitive processing steps and rapid constructions are additional advantages of 3D modular transfer printing. Thus, this method provides a fascinating route to generate flexible and stretchable 2D/3D metamaterials and metadevices with heterogeneous material components, complex device architectures, and diverse functionalities.
The OpenEarth Framework (OEF) for the 3D Visualization of Integrated Earth Science Data
NASA Astrophysics Data System (ADS)
Nadeau, David; Moreland, John; Baru, Chaitan; Crosby, Chris
2010-05-01
Data integration is increasingly important as we strive to combine data from disparate sources and assemble better models of the complex processes operating at the Earth's surface and within its interior. These data are often large, multi-dimensional, and subject to differing conventions for data structures, file formats, coordinate spaces, and units of measure. When visualized, these data require differing, and sometimes conflicting, conventions for visual representations, dimensionality, symbology, and interaction. All of this makes the visualization of integrated Earth science data particularly difficult. The OpenEarth Framework (OEF) is an open-source data integration and visualization suite of applications and libraries being developed by the GEON project at the University of California, San Diego, USA. Funded by the NSF, the project is leveraging virtual globe technology from NASA's WorldWind to create interactive 3D visualization tools that combine and layer data from a wide variety of sources to create a holistic view of features at, above, and beneath the Earth's surface. The OEF architecture is open, cross-platform, modular, and based upon Java. The OEF's modular approach to software architecture yields an array of mix-and-match software components for assembling custom applications. Available modules support file format handling, web service communications, data management, user interaction, and 3D visualization. File parsers handle a variety of formal and de facto standard file formats used in the field. Each one imports data into a general-purpose common data model supporting multidimensional regular and irregular grids, topography, feature geometry, and more. Data within these data models may be manipulated, combined, reprojected, and visualized. The OEF's visualization features support a variety of conventional and new visualization techniques for looking at topography, tomography, point clouds, imagery, maps, and feature geometry. 3D data such as
Rapid Assembly of Heterogeneous 3D Cell Microenvironments in a Microgel Array.
Li, Yiwei; Chen, Pu; Wang, Yachao; Yan, Shuangqian; Feng, Xiaojun; Du, Wei; Koehler, Stephan A; Demirci, Utkan; Liu, Bi-Feng
2016-05-01
Heterogeneous 3D cell microenvironment arrays are rapidly assembled by combining surface-wettability-guided assembly and microdroplet-array-based operations. This approach enables precise control over individual shapes, sizes, chemical concentrations, cell density, and 3D spatial distribution of multiple components. This technique provides a cost-effective solution to meet the increasing demand of stem cell research, tissue engineering, and drug screening. PMID:26991071
A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation
Riyanti, C.D. . E-mail: C.D.Riyanti@tudelft.nl; Kononov, A.; Erlangga, Y.A.; Vuik, C.; Oosterlee, C.W.; Plessix, R.-E.; Mulder, W.A.
2007-05-20
We investigate the parallel performance of an iterative solver for 3D heterogeneous Helmholtz problems related to applications in seismic wave propagation. For large 3D problems, the computation is no longer feasible on a single processor, and the memory requirements increase rapidly. Therefore, parallelization of the solver is needed. We employ a complex shifted-Laplace preconditioner combined with the Bi-CGSTAB iterative method and use a multigrid method to approximate the inverse of the resulting preconditioning operator. A 3D multigrid method with 2D semi-coarsening is employed. We show numerical results for large problems arising in geophysical applications.
3D Printing of Advanced Biocomposites on Earth and Beyond
NASA Technical Reports Server (NTRS)
Rothschild, Lynn J.; Gentry, Diana; Micks, Ashley
2015-01-01
Human exploration off planet is severely limited by the cost of launching materials into space and re-supply. Thus materials brought from earth must be light, stable and reliable at destination. Using traditional approaches a lunar or Mars base would require either transporting a hefty store of metals or heavy manufacturing equipment and construction materials for in situ extraction; both would severely limit any other mission objectives. Long-term human space presence requires periodic replenishment, adding a massive cost overhead. Even robotic missions often sacrifice science goals for heavy radiation and thermal protection. Biology has the potential to solve these problems because it can replicate and repair itself, and do a wide variety of chemical reactions including making food, fuel and materials. Synthetic biology enhances and expands life's evolved repertoire. Using organisms as feedstock, additive manufacturing could make possible the dream of producing bespoke tools, food, smart fabrics and even replacement organs on demand. Image what new products can be enabled by such a technology, on earth or beyond!
3D Printing of Advanced Biocomposites on Earth and Beyond
NASA Technical Reports Server (NTRS)
Rothschild, Lynn J.; Gentry, Diana M.; Micks, Ashley
2015-01-01
Human exploration off planet is severely limited by the cost of launching materials into space and re-supply. Thus materials brought from earth must be light, stable and reliable at destination. Using traditional approaches a lunar or Mars base would require either transporting a hefty store of metals or heavy manufacturing equipment and construction materials for in situ extraction; both would severely limit any other mission objectives. Long-term human space presence requires periodic replenishment, adding a massive cost overhead. Even robotic missions often sacrifice science goals for heavy radiation and thermal protection. Biology has the potential to solve these problems because it can replicate and repair itself, and do a wide variety of chemical reactions including making food, fuel and materials. Synthetic biology can greatly enhance and expand life's evolved repertoire. Using natural and synthetically altered organisms as the feedstock for additive manufacturing could one day make possible the dream of producing bespoke tools, food, smart fabrics and even replacement organs on demand. To this end our lab has produced a proof-of-concept bioprinter with nearly one-cell resolution. Genetically engineering yeast cells to secrete bioproducts subsequent to printing allows the potential to make biomaterials with a fine microstructure. Imagine a production system that, at a few micron scale resolution, can add mollusk shell for compressive strength per unit mass, spider silk or collagen for tensile strength per unit mass, and potentially biologically-deposited wires. Now imagine what new products can be enabled by such a technology, on earth or beyond
3D numerical analysis of crack propagation of heterogeneous notched rock under uniaxial tension
NASA Astrophysics Data System (ADS)
Wang, S. Y.; Sloan, S. W.; Sheng, D. C.; Tang, C. A.
2016-05-01
Macroscopic notches play an important role in evaluating the fracture process zone (FPZ) and the strengths of a heterogeneous rock mass. Crack initiation, propagation and coalescence for unnotched, single-notched and double-notched rock specimens are numerically simulated in a 3-D numerical model (RFPA3D). A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. For the unnotched case, special attention is given to the complete stress-strain curve and the corresponding AE events for the failure process of rock specimen. By comparing with published experimental results, the simulation results from RFPA3D are found to be satisfactory. For the single-notched case, the effect of the length and the depth of the single notch and the thickness of the specimen on the failure mode and peak stress are evaluated. The 3D FPZ is very different from that in two dimensions. For the double-notched case, the effects of the separation distance and overlap distance of the double notches, as well as influence of the homogeneity index (m) are also investigated. As the overlap distance increases, the direction of the principal tensile stress at each notch-end changes from a perpendicular direction (tensile stress field) to a nearly parallel direction (compressive stress field), which affects the evolution of the cracks from the two notches.
Scales of mantle heterogeneity emerging from 3-D models of advective stretching
NASA Astrophysics Data System (ADS)
Kellogg, L. H.; Conjeepuram, N.
2009-12-01
Heterogeneities are continually introduced into the mantle by subduction, and then are homogenized by stretching, folding, and finally diffusion. The stretching and folding components control the timescale of mixing in the mantle. Mixing has been studied in 2-D and to a lesser extent in 3-D models, often by using statistical analysis of separation of passive tracers. It has been proposed that mixing in 3-D time dependent convection may differ substantially from mixing in 2-D due to the different structure of the flow. To investigate the processes that determine the scales of heterogeneity in the mantle, we use a complementary method, computing the stretching experienced by passive, infinitesimal, ellipsoidal strain markers in 3-D models of mantle convection. This approach has an advantage over more commonly used methods of calculating separation of particles, because we obtain information about deformation (a mechanism to develop different scales of heterogeneity in the mantle) and about orientation of strain ellipsoids (which can result in fabrics that may lead to anisotropy). We investigate both kinematic and dynamic flows. In plate-driven kinematic flows, the toroidal component of the velocity field emerges as an important factor in mixing. Increasing the toroidal energy in the flow increases the complexity of the stretching patterns that develop and persist through time and homogenizes the stretching distribution. By computing the frequency size distribution of the strain ellipsoids we find that a marble cake upper mantle is a natural consequence of plate-driven flow. We also apply this method to evaluate the role of viscosity contrast in development of heterogeneity convection at different Rayleigh numbers. These models yield complex patterns in which tracers can separate or remain isolated, again leading to a marble-cake upper mantle. We use an innovative method of visualizing the distribution of stretching in 3-D to illustrate these results.
3D exploitation system for operational applications of Earth observation data
NASA Astrophysics Data System (ADS)
Bonavenia, Roberto; Colaiacomo, Lucio; Dooley, Colin; Menu, Eric; Palumbo, Gianluca; Prisco, Giulio; Valero, Juan L.; Villemaud, Olivier
1999-07-01
The operational applicability of Earth Observation data to facilitate decision making is demonstrated with examples taken from the work of the Western European Union Satellite Center in Madrid. Analysis and reporting techniques based on 3D representations of the surface of the Earth and Virtual Reality are described.
3D bioprinting of heterogeneous aortic valve conduits with alginate/gelatin hydrogels.
Duan, Bin; Hockaday, Laura A; Kang, Kevin H; Butcher, Jonathan T
2013-05-01
Heart valve disease is a serious and growing public health problem for which prosthetic replacement is most commonly indicated. Current prosthetic devices are inadequate for younger adults and growing children. Tissue engineered living aortic valve conduits have potential for remodeling, regeneration, and growth, but fabricating natural anatomical complexity with cellular heterogeneity remain challenging. In the current study, we implement 3D bioprinting to fabricate living alginate/gelatin hydrogel valve conduits with anatomical architecture and direct incorporation of dual cell types in a regionally constrained manner. Encapsulated aortic root sinus smooth muscle cells (SMC) and aortic valve leaflet interstitial cells (VIC) were viable within alginate/gelatin hydrogel discs over 7 days in culture. Acellular 3D printed hydrogels exhibited reduced modulus, ultimate strength, and peak strain reducing slightly over 7-day culture, while the tensile biomechanics of cell-laden hydrogels were maintained. Aortic valve conduits were successfully bioprinted with direct encapsulation of SMC in the valve root and VIC in the leaflets. Both cell types were viable (81.4 ± 3.4% for SMC and 83.2 ± 4.0% for VIC) within 3D printed tissues. Encapsulated SMC expressed elevated alpha-smooth muscle actin, while VIC expressed elevated vimentin. These results demonstrate that anatomically complex, heterogeneously encapsulated aortic valve hydrogel conduits can be fabricated with 3D bioprinting. PMID:23015540
A highly heterogeneous 3D PWR core benchmark: deterministic and Monte Carlo method comparison
NASA Astrophysics Data System (ADS)
Jaboulay, J.-C.; Damian, F.; Douce, S.; Lopez, F.; Guenaut, C.; Aggery, A.; Poinot-Salanon, C.
2014-06-01
Physical analyses of the LWR potential performances with regards to the fuel utilization require an important part of the work dedicated to the validation of the deterministic models used for theses analyses. Advances in both codes and computer technology give the opportunity to perform the validation of these models on complex 3D core configurations closed to the physical situations encountered (both steady-state and transient configurations). In this paper, we used the Monte Carlo Transport code TRIPOLI-4®; to describe a whole 3D large-scale and highly-heterogeneous LWR core. The aim of this study is to validate the deterministic CRONOS2 code to Monte Carlo code TRIPOLI-4®; in a relevant PWR core configuration. As a consequence, a 3D pin by pin model with a consistent number of volumes (4.3 millions) and media (around 23,000) is established to precisely characterize the core at equilibrium cycle, namely using a refined burn-up and moderator density maps. The configuration selected for this analysis is a very heterogeneous PWR high conversion core with fissile (MOX fuel) and fertile zones (depleted uranium). Furthermore, a tight pitch lattice is selcted (to increase conversion of 238U in 239Pu) that leads to harder neutron spectrum compared to standard PWR assembly. In these conditions two main subjects will be discussed: the Monte Carlo variance calculation and the assessment of the diffusion operator with two energy groups for the core calculation.
NASA Astrophysics Data System (ADS)
Wang, Minjie; Shen, Shuwei; Yang, Jie; Dong, Erbao; Xu, Ronald
2014-03-01
The performance of biomedical optical imaging devices heavily relies on appropriate calibration. However, many of existing calibration phantoms for biomedical optical devices are based on homogenous materials without considering the multi-layer heterogeneous structures observed in biological tissue. Using such a phantom for optical calibration may result in measurement bias. To overcome this problem, we propose a 3D printing method for freeform fabrication of tissue simulating phantoms with multilayer heterogeneous structure. The phantom simulates not only the morphologic characteristics of biological tissue but also absorption and scattering properties. The printing system is based on a 3D motion platform with coordinated control of the DC motors. A special jet nozzle is designed to mix base, scattering, and absorption materials at different ratios. 3D tissue structures are fabricated through layer-by-layer printing with selective deposition of phantom materials of different ingredients. Different mixed ratios of base, scattering and absorption materials have been tested in order to optimize the printing outcome. A spectrometer and a tissue spectrophotometer are used for characterizing phantom absorption and scattering properties. The goal of this project is to fabricate skin tissue simulating phantoms as a traceable standard for the calibration of biomedical optical spectral devices.
A hybrid-3D hillslope hydrological model for use in Earth system models
NASA Astrophysics Data System (ADS)
Hazenberg, P.; Fang, Y.; Broxton, P.; Gochis, D.; Niu, G.-Y.; Pelletier, J. D.; Troch, P. A.; Zeng, X.
2015-10-01
Hillslope-scale rainfall-runoff processes leading to a fast catchment response are not explicitly included in land surface models (LSMs) for use in earth system models (ESMs) due to computational constraints. This study presents a hybrid-3D hillslope hydrological model (h3D) that couples a 1-D vertical soil column model with a lateral pseudo-2D saturated zone and overland flow model for use in ESMs. By representing vertical and lateral responses separately at different spatial resolutions, h3D is computationally efficient. The h3D model was first tested for three different hillslope planforms (uniform, convergent and divergent). We then compared h3D (with single and multiple soil columns) with a complex physically based 3-D model and a simple 1-D soil moisture model coupled with an unconfined aquifer (as typically used in LSMs). It is found that simulations obtained by the simple 1-D model vary considerably from the complex 3-D model and are not able to represent hillslope-scale variations in the lateral flow response. In contrast, the single soil column h3D model shows a much better performance and saves computational time by 2-3 orders of magnitude compared with the complex 3-D model. When multiple vertical soil columns are implemented, the resulting hydrological responses (soil moisture, water table depth, and base flow along the hillslope) from h3D are nearly identical to those predicted by the complex 3-D model, but still saves computational time. As such, the computational efficiency of the h3D model provides a valuable and promising approach to incorporating hillslope-scale hydrological processes into continental and global-scale ESMs.
Zhu, Jiahua; Penfold, Scott N
2016-06-01
Correct modelling of the interaction parameters of patient tissues is of vital importance in proton therapy treatment planning because of the large dose gradients associated with the Bragg peak. Different 3D imaging techniques yield different information regarding these interaction parameters. Given the rapidly expanding interest in proton therapy, this review is written to make readers aware of the current challenges in accounting for tissue heterogeneities and the imaging systems that are proposed to tackle these challenges. A summary of the interaction parameters of interest in proton therapy and the current and developmental 3D imaging techniques used in proton therapy treatment planning is given. The different methods to translate the imaging data to the interaction parameters of interest are reviewed and a summary of the implementations in several commercial treatment planning systems is presented. PMID:27115163
Characterizing heterogeneity among virus particles by stochastic 3D signal reconstruction
NASA Astrophysics Data System (ADS)
Xu, Nan; Gong, Yunye; Wang, Qiu; Zheng, Yili; Doerschuk, Peter C.
2015-09-01
In single-particle cryo electron microscopy, many electron microscope images each of a single instance of a biological particle such as a virus or a ribosome are measured and the 3-D electron scattering intensity of the particle is reconstructed by computation. Because each instance of the particle is imaged separately, it should be possible to characterize the heterogeneity of the different instances of the particle as well as a nominal reconstruction of the particle. In this paper, such an algorithm is described and demonstrated on the bacteriophage Hong Kong 97. The algorithm is a statistical maximum likelihood estimator computed by an expectation maximization algorithm implemented in Matlab software.
NASA Astrophysics Data System (ADS)
Duclaux, Guillaume; Huismans, Ritske S.; May, Dave
2015-04-01
How and why do continents break? More than two decades of analogue and 2D plane-strain numerical experiments have shown that despite the origin of the forces driving extension, the geometry of continental rifts falls into three categories - or modes: narrow rift, wide rift, or core complex. The mode of extension itself is strongly influenced by the rheology (and rheological behaviour) of the modelled layered system. In every model, an initial thermal or mechanical heterogeneity, such as a weak seed or a notch, is imposed to help localise the deformation and avoid uniform stretching of the lithosphere by pure shear. While it is widely accepted that structural inheritance is a key parameter for controlling rift localisation - as implied by the Wilson Cycle - modelling the effect of lithospheric heterogeneities on the long-term tectonic evolution of an extending plate in full 3D remains challenging. Recent progress in finite-element methods applied to computational tectonics along with the improved accessibility to high performance computers, now enable to switch from plane strain thermo-mechanical experiments to full 3D high-resolution experiments. Here we investigate the role of mechanical heterogeneities on rift opening, linkage and propagation during extension of a layered lithospheric systems with pTatin3d, a geodynamics modeling package utilising the material-point-method for tracking material composition, combined with a multigrid finite-element method to solve heterogeneous, incompressible visco-plastic Stokes problems. The initial model setup consists in a box of 1200 km horizontally by 250 km deep. It includes a 35 km layer of continental crust, underlaid by 85 km of sub-continental lithospheric mantle, and an asthenospheric mantle. Crust and mantle have visco-plastic rheologies with a pressure dependent yielding, which includes strain weakening, and a temperature, stress, strain-rate-dependent viscosity based on wet quartzite rheology for the crust, and wet
Microdomain heterogeneity in 3D affects the mechanics of neonatal cardiac myocyte contraction.
Curtis, Matthew W; Budyn, Elisa; Desai, Tejal A; Samarel, Allen M; Russell, Brenda
2013-01-01
Cardiac muscle cells are known to adapt to their physical surroundings, optimizing intracellular organization and contractile function for a given culture environment. A previously developed in vitro model system has shown that the inclusion of discrete microscale domains (or microrods) in three dimensions (3D) can alter long-term growth responses of neonatal ventricular myocytes. The aim of this work was to understand how cellular contact with such a domain affects various mechanical changes involved in cardiac muscle cell remodeling. Myocytes were maintained in 3D gels over 5 days in the presence or absence of 100-μm-long microrods, and the effect of this local heterogeneity on cell behavior was analyzed via several imaging techniques. Microrod abutment resulted in approximately twofold increases in the maximum displacement of spontaneously beating myocytes, as based on confocal microscopy scans of the gel xy-plane or the myocyte long axis. In addition, microrods caused significant increases in the proportion of aligned myofibrils (≤20° deviation from long axis) in fixed myocytes. Microrod-related differences in axial contraction could be abrogated by long-term interruption of certain signals of the RhoA-/Rho-associated kinase (ROCK) or protein kinase C (PKC) pathway. Furthermore, microrod-induced increases in myocyte size and protein content were prevented by ROCK inhibition. In all, the data suggest that microdomain heterogeneity in 3D appears to promote the development of axially aligned contractile machinery in muscle cells, an observation that may have relevance to a number of cardiac tissue engineering interventions. PMID:22407215
Microdomain heterogeneity in 3D affects the mechanics of neonatal cardiac myocyte contraction
Curtis, Matthew W.; Budyn, Elisa; Desai, Tejal A.; Samarel, Allen M.
2012-01-01
Cardiac muscle cells are known to adapt to their physical surroundings, optimizing intracellular organization and contractile function for a given culture environment. A previously developed in vitro model system has shown that the inclusion of discrete microscale domains (or microrods) in three dimensions (3D) can alter long-term growth responses of neonatal ventricular myocytes. The aim of this work was to understand how cellular contact with such a domain affects various mechanical changes involved in cardiac muscle cell remodeling. Myocytes were maintained in 3D gels over 5 days in the presence or absence of 100 – μm-long microrods, and the effect of this local heterogeneity on cell behavior was analyzed via several imaging techniques. Microrod abutment resulted in approximately twofold increases in the maximum displacement of spontaneously beating myocytes, as based on confocal microscopy scans of the gel xy-plane or the myocyte long axis. In addition, microrods caused significant increases in the proportion of aligned myofibrils (≤20° deviation from long axis) in fixed myocytes. Microrod-related differences in axial contraction could be abrogated by long-term interruption of certain signals of the RhoA-/Rho-associated kinase (ROCK) or protein kinase C (PKC) pathway. Furthermore, microrod-induced increases in myocyte size and protein content were prevented by ROCK inhibition. In all, the data suggest that microdomain heterogeneity in 3D appears to promote the development of axially aligned contractile machinery in muscle cells, an observation that may have relevance to a number of cardiac tissue engineering interventions. PMID:22407215
NASA Astrophysics Data System (ADS)
Thesenga, David; Town, James
2014-05-01
In February 2000, the Space Shuttle Endeavour flew a specially modified radar system during an 11-day mission. The purpose of the multinational Shuttle Radar Topography Mission (SRTM) was to "obtain elevation data on a near-global scale to generate the most complete high-resolution digital topographic database of Earth" by using radar interferometry. The data and resulting products are now publicly available for download and give a view of the landscape removed of vegetation, buildings, and other structures. This new view of the Earth's topography allows us to see previously unmapped or poorly mapped regions of the Earth as well as providing a level of detail that was previously unknown using traditional topographic mapping techniques. Understanding and appreciating the geographic terrain is a complex but necessary requirement for middle school aged (11-14yo) students. Abstract in nature, topographic maps and other 2D renderings of the Earth's surface and features do not address the inherent spatial challenges of a concrete-learner and traditional methods of teaching can at times exacerbate the problem. Technological solutions such as 3D-imaging in programs like Google Earth are effective but lack the tactile realness that can make a large difference in learning comprehension and retention for these young students. First developed in the 1980's, 3D printers were not commercial reality until recently and the rapid rise in interest has driven down the cost. With the advent of sub US1500 3D printers, this technology has moved out of the high-end marketplace and into the local office supply store. Schools across the US and elsewhere in the world are adding 3D printers to their technological workspaces and students have begun rapid-prototyping and manufacturing a variety of projects. This project attempted to streamline the process of transforming SRTM data from a GeoTIFF format by way of Python code. The resulting data was then inputted into a CAD-based program for
Mach-wave coherence in 3D media with random heterogeneities
NASA Astrophysics Data System (ADS)
Vyas, Jagdish C.; Mai, P. Martin; Galis, Martin; Dunham, Eric M.; Imperatori, Walter
2016-04-01
We investigate Mach-waves coherence for complex super-shear ruptures embedded in 3D random media that lead to seismic scattering. We simulate Mach-wave using kinematic earthquake sources that include fault-regions over which the rupture propagates at super-shear speed. The local slip rate is modeled with the regularized Yoffe function. The medium heterogeneities are characterized by Von Karman correlation function. We consider various realizations of 3D random media from combinations of different values of correlation length (0.5 km, 2 km, 5 km), standard deviation (5%, 10%, 15%) and Hurst exponent (0.2). Simulations in a homogeneous medium serve as a reference case. The ground-motion simulations (maximum resolved frequency of 5 Hz) are conducted by solving the elasto-dynamic equations of motions using a generalized finite-difference method, assuming a vertical strike-slip fault. The seismic wavefield is sampled at numerous locations within the Mach-cone region to study the properties and evolution of the Mach-waves in scattering media. We find that the medium scattering from random heterogeneities significantly diminishes the coherence of Mach-wave in terms of both amplitude and frequencies. We observe that Mach-waves are considerably scattered at distances RJB > 20 km (and beyond) for random media with standard deviation 10%. The scattering efficiency of the medium for small Hurst exponents (H <= 0.2) is mainly controlled by the standard deviation of the velocity heterogeneities, rather than their correlation length, as both theoretical considerations and numerical experiments indicate. Based on our simulations, we propose that local super-shear ruptures may be more common in nature then reported, but are very difficult to detect due to the strong seismic scattering. We suggest that if an earthquake is recorded within 10-15 km fault perpendicular distance and has high PGA, then inversion should be carried out by allowing rupture speed variations from sub
A NURBS-based generalized finite element scheme for 3D simulation of heterogeneous materials
NASA Astrophysics Data System (ADS)
Safdari, Masoud; Najafi, Ahmad R.; Sottos, Nancy R.; Geubelle, Philippe H.
2016-08-01
A 3D NURBS-based interface-enriched generalized finite element method (NIGFEM) is introduced to solve problems with complex discontinuous gradient fields observed in the analysis of heterogeneous materials. The method utilizes simple structured meshes of hexahedral elements that do not necessarily conform to the material interfaces in heterogeneous materials. By avoiding the creation of conforming meshes used in conventional FEM, the NIGFEM leads to significant simplification of the mesh generation process. To achieve an accurate solution in elements that are crossed by material interfaces, the NIGFEM utilizes Non-Uniform Rational B-Splines (NURBS) to enrich the solution field locally. The accuracy and convergence of the NIGFEM are tested by solving a benchmark problem. We observe that the NIGFEM preserves an optimal rate of convergence, and provides additional advantages including the accurate capture of the solution fields in the vicinity of material interfaces and the built-in capability for hierarchical mesh refinement. Finally, the use of the NIGFEM in the computational analysis of heterogeneous materials is discussed.
What spherically symmetric viscosity structure produces the same PGR as a realistic 3D Earth?
NASA Astrophysics Data System (ADS)
Paulson, A.; Zhong, S.; Wahr, J.
2003-04-01
Observations of isostatic adjustment of the earth's surface due to transient loading provide important constraints on the mantle viscosity structure. However, most studies of this response have assumed a spherically symmetric (1D) earth. This study is motivated by the following question: when a one-dimensional viscosity model is derived from post-glacial rebound (PGR) observations, how does this 1D structure correspond to the three-dimensional structure of the earth? Using the 3D spherical finite element software CitcomSVE [Zhong et al., 2002], we are able to compute the earth's response to realistic glacial loading when the earth has a truly 3D viscosity structure. The loading is provided by the ICE-3G deglaciation history [Tushingham &Peltier, 1991]. The 3D viscosity structure is constructed by first selecting a priori a radial average viscosity (for example, ( 1021 \\: {Pa \\cdot s}) in the upper mantle and (2 × 1021 \\: {Pa \\cdot s}) in the lower mantle). The lateral variations about this radial structure are derived from seismic shear-velocity tomography models by converting velocities to temperature, then temperature to viscosity. The seismic tomography models used are S20RTS [Ritsema et al., 1999] and NA00 [Van der Lee, 2002]. From the computed isostatic response, we measure typical PGR observables: relative sea level change (RSLC) and (dot{J2}). These measurements are then treated as synthetic data, and we search for 1D (radially stratified) viscosity models, forced with the same glaciation history, that will best fit these synthetic PGR observations. We find that for sites near the center of a large glacial load (e.g., southern Hudson Bay), a local average of the 3D viscosity structure provides a reasonable 1D proxy. For sites along the periphery of the glacial load (e.g., Boston), it is much more difficult to find a 1D model that can reproduce the 3D observations. We also approach the problem by running an ensemble of 1D viscosity models, and finding
NASA Astrophysics Data System (ADS)
Coviello, Velio; Manconi, Andrea; Occhiena, Cristina; Arattano, Massimo; Scavia, Claudio
2013-04-01
reliability. Preliminary results show a relevant improvement by considering a more realistic 3D heterogeneous velocity model for the relocation of seismic recordings due to artificial sources with respect to a homogeneous one. Future work will focus on the evaluation of the efficiency of the here presented 3D velocity model for the location of natural microseismic sources. [1] Occhiena, C., Coviello, V., Arattano, M., Chiarle, M., Morra di Cella, U., Pirulli, M., Pogliotti, P., and Scavia, C.: Analysis of microseismic signals and temperature recordings for rock slope stability investigations in high mountain areas, Nat. Hazards Earth Syst. Sci., 12, 2283-2298, doi:10.5194/nhess-12-2283-2012, 2012.
Seo, Sujin; Zhou, Xiangfei; Liu, Gang Logan
2016-07-01
Plasmonic substrates have fixed sensitivity once the geometry of the structure is defined. In order to improve the sensitivity, significant research effort has been focused on designing new plasmonic structures, which involves high fabrication costs; however, a method is reported for improving sensitivity not by redesigning the structure but by simply assembling plasmonic nanoparticles (NPs) near the evanescent field of the underlying 3D plasmonic nanostructure. Here, a nanoscale Lycurgus cup array (nanoLCA) is employed as a base colorimetric plasmonic substrate and an assembly template. Compared to the nanoLCA, the NP assembled nanoLCA (NP-nanoLCA) exhibits much higher sensitivity for both bulk refractive index sensing and biotin-streptavidin binding detection. The limit of detection of the NP-nanoLCA is at least ten times smaller when detecting biotin-streptavidin conjugation. The numerical calculations confirm the importance of the additive plasmon coupling between the NPs and the nanoLCA for a denser and stronger electric field in the same 3D volumetric space. Tunable sensitivity is accomplished by controlling the number of NPs in each nanocup, or the number density of the hot spots. This simple yet scalable and cost-effective method of using additive heterogeneous plasmon coupling effects will benefit various chemical, medical, and environmental plasmon-based sensors. PMID:27206214
A 3D numerical study of antimicrobial persistence in heterogeneous multi-species biofilms.
Zhao, Jia; Shen, Ya; Haapasalo, Markus; Wang, Zhejun; Wang, Qi
2016-03-01
We develop a 3D hydrodynamic model to investigate the mechanism of antimicrobial persistence in a multi-species oral biofilm and its recovery after being treated by bisbiguanide chlorhexidine gluconate (CHX). In addition to the hydrodynamic transport in the spatially heterogeneous biofilm, the model also includes mechanisms of solvent-biomass interaction, bacterial phenotype conversion, and bacteria-drug interaction. A numerical solver for the model is developed using a second order numerical scheme in 3D space and time and implemented on GPUs for high-performance computing. The model is calibrated against a set of experimental data obtained using confocal laser scan microscopy (CLSM) on multi-species oral biofilms, where a quantitative agreement is reached. Our numerical results reveal that quorum sensing molecules and growth factors in this model are instrumental in biofilm formation and recovery after the antimicrobial treatment. In particular, we show that (i) young biofilms are more susceptible to the antimicrobial treatment than the mature ones, (ii) this phenomenon is strongly correlated with volume fractions of the persister and EPS in the biofilm being treated. This suggests that antimicrobial treatment should be best administered to biofilms earlier before they mature to produce a thick protective EPS layer. In addition, the numerical study also indicates that an antimicrobial effect can be achieved should a proper mechanism be devised to minimize the conversion of susceptible bacteria to persisters during and even after the treatment. PMID:26739374
Compilation of 3D global conductivity model of the Earth for space weather applications
NASA Astrophysics Data System (ADS)
Alekseev, Dmitry; Kuvshinov, Alexey; Palshin, Nikolay
2015-07-01
We have compiled a global three-dimensional (3D) conductivity model of the Earth with an ultimate goal to be used for realistic simulation of geomagnetically induced currents (GIC), posing a potential threat to man-made electric systems. Bearing in mind the intrinsic frequency range of the most intense disturbances (magnetospheric substorms) with typical periods ranging from a few minutes to a few hours, the compiled 3D model represents the structure in depth range of 0-100 km, including seawater, sediments, earth crust, and partly the lithosphere/asthenosphere. More explicitly, the model consists of a series of spherical layers, whose vertical and lateral boundaries are established based on available data. To compile a model, global maps of bathymetry, sediment thickness, and upper and lower crust thicknesses as well as lithosphere thickness are utilized. All maps are re-interpolated on a common grid of 0.25×0.25 degree lateral spacing. Once the geometry of different structures is specified, each element of the structure is assigned either a certain conductivity value or conductivity versus depth distribution, according to available laboratory data and conversion laws. A numerical formalism developed for compilation of the model, allows for its further refinement by incorporation of regional 3D conductivity distributions inferred from the real electromagnetic data. So far we included into our model four regional conductivity models, available from recent publications, namely, surface conductance model of Russia, and 3D conductivity models of Fennoscandia, Australia, and northwest of the United States.
Forcing Ferromagnetic Coupling Between Rare-Earth-Metal and 3d Ferromagnetic Films
NASA Astrophysics Data System (ADS)
Sanyal, Biplab; Antoniak, Carolin; Burkert, Till; Krumme, Bernhard; Warland, Anne; Stromberg, Frank; Praetorius, Christian; Fauth, Kai; Wende, Heiko; Eriksson, Olle
2010-04-01
Using density functional calculations, we have studied the magnetic properties of nanocomposites composed of rare-earth-metal elements in contact with 3d transition metals (Fe and Cr). We demonstrate the possibility to obtain huge magnetic moments in such nanocomposites, of order 10μB/rare-earth-metalatom, with a potential to reach the maximum magnetic moment of Fe-Co alloys at the top of the so-called Slater-Pauling curve. A first experimental proof of concept is given by thin-film synthesis of Fe/Gd and Fe/Cr/Gd nanocomposites, in combination with x-ray magnetic circular dichroism.
Development of hybrid 3-D hydrological modeling for the NCAR Community Earth System Model (CESM)
Zeng, Xubin; Troch, Peter; Pelletier, Jon; Niu, Guo-Yue; Gochis, David
2015-11-15
This is the Final Report of our four-year (3-year plus one-year no cost extension) collaborative project between the University of Arizona (UA) and the National Center for Atmospheric Research (NCAR). The overall objective of our project is to develop and evaluate the first hybrid 3-D hydrological model with a horizontal grid spacing of 1 km for the NCAR Community Earth System Model (CESM).
Quasi-heterogeneous efficient 3-D discrete ordinates CANDU calculations using Attila
Preeti, T.; Rulko, R.
2012-07-01
In this paper, 3-D quasi-heterogeneous large scale parallel Attila calculations of a generic CANDU test problem consisting of 42 complete fuel channels and a perpendicular to fuel reactivity device are presented. The solution method is that of discrete ordinates SN and the computational model is quasi-heterogeneous, i.e. fuel bundle is partially homogenized into five homogeneous rings consistently with the DRAGON code model used by the industry for the incremental cross-section generation. In calculations, the HELIOS-generated 45 macroscopic cross-sections library was used. This approach to CANDU calculations has the following advantages: 1) it allows detailed bundle (and eventually channel) power calculations for each fuel ring in a bundle, 2) it allows the exact reactivity device representation for its precise reactivity worth calculation, and 3) it eliminates the need for incremental cross-sections. Our results are compared to the reference Monte Carlo MCNP solution. In addition, the Attila SN method performance in CANDU calculations characterized by significant up scattering is discussed. (authors)
3D Observations of Dispersion, Mixing and Reaction in Heterogeneous Rocks
NASA Astrophysics Data System (ADS)
Boon, M.; Bijeljic, B.; Krevor, S. C.
2015-12-01
Rock structure heterogeneity can have a significant effect on dispersion, mixing and reaction of aqueous components in porous media. To observe the effect of pore structure heterogeneity on reactive transport, core flooding experiments were carried out for a sandstone and two carbonate rocks of different heterogeneity for eight different Peclet numbers ranging from 0.5 to 100. The rock cores were 20cm long and had a diameter of 7.62cm. A device consisting out of three annular regions was used for injection (fig. 1A). An analytical solution to the flow and transport equations for this new inlet configuration was derived to design the experiments (fig. 1B).The dispersion, mixing and reaction of chemical components were visualised in 3D with the use of chemical dopants and a medical CT scanner (fig. 1C). Local transverse dispersion coefficients (Dt) were calculated from the change in variance of the transverse distance travelled by the chemical dopant along the core. The change in variance along the core showed a characteristic pattern for each rock that was independent of spatial location. Heterogeneity was characterized by the spread in local transverse dispersion coefficients. For Peclet number 2, for the homogenous rock the local transverse dispersion coefficients ranged from 4.1x10-4 cm2 min-1 to 5.9x10-4 cm2 min-1 and for the most heterogeneous rock from 2.5x10-3 cm2 min-1 to 7.2x10-3 cm2 min-1. For the reactive transport experiments, an ICP-MS was used to measure the effluent. The core flooding experiments were modelled using both, the CrunchFlow and ToughReact reactive transport codes. High quality data sets of the space and time evolution of the concentration in non-reactive and reactive core-flooding experiments like these can be used as future benchmark test for numerical models. Furthermore, these observations can be used in the development of upscaling techniques for accurate and efficient modelling of chemical processes during flow in porous media.
NASA Astrophysics Data System (ADS)
Jung, Jin Woo; Lee, Jung-Seob; Cho, Dong-Woo
2016-02-01
Recently, much attention has focused on replacement or/and enhancement of biological tissues via the use of cell-laden hydrogel scaffolds with an architecture that mimics the tissue matrix, and with the desired three-dimensional (3D) external geometry. However, mimicking the heterogeneous tissues that most organs and tissues are formed of is challenging. Although multiple-head 3D printing systems have been proposed for fabricating heterogeneous cell-laden hydrogel scaffolds, to date only the simple exterior form has been realized. Here we describe a computer-aided design and manufacturing (CAD/CAM) system for this application. We aim to develop an algorithm to enable easy, intuitive design and fabrication of a heterogeneous cell-laden hydrogel scaffolds with a free-form 3D geometry. The printing paths of the scaffold are automatically generated from the 3D CAD model, and the scaffold is then printed by dispensing four materials; i.e., a frame, two kinds of cell-laden hydrogel and a support. We demonstrated printing of heterogeneous tissue models formed of hydrogel scaffolds using this approach, including the outer ear, kidney and tooth tissue. These results indicate that this approach is particularly promising for tissue engineering and 3D printing applications to regenerate heterogeneous organs and tissues with tailored geometries to treat specific defects or injuries.
Jung, Jin Woo; Lee, Jung-Seob; Cho, Dong-Woo
2016-01-01
Recently, much attention has focused on replacement or/and enhancement of biological tissues via the use of cell-laden hydrogel scaffolds with an architecture that mimics the tissue matrix, and with the desired three-dimensional (3D) external geometry. However, mimicking the heterogeneous tissues that most organs and tissues are formed of is challenging. Although multiple-head 3D printing systems have been proposed for fabricating heterogeneous cell-laden hydrogel scaffolds, to date only the simple exterior form has been realized. Here we describe a computer-aided design and manufacturing (CAD/CAM) system for this application. We aim to develop an algorithm to enable easy, intuitive design and fabrication of a heterogeneous cell-laden hydrogel scaffolds with a free-form 3D geometry. The printing paths of the scaffold are automatically generated from the 3D CAD model, and the scaffold is then printed by dispensing four materials; i.e., a frame, two kinds of cell-laden hydrogel and a support. We demonstrated printing of heterogeneous tissue models formed of hydrogel scaffolds using this approach, including the outer ear, kidney and tooth tissue. These results indicate that this approach is particularly promising for tissue engineering and 3D printing applications to regenerate heterogeneous organs and tissues with tailored geometries to treat specific defects or injuries. PMID:26899876
Jung, Jin Woo; Lee, Jung-Seob; Cho, Dong-Woo
2016-01-01
Recently, much attention has focused on replacement or/and enhancement of biological tissues via the use of cell-laden hydrogel scaffolds with an architecture that mimics the tissue matrix, and with the desired three-dimensional (3D) external geometry. However, mimicking the heterogeneous tissues that most organs and tissues are formed of is challenging. Although multiple-head 3D printing systems have been proposed for fabricating heterogeneous cell-laden hydrogel scaffolds, to date only the simple exterior form has been realized. Here we describe a computer-aided design and manufacturing (CAD/CAM) system for this application. We aim to develop an algorithm to enable easy, intuitive design and fabrication of a heterogeneous cell-laden hydrogel scaffolds with a free-form 3D geometry. The printing paths of the scaffold are automatically generated from the 3D CAD model, and the scaffold is then printed by dispensing four materials; i.e., a frame, two kinds of cell-laden hydrogel and a support. We demonstrated printing of heterogeneous tissue models formed of hydrogel scaffolds using this approach, including the outer ear, kidney and tooth tissue. These results indicate that this approach is particularly promising for tissue engineering and 3D printing applications to regenerate heterogeneous organs and tissues with tailored geometries to treat specific defects or injuries. PMID:26899876
On the computation of long period seismograms in a 3-D earth using normal mode based approximations
NASA Astrophysics Data System (ADS)
Romanowicz, Barbara A.; Panning, Mark P.; Gung, Yuancheng; Capdeville, Yann
2008-11-01
Tomographic inversions for large-scale structure of the earth's mantle involve a forward modelling step of wave propagation through 3-D heterogeneity. Until now, most investigators have worked in the framework of the simplest theoretical assumptions, namely the infinite frequency `ray theory' in the case of body wave traveltime inversions, or the `path-average' approximation (PAVA) to normal mode perturbation theory, in the case of surface waves and long-period waveforms. As interest is shifting to mapping shorter wavelength structures, the need for a more accurate theoretical account of the interaction of seismic waves with mantle heterogeneity, coupled with improvements in path coverage, has been realized. Here we discuss different levels of approximations used in the context of normal mode perturbation theory, when modelling time domain seismic waveforms. We compare the performance of asymptotic approximations, which collapse the effects of 3-D structure onto the great circle vertical plane: the 1-D PAVA and a 2-D approximation called non-linear asymptotic coupling theory (NACT), which both are zeroth order asymptotic approximations. We then discuss how off-vertical plane effects can be introduced using higher order asymptotics. These computationally efficient approximations are compared to the linear Born formalism (BORN), which computes scattering integrals over the entire surface of the sphere. We point out some limitations of this linear formalism in the case of spatially extended anomalies, and show how that can be remedied through the introduction of a non-linear term (NBORN). All these approximations are referenced to a precise 3-D numerical computation afforded by the spectral element method. We discuss simple geometries, and explore a range of sizes of anomalies compared to the wavelength of the seismic waves considered, thus illustrating the range of validity and limitations of the various approximations considered.
Synthetic Seismograms for Realistic 3D Earth Model with Anisotropic Inner Core
NASA Astrophysics Data System (ADS)
Tsuboi, S.; Tono, Y.
2006-12-01
We have demonstrated that we can calculate global theoretical seismograms for realistic 3D Earth models based upon the combination of a precise numerical technique (the spectral-element method) and a sufficiently fast supercomputer (the Earth Simulator) [Tsuboi et al, 2003]. Here we have calculated synthetic seismograms by using model S20RTS of the mantle (Ritsema et al., 1999), model CRUST2.0 of the crust (Basin et al., 2000), topography and bathymetry model ETOPO5, and anisotropic inner core model (Ishii 2002). The calculations are performed on 4056 processors, which require 507 out of 640 nodes of the Earth Simulator. These synthetics are computed by using SPECFEM3D(Komatitsch and Tromp, 2002) and are accurate up to 3.5 seconds. We have calculated these synthetics with aisotropic inner core model for several earthquakes and compared with the synthetics which are calculated for isotropic inner core model. Preliminary comparison shows that the travel time differences between anisotropic inner core model and isotropic core model for PKPab phases are at most a few seconds. There seems to be no significant differences in waveforms of PKP phases. These differences in travel times may help us to improve inner core fine structure by comparing these synthetics with observation.
Glacial isostatic adjustment model with composite 3-D Earth rheology for Fennoscandia
NASA Astrophysics Data System (ADS)
van der Wal, Wouter; Barnhoorn, Auke; Stocchi, Paolo; Gradmann, Sofie; Wu, Patrick; Drury, Martyn; Vermeersen, Bert
2013-07-01
Models for glacial isostatic adjustment (GIA) can provide constraints on rheology of the mantle if past ice thickness variations are assumed to be known. The Pleistocene ice loading histories that are used to obtain such constraints are based on an a priori 1-D mantle viscosity profile that assumes a single deformation mechanism for mantle rocks. Such a simplified viscosity profile makes it hard to compare the inferred mantle rheology to inferences from seismology and laboratory experiments. It is unknown what constraints GIA observations can provide on more realistic mantle rheology with an ice history that is not based on an a priori mantle viscosity profile. This paper investigates a model for GIA with a new ice history for Fennoscandia that is constrained by palaeoclimate proxies and glacial sediments. Diffusion and dislocation creep flow law data are taken from a compilation of laboratory measurements on olivine. Upper-mantle temperature data sets down to 400 km depth are derived from surface heatflow measurements, a petrochemical model for Fennoscandia and seismic velocity anomalies. Creep parameters below 400 km are taken from an earlier study and are only varying with depth. The olivine grain size and water content (a wet state, or a dry state) are used as free parameters. The solid Earth response is computed with a global spherical 3-D finite-element model for an incompressible, self-gravitating Earth. We compare predictions to sea level data and GPS uplift rates in Fennoscandia. The objective is to see if the mantle rheology and the ice model is consistent with GIA observations. We also test if the inclusion of dislocation creep gives any improvements over predictions with diffusion creep only, and whether the laterally varying temperatures result in an improved fit compared to a widely used 1-D viscosity profile (VM2). We find that sea level data can be explained with our ice model and with information on mantle rheology from laboratory experiments
3D heterogeneous sensor system on a chip for defense and security applications
NASA Astrophysics Data System (ADS)
Bhansali, Shekhar; Chapman, Glenn H.; Friedman, Eby G.; Ismail, Yehea; Mukund, P. R.; Tebbe, Dennis; Jain, Vijay K.
2004-09-01
This paper describes a new concept for ultra-small, ultra-compact, unattended multi-phenomenological sensor systems for rapid deployment, with integrated classification-and-decision-information extraction capability from a sensed environment. We discuss a unique approach, namely a 3-D Heterogeneous System on a Chip (HSoC) in order to achieve a minimum 10X reduction in weight, volume, and power and a 10X or greater increase in capability and reliability -- over the alternative planar approaches. These gains will accrue from (a) the avoidance of long on-chip interconnects and chip-to-chip bonding wires, and (b) the cohabitation of sensors, preprocessing analog circuitry, digital logic and signal processing, and RF devices in the same compact volume. A specific scenario is discussed in detail wherein a set of four types of sensors, namely an array of acoustic and seismic sensors, an active pixel sensor array, and an uncooled IR imaging array are placed on a common sensor plane. The other planes include an analog plane consisting of transductors and A/D converters. The digital processing planes provide the necessary processing and intelligence capability. The remaining planes provide for wireless communications/networking capability. When appropriate, this processing and decision-making will be accomplished on a collaborative basis among the distributed sensor nodes through a wireless network.
NASA Astrophysics Data System (ADS)
Arroucau, Pierre; Custódio, Susana
2015-04-01
Solving inverse problems requires an estimate of data uncertainties. This usually takes the form of a data covariance matrix, which determines the shape of the model posterior distribution. Those uncertainties are yet not always known precisely and it is common practice to simply set them to a fixed, reasonable value. In the case of earthquake location, the hypocentral parameters (longitude, latitude, depth and origin time) are typically inverted for using seismic phase arrival times. But quantitative data variance estimates are rarely provided. Instead, arrival time catalogs usually associate phase picks with a quality factor, which is subsequently interpreted more or less arbitrarily in terms of data uncertainty in the location procedure. Here, we present a hierarchical Bayesian algorithm for earthquake location in 3D heterogeneous media, in which not only the earthquake hypocentral parameters, but also the P- and S-wave arrival time uncertainties, are inverted for, hence allowing more realistic posterior model covariance estimates. Forward modeling is achieved by means of the Fast Marching Method (FMM), an eikonal solver which has the ability to take interfaces into account, so direct, reflected and refracted phases can be used in the inversion. We illustrate the ability of our algorithm to retrieve earthquake hypocentral parameters as well as data uncertainties through synthetic examples and using a subset of arrival time catalogs for mainland Portugal and its Atlantic margin.
NASA Astrophysics Data System (ADS)
Arroucau, P.; Custodio, S.
2014-12-01
Solving inverse problems requires an estimate of data uncertainties. This usually takes the form of a data covariance matrix, which determines the shape of the model posterior distribution. Those uncertainties are yet not always known precisely and it is common practice to simply set them to a fixed, reasonable value. In the case of earthquake location, the hypocentral parameters (longitude, latitude, depth and origin time) are typically inverted for using seismic phase arrival times. But quantitative data variance estimates are rarely provided. Instead, arrival time catalogs usually associate phase picks with a quality factor, which is subsequently interpreted more or less arbitrarily in terms of data uncertainty in the location procedure. Here, we present a hierarchical Bayesian algorithm for earthquake location in 3D heterogeneous media, in which not only the earthquake hypocentral parameters, but also the P- and S-wave arrival time uncertainties, are inverted for, hence allowing more realistic posterior model covariance estimates. Forward modeling is achieved by means of the Fast Marching Method (FMM), an eikonal solver which has the ability to take interfaces into account, so direct, reflected and refracted phases can be used in the inversion. We illustrate the ability of our algorithm to retrieve earthquake hypocentral parameters as well as data uncertainties through synthetic examples and using a subset of arrival time catalogs for mainland Portugal and its Atlantic margin.
NASA Astrophysics Data System (ADS)
Tandon, K.; Egbert, G.; Siripunvaraporn, W.
2003-12-01
We are developing a modular system for three-dimensional inversion of electromagnetic (EM) induction data, using an object oriented programming approach. This approach allows us to modify the individual components of the inversion scheme proposed, and also reuse the components for variety of problems in earth science computing howsoever diverse they might be. In particular, the modularity allows us to (a) change modeling codes independently of inversion algorithm details; (b) experiment with new inversion algorithms; and (c) modify the way prior information is imposed in the inversion to test competing hypothesis and techniques required to solve an earth science problem. Our initial code development is for EM induction equations on a staggered grid, using iterative solution techniques in 3D. An example illustrated here is an experiment with the sensitivity of 3D magnetotelluric inversion to uncertainties in the boundary conditions required for regional induction problems. These boundary conditions should reflect the large-scale geoelectric structure of the study area, which is usually poorly constrained. In general for inversion of MT data, one fixes boundary conditions at the edge of the model domain, and adjusts the earth?s conductivity structure within the modeling domain. Allowing for errors in specification of the open boundary values is simple in principle, but no existing inversion codes that we are aware of have this feature. Adding a feature such as this is straightforward within the context of the modular approach. More generally, a modular approach provides an efficient methodology for setting up earth science computing problems to test various ideas. As a concrete illustration relevant to EM induction problems, we investigate the sensitivity of MT data near San Andreas Fault at Parkfield (California) to uncertainties in the regional geoelectric structure.
Earth Science Research Discovery, Integration, 3D Visualization and Analysis using NASA World Wind
NASA Astrophysics Data System (ADS)
Alameh, N.; Hogan, P.
2008-12-01
NASA plays a leadership role in the world of Advanced Information Technologies. Part of our mission is to leverage those technologies to increase the usability of the growing amount of earth observation produced by the science community. NASA World Wind open source technology provides a complete 3D visualization platform that is being continually advanced by NASA, its partners and the open source community. The technology makes scientific data and observations more accessible to Earth scientists and offers them a standards-based extensible platform to manipulate and analyze that data. The API-centric architecture of World Wind's SDK allows others to readily extend or embed this technology (including in web pages). Such multiple approaches to using the technology accelerate opportunities for the research community to provide "advances in fundamental understanding of the Earth system and increased application of this understanding to serve the nation and the people of the world" (NRC Decadal Survey). The opportunities to advance this NASA Open Source Agreement (NOSA) technology by leveraging advances in web services, interoperability, data discovery mechanisms, and Sensor Web are unencumbered by proprietary constraints and therefore provide the basis for an evolving platform that can reliably service the needs of the Earth Science, Sensor Web and GEOSS communities. The ability for these communities to not only use this technology in an unrestricted manner but to also participate in advancing it leads to accelerated innovation and maximum exchange of information. 3 characteristics enable World Wind to push the frontier in Advanced Information Systems: 1- World Wind provides a unifying information browser to enable a variety of 3D geospatial applications. World Wind consists of a coherent suite of modular components to be used selectively or in concert with any number of programs. 2- World Wind technology can be embedded as part of any application and hence makes it
Sridharan, Anush; Eisenbrey, John R.; Machado, Priscilla; Ojeda-Fournier, Haydee; Wilkes, Annina; Sevrukov, Alexander; Mattrey, Robert F.; Wallace, Kirk; Chalek, Carl L.; Thomenius, Kai E.; Forsberg, Flemming
2015-01-01
Ability to visualize breast lesion vascularity and quantify the vascular heterogeneity using contrast-enhanced 3-D harmonic (HI) and subharmonic (SHI) ultrasound imaging was investigated in a clinical population. Patients (n = 134) identified with breast lesions on mammography were scanned using power Doppler imaging, contrast-enhanced 3-D HI, and 3-D SHI on a modified Logiq 9 scanner (GE Healthcare). A region of interest corresponding to ultrasound contrast agent flow was identified in 4D View (GE Medical Systems) and mapped to raw slice data to generate a map of time-intensity curves for the lesion volume. Time points corresponding to baseline, peak intensity, and washout of ultrasound contrast agent were identified and used to generate and compare vascular heterogeneity plots for malignant and benign lesions. Vascularity was observed with power Doppler imaging in 84 lesions (63 benign and 21 malignant). The 3-D HI showed flow in 8 lesions (5 benign and 3 malignant), whereas 3-D SHI visualized flow in 68 lesions (49 benign and 19 malignant). Analysis of vascular heterogeneity in the 3-D SHI volumes found benign lesions having a significant difference in vascularity between central and peripheral sections (1.71 ± 0.96 vs. 1.13 ± 0.79 dB, p < 0.001, respectively), whereas malignant lesions showed no difference (1.66 ± 1.39 vs. 1.24 ± 1.14 dB, p = 0.24), indicative of more vascular coverage. These preliminary results suggest quantitative evaluation of vascular heterogeneity in breast lesions using contrast-enhanced 3-D SHI is feasible and able to detect variations in vascularity between central and peripheral sections for benign and malignant lesions. PMID:25935933
Moving beyond flat earth: dense 3D scene reconstruction from a single FL-LWIR camera
NASA Astrophysics Data System (ADS)
Stone, K.; Keller, J. M.; Anderson, D. T.
2013-06-01
In previous work an automatic detection system for locating buried explosive hazards in forward-looking longwave infrared (FL-LWIR) and forward-looking ground penetrating radar (FL-GPR) data was presented. This system consists of an ensemble of trainable size-contrast filters prescreener coupled with a secondary classification step which extracts cell-structured image space features, such as local binary patterns (LBP), histogram of oriented gradients (HOG), and edge histogram descriptors (EHD), from multiple looks and classifies the resulting feature vectors using a support vector machine. Previously, this system performed image space to UTM coordinate mapping under a flat earth assumption. This limited its applicability to flat terrain and short standoff distances. This paper demonstrates a technique for dense 3D scene reconstruction from a single vehicle mounted FL-LWIR camera. This technique utilizes multiple views and standard stereo vision algorithms such as polar rectification and optimal correction. Results for the detection algorithm using this 3D scene reconstruction approach on data from recent collections at an arid US Army test site are presented. These results are compared to those obtained under the flat earth assumption, with special focus on rougher terrain and longer standoff distance than in previous experiments. The most recent collection also allowed comparison between uncooled and cooled FL-LWIR cameras for buried explosive hazard detection.
3D Global Magnetohydrodynamic Simulations of the Solar Wind/Earth's Magnetosphere Interaction
NASA Astrophysics Data System (ADS)
Yalim, M. S.; Poedts, S.
2014-09-01
In this paper, we present results of real-time 3D global magnetohydrodynamic (MHD) simulations of the solar wind interaction with the Earth's magnetosphere using time-varying data from the NASA Advanced Composition Explorer (ACE) satellite during a few big magnetic storm events of the previous and current solar cycles, namely the 06 April 2000, 20 November 2003 and 05 April 2010 storms. We introduce a numerical magnetic storm index and compare the geo-effectiveness of these events in terms of this storm index which is a measure for the resulting global perturbation of the Earth's magnetic field. Steady simulations show that the upstream solar wind plasma parameters enter the low-β switch-on regime for some time intervals during a magnetic storm causing a complex dimpled bow shock structure. We also investigate the traces of such bow shock structures during time-dependent simulations of the events. We utilize a 3D, implicit, parallel, unstructured grid, compressible finite volume ideal MHD solver with an anisotropic grid adaptation technique for the computer simulations.
NASA Astrophysics Data System (ADS)
Baraka, S. M.; Ben-Jaffel, L. B.
2014-12-01
We use particle-in-cell PIC 3D Electromagnetic, relativistic global code to address large-scale problems in magnetosphere electrodynamics. Terrestrial bow shock is simulated as an example. 3D Magnetohydrodynamics model ,MHD GUMICS in CCMC project, have been used in parallel with PIC under same scaled Solar wind (SW) and IMF conditions. We report new results from the coupling between the two models. Further investigations are required for confirmations of these results. In both codes the Earth's bow shock position is found at ~14.8 RE along the Sun-Earth line, and ~29 RE on the dusk side which is consistent with past in situ observation. Both simulations reproduce the theoretical jump conditions at the shock. However, PIC code density and temperature distributions are inflated and slightly shifted sunward when compared to MHD results. Reflected ions upstream of the bow shock may cause this sunward shift for density and temperature. Distribution of reflected ions and electrons are shown in the foreshock region, within the transition of the shock and in the downstream. The current version of PIC code can be run under modest computing facilities and resources. Additionally, existing MHD simulations should be useful to calibrate scaled properties of plasma resulting from PIC simulations for comparison with observations. Similarities and drawbacks of the results obtained by the two models are listed. The ultimate goal of using these different models in a complimentary manner rather than competitive is to better understand the macrostructure of the magnetosphere
NASA Astrophysics Data System (ADS)
Duru, Kenneth; Dunham, Eric M.
2016-01-01
Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration; b) dynamic propagation of earthquake ruptures along nonplanar faults; and c) accurate propagation of seismic waves in heterogeneous media with free surface topography. We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts (SBP) finite difference operators in space. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. The finite difference stencils used in this paper are sixth order accurate in the interior and third order accurate close to the boundaries. However, the method is applicable to any spatial operator with a diagonal norm satisfying the SBP property. Time stepping is performed with a 4th order accurate explicit low storage Runge-Kutta scheme, thus yielding a globally fourth order accurate method in both space and time. We show numerical simulations on band limited self-similar fractal faults revealing the complexity of rupture dynamics
Cooperative 3D and 2D mapping with heterogenous ground robots
NASA Astrophysics Data System (ADS)
Rogers, John G., III; Baran, David; Stump, Ethan; Young, Stuart; Christensen, Henrik I.
2012-06-01
Efficient and accurate 3D mapping is desirable in disaster recovery as well as urban warfare situations. The speed with which these maps can be generated is vital to provide situational awareness in these situations. A team of mobile robots can work together to build maps more quickly. We present an algorithm by which a team of mobile robots can merge 2D and 3D measurements to build a 3D map, together with experiments performed at a military test facility.
A 3D Earth orbit model; visualization and analysis of Milankovitch cycles and insolation
NASA Astrophysics Data System (ADS)
Gilb, R. D.; Kostadinov, T. S.
2012-12-01
An astronomically precise and accurate Earth orbit graphical model, Earth orbit v2.0, is presented. The model offers 3D visualizations of Earth's orbital geometry, Milankovitch parameters and the ensuing insolation forcings. Prevalent paleoclimatic theories invoke Milankovitch cycles as a major forcing mechanism capable of shifting Earth's climate regimes on time scales of tens to hundreds of thousands of years. Variability of eccentricity (ellipticity of orbit), precession (longitude of perihelion) and obliquity (Earth's axial tilt) changes parameters such as amplitude of seasonal insolation, timing of seasons with respect to perihelion, and total annual insolation. Hays et al. (1976) demonstrated a strong link between Milankovitch cycles and paleoclimatological records, which has been confirmed and expanded many times since (e.g. Berger et al., 1994; Berger et al., 2010). The complex interplay of several orbital parameters on various time scales makes assessment and visualization of Earth's orbit and spatio-temporal insolation variability challenging. It is difficult to appreciate the pivotal importance of Kepler's laws of planetary motion in controlling the effects of Milankovitch cycles on insolation patterns on various spatio-temporal scales. These factors also make Milankovitch theory difficult to teach effectively. The model allows substantial user control in a robust, yet intuitive and user-friendly graphical user interface (GUI) developed in Matlab. We present the user with a choice between Berger et al. (1978) and Laskar et al. (2004) astronomical solutions for eccentricity, obliquity and precession. Berger solutions span from -1 Myr to +1 Myr, while Laskar provides solutions from -101 Myr to +21 Myr since J2000. Users can also choose a "demo" mode which allows the three Milankovitch parameters to be varied independently of each other, so the user can isolate the effects of each on orbital geometry and insolation. For example, extreme eccentricity can be
3D Elastic Seismic Wave Propagation Code
Energy Science and Technology Software Center (ESTSC)
1998-09-23
E3D is capable of simulating seismic wave propagation in a 3D heterogeneous earth. Seismic waves are initiated by earthquake, explosive, and/or other sources. These waves propagate through a 3D geologic model, and are simulated as synthetic seismograms or other graphical output.
NASA Astrophysics Data System (ADS)
Parisi, Laura; Ferreira, Ana M. G.
2016-04-01
The surface wave full ray theory (FRT) is an efficient tool to calculate synthetic waveforms of surface waves. It combines the concept of local modes with exact ray tracing as a function of frequency, providing a more complete description of surface wave propagation than the widely used great circle approximation (GCA). The purpose of this study is to evaluate the ability of the FRT approach to model teleseismic long-period surface waveforms (T ˜ 45-150 s) in the context of current 3-D Earth models to empirically assess its validity domain and its scope for future studies in seismic tomography. To achieve this goal, we compute vertical and horizontal component fundamental mode synthetic Rayleigh waveforms using the FRT, which are compared with calculations using the highly accurate spectral element method. We use 13 global earth models including 3-D crustal and mantle structure, which are derived by successively varying the strength and lengthscale of heterogeneity in current tomographic models. For completeness, GCA waveforms are also compared with the spectral element method. We find that the FRT accurately predicts the phase and amplitude of long-period Rayleigh waves (T ˜ 45-150 s) for almost all the models considered, with errors in the modelling of the phase (amplitude) of Rayleigh waves being smaller than 5 per cent (10 per cent) in most cases. The largest errors in phase and amplitude are observed for T ˜ 45 s and for the three roughest earth models considered that exhibit shear wave anomalies of up to ˜20 per cent, which is much larger than in current global tomographic models. In addition, we find that overall the GCA does not predict Rayleigh wave amplitudes well, except for the longest wave periods (T ˜ 150 s) and the smoothest models considered. Although the GCA accurately predicts Rayleigh wave phase for current earth models such as S20RTS and S40RTS, FRT's phase errors are smaller, notably for the shortest wave periods considered (T ˜ 45 s and
KAGLVis - On-line 3D Visualisation of Earth-observing-satellite Data
NASA Astrophysics Data System (ADS)
Szuba, Marek; Ameri, Parinaz; Grabowski, Udo; Maatouki, Ahmad; Meyer, Jörg
2015-04-01
One of the goals of the Large-Scale Data Management and Analysis project is to provide a high-performance framework facilitating management of data acquired by Earth-observing satellites such as Envisat. On the client-facing facet of this framework, we strive to provide visualisation and basic analysis tool which could be used by scientists with minimal to no knowledge of the underlying infrastructure. Our tool, KAGLVis, is a JavaScript client-server Web application which leverages modern Web technologies to provide three-dimensional visualisation of satellite observables on a wide range of client systems. It takes advantage of the WebGL API to employ locally available GPU power for 3D rendering; this approach has been demonstrated to perform well even on relatively weak hardware such as integrated graphics chipsets found in modern laptop computers and with some user-interface tuning could even be usable on embedded devices such as smartphones or tablets. Data is fetched from the database back-end using a ReST API and cached locally, both in memory and using HTML5 Web Storage, to minimise network use. Computations, calculation of cloud altitude from cloud-index measurements for instance, can depending on configuration be performed on either the client or the server side. Keywords: satellite data, Envisat, visualisation, 3D graphics, Web application, WebGL, MEAN stack.
3D time-domain airborne EM modeling for an arbitrarily anisotropic earth
NASA Astrophysics Data System (ADS)
Yin, Changchun; Qi, Yanfu; Liu, Yunhe
2016-08-01
Time-domain airborne EM data is currently interpreted based on an isotropic model. Sometimes, it can be problematic when working in the region with distinct dipping stratifications. In this paper, we simulate the 3D time-domain airborne EM responses over an arbitrarily anisotropic earth with topography by edge-based finite-element method. Tetrahedral meshes are used to describe the abnormal bodies with complicated shapes. We further adopt the Backward Euler scheme to discretize the time-domain diffusion equation for electric field, obtaining an unconditionally stable linear equations system. We verify the accuracy of our 3D algorithm by comparing with 1D solutions for an anisotropic half-space. Then, we switch attentions to effects of anisotropic media on the strengths and the diffusion patterns of time-domain airborne EM responses. For numerical experiments, we adopt three typical anisotropic models: 1) an anisotropic anomalous body embedded in an isotropic half-space; 2) an isotropic anomalous body embedded in an anisotropic half-space; 3) an anisotropic half-space with topography. The modeling results show that the electric anisotropy of the subsurface media has big effects on both the strengths and the distribution patterns of time-domain airborne EM responses; this effect needs to be taken into account when interpreting ATEM data in areas with distinct anisotropy.
A 3D Visualization and Analysis Model of the Earth Orbit, Milankovitch Cycles and Insolation.
NASA Astrophysics Data System (ADS)
Kostadinov, Tihomir; Gilb, Roy
2013-04-01
Milankovitch theory postulates that periodic variability of Earth's orbital elements is a major climate forcing mechanism. Although controversies remain, ample geologic evidence supports the major role of the Milankovitch cycles in climate, e.g. glacial-interglacial cycles. There are three Milankovitch orbital parameters: orbital eccentricity (main periodicities of ~100,000 and ~400,000 years), precession (quantified as the longitude of perihelion, main periodicities 19,000-24,000 years) and obliquity of the ecliptic (Earth's axial tilt, main periodicity 41,000 years). The combination of these parameters controls the spatio-temporal patterns of incoming solar radiation (insolation) and the timing of the seasons with respect to perihelion, as well as season duration. The complex interplay of the Milankovitch orbital parameters on various time scales makes assessment and visualization of Earth's orbit and insolation variability challenging. It is difficult to appreciate the pivotal importance of Kepler's laws of planetary motion in controlling the effects of Milankovitch cycles on insolation patterns. These factors also make Earth-Sun geometry and Milankovitch theory difficult to teach effectively. Here, an astronomically precise and accurate Earth orbit visualization model is presented. The model offers 3D visualizations of Earth's orbital geometry, Milankovitch parameters and the ensuing insolation forcings. Both research and educational uses are envisioned for the model, which is developed in Matlab® as a user-friendly graphical user interface (GUI). We present the user with a choice between the Berger et al. (1978) and Laskar et al. (2004) astronomical solutions for eccentricity, obliquity and precession. A "demo" mode is also available, which allows the three Milankovitch parameters to be varied independently of each other (and over much larger ranges than the naturally occurring ones), so the user can isolate the effects of each parameter on orbital geometry
Early Earth tectonics: A high-resolution 3D numerical modelling approach
NASA Astrophysics Data System (ADS)
Fischer, R.; Gerya, T.
2015-12-01
Early Earth had a higher amount of radiogenic elements as well as a higher amount of leftover primordial heat. Both contribute to the increased temperature in the Earth's interior and it is mainly this increased mantle potential temperature Tp that controls the dynamics of the crust and upper mantle and the predominant style of tectonics in the Early Earth. We conduct 3D petrological-magmatic-thermomechanical numerical modelling experiments of the crust and upper mantle under Early Earth conditions using a plume tectonics model setup. For varying crustal structures and a mantle potential temperature increase (ΔTp, compared to present day conditions), a hot lower thermal boundary layer introduces spontaneously developing mantle plumes and after repeated melt removal, depleted mantle lithosphere is formed self-consistently. New crust is produced in the form of both volcanics and plutonics. For an increase in mantle potential temperature ΔTp= 250 K, presumably corresponding to an Archean mantle, models show large amounts of subcrustal decompression melting and consequently large amounts of magmatism, which in turn influence the dynamics. In a first active phase (10-20 Ma) rising diapirs within the crust lead to the formation of the typical dome and keel pattern (e.g. Kaapvaal craton in South Africa, Pilbara craton in northwest Australia). A long passive phase follows with strong growth of crust and lithosphere. Both crust and lithosphere thickness are regulated by thermal-chemical instabilities assisted by lower crust eclogitization. Eclogitization depth is reached after ~80 Ma and linear or cylindrical drips originate at the crust or lithosphere bottom. Delamination of lower crust and lithosphere then occurs as one 'catastrophic' event within the next 20 Ma.
Early Earth tectonics: A high-resolution 3D numerical modelling approach
NASA Astrophysics Data System (ADS)
Fischer, R.; Gerya, T.
2014-12-01
Early Earth had a higher amount of remaining radiogenic elements as well as a higher amount of leftover primordial heat. Both contributed to the increased temperature in the Earth's interior and it is mainly this increased mantle potential temperature ΔTp that controls the dynamics of the crust and upper mantle and the style of Early Earth tectonics. For a minor increase in temperature ΔTp < 175 K a subduction-collision style ensues which is largely similar to present day plate tectonics. For a moderate increase in ΔTp = 175-250 K subduction can still occur, however plates are strongly weakened and buckling, delamination and Rayleigh-Taylor style dripping of the plate is observed in addition. For higher temperatures ΔTp > 250 K no subduction can be observed anymore and tectonics is dominated by delamination and Rayleigh-Taylor instabilities. We conduct 3D petrological-thermomechanical numerical modelling experiments of the crust and upper mantle under Early Earth conditions and a plume tectonics model setup. For varying crustal structures and an increased mantle potential temperature ΔTp, a thermal anomaly in the bottom temperature boundary introduces a plume. The model is able to self-sufficiently form depleted mantle lithosphere after repeated melt removal. New crust can be produced in the form of volcanics or plutonics. To simulate differentiation the newly formed crust can have a range in composition from basaltic over dacitic to granitic depending on its source rock. Models show large amounts of subcrustal decompression melting and consequently large amounts of new formed crust which in turn influences the dynamics. Mantle and crust are convecting separately. Dome-shaped plutons of mafic or felsic composition can be observed in the crust. Between these domes elongated belts of upper crust, volcanics and sediments are formed. These structures look similar to, for example, the Kaapvaal craton in South Africa where the elongated shape of the Barberton
Heterogeneous Differentiation of Human Mesenchymal Stem Cells in 3D Extracellular Matrix Composites
Jung, Jangwook P.; Bache-Wiig, Meredith K.; Provenzano, Paolo P.; Ogle, Brenda M.
2016-01-01
Abstract Extracellular matrix (ECM) proteins are structural elements of tissue and also potent signaling molecules. Previously, our laboratory showed that ECM of 2D coatings can trigger differentiation of bone marrow-derived mesenchymal stem cells (MSCs) into mesodermal lineages in an ECM-specific manner over 14 days, in some cases comparable to chemical induction. To test whether a similar effect was possible in a 3D, tissue-like environment, we designed a synthetic-natural biomaterial composite. The composite can present whole-molecule ECM proteins to cells, even those that do not spontaneously form hydrogels ex vivo, in 3D. To this end, we entrapped collagen type I, laminin-111, or fibronectin in ECM composites with MSCs and directly compared markers of mesodermal differentiation including cardiomyogenic (ACTC1), osteogenic (SPP1), adipogenic (PPARG), and chondrogenic (SOX9) in 2D versus 3D. We found the 3D condition largely mimicked the 2D condition such that the addition of type I collagen was the most potent inducer of differentiation to all lineages tested. One notable difference between 2D and 3D was pronounced adipogenic differentiation in 3D especially in the presence of exogenous collagen type I. In particular, PPARG gene expression was significantly increased ∼16-fold relative to chemical induction, in 3D and not in 2D. Unexpectedly, 3D engagement of ECM proteins also altered immunomodulatory function of MSCs in that expression of IL-6 gene was elevated relative to basal levels in 2D. In fact, levels of IL-6 gene expression in 3D composites containing exogenously supplied collagen type I or fibronectin were statistically similar to levels attained in 2D with tumor necrosis factor-α (TNF-α) stimulation and these levels were sustained over a 2-week period. Thus, this novel biomaterial platform allowed us to compare the biochemical impact of whole-molecule ECM proteins in 2D versus 3D indicating enhanced adipogenic differentiation and IL-6 expression
Heterogeneous Differentiation of Human Mesenchymal Stem Cells in 3D Extracellular Matrix Composites.
Jung, Jangwook P; Bache-Wiig, Meredith K; Provenzano, Paolo P; Ogle, Brenda M
2016-01-01
Extracellular matrix (ECM) proteins are structural elements of tissue and also potent signaling molecules. Previously, our laboratory showed that ECM of 2D coatings can trigger differentiation of bone marrow-derived mesenchymal stem cells (MSCs) into mesodermal lineages in an ECM-specific manner over 14 days, in some cases comparable to chemical induction. To test whether a similar effect was possible in a 3D, tissue-like environment, we designed a synthetic-natural biomaterial composite. The composite can present whole-molecule ECM proteins to cells, even those that do not spontaneously form hydrogels ex vivo, in 3D. To this end, we entrapped collagen type I, laminin-111, or fibronectin in ECM composites with MSCs and directly compared markers of mesodermal differentiation including cardiomyogenic (ACTC1), osteogenic (SPP1), adipogenic (PPARG), and chondrogenic (SOX9) in 2D versus 3D. We found the 3D condition largely mimicked the 2D condition such that the addition of type I collagen was the most potent inducer of differentiation to all lineages tested. One notable difference between 2D and 3D was pronounced adipogenic differentiation in 3D especially in the presence of exogenous collagen type I. In particular, PPARG gene expression was significantly increased ∼16-fold relative to chemical induction, in 3D and not in 2D. Unexpectedly, 3D engagement of ECM proteins also altered immunomodulatory function of MSCs in that expression of IL-6 gene was elevated relative to basal levels in 2D. In fact, levels of IL-6 gene expression in 3D composites containing exogenously supplied collagen type I or fibronectin were statistically similar to levels attained in 2D with tumor necrosis factor-α (TNF-α) stimulation and these levels were sustained over a 2-week period. Thus, this novel biomaterial platform allowed us to compare the biochemical impact of whole-molecule ECM proteins in 2D versus 3D indicating enhanced adipogenic differentiation and IL-6 expression of MSC in
3D Sun-to-Earth Solar Wind Modeling by SIP-CESE-MHD Model
NASA Astrophysics Data System (ADS)
Feng, Xueshang
2012-07-01
3D Sun-to-Earth Solar Wind Modeling by SIP-CESE-MHD Model Xueshang Feng, Xiang Changqing, Jiang Chaowei State Key Lab of Space Weather/CSSAR, CAS Beijing 100190 The objective is to present our solar-interplanetary space-time conservation element and solution element (CESE) model (SIP-CESE MHD model) (Feng, Zhou and Wu, Astrophys. J. 655, 1110, 2007; Feng et al., Astrophys. J. 723, 300, 2010; Feng et al., ApJ, 734, 50, 2011) with adaptive mesh refinement (AMR) implementation under six-component/Yin-Yang grid system. In this talk, we present the results of applying the SIP-AMR-CESE MHD model for modeling the coronal mass ejection and the solar wind background of different solar activity phases by comparison with SOHO observations and other spacecraft data from OMNI. Our numerical results show overall good agreements in the solar corona and in interplanetary space with these multiple spacecraft observations. From the modeler's experience, some limitations are addressed for this kind of initial-value boundary problems in Sun-to-Earth MHD modeling, and future out-of-ecliptic and in situ observations of the Sun and solar wind will be a solution to these limitations.
Early Earth plume-lid tectonics: A high-resolution 3D numerical modelling approach
NASA Astrophysics Data System (ADS)
Fischer, Ria; Gerya, Taras
2016-04-01
Early Earth had a higher amount of radiogenic elements as well as a higher amount of leftover primordial heat. Both contribute to the increased temperature in the Earth's interior and it is mainly this increased mantle potential temperature Tp that controls the dynamics of the crust and upper mantle and the predominant style of tectonics in the Archean Earth. We conduct 3D petrological-magmatic-thermomechanical numerical modelling experiments of the crust and upper mantle under Archean conditions using a plume-lid tectonics model setup. For varying crustal compositions and a mantle potential temperature increase ΔTp = 250K (compared to present day conditions), a hot lower thermal boundary layer introduces spontaneously developing mantle plumes and after repeated melt removal, depleted mantle lithosphere is formed self-consistently. New crust is produced in the form of both volcanic and plutonic magmatism. Models show large amounts of subcrustal decompression melting and production of new crust which in turn influences the dynamics. On short-term (10 ‑ 20Myr) rising diapirs and sinking basaltic crust lead to crustal overturn and to the formation of the typical Archean dome-and-keel pattern. On long-term a long (˜ 80Myr) passive 'growth phase' with strong growth of crust and lithosphere is observed. Both crust and lithosphere thickness are regulated by thermochemical instabilities assisted by lower crustal eclogitisation and a subcrustal small-scale convection area. Delamination of lower crust and lithosphere is initiated by linear or cylindrical eclogite drips and occurs as one 'catastrophic' event within a 20Myr 'removal phase'.
Heterogeneity and Flow in the Deep Earth
NASA Astrophysics Data System (ADS)
Cottaar, Sanne
Over the past half century, study of deep regions in the Earth has revealed them to be complex and dynamic. Much of our knowledge comes from seismological data, and ultimately these observations need to be linked to results from geodynamics and mineral physics in order to make inferences about compositional heterogeneities and flow. One example of a strongly-heterogeneous region is the lower thermal boundary layer of the mantle, i.e. the layer of several hundred kilometers thickness above the core-mantle boundary, commonly referred to as D″. This region appears to be characterized by two large provinces of distinctive slow shear velocities, 4000-5000 km across, one beneath the Pacific and one beneath Africa. Surrounding these regions, seismic velocities are faster, and often interpreted as corresponding to a graveyard of slabs. A recently discovered phase transition from perovskite to postperovskite may also occur in this depth range and has been associated with an intermittently observed seismic discontinuity at the top of D″. This study adds to the seismological evidence for complexities both in isotropic seismic velocities as well as in anisotropic velocities and how those can be linked to flow of material. We map a distinctive small "pile" of slow shear velocity beneath Russia through direct waveform evidence. This "pile" is less than 1000 km across, and thus much smaller than the Pacific and African provinces. Its height is several hundreds of kilometers and its velocity reduction suggests it is composed of the same material as the large provinces of slow shear velocity. Beneath Hawaii, at the northern edge of the Pacific province, we find an extended thin zone of ultra-low velocities. This is the first time the three-dimensional extent of such a zone is constrained with some accuracy. The constraints on its morphology come from the presence of strong postcursors to shear waves diffracted along the core-mantle boundary, delayed by 30 to 50 seconds with
3-D density models within an ellipsoidal-Earth from inversion of geoid anomalies
NASA Astrophysics Data System (ADS)
Chaves, C. M.; Ussami, N.
2013-12-01
Modeling density perturbations is very important to understand geodynamic processes which occur within the Earth's mantle. Commonly, the Earth's density is predicted by converting a velocity model into a density model using either a constant scaling factor or a relationship provided by mineral physics. Nonetheless, several factors such as temperature, composition and melting can affect the wave propagation speed so that a seismically converted density model may not retrieve the actual density distribution. This limitation may hamper the modeling the geodynamic processes. Due to advances in satellite-derived gravity data acquisition (e.g. GRACE, GOCE), the gravity field is now obtained with an unprecedented accuracy and resolution allowing us to estimate more uniformly the 3-D density distribution for the whole Earth. Here we present a computational algorithm to invert geoid anomalies in order to estimate density variations in the mantle. Using an ellipsoidal-Earth approximation, the model space is represented by a set of tesseroids. From a synthetic geoid anomaly caused by a plume tail ascending through the mantle with Gaussian noise added, the inversion code is capable to recover with good accuracy the density contrast and the body geometry when compared to the synthetic model. This algorithm was also tested in a natural case study, where geoid anomalies from the Yellowstone Province (YP) were inverted. The estimated density model (EDM) has a predominantly negative density contrast (~ -50 kg/m3) relative to the surrounding upper mantle and extends to the depth of 1000 km. The EDM exhibits an anti-correlation of up to -0.7 with one of the most recent S-velocity model for the western United States. The predicted dynamic topography from the EDM explains almost 80 % of the observed dynamic topography in the YP. From our results, we conclude that a joint-interpretation of density anomalies derived from geoid and velocity perturbations from seismic tomography models
NASA Astrophysics Data System (ADS)
El Said, Bassam; Ivanov, Dmitry; Long, Andrew C.; Hallett, Stephen R.
2016-03-01
3D composite materials are characterized by complex internal yarn architectures, leading to complex deformation and failure development mechanisms. Net-shaped preforms, which are originally periodic in nature, lose their periodicity when the fabric is draped, deformed on a tool, and consolidated to create geometrically complex composite components. As a result, the internal yarn architecture, which dominates the mechanical behaviour, becomes dependent on the structural geometry. Hence, predicting the mechanical behaviour of 3D composites requires an accurate representation of the yarn architecture within structural scale models. When applied to 3D composites, conventional finite element modelling techniques are limited to either homogenised properties at the structural scale, or the unit cell scale for a more detailed material property definition. Consequently, these models fail to capture the complex phenomena occurring across multiple length scales and their effects on a 3D composite's mechanical response. Here a multi-scale modelling approach based on a 3D spatial Voronoi tessellation is proposed. The model creates an intermediate length scale suitable for homogenisation to deal with the non-periodic nature of the final material. Information is passed between the different length scales to allow for the effect of the structural geometry to be taken into account on the smaller scales. The stiffness and surface strain predictions from the proposed model have been found to be in good agreement with experimental results. The proposed modelling framework has been used to gain important insight into the behaviour of this category of materials. It has been observed that the strain and stress distributions are strongly dependent on the internal yarn architecture and consequently on the final component geometry. Even for simple coupon tests, the internal architecture and geometric effects dominate the mechanical response. Consequently, the behaviour of 3D woven
Microfluidic Bioprinting of Heterogeneous 3D Tissue Constructs Using Low-Viscosity Bioink.
Colosi, Cristina; Shin, Su Ryon; Manoharan, Vijayan; Massa, Solange; Costantini, Marco; Barbetta, Andrea; Dokmeci, Mehmet Remzi; Dentini, Mariella; Khademhosseini, Ali
2016-01-27
A novel bioink and a dispensing technique for 3D tissue-engineering applications are presented. The technique incorporates a coaxial extrusion needle using a low-viscosity cell-laden bioink to produce highly defined 3D biostructures. The extrusion system is then coupled to a microfluidic device to control the bioink arrangement deposition, demonstrating the versatility of the bioprinting technique. This low-viscosity cell-responsive bioink promotes cell migration and alignment within each fiber organizing the encapsulated cells. PMID:26606883
Rapid probabilistic source characterisation in 3D earth models using learning algorithms
NASA Astrophysics Data System (ADS)
Valentine, A. P.; Kaeufl, P.; Trampert, J.
2015-12-01
triggers; if longer windows are used, constraints on focal mechanism can also be obtained. We demonstrate that the use of 3D wave propagation allows results to be constrained better than is possible when only 1D earth models are used.
NASA Astrophysics Data System (ADS)
Begnaud, M. L.; Ballard, S.; Young, C. J.; Hipp, J. R.; Encarnacao, A.; Phillips, W. S.; Chael, E. P.; Rowe, C. A.
2012-12-01
We are developing a global 3D P wave velocity model of the Earth's crust and mantle using seismic tomography to assess improvement to seismic event locations obtained using high quality 3D Earth models in lieu of 1D and 2/2.5D models. We present the most recent version of SALSA3D (SAndia LoS Alamos 3D) version 1.9, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth (GT) events. Our model is derived from the latest version of the GT catalog of P/Pn travel-time picks assembled by Los Alamos National Laboratory. For this current version, we employ more robust data quality control measures than previously used, as well as additional global GT data sources. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays into representative rays. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions. For our starting model, we use a simplified layer crustal model derived from the NNSA Unified model in Eurasia and Crust 2.0 model everywhere else, overlying a uniform ak135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only in areas where the data warrant such a refinement. In previous versions, we based this refinement on velocity changes from previous model iterations. For the current version, we utilize the diagonal of the model resolution matrix to control where grid refinement occurs, resulting in more consistent and continuous areas of refinement than before. In addition to the changes in grid refinement, we also employ a more robust convergence criterion between successive grid refinements, allowing a better fit to first broader
Early Earth tectonics: A high-resolution 3D numerical modelling approach
NASA Astrophysics Data System (ADS)
Fischer, Ria; Gerya, Taras
2015-04-01
Early Earth had a higher amount of remaining radiogenic elements as well as a higher amount of leftover primordial heat. Both contributed to the increased temperature in the Earth's interior and it is mainly this increased mantle potential temperature ΔTp that controls the dynamics of the crust and upper mantle and the style of Early Earth tectonics. We conduct 3D petrological-thermomechanical numerical modelling experiments of the crust and upper mantle under Early Earth conditions using a plume tectonics model setup. For varying crustal structures and an increased mantle potential temperature ΔTp, a hot lower thermal boundary layer is used to introduce spontaneously developing mantle plumes. The model is able to self-sufficiently form depleted mantle lithosphere after repeated melt removal. New crust can be produced in the form of volcanics and/or plutonics. To simulate differentiation the newly formed crust can have a range in composition from basaltic to granitic depending on its source rock. For a major increase in the mantle temperature, presumably corresponding to an Archean mantle (ΔTp = 200 - 300K compared to present day conditions), models show large amounts of subcrustal decompression melting and consequently large amounts of volcanics, which in turn influence the dynamics. Mantle and crust are convecting separately. Dome-shaped felsic plutons can be observed in the crust. Between these domes elongated belts of downwelling basalt and sediments are formed. Both crust and lithosphere thickness are regulated by thermo-chemical instabilities assisted by lower crust eclogitization: linear or cylindrical drips originating at the crust or lithosphere bottom or delamination of lower crust or lithosphere. Very similar examples of dome and belt structures are still preserved in Archean cratons. One example is the Kaapvaal craton is South Africa where the elongated shape of the Barberton Greenstone Belt, mainly built from mafic rocks and sediments, is surrounded
A statistical measure of tissue heterogeneity with application to 3D PET sarcoma data.
O'Sullivan, Finbarr; Roy, Supratik; Eary, Janet
2003-07-01
In vivo measurement of local tissue characteristics by modern bioimaging techniques such as positron emission tomography (PET) provides the opportunity to analyze quantitatively the role that tissue heterogeneity may play in understanding biological function. This paper develops a statistical measure of the heterogeneity of a tissue characteristic that is based on the deviation of the distribution of the tissue characteristic from a unimodal elliptically contoured spatial pattern. An efficient algorithm is developed for computation of the measure based on volumetric region of interest data. The technique is illustrated by application to data from PET imaging studies of fluorodeoxyglucose utilization in human sarcomas. A set of 74 sarcoma patients (with five-year follow-up survival information) were evaluated for heterogeneity as well as a number of other potential prognostic indicators of survival. A Cox proportional hazards analysis of these data shows that the degree of heterogeneity of the sarcoma is the major risk factor associated with patient death. Some theory is developed to analyze the asymptotic statistical behavior of the heterogeneity estimator. In the context of data arising from Poisson deconvolution (PET being the prime example), the heterogeneity estimator, which is a non-linear functional of the PET image data, is consistent and converges at a rate that is parametric in the injected dose. PMID:12925510
NASA Astrophysics Data System (ADS)
Petersson, Anders; Rodgers, Arthur
2010-05-01
The finite difference method on a uniform Cartesian grid is a highly efficient and easy to implement technique for solving the elastic wave equation in seismic applications. However, the spacing in a uniform Cartesian grid is fixed throughout the computational domain, whereas the resolution requirements in realistic seismic simulations usually are higher near the surface than at depth. This can be seen from the well-known formula h ≤ L-P which relates the grid spacing h to the wave length L, and the required number of grid points per wavelength P for obtaining an accurate solution. The compressional and shear wave lengths in the earth generally increase with depth and are often a factor of ten larger below the Moho discontinuity (at about 30 km depth), than in sedimentary basins near the surface. A uniform grid must have a grid spacing based on the small wave lengths near the surface, which results in over-resolving the solution at depth. As a result, the number of points in a uniform grid is unnecessarily large. In the wave propagation project (WPP) code, we address the over-resolution-at-depth issue by generalizing our previously developed single grid finite difference scheme to work on a composite grid consisting of a set of structured rectangular grids of different spacings, with hanging nodes on the grid refinement interfaces. The computational domain in a regional seismic simulation often extends to depth 40-50 km. Hence, using a refinement ratio of two, we need about three grid refinements from the bottom of the computational domain to the surface, to keep the local grid size in approximate parity with the local wave lengths. The challenge of the composite grid approach is to find a stable and accurate method for coupling the solution across the grid refinement interface. Of particular importance is the treatment of the solution at the hanging nodes, i.e., the fine grid points which are located in between coarse grid points. WPP implements a new, energy
DynEarthSol3D: numerical studies of basal crevasses and calving blocks
NASA Astrophysics Data System (ADS)
Logan, E.; Lavier, L. L.; Choi, E.; Tan, E.; Catania, G. A.
2014-12-01
DynEarthSol3D (DES) is a thermomechanical model for the simulation of dynamic ice flow. We present the application of DES toward two case studies - basal crevasses and calving blocks - to illustrate the potential of the model to aid in understanding calving processes. Among the advantages of using DES are: its unstructured meshes which adaptively resolve zones of high interest; its use of multiple rheologies to simulate different types of dynamic behavior; and its explicit and parallel numerical core which both make the implementation of different boundary conditions easy and the model highly scalable. We examine the initiation and development of both basal crevasses and calving blocks through time using visco-elasto-plastic rheology. Employing a brittle-to-ductile transition zone (BDTZ) based on local strain rate shows that the style and development of brittle features like crevasses differs markedly on the rheological parameters. Brittle and ductile behavior are captured by Mohr-Coulomb elastoplasticity and Maxwell viscoelasticity, respectively. We explore the parameter spaces which define these rheologies (including temperature) as well as the BDTZ threshold (shown in the literature as 10-7 Pa s), using time-to-failure as a metric for accuracy within the model. As the time it takes for a block of ice to fail can determine an iceberg's size, this work has implications for calving laws.
Estimation of Hydraulic Fracturing in the Earth Fill Dam by 3-D Analysis
NASA Astrophysics Data System (ADS)
Nishimura, Shin-Ichi
It is necessary to calculate strength and strain for estimation of hydraulic fracturing in the earth fill dam, and to which the FEM is effective. 2-D analysis can produce good results to some extent if an embankment is linear and the plain strain condition can be set to the cross section. However, there may be some conditions not possible to express in the 2-D plain because the actual embankment of agricultural reservoirs is formed by straight and curved lines. Moreover, it may not be possible to precisely calculate strain in the direction of dam axis because the 2-D analysis in the cross section cannot take the shape in the vertical section into consideration. Therefore, we performed 3-D built up analysis targeting the actually-leaked agricultural reservoir to examine hazards of hydraulic fracturing based on the shape of an embankment and by rapid impoundment of water. It resulted in the occurrence of hydraulic fracturing to develop by water pressure due to the vertical cracks caused by tensile strain in the valley and refractive section of the foundation.
Development of hybrid 3-D hydrological modeling for the NCAR Community Earth System Model (CESM)
Zeng, Xubin; Troch, Peter; Pelletier, Jon; Niu, Guo-Yue; Gochis, David
2015-11-15
This is the Final Report of our four-year (3-year plus one-year no cost extension) collaborative project between the University of Arizona (UA) and the National Center for Atmospheric Research (NCAR). The overall objective of our project is to develop and evaluate the first hybrid 3-D hydrological model with a horizontal grid spacing of 1 km for the NCAR Community Earth System Model (CESM). We have made substantial progress in model development and evaluation, computational efficiencies and software engineering, and data development and evaluation, as discussed in Sections 2-4. Section 5 presents our success in data dissemination, while Section 6 discusses the scientific impacts of our work. Section 7 discusses education and mentoring success of our project, while Section 8 lists our relevant DOE services. All peer-reviewed papers that acknowledged this project are listed in Section 9. Highlights of our achievements include: • We have finished 20 papers (most published already) on model development and evaluation, computational efficiencies and software engineering, and data development and evaluation • The global datasets developed under this project have been permanently archived and publicly available • Some of our research results have already been implemented in WRF and CLM • Patrick Broxton and Michael Brunke have received their Ph.D. • PI Zeng has served on DOE proposal review panels and DOE lab scientific focus area (SFA) review panels
NASA Astrophysics Data System (ADS)
Begnaud, M. L.; Ballard, S.; Young, C. J.; Hipp, J. R.; Chang, M.; Encarnacao, A.; Rowe, C. A.; Phillips, W. S.; Steck, L.
2011-12-01
To test the hypothesis that high quality 3D Earth models will produce seismic event locations that are more accurate and more precise than currently used 1D and 2/2.5D models, we are developing a global 3D P wave velocity model of the Earth's crust and mantle using seismic tomography. In this paper, we present the most recent version of our model, SALSA3D (SAndia LoS Alamos 3D) version 1.7, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth (GT) events, compared to existing models and/or systems. Our model is derived from the latest version of the GT catalog of P and Pn travel time picks assembled by Los Alamos National Laboratory. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays to produce representative rays. Reduction in the total number of ray paths is ~50%. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions. For our starting model, we use a simplified layer crustal model derived from the NNSA Unified model in Eurasia and Crust 2.0 model elsewhere, over a uniform ak135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only in areas where the data warrant it. In previous versions of SALSA3D, we based this refinement on velocity changes from previous model iterations. For version 1.7, we utilize the diagonal of the model resolution matrix to control where grid refinement occurs, resulting in more consistent and continuous areas of refinement than before. Our approach produces a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. We
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
NASA Astrophysics Data System (ADS)
Zhang, Mian; Huang, Cheng-li
2012-08-01
Generalized spherical harmonics (GSH) are usually applied on the problems where the Earth model is elliptical and elastic stress tensor is involved in, as stress tensor can’t be represented in vector spherical harmonics. However, the divergence of the te ns or and a vector dot - product with the tensor are only needed on computation rotation modes of the Earth which can be written in the vector spherical harmonics. We extend the equations on the spherical Earth to asymmetric 3D model by means of linear operator method. This method doesn’t use the complicated generalized spherical harmonics nor Wigner 3 - j symbol. As a validation of this method, the practical calculation of rotational modes of 3D Earth will be made and discussed.
Zhang, Yan; Guo, Zhen; Xie, Shuang; Li, Hui-Li; Zhu, Wen-Hua; Liu, Li; Dong, Xun-Qing; He, Wei-Xun; Ren, Jin-Chao; Liu, Ling-Zhi; Powell, Annie K
2015-11-01
Three isostructural cyano-bridged 3d-4f compounds, [YFe(CN)6(hep)2(H2O)4] (1), [DyFe(CN)6(hep)2(H2O)4] (2), and [DyCo(CN)6(hep)2(H2O)4] (3), were successfully assembled by site-targeted substitution of the 3d or rare-earth ions. All compounds have been structurally characterized to display slightly distorted pentagonal-bipyramidal local coordination geometry around the rare-earth ions. Magnetic analyses revealed negligible magnetic coupling in compound 1, antiferromagnetic intradimer interaction in 2, and weak ferromagnetic coupling through dipolar-dipolar interaction in 3. Under an applied direct-current (dc) field, 1 (Hdc = 2.5 kOe, τ0 = 1.3 × 10(-7) s, and Ueff/kB = 23 K) and 3 (Hdc = 2.0 kOe, τ0 = 7.1 × 10(-11) s, and Ueff/kB = 63 K) respectively indicated magnetic relaxation behavior based on a single [Fe(III)]LS ion and a Dy(III) ion; nevertheless, 2 (Hdc = 2.0 kOe, τ0 = 9.7 × 10(-8) s, and Ueff/kB = 23 K) appeared to be a single-molecule magnet based on a cyano-bridged DyFe dimer. Compound 1, which can be regarded as a single-ion magnet of the [Fe(III)]LS ion linked to a diamagnetic Y(III) ion in a cyano-bridged heterodimer, represents one of the rarely investigated examples based on a single Fe(III) ion explored in magnetic relaxation behavior. It demonstrated that the introduction of intradimer magnetic interaction of 2 through a cyano bridge between Dy(III) and [Fe(III)]LS ions negatively affects the energy barrier and χ″(T) peak temperature compared to 3. PMID:26473654
Visualising, segmenting and analysing heterogenous glacigenic sediments using 3D x-ray CT.
NASA Astrophysics Data System (ADS)
Carr, Simon; Diggens, Lucy; Groves, John; O'Sullivan, Catherine; Marsland, Rhona
2015-04-01
Whilst there has been significant application of 3D x-ray CT to geological contexts, much of this work has focused on examining properties such as porosity, which are important in reservoir assessment and hydrological evaluations. There has been considerably less attention given to the analysis of the properties of sediments themselves. One particular challenge in CT analysis is to effectively observe and discriminate the relationships between the skeleton and matrix of a sediment. This is particularly challenging in glacial sediments, which comprise an admixture of particles of a wide range of size, morphology and composition within a variably-consolidated sediment body. A key sedimentological component of glacial sediments is their fabric properties. Till fabric data has long been applied to the analysis of the coupling between glaciers and their deformable substrates. This work has typically focused on identifying former ice-flow directions, processes of till deformation and emplacement, and such data is often used to reconcile the sedimentary evidence of former glaciation with the predicted glacier and ice-sheet dynamics derived from numerical models. The collection and interpretation of till fabric data has received significant criticism in recent years, with issues such as low sample populations (typically ~50 grains per sample), small-scale spatial variation in till fabric and operator bias during data collection, all of which compromise the reliability of macro-scale till fabric analysis. Recent studies of micro-scale till fabrics have substantially added to our understanding, and suggest there is systematic variation in particle fabric as a function of particle size. However, these findings are compromised by the 2D nature of the samples (derived from thin sections) capturing only apparent orientations of particles, and are again limited to relatively small datasets. As such, there are fundamental limitations in the quality and application of till fabric
A 3D pencil-beam-based superposition algorithm for photon dose calculation in heterogeneous media
NASA Astrophysics Data System (ADS)
Tillikainen, L.; Helminen, H.; Torsti, T.; Siljamäki, S.; Alakuijala, J.; Pyyry, J.; Ulmer, W.
2008-07-01
In this work, a novel three-dimensional superposition algorithm for photon dose calculation is presented. The dose calculation is performed as a superposition of pencil beams, which are modified based on tissue electron densities. The pencil beams have been derived from Monte Carlo simulations, and are separated into lateral and depth-directed components. The lateral component is modeled using exponential functions, which allows accurate modeling of lateral scatter in heterogeneous tissues. The depth-directed component represents the total energy deposited on each plane, which is spread out using the lateral scatter functions. Finally, convolution in the depth direction is applied to account for tissue interface effects. The method can be used with the previously introduced multiple-source model for clinical settings. The method was compared against Monte Carlo simulations in several phantoms including lung- and bone-type heterogeneities. Comparisons were made for several field sizes for 6 and 18 MV energies. The deviations were generally within (2%, 2 mm) of the field central axis dmax. Significantly larger deviations (up to 8%) were found only for the smallest field in the lung slab phantom for 18 MV. The presented method was found to be accurate in a wide range of conditions making it suitable for clinical planning purposes.
A New Global Model for 3-D variations in P Wave Speed in Earth's Mantle
NASA Astrophysics Data System (ADS)
Karason, H.; van der Hilst, R. D.; Li, C.
2003-12-01
In an effort to improve the resolution of mantle structure we have combined complementary data sets of short- and long period (absolute and differential) travel time residuals. Our new model is based on short period P (N\\~7.7x10**6), pP (N\\~2.3x10**5), and PKP (N\\~16x10**4) data from the catalog by Engdahl et al (BSSA, 1998), short-period PKP differential times (N\\~1600) measured by McSweeney & Creager, and long-period differential PP-P times - N\\~20,000 measured by Bolton & Masters and N\\~18,000 by Ritsema - and Pdiff-PKP (N\\~560) measured by Wysession. Inversion tests, spectral analysis, and comparison with geology indicate that the large-scale upper mantle structure is better constrained with the addition of PP-P, whereas the Pdiff and PKP data help constrain deep mantle structure (Karason & Van der Hilst, JGR, 2001). The long period data were measured by cross-correlation. We solved the system of equations using 400 iterations of the iterative algorithm LSQR For the short period (1 Hz) data we use a high frequency approximation and trace rays through a fine grid of constant slowness cells to invert for mantle structure. For low frequency Pdiff and PP data we account for sensitivity to structure away from the optical ray path with 3-D Frechet derivatives (sensitivity kernels) estimated from single forward scattering and projected onto basis functions (constant slowness blocks) used for model parameterization. With such kernels the low frequency data can constrain long wavelength heterogeneity without keeping the short period data from mapping details in densely sampled regions. In addition to finite frequency sensitivity kernels we optimized the localization by using a parameterization that adapts to spatial resolution, with small cells in regions of dense sampling and larger cells in regions where sampling is more sparse (the total number of cells was \\~ 350,000). Finally, we corrected all travel times and surface reflections for lateral variations in
NASA Astrophysics Data System (ADS)
Prokein, P.; Nolan, M.
2004-12-01
In spring of 2004 we received a Small Grant for Exploratory Research (SGER) from the NSF's Office of Polar Program's Arctic Logistics and Research Support program to create an internet-based, interactive, 3D terrain and data visualization system of the Arctic. A preliminary version of this application, called EarthSLOT, can now be found on-line at www.earthslot.org. EarthSLOT allows users to visualize the earth, whether as a spinning globe from space or from the sea-floor looking up at a mid-ocean ridge or anywhere in between. Flight controls range from completely interactive to following pre-planned routes, and the visualizations can be done real-time over the internet or saved as screen shots or MPG movies. The purpose of this project is to put easy-to-use 3D terrain and visualization tools into the hands of many users at little or no cost to them, by taking care of the complicated and expensive work ourselves. Therefore EarthSLOT may be an ideal tool for scientists with low outreach budgets to share their research with other scientists or the public. Those on a very low budget can use EarthSLOT for free, as can any ordinary user, without modifying it or adding their own data. Example uses would be analyzing the terrain surrounding a field site, adding a 3D visualization of a study area to a presentation, or exploring the vector data added by others to their study areas. Those with about \\$1200 to spend on the necessary commercial software can add their own content to the existing application. For example, an ecologist can add or create shapefiles that outlines their study plots and then link those outlines to web pages containing data or further information. Or a glaciologist can superimpose the locations of mass balance stakes and weather stations on the surface of a glacier, then create an mpg movie that starts in space and zooms down to the stake level to visualize how weather systems on a planetary scale relate to the local measurements. Or scientists
Performance Modeling for 3D Visualization in a Heterogeneous Computing Environment
Bowman, Ian; Shalf, John; Ma, Kwan-Liu; Bethel, Wes
2004-06-30
The visualization of large, remotely located data sets necessitates the development of a distributed computing pipeline in order to reduce the data, in stages, to a manageable size. The required baseline infrastructure for launching such a distributed pipeline is becoming available, but few services support even marginally optimal resource selection and partitioning of the data analysis workflow. We explore a methodology for building a model of overall application performance using a composition of the analytic models of individual components that comprise the pipeline. The analytic models are shown to be accurate on a testbed of distributed heterogeneous systems. The prediction methodology will form the foundation of a more robust resource management service for future Grid-based visualization applications.
NASA Astrophysics Data System (ADS)
Sidler, Rolf; Carcione, José M.; Holliger, Klaus
2014-02-01
We present a novel approach for the comprehensive, flexible and accurate simulation of poroelastic wave propagation in 3-D cylindrical coordinates. An important application of this method is the realistic modelling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, as of yet largely unresolved, problem in exploration geophysics. To this end, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the centre followed by the mudcake and/or casing, and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction and Fourier expansions in the vertical and azimuthal directions as well as a Runge-Kutta integration scheme for the time evolution. Trigonometric interpolation and a domain decomposition method based on the method of characteristics are used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces as well as to reduce the number of gridpoints in the innermost domain for computational efficiency. We apply this novel modelling approach to the particularly challenging scenario of near-surface borehole environments. To this end, we compare 3-D heterogeneous and corresponding rotationally invariant simulations, assess the sensitivity of Stoneley waves to formation permeability in the presence of a casing and evaluate the effects of an excavation damage zone behind a casing on sonic log recordings. Our results indicate that only first arrival times of fast modes are reasonably well described by rotationally invariant approximations of 3-D heterogenous media. We also find that Stoneley waves are indeed remarkably sensitive to the average permeability behind a perforated PVC casing, and that the presence of an excavation damage zone behind a casing tends to dominate the overall signature of recorded seismograms.
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-28
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
NASA Astrophysics Data System (ADS)
Luquet, David; Marchiano, Régis; Coulouvrat, François
2015-10-01
Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D
A 3D Cloud-Construction Algorithm for the EarthCARE Satellite Mission
NASA Technical Reports Server (NTRS)
Barker, H. W.; Jerg, M. P.; Wehr, T.; Kato, S.; Donovan, D. P.; Hogan, R. J.
2011-01-01
This article presents and assesses an algorithm that constructs 3D distributions of cloud from passive satellite imagery and collocated 2D nadir profiles of cloud properties inferred synergistically from lidar, cloud radar and imager data.
SALSA3D - A Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location
NASA Astrophysics Data System (ADS)
Ballard, S.; Begnaud, M. L.; Young, C. J.; Hipp, J. R.; Chang, M.; Encarnacao, A. V.; Rowe, C. A.; Phillips, W. S.; Steck, L.
2010-12-01
To test the hypothesis that high quality 3D Earth models will produce seismic event locations which are more accurate and more precise, we are developing a global 3D P wave velocity model of the Earth’s crust and mantle using seismic tomography. In this paper, we present the most recent version of our model, SALSA3D version 1.5, and demonstrate its ability to reduce mislocations for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. Our model is derived from the latest version of the Ground Truth (GT) catalog of P and Pn travel time picks assembled by Los Alamos National Laboratory. To prevent over-weighting due to ray path redundancy and to reduce the computational burden, we cluster rays to produce representative rays. Reduction in the total number of ray paths is ~50%. The model is represented using the triangular tessellation system described by Ballard et al. (2009), which incorporates variable resolution in both the geographic and radial dimensions.. For our starting model, we use a simplified two layer crustal model derived from the Crust 2.0 model over a uniform AK135 mantle. Sufficient damping is used to reduce velocity adjustments so that ray path changes between iterations are small. We obtain proper model smoothness by using progressive grid refinement, refining the grid only around areas with significant velocity changes from the starting model. At each grid refinement level except the last one we limit the number of iterations to prevent convergence thereby preserving aspects of broad features resolved at coarser resolutions. Our approach produces a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. This scheme is computationally expensive, so we use a distributed computing framework based on the Java Parallel Processing Framework, providing us with ~400 processors. Resolution of our model is assessed using a
Curtis, Matthew W.; Sharma, Sadhana; Desai, Tejal A.
2011-01-01
Cardiac myocytes are known to be influenced by the rigidity and topography of their physical microenvironment. It was hypothesized that 3D heterogeneity introduced by purely physical microdomains regulates cardiac myocyte size and contraction. This was tested in vitro using polymeric microstructures (G′=1.66 GPa) suspended with random orientation in 3D by a soft Matrigel matrix (G′=22.9 Pa). After 10 days of culture, the presence of 100 μm-long microstructures in 3D gels induced fold increases in neonatal rat ventricular myocyte size (1.61±0.06, p<0.01) and total protein/cell ratios (1.43± 0.08, p<0.05) that were comparable to those induced chemically by 50 μM phenylephrine treatment. Upon attachment to microstructures, individual myocytes also had larger cross-sectional areas (1.57±0.05, p<0.01) and higher average rates of spontaneous contraction (2.01±0.08, p<0.01) than unattached myocytes. Furthermore, the inclusion of microstructures in myocyte-seeded gels caused significant increases in the expression of beta-1 adrenergic receptor (β1-AR, 1.19±0.01), cardiac ankyrin repeat protein (CARP, 1.26±0.02), and sarcoplasmic reticulum calcium-ATPase (SERCA2, 1.59±0.12, p<0.05), genes implicated in hypertrophy and contractile activity. Together, the results demonstrate that cardiac myocyte behavior can be controlled through local 3D microdomains alone. This approach of defining physical cues as independent features may help to advance the elemental design considerations for scaffolds in cardiac tissue engineering and therapeutic microdevices. PMID:20668947
Hwang, In Hong; Bae, Jeong Mi; Kim, Wan-Seok; Jo, Young Dan; Kim, Cheal; Kim, Youngmee; Kim, Sung-Jin; Huh, Seong
2012-11-01
We report bifunctional three-dimensional (3D) Cu-MOFs with high selectivity of CO(2) over N(2) and H(2) as well as high catalytic activity for transesterification of esters. The Cu-MOFs containing Cu(2) dinuclear units connected by glutarates and bipyridyl ligands are formulated as [{Cu(2)(Glu)(2)(μ-bpa)}·(CH(3)CN)](n) (1) and [{Cu(2)(Glu)(2)(μ-bpp)}·(C(3)H(6)O)](n) (2) (Glu = glutarate, bpa = 1,2-bis(4-pyridyl)ethane, bpp = 1,3-bis(4-pyridyl)propane). These two new bifunctional 3D Cu-MOFs possess very similar pore shape with different pore dimensions. Their gas sorption behaviors were investigated by using CO(2), N(2) and H(2) at suitable temperatures. Both MOFs exhibited good CO(2) selectivity over N(2) and H(2). MOF 1 having a smaller pore dimension exhibited much higher CO(2) adsorption enthalpy than MOF 2 having a larger pore dimension. However, MOF 2 exhibited more enhanced CO(2) uptake ability than MOF 1. A subtle variation of pore dimension indeed influenced the CO(2) uptake ability somewhat significantly especially at higher temperatures such as 273 K and 298 K. Heterogeneous catalytic activities of the MOFs were also investigated in detail. Only MOF 1 appeared to be an efficient, mild, and easily recyclable heterogeneous catalyst for the transesterification of esters and constitutes a promising class of heterogeneous catalysts that allowed reuse without a significant loss of activity through twenty runs with ester. PMID:22968940
Assessing the habitability of planets with Earth-like atmospheres with 1D and 3D climate modeling
NASA Astrophysics Data System (ADS)
Godolt, M.; Grenfell, J. L.; Kitzmann, D.; Kunze, M.; Langematz, U.; Patzer, A. B. C.; Rauer, H.; Stracke, B.
2016-07-01
Context. The habitable zone (HZ) describes the range of orbital distances around a star where the existence of liquid water on the surface of an Earth-like planet is in principle possible. The applicability of one-dimensional (1D) climate models for the estimation of the HZ boundaries has been questioned by recent three-dimensional (3D) climate studies. While 3D studies can calculate the water vapor, ice albedo, and cloud feedback self-consistently and therefore allow for a deeper understanding and the identification of relevant climate processes, 1D model studies rely on fewer model assumptions and can be more easily applied to the large parameter space possible for extrasolar planets. Aims: We evaluate the applicability of 1D climate models to estimate the potential habitability of Earth-like extrasolar planets by comparing our 1D model results to those of 3D climate studies in the literature. We vary the two important planetary properties, surface albedo and relative humidity, in the 1D model. These depend on climate feedbacks that are not treated self-consistently in most 1D models. Methods: We applied a cloud-free 1D radiative-convective climate model to calculate the climate of Earth-like planets around different types of main-sequence stars with varying surface albedo and relative humidity profile. We compared the results to those of 3D model calculations available in the literature and investigated to what extent the 1D model can approximate the surface temperatures calculated by the 3D models. Results: The 1D parameter study results in a large range of climates possible for an Earth-sized planet with an Earth-like atmosphere and water reservoir at a certain stellar insolation. At some stellar insolations the full spectrum of climate states could be realized, i.e., uninhabitable conditions due to surface temperatures that are too high or too low as well as habitable surface conditions, depending only on the relative humidity and surface albedo assumed. When
NASA Astrophysics Data System (ADS)
Huber, E.; Huggenberger, P.
2012-12-01
Accurate predictions on groundwater flow and transport behavior within fluvial and glaciofluvial sediments, but also interaction with surface water bodies, rely on knowledge of distributed aquifer properties. The complexity of the depositional and erosional processes in fluvial systems leads to highly heterogeneous distributions of hydrogeological parameters. The system dynamics, such as aggradation rates and channel mobility of alluvial systems; its influence on the preservation potential of the key depositional elements in the geological record; and its influence on the heterogeneity scales and the relevance for groundwater hydraulics is topic of the presentation. The aims of our work are to find a relation between surface morphological structures and the sedimentary structures in vertical profiles (i.e. gravel pits or GPR sections) and to derive rules for the interpretation of horizontal time-slices from 3D GPR data. Based on these data we set-up conceptual models of the structures of coarse alluvial systems at different scales which can be tested by stochastic methods. Relevant depositional elements and a hierarchy or genetic relationship of such elements will be defined based on the knowledge of depositional processes in alluvial systems inferred from: field observations after major flood events; 2D and 3D GPR data; and from existing data derived from laboratory flumes. Extensive geophysical field experiments within the Tagliamento alluvial system gave new insights to the sedimentary structures developing at high flows. Owing to the fact that rivers often destroy at least part of their bed during or shortly after large floods and subsequently rebuild, it is not easy to establish a simple relationship between surface morphology and the sedimentary structures found in vertical sections of many alluvial outcrops. According to these findings we suppose that surface or near-surface structures will not catch the essence of heterogeneity of alluvial aquifers
NASA Astrophysics Data System (ADS)
Meheust, Y.; De Dreuzy, J.; Pichot, G.
2011-12-01
Flow channeling and permeability scaling in fractured media have been classically addressed either at the fracture- or at the network- scales. In the latter case they are linked to the topological structure of the network, while at the fracture scale they are controlled by the variability of the local aperture distribution inside individual fractures. In this study we analyze these two combined effects, investigating how flow localization below the scale of individual fractures influences that at the network scale and the resulting medium permeability. This is done by use of a new highly-resolved 3D discrete fracture network model (DFN). The local apertures of individual fractures are distributed according to a truncated Gaussian law, and exhibit self-affine spatial correlations that are bounded by an upper cutoff scale Lc; Lc and the fracture closure, defined as the ratio of the aperture fluctuations at scale Lc to the mean aperture, are considered homogeneous over the DFN. The network topology is controlled by a homogeneous scalar fracture density and a power law fracture length distribution. We have varied these features to investigate a large variety of DFN topologies, from sparse networks with varying degrees of fracture interconnections, flow bottlenecks and dead-ends (Fig. 1a), to dense well-connected networks (Fig. 1b). We have also investigated a large range of fracture closures, performing extensive simulations of about 105 different DFN realizations. At the fracture scale, accounting for local aperture fluctuations leads to a monotical deviation (which can exceed 50%) of the equivalent fracture transmissivity from the parallel plate behavior. At the network scale we observe a complex interaction between flow channeling within fracture planes and flow localization in the network. This interaction is controlled by the location of fracture interactions with respect to that of low local transmissivity zones (particularly the closed zones), in the fracture
McKirahan, J.R.; Goldstein, R.H.; Franseen, E.K.
2005-01-01
This study analyzes the three-dimensional variability of a 20-meter-thick section of Pennsylvanian (Missourian) strata over a 600 km2 area of northeastern Kansas, USA. It hypothesizes that sea-level changes interact with subtle variations in paleotopography to influence the heterogeneity of potential reservoir systems in mixed carbonate-silidclastic systems, commonly produdng build-and-fill sequences. For this analysis, ten lithofacies were identified: (1) phylloid algal boundstone-packstone, (2) skeletal wackestone-packstone, (3) peloidal, skeletal packstone, (4) sandy, skeletal grainstone-packstone, (5) oolite grainstone-packstone, (6) Osagia-brachiopod packstone, (7) fossiliferous siltstone, (8) lenticular bedded-laminated siltstone and fine sandstone, (9) organic-rich mudstone and coal, and (10) massive mudstone. Each facies can be related to depositional environment and base-level changes to develop a sequence stratigraphy consisting of three sequence boundaries and two flooding surfaces. Within this framework, eighteen localities are used to develop a threedimensional framework of the stratigraphy and paleotopography. The studied strata illustrate the model of "build-and-fill". In this example, phylloid algal mounds produce initial relief, and many of the later carbonate and silidclastic deposits are focused into subtle paleotopographic lows, responding to factors related to energy, source, and accommodation, eventually filling the paleotopography. After initial buildup of the phylloid algal mounds, marine and nonmarine siliciclastics, with characteristics of both deltaic lobes and valley fills, were focused into low areas between mounds. After a sea-level rise, oolitic carbonates formed on highs and phylloid algal facies accumulated in lows. A shift in the source direction of siliciclastics resulted from flooding or filling of preexisting paleotopographic lows. Fine-grained silidclastics were concentrated in paleotopographic low areas and resulted in clay
Opp, Daniel; Nelms, Benjamin E; Zhang, Geoffrey; Stevens, Craig; Feygelman, Vladimir
2013-01-01
3DVH software (Sun Nuclear Corp., Melbourne, FL) is capable of generating a volumetric patient VMAT dose by applying a volumetric perturbation algorithm based on comparing measurement-guided dose reconstruction and TPS-calculated dose to a cylindrical phantom. The primary purpose of this paper is to validate this dose reconstruction on an anthropomorphic heterogeneous thoracic phantom by direct comparison to independent measurements. The dosimetric insert to the phantom is novel, and thus the secondary goal is to demonstrate how it can be used for the hidden target end-to-end testing of VMAT treatments in lung. A dosimetric insert contains a 4 cm diameter unit-density spherical target located inside the right lung (0.21 g/cm(3) density). It has 26 slots arranged in two orthogonal directions, milled to hold optically stimulated luminescent dosimeters (OSLDs). Dose profiles in three cardinal orthogonal directions were obtained for five VMAT plans with varying degrees of modulation. After appropriate OSLD corrections were applied, 3DVH measurement-guided VMAT dose reconstruction agreed 100% with the measurements in the unit density target sphere at 3%/3 mm level (composite analysis) for all profile points for the four less-modulated VMAT plans, and for 96% of the points in the highly modulated C-shape plan (from TG-119). For this latter plan, while 3DVH shows acceptable agreement with independent measurements in the unit density target, in the lung disagreement with experiment is relatively high for both the TPS calculation and 3DVH reconstruction. For the four plans excluding the C-shape, 3%/3 mm overall composite analysis passing rates for 3DVH against independent measurement ranged from 93% to 100%. The C-shape plan was deliberately chosen as a stress test of the algorithm. The dosimetric spatial alignment hidden target test demonstrated the average distance to agreement between the measured and TPS profiles in the steep dose gradient area at the edge of the 2 cm
NASA Astrophysics Data System (ADS)
Obermann, Anne; Planès, Thomas; Hadziioannou, Céline; Campillo, Michel
2016-07-01
In the context of seismic monitoring, recent studies made successful use of seismic coda waves to locate medium changes on the horizontal plane. Locating the depth of the changes, however, remains a challenge. In this paper, we use 3-D wavefield simulations to address two problems: firstly, we evaluate the contribution of surface and body wave sensitivity to a change at depth. We introduce a thin layer with a perturbed velocity at different depths and measure the apparent relative velocity changes due to this layer at different times in the coda and for different degrees of heterogeneity of the model. We show that the depth sensitivity can be modelled as a linear combination of body- and surface-wave sensitivity. The lapse-time dependent sensitivity ratio of body waves and surface waves can be used to build 3-D sensitivity kernels for imaging purposes. Secondly, we compare the lapse-time behavior in the presence of a perturbation in horizontal and vertical slabs to address, for instance, the origin of the velocity changes detected after large earthquakes.
EarthServer - an FP7 project to enable the web delivery and analysis of 3D/4D models
NASA Astrophysics Data System (ADS)
Laxton, John; Sen, Marcus; Passmore, James
2013-04-01
EarthServer aims at open access and ad-hoc analytics on big Earth Science data, based on the OGC geoservice standards Web Coverage Service (WCS) and Web Coverage Processing Service (WCPS). The WCS model defines "coverages" as a unifying paradigm for multi-dimensional raster data, point clouds, meshes, etc., thereby addressing a wide range of Earth Science data including 3D/4D models. WCPS allows declarative SQL-style queries on coverages. The project is developing a pilot implementing these standards, and will also investigate the use of GeoSciML to describe coverages. Integration of WCPS with XQuery will in turn allow coverages to be queried in combination with their metadata and GeoSciML description. The unified service will support navigation, extraction, aggregation, and ad-hoc analysis on coverage data from SQL. Clients will range from mobile devices to high-end immersive virtual reality, and will enable 3D model visualisation using web browser technology coupled with developing web standards. EarthServer is establishing open-source client and server technology intended to be scalable to Petabyte/Exabyte volumes, based on distributed processing, supercomputing, and cloud virtualization. Implementation will be based on the existing rasdaman server technology developed. Services using rasdaman technology are being installed serving the atmospheric, oceanographic, geological, cryospheric, planetary and general earth observation communities. The geology service (http://earthserver.bgs.ac.uk/) is being provided by BGS and at present includes satellite imagery, superficial thickness data, onshore DTMs and 3D models for the Glasgow area. It is intended to extend the data sets available to include 3D voxel models. Use of the WCPS standard allows queries to be constructed against single or multiple coverages. For example on a single coverage data for a particular area can be selected or data with a particular range of pixel values. Queries on multiple surfaces can be
NASA Astrophysics Data System (ADS)
Chen, Lin; Gerya, Taras V.
2016-04-01
Preexisting lateral variations in crustal thickness and lithospheric thermal state are documented for the formation of some orogenic plateaux. Here we use high-resolution 3-D thermo-mechanical simulations to investigate the influence of preexisting lateral lithospheric strength heterogeneity on the growth of orogenic plateau. The modeling results illustrate an episodic scenario for plateau growth: (1) an early rapid growth stage, characterized by rapid surface uplift and intensive crustal buckling and thickening; (2) an outward spreading stage, characterized by significant lateral expansion of the plateau edges; and (3) a mature stage, characterized by the development of the intracrustal partial melting and subduction of the surrounding lithosphere under the plateau. Sensitivity analyses indicate that lateral variation in crustal thickness favors outward spreading of orogenic plateau, while lateral variation in geothermal gradient favors crustal buckling. The model in absence of lateral strength heterogeneity leads to progressive migration of orogenic belt. Our models show that the plateau's lower crust is largely coupled with underlying lithospheric mantle and does not flow into the surrounding lithospheres, casting doubt on the lower crust flow model. We suggest that the Himalayan-Tibetan orogenic system can be best understood within the framework that the proto-southern Asian margin was fairly weak prior to the India-Asia collision to steer the formation of a large hot orogenic plateau there.
NASA Astrophysics Data System (ADS)
Boisvert, Alexandre
sont par consequent bases sur des resultats de simulations avec le logiciel Cadence. Mots-cles : Circuit d'etouffement, Photodiodes a avalanche monophotoniques (PAMP), Single Photon Avalanche Diode (SPAD), Integration 3D heterogene, Drain-Extended MOS (DEMOS), CMOS 130 nm 3D Tezzaron/Chartered, Tomographie d'emission par positrons (TEP)
NASA Astrophysics Data System (ADS)
Dick, Jonathan; Tetzlaff, Doerthe; Bradford, John; Soulsby, Chris
2016-04-01
Soil moisture is a fundamental component of the water cycle that influences many hydrological processes, such as flooding, solute transport, biogeochemical processes, and land-atmosphere interactions. The relationship between vegetation and soil moisture is complex and reciprocal. Soil moisture may affect vegetation distribution due to its function as the primary source of water, in turn the structure of vegetation canopies regulate water partitioning into interception, throughfall and steam flow. Such spatial differences in inputs, together with complex patterns of water uptake from distributed root networks can create marked heterogeneity in soil moisture dynamics at small scales. Traditional methods of monitoring soil moisture have revolved around limited point measurements, but improved geophysical techniques have facilitated a trend towards more spatially distributed measurements to help understand this heterogeneity. Here, we present a study using 3D ERT surveys in a 3.2km upland catchment in the Scottish Highlands where increasing afforestation (for climate change adaptation, biofuels and conservation) has the potential to increase interception losses and reduce soil moisture storage. The study combined 3D surveys, traditional point measurements and laboratory analysis of soil cores to assess the plot scale soil moisture dynamics in podzolic soils under forest stands of 15m high Scots pine (Pinus sylvestris) and adjacent non-forest plots dominated by heather (Calluna vulgaris) shrubs (<0.5m high). These dominant species are typical of forest and non-forest vegetation communities the Scottish Highlands. Results showed differences in the soil moisture dynamics under the different vegetation types, with heterogeneous patterns in the forested site mainly correlated with canopy cover and mirroring interception losses. Temporal variability in the forested site was greater, probably due to the interception, and increased evapotranspiration losses relative to the
Earth abundant bimetallic nanoparticles for heterogeneous catalysis
NASA Astrophysics Data System (ADS)
Senn, Jonathan F., Jr.
Polymer exchange membrane fuel cells have the potential to replace current fossil fuel-based technologies in terms of emissions and efficiency, but CO contamination of H2 fuel, which is derived from steam methane reforming, leads to system inefficiency or failure. Solutions currently under development are bimetallic nanoparticles comprised of earth-abundant metals in different architectures to reduce the concentration of CO by PROX during fuel cell operation. Chapter One introduces the Pt-Sn and Co-Ni bimetallic nanoparticle systems, and the intermetallic and core-shell architectures of interest for catalytic evaluation. Application, theory, and studies associated with the efficacy of these nanoparticles are briefly reviewed. Chapter Two describes the concepts of the synthetic and characterization methods used in this work. Chapter Three presents the synthetic, characterization, and catalytic findings of this research. Pt, PtSn, PtSn2, and Pt 3Sn nanoparticles have been synthesized and supported on gamma-Al2O3. Pt3Sn was shown to be an effective PROX catalyst in various gas feed conditions, such as the gas mixture incorporating 0.1% CO, which displayed a light-off temperatures of ˜95°C. Co and Ni monometallic and CoNi bimetallic nanoparticles have been synthesized and characterized, ultimately leading to the development of target Co Ni core-shell nanoparticles. Proposed studies of catalytic properties of these nanoparticles in preferential oxidation of CO (PROX) reactions will further elucidate the effects of different crystallographic phases, nanoparticle-support interactions, and architecture on catalysis, and provide fundamental understanding of catalysis with nanoparticles composed of earth abundant metals in different architectures.
3D Online Visualization and Synergy of NASA A-Train Data Using Google Earth
NASA Technical Reports Server (NTRS)
Chen, Aijun; Kempler, Steven; Leptoukh, Gregory; Smith, Peter
2010-01-01
This poster presentation reviews the use of Google Earth to assist in three dimensional online visualization of NASA Earth science and geospatial data. The NASA A-Train satellite constellation is a succession of seven sun-synchronous orbit satellites: (1) OCO-2 (Orbiting Carbon Observatory) (will launch in Feb. 2013), (2) GCOM-W1 (Global Change Observation Mission), (3) Aqua, (4) CloudSat, (5) CALIPSO (Cloud-Aerosol Lidar & Infrared Pathfinder Satellite Observations), (6) Glory, (7) Aura. The A-Train makes possible synergy of information from multiple resources, so more information about earth condition is obtained from the combined observations than would be possible from the sum of the observations taken independently
NASA Astrophysics Data System (ADS)
Wehrer, Markus; Slater, Lee
2015-04-01
flow fraction was observed to be independent of precipitation rate. This suggests the presence of a fingering process driven by textural heterogeneities. As a consequence, preferential transport of the conservative and the reactive tracer also occurred. We found that 3D ERT can serve to quantitatively characterize shape measures of both tracer breakthroughs and water content dynamics. In particular, shape measures influenced by the advective propagation of the tracer peak, like mean velocity and normalized first central moment, are highly correlated between ERT data and validation data (consisting of tracer measurements in seepage water samples). Using shape measures proved to be advantageous over interpretation of ERT data with spatially uncertain petrophysical functions for the characterization of heterogeneous flow and transport. Consequently, for future applications of ERT in soil hydrological modeling, the use of temporal moments is recommended.
Sunlight effects on the 3D polar current system determined from low Earth orbit measurements
NASA Astrophysics Data System (ADS)
Laundal, Karl M.; Finlay, Christopher C.; Olsen, Nils
2016-08-01
Interaction between the solar wind and the Earth's magnetosphere is associated with large-scale currents in the ionosphere at polar latitudes that flow along magnetic field lines (Birkeland currents) and horizontally. These current systems are tightly linked, but their global behaviors are rarely analyzed together. In this paper, we present estimates of the average global Birkeland currents and horizontal ionospheric currents from the same set of magnetic field measurements. The magnetic field measurements, from the low Earth orbiting Swarm and CHAMP satellites, are used to co-estimate poloidal and toroidal parts of the magnetic disturbance field, represented in magnetic apex coordinates. The use of apex coordinates reduces effects of longitudinal and hemispheric variations in the Earth's main field. We present global currents from both hemispheres during different sunlight conditions. The results show that the Birkeland currents vary with the conductivity, which depends most strongly on solar EUV emissions on the dayside and on particle precipitation at pre-midnight magnetic local times. In sunlight, the horizontal equivalent current flows in two cells, resembling an opposite ionospheric convection pattern, which implies that it is dominated by Hall currents. By combining the Birkeland current maps and the equivalent current, we are able to calculate the total horizontal current, without any assumptions about the conductivity. We show that the total horizontal current is close to zero in the polar cap when it is dark. That implies that the equivalent current, which is sensed by ground magnetometers, is largely canceled by the horizontal closure of the Birkeland currents.
Exploring the possible climates of the Archean Earth with a 3D GCM
NASA Astrophysics Data System (ADS)
Charnay, B.; Forget, F.; Wordsworth, R.; Leconte, J.; Millour, E.; Codron, F.
2012-12-01
The climate of the early Earth and the faint young Sun paradox have been hitherto essentially studied with 1D models. Different hypothesis have been suggest to solve the faint young Sun paradox as: a higher greenhouse effect by CO2 and CH4 [1], a reduced planetary albedo by a reduced continental area and thinner archean clouds [2] or a higher atmospheric pressure [3]. The problem of 1D models is that they do not simulate the oceanic ice and cloud feedback well. A new global climate model has been developped by our team to study the climates of exoplanets [4] and primitive atmospheres [5]. A correlated-k radiative transfert model is use to simulate atmospheres of N2, CO2 and CH4. Oceanic transport and oceanic ice formation is computed through a slab-ocean module [6], cloud formation and precipitations are based on physical and robust parametrizations. First, we validate our model by simulating the modern Earth. Then we apply it to simulate the climates during the Archean and to test different hypothesis supposed to solve the faint young Sun paradox. We quantify the effect of the lack of ozone, the faster rotation rate and the reduced continental area of the early Earth on the climate. We simulate climates for low, moderate and high CO2/CH4 concentration. We test the effect the radius of cloud droplets and the effect of a higher atmospheric pressure. In particular, by fixing the radius of cloud droplets or the amount of CCN, we show that thinner clouds with a low CO2 and CH4 partial pressure [2] is not a solution over all the Archean, higher greenhouse gas concentrations are required. References [1] Haqq-Misra et al.: A revised, hazy methane greenhouse for the Archean Earth, Astrobiology,8,1127-1137 (2008) [2] Rosing et al.: No climate paradox under the faint early Sun, Nature, 464, 744-747 (2010) [3] Goldblatt et al.: Nitrogen-enhanced greenhouse warming on early Earth, Nature Geoscience, 2, 891-896 (2009) [4] Wordsworth et al.: Gliese 581d is the First Discovered
Strong, Multi-Scale Heterogeneity in Earth's Lowermost Mantle.
Tkalčić, Hrvoje; Young, Mallory; Muir, Jack B; Davies, D Rhodri; Mattesini, Maurizio
2015-01-01
The core mantle boundary (CMB) separates Earth's liquid iron outer core from the solid but slowly convecting mantle. The detailed structure and dynamics of the mantle within ~300 km of this interface remain enigmatic: it is a complex region, which exhibits thermal, compositional and phase-related heterogeneity, isolated pockets of partial melt and strong variations in seismic velocity and anisotropy. Nonetheless, characterising the structure of this region is crucial to a better understanding of the mantle's thermo-chemical evolution and the nature of core-mantle interactions. In this study, we examine the heterogeneity spectrum from a recent P-wave tomographic model, which is based upon trans-dimensional and hierarchical Bayesian imaging. Our tomographic technique avoids explicit model parameterization, smoothing and damping. Spectral analyses reveal a multi-scale wavelength content and a power of heterogeneity that is three times larger than previous estimates. Inter alia, the resulting heterogeneity spectrum gives a more complete picture of the lowermost mantle and provides a bridge between the long-wavelength features obtained in global S-wave models and the short-scale dimensions of seismic scatterers. The evidence that we present for strong, multi-scale lowermost mantle heterogeneity has important implications for the nature of lower mantle dynamics and prescribes complex boundary conditions for Earth's geodynamo. PMID:26674394
NASA Astrophysics Data System (ADS)
Ryu, Dongok; Kim, Sug-Whan; Kim, Dae Wook; Lee, Jae-Min; Lee, Hanshin; Park, Won Hyun; Seong, Sehyun; Ham, Sun-Jeong
2010-09-01
Understanding the Earth spectral bio-signatures provides an important reference datum for accurate de-convolution of collapsed spectral signals from potential earth-like planets of other star systems. This study presents a new ray tracing computation method including an improved 3D optical earth model constructed with the coastal line and vegetation distribution data from the Global Ecological Zone (GEZ) map. Using non-Lambertian bidirectional scattering distribution function (BSDF) models, the input earth surface model is characterized with three different scattering properties and their annual variations depending on monthly changes in vegetation distribution, sea ice coverage and illumination angle. The input atmosphere model consists of one layer with Rayleigh scattering model from the sea level to 100 km in altitude and its radiative transfer characteristics is computed for four seasons using the SMART codes. The ocean scattering model is a combination of sun-glint scattering and Lambertian scattering models. The land surface scattering is defined with the semi empirical parametric kernel method used for MODIS and POLDER missions. These three component models were integrated into the final Earth model that was then incorporated into the in-house built integrated ray tracing (IRT) model capable of computing both spectral imaging and radiative transfer performance of a hypothetical space instrument as it observes the Earth from its designated orbit. The IRT model simulation inputs include variation in earth orientation, illuminated phases, and seasonal sea ice and vegetation distribution. The trial simulation runs result in the annual variations in phase dependent disk averaged spectra (DAS) and its associated bio-signatures such as NDVI. The full computational details are presented together with the resulting annual variation in DAS and its associated bio-signatures.
NASA Astrophysics Data System (ADS)
Macquet, M.; Paul, A.; Pedersen, H.
2013-12-01
A 3-D Vs model is the output of ambient noise tomography with the greatest add-on value. Computing this model is however difficult, especially in regions with very heterogeneous crustal structure. We propose here a new approach to invert group velocity data from ambient noise tomography for 3-D Vs structure, which gives good result even with strong lateral variations of the seismic structure. The study region is the Pyrenees mountain range s.l., at the border between France and Spain, as it includes 2 thick sedimentary basins surrounded by crystalline outcrops in the Pyrenees and the Massif Central. We use data of a dense temporary seismic broadband array of 49 stations installed in southwestern France at the end of 2010 in the framework of the PYROPE (PYRenean Observational Portable Experiment) project. The dataset also includes records of 70 broadband stations of the third leg of the IberArray project installed at the same time period in northern Spain. The two dense arrays with 60 km average inter-station distance make it possible to obtain high-resolution images of the lithosphere of the Pyrenees, the Ebro and Aquitaine basins and the Bay of Biscay. To complement the two temporary arrays and avoid smearing effects along their edges, we also used records of the French and Catalan permanent broadband networks. As a whole, the dataset includes records of 158 broadband stations. Using one year of data, we computed 12324 two-station correlations of ambient seismic noise records. We applied the processing in two overlapping period intervals: 5s-50s and 30s-75s, the latter to extend the subsequent measurement of group velocities to periods outside the first and secondary microseismic peak. We merged the common part of the 2 dispersion curves using a linear weighting to obtain, for each station pair, a single dispersion curve in the 5-75s period range. These group velocity curves are subsequently used to obtain 2-D group velocity maps using the tomographic method of
NASA Astrophysics Data System (ADS)
Green, A.; Gribenko, A.; Cuma, M.; Zhdanov, M. S.
2008-12-01
In this paper we apply 3D inversion to MT data collected in Oregon as a part of the EarthScope project. We use the integral equation method as a forward modeling engine. Quasi-analytical approximation with a variable background (QAVB) method of Frechet derivative calculation is applied. This technique allows us to simplify the inversion algorithm and to use just one forward modeling on every iteration step. The receiver footprint approach considerably reduces the computational resources needed to invert the large volumes of data covering vast areas. The data set, which was used in the inversion, was obtained through the Incorporated Research Institutions for Seismology (IRIS). The long-period MT data was collected in Eastern Oregon in 2006. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region as well as with 3D MT inversion results obtained by other researchers. The geoelectrical model of the Oregon deep interior produced by 3D inversion indicates several lithospheres' electrical conductivity anomalies, including a linear zone marked by low-high conductivity transition along the Klamath Blue Mountain Lineament associated with a linear trend of gravity minima. High electrical conductivity values occur in the upper crust under the accreted terrains in the Blue Mountains region.
Recipe for High Moment Materials with Rare-earth and 3d Transition Metal Composites
Autieri, Carmine; Kumar, P. Anil; Walecki, Dirk; Webers, Samira; Gubbins, Mark A.; Wende, Heiko; Sanyal, Biplab
2016-01-01
Materials with high volume magnetization are perpetually needed for the generation of sufficiently large magnetic fields by writer pole of magnetic hard disks, especially for achieving increased areal density in storage media. In search of suitable materials combinations for this purpose, we have employed density functional theory to predict the magnetic coupling between iron and gadolinium layers separated by one to several monolayers of 3d transition metals (Sc-Zn). We demonstrate that it is possible to find ferromagnetic coupling for many of them and in particular for the early transition metals giving rise to high moment. Cr and Mn are the only elements able to produce a significant ferromagnetic coupling for thicker spacer layers. We also present experimental results on two trilayer systems Fe/Sc/Gd and Fe/Mn/Gd. From the experiments, we confirm a ferromagnetic coupling between Fe and Gd across a 3 monolayers Sc spacer or a Mn spacer thicker than 1 monolayer. In addition, we observe a peculiar dependence of Fe/Gd magnetic coupling on the Mn spacer thickness. PMID:27381456
Recipe for High Moment Materials with Rare-earth and 3d Transition Metal Composites
NASA Astrophysics Data System (ADS)
Autieri, Carmine; Kumar, P. Anil; Walecki, Dirk; Webers, Samira; Gubbins, Mark A.; Wende, Heiko; Sanyal, Biplab
2016-07-01
Materials with high volume magnetization are perpetually needed for the generation of sufficiently large magnetic fields by writer pole of magnetic hard disks, especially for achieving increased areal density in storage media. In search of suitable materials combinations for this purpose, we have employed density functional theory to predict the magnetic coupling between iron and gadolinium layers separated by one to several monolayers of 3d transition metals (Sc-Zn). We demonstrate that it is possible to find ferromagnetic coupling for many of them and in particular for the early transition metals giving rise to high moment. Cr and Mn are the only elements able to produce a significant ferromagnetic coupling for thicker spacer layers. We also present experimental results on two trilayer systems Fe/Sc/Gd and Fe/Mn/Gd. From the experiments, we confirm a ferromagnetic coupling between Fe and Gd across a 3 monolayers Sc spacer or a Mn spacer thicker than 1 monolayer. In addition, we observe a peculiar dependence of Fe/Gd magnetic coupling on the Mn spacer thickness.
Recipe for High Moment Materials with Rare-earth and 3d Transition Metal Composites.
Autieri, Carmine; Kumar, P Anil; Walecki, Dirk; Webers, Samira; Gubbins, Mark A; Wende, Heiko; Sanyal, Biplab
2016-01-01
Materials with high volume magnetization are perpetually needed for the generation of sufficiently large magnetic fields by writer pole of magnetic hard disks, especially for achieving increased areal density in storage media. In search of suitable materials combinations for this purpose, we have employed density functional theory to predict the magnetic coupling between iron and gadolinium layers separated by one to several monolayers of 3d transition metals (Sc-Zn). We demonstrate that it is possible to find ferromagnetic coupling for many of them and in particular for the early transition metals giving rise to high moment. Cr and Mn are the only elements able to produce a significant ferromagnetic coupling for thicker spacer layers. We also present experimental results on two trilayer systems Fe/Sc/Gd and Fe/Mn/Gd. From the experiments, we confirm a ferromagnetic coupling between Fe and Gd across a 3 monolayers Sc spacer or a Mn spacer thicker than 1 monolayer. In addition, we observe a peculiar dependence of Fe/Gd magnetic coupling on the Mn spacer thickness. PMID:27381456
NASA Astrophysics Data System (ADS)
Huang, Guo-Jiao; Bai, Chao-Ying; Greenhalgh, Stewart
2013-09-01
The traditional grid/cell-based wavefront expansion algorithms, such as the shortest path algorithm, can only find the first arrivals or multiply reflected (or mode converted) waves transmitted from subsurface interfaces, but cannot calculate the other later reflections/conversions having a minimax time path. In order to overcome the above limitations, we introduce the concept of a stationary minimax time path of Fermat's Principle into the multistage irregular shortest path method. Here we extend it from Cartesian coordinates for a flat earth model to global ray tracing of multiple phases in a 3-D complex spherical earth model. The ray tracing results for 49 different kinds of crustal, mantle and core phases show that the maximum absolute traveltime error is less than 0.12 s and the average absolute traveltime error is within 0.09 s when compared with the AK135 theoretical traveltime tables for a 1-D reference model. Numerical tests in terms of computational accuracy and CPU time consumption indicate that the new scheme is an accurate, efficient and a practical way to perform 3-D multiphase arrival tracking in regional or global traveltime tomography.
NASA Astrophysics Data System (ADS)
Chen, P.; Lee, E.; Jordan, T. H.; Maechling, P. J.
2009-12-01
Accurate and rapid CMT inversion is important for seismic hazard analysis. We have developed an algorithm for very rapid CMT inversions in a 3D Earth structure model and applied it on small to medium-sized earthquakes recorded by the Southern California Seismic Network (SCSN). Our CMT inversion algorithm is an integral component of the scattering-integral (SI) method for full-3D waveform tomography (F3DT). In the SI method for F3DT, the sensitivity (Fréchet) kernels are constructed through the temporal convolution between the earthquake wavefield (EWF) and the receiver Green tensor (RGT), which is the wavefield generated by 3 orthogonal unit impulsive body forces acting at the receiver location. The RGTs are also the partial derivatives of the waveform with respect to the moment tensors. In this study, our RGTs are computed in a 3D seismic structure model for Southern California (CVM4SI1) using the finite-difference method, which allows us to account for 3D path effects in our source inversion. We used three component broadband waveforms below 0.2 Hz. An automated waveform-picking algorithm based on continuous wavelet transform is applied on observed waveforms to pick P, S and surface waves. A multi-scale grid-searching algorithm is then applied on the picked waveforms to find the optimal strike, dip and rake values that minimize the amplitude misfit and maximize the correlation coefficient. In general, our CMT solutions agree with solutions inverted using other methods and provide better fit to the observed waveforms.
NASA Astrophysics Data System (ADS)
Mu, D.; Lee, E.; Chen, P.; Jordan, T. H.; Maechling, P. J.
2010-12-01
Accurate and rapid CMT inversion is important for seismic hazard analysis. We have developed an algorithm for very rapid CMT inversions in a 3D Earth structure model and applied it on small to medium-sized earthquakes recorded by the Southern California Seismic Network (SCSN). Our CMT inversion algorithm is an integral component of the scattering-integral (SI) method for full-3D waveform tomography (F3DT). In the SI method for F3DT, the sensitivity (Fréchet) kernels are constructed through the temporal convolution between the earthquake wavefield (EWF) from the source and the receiver Green tensor (RGT) from the receiver. In this study, our RGTs were computed in a 3D seismic structure model for Southern California (CVM4SI1) using the finite-difference method, which allows us to account for 3D path effects in our source inversion. By storing the RGTs, synthetic seismograms for any source in our modeling volume could be generated rapidly by applying the reciprocity principle. An automated waveform-picking algorithm based on continuous wavelet transform is applied on observed waveforms to pick P, S and surface waves. A grid-searching algorithm is then applied on the picked waveforms to find an optimal focal mechanism that minimizes the amplitude misfit and maximize the weighted correlation coefficient. The grid-search result is then used as the initial solution in a gradient-based optimization algorithm that minimizes the L2 norm of the generalized seismological data functionals (GSDF), which quantifies waveform differences between observed and synthetic seismograms using frequencies-dependent phase-delay and amplitude anomalies. In general, our CMT solutions agree with solutions inverted using other methods and provide better fit to the observed waveforms.
Adsorption of alkali, alkaline-earth, and 3d transition metal atoms on silicene
NASA Astrophysics Data System (ADS)
Sahin, H.; Peeters, F. M.
2013-02-01
The adsorption characteristics of alkali, alkaline-earth, and transition metal adatoms on silicene, a graphene-like monolayer structure of silicon are analyzed by means of first-principles calculations. In contrast to graphene, interaction between the metal atoms and the silicene surface is quite strong due to its highly reactive buckled hexagonal structure. In addition to structural properties, we also calculate the electronic band dispersion, net magnetic moment, charge transfer, work function, and dipole moment of the metal adsorbed silicene sheets. Alkali metals, Li, Na, and K, adsorb to hollow sites without any lattice distortion. As a consequence of the significant charge transfer from alkalis to silicene, metalization of silicene takes place. Trends directly related to atomic size, adsorption height, work function, and dipole moment of the silicene/alkali adatom system are also revealed. We found that the adsorption of alkaline-earth metals on silicene is entirely different from their adsorption on graphene. The adsorption of Be, Mg, and Ca turns silicene into a narrow gap semiconductor. Adsorption characteristics of eight transition metals Ti, V, Cr, Mn, Fe, Co, Mo, and W are also investigated. As a result of their partially occupied d orbital, transition metals show diverse structural, electronic, and magnetic properties. Upon the adsorption of transition metals, depending on the adatom type and atomic radius, the system can exhibit metal, half-metal, and semiconducting behavior. For all metal adsorbates, the direction of the charge transfer is from adsorbate to silicene, because of its high surface reactivity. Our results indicate that the reactive crystal structure of silicene provides a rich playground for functionalization at nanoscale.
NASA Astrophysics Data System (ADS)
Schuberth, Bernhard; Zaroli, Christophe; Nolet, Guust
2015-04-01
Of particular interest for the tectonic evolution of the Atlantic region is the influence of lower mantle structure under Africa on flow in the upper mantle beneath the ocean basin. Along with its Pacific counterpart, the large African anomaly in the lowermost mantle with strongly reduced seismic velocities has received considerable attention in seismological and geodynamic studies. Several seismological observations are typically taken as an indication that these two anomalies are being caused by large-scale compositional variations and that they are piles of material with higher density than normal mantle rock. This would imply negative buoyancy in the lowermost mantle under Africa, which has important implications for the flow at shallower depth and inferences on the processes that led to the formation of the Atlantic Ocean basin. However, a large number of recent studies argue for a strong thermal gradient across the core-mantle boundary that might provide an alternative explanation for the lower mantle anomaly through the resulting large lateral temperature variations. Recently, we developed a new joint forward modeling approach to test such geodynamic hypotheses directly against the seismic observations: Seismic heterogeneity is predicted by converting the temperature field of a high-resolution 3-D mantle circulation model into seismic velocities using thermodynamic models of mantle mineralogy. 3-D global wave propagation in the synthetic elastic structures is then simulated using a spectral element method. Being based on forward modelling only, this approach allows us to generate synthetic wavefields and seismograms independently of seismic observations. The statistics of observed long-period body wave traveltime variations show a markedly different behaviour for P- and S-waves: the standard deviation of P-wave delay times stays almost constant with ray turning depth, while that of the S-wave delay times increases strongly throughout the mantle. In an
NASA Astrophysics Data System (ADS)
Greenhalgh, Stewart; Zhou, Bing; Maurer, Hansruedi
2010-05-01
We have developed a modified version of the spectral element method (SEM), called the Gaussian Quadrature Grid (GQG) approach, for frequency domain 3D seismic modelling in arbitrary heterogeneous, anisotropic media. The model may incorporate an arbitrary free-surface topography and irregular subsurface interfaces. Unlike the SEM ,it does not require a powerful mesh generator such as the Delauney Triangular or TetGen. Rather, the GQG approach replaces the element mesh with Gaussian quadrature abscissae to directly sample the physical properties of the model parameters and compute the weighted residual or variational integral. This renders the model discretisation simple and easily matched to the model topography, as well as direct control of the model paramterisation for subsequent inversion. In addition, it offers high accuracy in numerical modelling provided that an appropriate density of the Gaussian quadrature abscissae is employed. The second innovation of the GQG is the incorporation of a new implementation of perfectly matched layers to suppress artificial reflections from the domain margins. We employ PML model parameters (specified complex valued density and elastic moduli) rather than explicitly solving the governing wave equation with a complex co-ordinate system as in conventional approaches. Such an implementation is simple, general, effective and easily extendable to any class of anisotropy and other numerical modelling methods. The accuracy of the GQG approach is controlled by the number of Gaussian quadrature points per minimum wavelength, the so-called sampling density. The optimal sampling density should be the one which enables high definition of geological characteristics and high precision of the variational integral evaluation and spatial differentiation. Our experiments show that satisfactory results can be obtained using sampling densities of 5 points per minimum wavelength. Efficiency of the GQG approach mainly depends on the linear
NASA Astrophysics Data System (ADS)
Bianco, T. A.; Ito, G.; van Hunen, J.; Ballmer, M.; Mahoney, J. J.
2006-12-01
Spatial variations in magma geochemistry among hotspot volcanoes hold clues to the dynamics and composition of the mantle feeding hotspot volcanism. We use a 3D geodynamic model of plume-lithosphere interaction to explore the causes of spatial patterns of magmatic volumes and compositions at intraplate hotspots. This study focuses on coupling between upper mantle flow, heat transfer, and melting of a heterogeneous (veined) plume. We assume multiple lithologies have different solidi, trace-element, and isotope composition. We use the Cartesian finite-element code, CITCOM, (Zhong and Watts, 2002) to simulate mantle convection with the extended Boussinesq approximation in a volume of upper mantle 400 km in thickness. A parameterized melting model is used to simulate melting of materials with different water contents (Katz et al., 2003). Melt depletion (F) for each lithology is calculated at finite element nodes as a function of temperature, pressure, and water content and is advected using particle tracers. We quantify the response of the geographic pattern of the volume and composition of magmas to different lithospheric thicknesses, and plume temperatures and viscosities, which together control the melting rates and sizes of the melting zones for the different lithologies. In the case of two-lithologies, preliminary results of a sluggishly convecting plume rising beneath thick lithosphere (60-100 km) predict that the melting zone of the least refractory "lithology 1" is wider than that of the more refractory "lithology 2". This leads to the prediction that on the surface, the isotope signature of lithology 1 is most prominent at the leading edge (i.e., upwind edge of plate motion) of the hotspot, whereas the isotope signature of lithology 2 is strongest at the hotspot center. This pattern will likely change for plumes convecting more vigorously or thinner lithosphere.
NASA Astrophysics Data System (ADS)
Kreuter, A.; Buras, R.; Mayer, B.; Webb, A.; Kift, R.; Bais, A.; Kouremeti, N.; Blumthaler, M.
2014-06-01
We present a unique case study of the solar global irradiance in a highly heterogeneous albedo environment at the Arctic coast. Diodearray spectroradiometers were deployed at three sites around Ny Ålesund, Svalbard, and spectral irradiances were simultaneously measured under clear-sky conditions during a 24 h period. The 3-D radiative transfer model MYSTIC is applied to simulate the measurements in various model scenarios. First, we model the effective albedos of ocean and snow and consequently around each measurement site. The effective albedos at 340 nm increase from 0.57 to 0.75, from the coastal site in the west towards the site 20 km east, away from the coast. The observed ratios of the global irradiance indicate a 15% higher average irradiance, at 340 nm east relative to west, due to the higher albedo. The comparison of our model scenarios suggest a snow albedo of > 0.9 and confirm the observation that drift ice has moved into the Fjord during the day. The local time shift between the locations causes a hysteresis-like behavior of these east-west ratios with solar zenith angle (SZA). The observed hysteresis, however, is larger and, at 340 nm, can be explained by the drift ice. At 500 nm, a plausible explanation is a detector tilt of about 1°. The ratios between afternoon and morning irradiances at the same SZA are investigated, which confirm the above conclusions. At the coastal site, the measured irradiance is significantly higher in the afternoon than in the morning. Besides the effect of changing drift ice and detector tilt, the small variations of the aerosol optical depth have to be considered also at the other stations to reduce the discrepancies between model and observations. Remaining discrepancies are possibly due to distant high clouds.
NASA Astrophysics Data System (ADS)
Kreuter, A.; Buras, R.; Mayer, B.; Webb, A.; Kift, R.; Bais, A.; Kouremeti, N.; Blumthaler, M.
2014-02-01
We present a unique case study of the solar global irradiance in a highly heterogeneous albedo environment at the arctic coast. Diodearray spectroradiometers were deployed at three sites around Ny Ålesund, Svalbard, and spectral irradiances were simultaneously measured under clear sky conditions during a 24 h period. The 3-D radiative transfer model MYSTIC is applied to simulate the measurements in various model scenarios. First, we model the effective albedos of ocean and snow and consequently around each measurement site. The effective albedos at 340 nm increase from 0.57 to 0.75, from the coastal site in the west towards the site 20 km east, away from the coast. The observed ratios of the global irradiance indicate a 15% higher average irradiance at east relative to west at 340 nm due to the higher albedo. The comparison to our model scenarios suggest a snow albedo of >0.9 and confirm the observation that drift ice has moved into the Fjord during the day. The local time shift between the locations causes a hysteresis-like behavior of these east-west ratios with solar zenith angle (SZA). The observed hysteresis, however, is larger and, at 340 nm, can be explained by the drift ice. At 500 nm, a plausible explanation is a detector tilt of about 1°. The ratios between afternoon and morning irradiances at the same SZA are investigated, which confirm the above conclusions. At the coastal site, the measured irradiance is significantly higher in the afternoon than in the morning. Besides the effect of changing drift ice and detector tilt, the small variations of the aerosol optical depth have to be considered also at the other stations to reduce the discrepancies between model and observations. Remaining discrepancies are possibly due to high thin clouds.
A study of the earth radiation budget using a 3D Monte-Carlo radiative transer code
NASA Astrophysics Data System (ADS)
Okata, M.; Nakajima, T.; Sato, Y.; Inoue, T.; Donovan, D. P.
2013-12-01
The purpose of this study is to evaluate the earth's radiation budget when data are available from satellite-borne active sensors, i.e. cloud profiling radar (CPR) and lidar, and a multi-spectral imager (MSI) in the project of the Earth Explorer/EarthCARE mission. For this purpose, we first developed forward and backward 3D Monte Carlo radiative transfer codes that can treat a broadband solar flux calculation including thermal infrared emission calculation by k-distribution parameters of Sekiguchi and Nakajima (2008). In order to construct the 3D cloud field, we tried the following three methods: 1) stochastic cloud generated by randomized optical thickness each layer distribution and regularly-distributed tilted clouds, 2) numerical simulations by a non-hydrostatic model with bin cloud microphysics model and 3) Minimum cloud Information Deviation Profiling Method (MIDPM) as explained later. As for the method-2 (numerical modeling method), we employed numerical simulation results of Californian summer stratus clouds simulated by a non-hydrostatic atmospheric model with a bin-type cloud microphysics model based on the JMA NHM model (Iguchi et al., 2008; Sato et al., 2009, 2012) with horizontal (vertical) grid spacing of 100m (20m) and 300m (20m) in a domain of 30km (x), 30km (y), 1.5km (z) and with a horizontally periodic lateral boundary condition. Two different cell systems were simulated depending on the cloud condensation nuclei (CCN) concentration. In the case of horizontal resolution of 100m, regionally averaged cloud optical thickness,
Incorporation of 3-D Scanning Lidar Data into Google Earth for Real-time Air Pollution Observation
NASA Astrophysics Data System (ADS)
Chiang, C.; Nee, J.; Das, S.; Sun, S.; Hsu, Y.; Chiang, H.; Chen, S.; Lin, P.; Chu, J.; Su, C.; Lee, W.; Su, L.; Chen, C.
2011-12-01
3-D Differential Absorption Scanning Lidar (DIASL) system has been designed with small size, light weight, and suitable for installation in various vehicles and places for monitoring of air pollutants and displays a detailed real-time temporal and spatial variability of trace gases via the Google Earth. The fast scanning techniques and visual information can rapidly identify the locations and sources of the polluted gases and assess the most affected areas. It is helpful for Environmental Protection Agency (EPA) to protect the people's health and abate the air pollution as quickly as possible. The distributions of the atmospheric pollutants and their relationship with local metrological parameters measured with ground based instruments will also be discussed. Details will be presented in the upcoming symposium.
NASA Astrophysics Data System (ADS)
Luo, Cong; Friederich, Wolfgang
2016-04-01
Realistic shallow seismic wave propagation simulation is an important tool for studying induced seismicity (e.g., during geothermal energy development). However over a long time, there is a significant problem which constrains computational seismologists from performing a successful simulation conveniently: pre-processing. Conventional pre-processing has often turned out to be inefficient and unrobust because of the miscellaneous operations, considerable complexity and insufficiency of available tools. An integrated web-based platform for shallow seismic wave propagation simulation has been built. It is aiming at providing a user-friendly pre-processing solution, and cloud-based simulation abilities. The main features of the platform for the user include: revised digital elevation model (DEM) retrieving and processing mechanism; generation of multi-layered 3D shallow Earth model geometry (the computational domain) with user specified surface topography based on the DEM; visualization of the geometry before the simulation; a pipeline from geometry to fully customizable hexahedral element mesh generation; customization and running the simulation on our HPC; post-processing and retrieval of the results over cloud. Regarding the computational aspect, currently the widely accepted specfem3D is chosen as the computational package; packages using different types of elements can be integrated as well in the future. According to our trial simulation experiments, this web-based platform has produced accurate waveforms while significantly simplifying and enhancing the pre-processing and improving the simulation success rate.
NASA Astrophysics Data System (ADS)
Capdeville, Y.; Gung, Y.; Romanowicz, B.
2002-12-01
The spectral element method (SEM) has recently been adapted successfully for global spherical earth wave propagation applications. Its advantage is that it provides a way to compute exact seismograms in a 3D earth, without restrictions on the size or wavelength of lateral heterogeneity at any depth, and can handle diffraction and other interactions with major structural boundaries. Its disadvantage is that it is computationally heavy. In order to partly address this drawback, a coupled SEM/normal mode method was developed (Capdeville et al., 2000). This enables us to more efficiently compute bodywave seismograms to realistically short periods (10s or less). In particular, the coupled SEM/normal mode method is a powerful tool to test the validity of some analytical approximations that are currently used in global waveform tomography, and that are considerably faster computationally. Here, we focus on several approximations based on normal mode perturbation theory: the classical "path-average approximation" (PAVA) introduced by Woodhouse and Dziewonski (1984) and well suited for fundamental mode surface waves (1D sensitivity kernels); the non-linear asymptotic coupling theory (NACT), which introduces coupling between mode branches and 2D kernels in the vertical plane containing the source and the receiver (Li and Tanimoto, 1993; Li and Romanowicz, 1995); an extension of NACT which includes out of plane focusing terms computed asymptotically (e.g. Romanowicz, 1987) and introduces 3D kernels; we also consider first order perturbation theory without asymptotic approximations, such as developed for example by Dahlen et al. (2000). We present the results of comparisons of realistic seismograms for different models of heterogeneity, varying the strength and sharpness of the heterogeneity and its location in depth in the mantle. We discuss the consequences of different levels of approximations on our ability to resolve 3D heterogeneity in the earth's mantle.
Spectral element modeling of 3D wave propagation in the Earth: the graver part of the spectrum
NASA Astrophysics Data System (ADS)
Chaljub, E.; Valette, B.
2003-04-01
The Spectral Element Method (SEM) has been recently established as a new reference to compute synthetic seismograms in 3D models of the Earth. So far, all the studies involving the SEM have been performed within the Cowling approximation, i. e. neglecting the variations of the gravity field during wave propagation. For low-frequency studies (typically less than 5 mHz) the previous assumption fails and the complete treatment of self-gravitation has to be considered. This requires the introduction of the mass redistribution potential (MRP) which has to satisfy Poisson's equation everywhere in space. Unlike spherical harmonics approaches, the use of a grid based method does not provide a natural framework for the resolution of the exterior problem. However, we show that a Dirichlet-to-Neumann operator at the surface of the Earth provides a simple and efficient solution to this problem. A special attention is needed for the fluid parts to avoid spurious oscillations. To this end, we introduce a general two-potentials formulation which allows to take any density stratification into account. Contrary to other studies that considered the velocity potential, our decomposition is applied to the displacement field in order to obtain natural solid-fluid boundary conditions for the MRP. At each time step, the MRP is computed from the displacement field through a preconditioned conjugate gradient algorithm. This procedure has to be accurate enough in order to ensure a stable calculation on long time series. Some examples of synthetic seismograms computed in spherical Earth models illustrate the accuracy of our approach.
NASA Astrophysics Data System (ADS)
Yanagisawa, Takatoshi; Yamagishi, Yasuko; Hamano, Yozo; Stegman, Dave R.; Suetsugu, Daisuke; Bina, Craig; Inoue, Toru; Wiens, Douglas; Jellinek, Mark
2010-11-01
Seismic tomography reveals the natural mode of convection in the Earth is whole mantle with subducted slabs clearly seen as continuous features into the lower mantle. However, simultaneously existing alongside these deep slabs are stagnant slabs which are, if only temporarily, trapped in the upper mantle. Previous numerical models of mantle convection have observed a range of behavior for slabs in the transition zone depending on viscosity stratification and mineral phase transitions, but typically only exhibit flat-lying slabs when mantle convection is layered or trench migration is imposed. We use 3-D spherical models of mantle convection which range up to Earth-like conditions in Rayleigh number to systematically investigate three effects on mantle dynamics: (1) the mineral phase transitions, (2) a strongly temperature-dependent viscosity with plastic yielding at shallow depth, and (3) a viscosity increase in the lower mantle. First a regime diagram is constructed for isoviscous models over a wide range of Rayleigh number and Clapeyron slope for which the convective mode is determined. It agrees very well with previous results from 2-D simulations by Christensen and Yuen (1985), suggesting present-day Earth is in the intermittent convection mode rather than layered or strictly whole mantle. Two calculations at Earth-like conditions (Ra and RaH = 2 í 107 and 5 í 108, respectively) which include effects (2) and (3) are produced with and without the effect of the mineral phase transitions. The first calculation (without the phase transition) successfully produces plate-like behavior with a long wavelength structure and surface heat flow similar to Earth's value. While the observed convective flow pattern in the lower mantle is broader compared to isoviscous models, it basically shows the behavior of whole mantle convection, and does not exhibit any slab flattening at the viscosity increase at 660 km depth. The second calculation which includes the phase
Strong, Multi-Scale Heterogeneity in Earth's Lowermost Mantle
NASA Astrophysics Data System (ADS)
Tkalčić, Hrvoje; Young, Mallory
2014-05-01
The ~300 km thick layer above the Earth's core mantle boundary remains largely an enigma and has proven to be far more than a simple dividing line; rather it is a complex region with a range of proposed phenomena such as thermal and compositional heterogeneity, partial melting and anisotropy. Characterizing the heterogeneity in the lowermost mantle will prove crucial to accurately understanding key geodynamical processes within our planet. Here we obtain compressional wave (P-wave) velocity images and uncertainty estimates for the lowermost mantle using old and newly collected travel time data sensitive to the lowermost mantle and core and collected by waveform cross-correlation. The images obtained by the inversion technique are void of explicit model parameterization and smoothing. To attest to the impressive capabilities of the transdimensional and hierarchical Bayesian inversion scheme, we design a comprehensive, all-embracing synthetic resolution test demonstrating the retrieval of velocity discontinuities, smooth velocity transitions, structures of varying scales and strengths. Subsequent spectral analyses reveal a power of heterogeneity three times larger than previous estimates and a multi-scale wavelength content in the P-wave velocity field of the lowermost mantle. The newly obtained P-wave tomographic images of the lowermost mantle are not dominated by harmonic degree 2 structure as is the case for tomographic images derived from S-wave data. Instead, the heterogeneity size is more uniformly distributed between about 500 and 6000 km. Inter alia, the resulting heterogeneity spectrum provides a bridge between the long-wavelength features of previous global models and the very short-scale dimensions of scatterers mapped in independent studies. Because the long scale features are less dominant in our model than in S-wave velocity maps, we cannot reasonably determine a correlation between them and the position of detected ultra-low velocity zones.
Constraints on the heterogeneity spectrum of Earth's upper mantle
NASA Astrophysics Data System (ADS)
Mancinelli, Nicholas; Shearer, Peter; Liu, Qinya
2016-05-01
We constrain the heterogeneity spectrum of Earth's upper mantle at scales from a few kilometers to tens of thousands of kilometers using observations from high-frequency scattering, long-period scattering, and tomography. Tomography and high-frequency scattering constraints are drawn from previous studies, but constraints on mantle heterogeneity at intermediate scales (5-500 km) are lacking. To address this, we stack ˜15,000 long-period P coda envelopes to characterize the globally averaged scattered wavefield at periods from 5 to 60 s and at ranges from 50 to 98°. To fit these observations, we consider models of random mantle heterogeneity and compute the corresponding global wavefield using both a ray theoretical "seismic particle" approach and full spectral element simulations. Von Kármán random media distributed throughout the uppermost 600 km of the mantle with a = 2000 km, ɛ = 10%, and κ = 0.05 provide a good fit to the time, range, and frequency dependence of the stacks, although there is a trade-off between ɛ and the thickness of the assumed scattering layer. This random media model also fits previously published 1 Hz stacks of P coda and agrees with constraints on long-wavelength structure from tomography. Finally, we explore geodynamically plausible scenarios that might be responsible for the RMS and falloff rate of the proposed spectrum, including a self-similar mixture of basalt and harzburgite.
Senthil Raja, Duraisamy; Lin, Pin-Chun; Liu, Wei-Ren; Zhan, Jun-Xiang; Fu, Xin-Yi; Lin, Chia-Her
2015-05-01
A series of new alkaline-earth metal diphosphonate frameworks were successfully synthesized under solvothermal reaction condition (160 °C, 3 d) using 1-hydroxyethylidene-1,1-diphosphonic acid (CH3C(OH)(H2PO3)2, hedpH4) as a diphosphonate building block and Mg(II), Ca(II), Sr(II), or Ba(II) ions as alkaline-earth metal ion centers in water, dimethylformamide, and/or EtOH media. These diphosphonate frameworks, (H2NMe2)4[Mg(hedpH2)3]·3H2O (1), (H2NMe2)2[Ca(hedpH2)2] (2), (H2NMe2)2[Sr3(hedpH2)4(H2O)2] (3), and [Ba3(hedpH2)3]·H2O (4) exhibited interesting structural topologies (zero-, one-, two-, and three-dimensional (0D, 1D, 2D, and 3D, respectively)), which are mainly depending on the metal ions and the solvents used in the synthesis. The single-crystal analysis of these newly synthesized compounds revealed that 1 was a 0D molecule, 2 has 1D chains, 3 was a 3D molecule, and 4 has 2D layers. All compounds were further characterized using thermogravimetric analysis, solid-state (31)P NMR, powder X-ray diffraction analysis, UV-vis spectra, and infrared spectroscopy. In addition, Eu(III)- and Tb(III)-doped compounds of 1-4, namely, (H2NMe2)4[Ln(x)Mg(1-x)(hedpH2)2(hedpH(2-x))]·3H2O (1Ln), (H2NMe2)2[Ln(x)Ca(1-x)(hedpH2)(hedpH(2-x))] (2Ln), (H2NMe2)2[Ln(x)Sr(3-x)(hedpH2)3(hedpH(2-x))(H2O)2] (3Ln), and [Ln(x)Ba(3-x)(hedpH2)2(hedpH(2-x))]·H2O (4Ln) (where Ln = Eu, Tb), were synthesized, and their photoluminescence properties were studied. The quantum yield of 1Eu-4Eu was measured with reference to commercial red phosphor, Y2O2S:Eu(3+) (YE), and the quantum yield of terbium-doped compounds 1Tb-4Tb was measured with reference to commercial green-emitting phosphor CeMgAl10O17:Tb(3+). Interestingly, the compound 2Eu showed very high quantum yield of 92.2%, which is better than that of the reference commercial red phosphor, YE (90.8%). PMID:25871285
Fine scale heterogeneity in the Earth's mantle - observation and interpretation
NASA Astrophysics Data System (ADS)
Thybo, H.
2012-12-01
High resolution seismic data has over the last decade provided significant evidence for pronounced fine scale heterogeneity in the Earth's mantle at an unprecedented detail. Seismic tomography developed tremendously during the last 20-30 years. The results show overall structure in the mantle which can be correlated to main plate tectonic features, such as oceanic spreading centres, continental rift zones and subducting slabs. Much seismological mantle research is now concentrated on imaging fine scale heterogeneity, which may be detected and imaged with high-resolution seismic data with dense station spacing and at high frequency, e.g. from the Russian Peaceful Nuclear Explosion (PNE) data set and array recordings of waves from natural seismic sources. Mantle body waves indicate pronounced heterogeneity at three depth levels whereas other depth intervals appear transparent, at least in the frequency band of 0.5-15 Hz: (1) The Mantle Low-Velocity Zone (LVZ) is a global feature which has been detected in more than 50 long-range seismic profiles (Thybo and Perchuc, Science, 1997). Since then numerous studies based on receiver functions, surface waves, and controlled source seismology have confirmed the presence of this zone. The data demonstrates that the top of the LVZ everywhere is at a depth of 100±20 km. A pronounced coda shows that the zone is highly heterogeneous at characteristic scale lengths of 5-15 by 2-6 km. We interpret that the rocks in the LVZ have a temperature close to the solidus or even may contain small fractions of partial melt. The solidus of mantle rocks is very low below a depth of ca. 90 km if volatiles are present due to a characteristic kink in the solidus which is much lower than for dry mantle rocks. We suggest that the rocks are in a totally solid state below the LVZ and that the depth to the interface to fully solid rocks is an indicator of the thermal state of the upper mantle. (2) Significant scattering from around the top of the
Nakaguchi, Yuji; Ono, Takeshi; Maruyama, Masato; Nagasue, Nozomu; Shimohigashi, Yoshinobu; Kai, Yudai
2015-01-01
In this study, we evaluated the performance of a three-dimensional (3D) dose verification system, COMPASS version 3, which has a dedicated beam models and dose calculation engine. It was possible to reconstruct the 3D dose distributions in patient anatomy based on the measured fluence using the MatriXX 2D array. The COMPASS system was compared with Monte Carlo simulation (MC), glass rod dosimeter (GRD), and 3DVH, using an anthropomorphic phantom for intensity-modulated radiation therapy (IMRT) dose verification in clinical neck cases. The GRD measurements agreed with the MC within 5% at most measurement points. In addition, most points for COMPASS and 3DVH also agreed with the MC within 5%. The COMPASS system showed better results than 3DVH for dose profiles due to individual adjustments, such as beam modeling for each linac. Regarding the dose-volume histograms, there were no large differences between MC, analytical anisotropic algorithm (AAA) in Eclipse treatment planning system (TPS), 3DVH, and the COMPASS system. However, AAA underestimated the dose to the clinical target volume and Rt-Parotid slightly. This is because AAA has some problems with dose calculation accuracy. Our results indicated that the COMPASS system offers highly accurate 3D dose calculation for clinical IMRT quality assurance. Also, the COMPASS system will be useful as a commissioning tool in routine clinical practice for TPS. PMID:25679177
NASA Astrophysics Data System (ADS)
Toledo-Redondo, S.; Salinas, A.; Fornieles, J.; Portí, J.; Lichtenegger, H. I. M.
2016-06-01
Schumann resonances can be found in planetary atmospheres, inside the cavity formed by the conducting surface of the planet and the lower ionosphere. They are a powerful tool to investigate both the electric processes that occur in the atmosphere and the characteristics of the surface and the lower ionosphere. Results from a full 3-D model of the Earth-ionosphere electromagnetic cavity based on the Transmission-Line Modeling (TLM) method are presented. A Cartesian scheme with homogeneous cell size of 10 km is used to minimize numerical dispersion present in spherical schemes. Time and frequency domain results have been obtained to study the resonance phenomenon. The effect of conductivity on the Schumann resonances in the cavity is investigated by means of numerical simulations, studying the transition from resonant to nonresonant response and setting the conductivity limit for the resonances to develop inside the cavity. It is found that the transition from resonant to nonresonant behavior occurs for conductivity values above roughly 10-9 S/m. For large losses in the cavity, the resonances are damped, but, in addition, the peak frequencies change according to the local distance to the source and with the particular electromagnetic field component. These spatial variations present steep variations around each mode's nodal position, covering distances around 1/4 of the mode wavelength, the higher modes being more sensitive to this effect than the lower ones. The dependence of the measured frequency on the distance to the source and particular component of the electric field offers information on the source generating these resonances.
NASA Astrophysics Data System (ADS)
Hazenberg, P.; Broxton, P. D.; Brunke, M.; Gochis, D. J.; Lawrence, D. M.; Niu, G. Y.; Pelletier, J. D.; Troch, P. A. A.; Zeng, X.
2014-12-01
Traditionally, LSMs for use in Earth system models (ESMs) only account for one-dimensional (1-D) vertical hydrological processes at very coarse resolutions (~100 km). Recently, there has been interest in simulating lateral exchange of surface and subsurface water, as the grid resolution of ESMs increases (currently towards ~20 km) due to advances in computational power.In the current work, we present a new physically-based hydrological model capable of simulating lateral flow at a hyperresolution (1 km pixels) coupled with the vertical soil column of the Community Land Model (CLM), which has a much coarser resolution. Our new approach also uses sub-pixel topographic information to represent small scale lateral-flow processes. As input to our model, we use 1 km data of required surface and subsurface information, such depth to bedrock, hillslope width functions, fractional coverage of hillslopes/wetlands/riparian zones within each 1 km pixel, and a river network. Such data are largely derived from higher-resolution (30 m) topographic data, and will be made available to the community as global products.We have tested our model against measurements over a well instrumented, artificial hillslope at the University of Arizona's Biosphere 2 and found that both our approach and a full 3-D physically-based hydrological model can realistically simulate hydrological states and fluxes. However, computationally, our approach is 2 - 3 orders of magnitude faster than the latter. We are currently testing it for an ESM grid box that covers much of southeastern New York State, which includes the headwaters of the Delaware River, the Susquehana River, and parts of the Mohawk River Basin. Our approach improves upon the 1-D hydrological representation in CLM by representing the subgrid topographic and geomorphological variability in the grid box. The impact of representing such subgrid variability on surface water and energy fluxes will be discussed in our presentation.
Predicting lower mantle heterogeneity from 4-D Earth models
NASA Astrophysics Data System (ADS)
Flament, Nicolas; Williams, Simon; Müller, Dietmar; Gurnis, Michael; Bower, Dan J.
2016-04-01
basal layer ˜ 4% denser than ambient mantle. Increasing convective vigour (Ra ≈ 5 x 108) or decreasing the density of the basal layer decreases both the accuracy and sensitivity of the predicted lower mantle structure. References: D. J. Bower, M. Gurnis, N. Flament, Assimilating lithosphere and slab history in 4-D Earth models. Phys. Earth Planet. Inter. 238, 8-22 (2015). V. Lekic, S. Cottaar, A. Dziewonski, B. Romanowicz, Cluster analysis of global lower mantle tomography: A new class of structure and implications for chemical heterogeneity. Earth Planet. Sci. Lett. 357, 68-77 (2012).
NASA Astrophysics Data System (ADS)
Matos, J. R.; Welty, C.; Packman, A.
2005-12-01
The main purpose of the simulations in this research is the analysis of three-dimensional surface-groundwater interchange in heterogeneous systems. The effects of channel pattern, bed forms and aquifer heterogeneity on flow interactions between stream and groundwater systems are examined in order to contribute for a better understanding of the hyporheic process. A two-dimensional approach was also adopted to allow comparisons with the three-dimensional results. The grid was designed using the correlation scales of the heterogeneous fields and the scale of the stream meanders. MODFLOW and MODPATH were used to evaluate magnitude, direction and spatial distribution of the exchange flow. PMWIN and PMPATH were used as pre and post-processors during the construction of the models and analysis of results. Gaining and losing streams as well as parallel flow and flow across streams were simulated as idealized cases intended to describe how properties of the streambed and aquifer in low-gradient lowland streams contribute to hyporheic exchange. At first a straight river was analyzed then meandering streams were created with a sine curve and variations on wavelength and amplitude. Bed forms were simulated assuming a sinusoidal distribution of pressure head in the bed surface. Aspects of the influence of bedforms on mechanisms such as "pumping" and "turnover" are expected to be addressed with simulations. Flow velocities between 20 and 40 cm/s in the channel were tested with the objective of showing the influence of river morphology and natural bed forms on the flow exchange in the hyporheic zone. Several meander cycles and four levels of hydraulic conductivity variance were analyzed. Results of flow variances along the cross-sections and wetted perimeter show the increasing on hyporheic exchange as the degree of heterogeneity increases. Particle tracking was performed to define hyporheic residence time distributions. When comparing the homogeneous fields with all degrees of
Smith, Stuart J.; Wilson, Martin; Ward, Jennifer H.; Rahman, Cheryl V.; Peet, Andrew C.; Macarthur, Donald C.; Rose, Felicity R. A. J.; Grundy, Richard G.; Rahman, Ruman
2012-01-01
Introduction Physiologically relevant pre-clinical ex vivo models recapitulating CNS tumor micro-environmental complexity will aid development of biologically-targeted agents. We present comprehensive characterization of tumor aggregates generated using the 3D Rotary Cell Culture System (RCCS). Methods CNS cancer cell lines were grown in conventional 2D cultures and the RCCS and comparison with a cohort of 53 pediatric high grade gliomas conducted by genome wide gene expression and microRNA arrays, coupled with immunohistochemistry, ex vivo magnetic resonance spectroscopy and drug sensitivity evaluation using the histone deacetylase inhibitor, Vorinostat. Results Macroscopic RCCS aggregates recapitulated the heterogeneous morphology of brain tumors with a distinct proliferating rim, necrotic core and oxygen tension gradient. Gene expression and microRNA analyses revealed significant differences with 3D expression intermediate to 2D cultures and primary brain tumors. Metabolic profiling revealed differential profiles, with an increase in tumor specific metabolites in 3D. To evaluate the potential of the RCCS as a drug testing tool, we determined the efficacy of Vorinostat against aggregates of U87 and KNS42 glioblastoma cells. Both lines demonstrated markedly reduced sensitivity when assaying in 3D culture conditions compared to classical 2D drug screen approaches. Conclusions Our comprehensive characterization demonstrates that 3D RCCS culture of high grade brain tumor cells has profound effects on the genetic, epigenetic and metabolic profiles of cultured cells, with these cells residing as an intermediate phenotype between that of 2D cultures and primary tumors. There is a discrepancy between 2D culture and tumor molecular profiles, and RCCS partially re-capitulates tissue specific features, allowing drug testing in a more relevant ex vivo system. PMID:23272238
NASA Astrophysics Data System (ADS)
Mousavi, Sima; Bauer, Klaus; Korn, Michael
2014-05-01
3-D Vp and Vp/Vs structure of the geodynamically active NW Bohemia/Vogtland area, located at the border region between Germany and Czech republic, has been determined from local earthquake tomography using 543 earthquakes which have been recorded during 2000 to 2010. This region is known for the occurrence of earthquake swarms that are supposed to be triggered by fluid upwelling in the crust, although fluid behaviour and migration paths in the subsurface of NW Bohemia is still poorly known. The events used in this study were selected based on a minimum 12 P and S phase observations and an azimuthal gap less than 160º. This data set is employed to derive a minimum 1-D velocity model and to relocate the hypocenters. The minimum 1-D velocity model is then used as an initial model in non-linear inversion to derive 3-D P-velocity and Vp/Vs ratio. Using synthetic tests, it can be shown that a high resolution is obtained in the central part of the studied region with the given source and receiver configuration. Two branches of high Vp/Vs ratio anomalies have been detected above the swarm quakes' focal zone. These anomalies support the existence of two main fluid passages toward Bad Brambach and Bublak moffette. Another interesting result is a high Vp/Vs line-like anomaly along Mariánské Lázně fault where most of the swarm quakes occur, which could be due to a fluid saturated area around the cracked zone of the fault plain. Hypocenters in the swarm region are located in a low Vp and Vp/Vs anomaly. The correlation between the detected Vp and Vp/Vs anomalies and the location of earthquake swarm suggests a model in which CO2 as part of magmatic fluids exist in a vast area beneath NW Bohemia and frequently migrate up to the surface.
Guerin, P.; Baudron, A. M.; Lautard, J. J.
2006-07-01
This paper describes a new technique for determining the pin power in heterogeneous core calculations. It is based on a domain decomposition with overlapping sub-domains and a component mode synthesis technique for the global flux determination. Local basis functions are used to span a discrete space that allows fundamental global mode approximation through a Galerkin technique. Two approaches are given to obtain these local basis functions: in the first one (Component Mode Synthesis method), the first few spatial eigenfunctions are computed on each sub-domain, using periodic boundary conditions. In the second one (Factorized Component Mode Synthesis method), only the fundamental mode is computed, and we use a factorization principle for the flux in order to replace the higher order Eigenmodes. These different local spatial functions are extended to the global domain by defining them as zero outside the sub-domain. These methods are well-fitted for heterogeneous core calculations because the spatial interface modes are taken into account in the domain decomposition. Although these methods could be applied to higher order angular approximations - particularly easily to a SPN approximation - the numerical results we provide are obtained using a diffusion model. We show the methods' accuracy for reactor cores loaded with UOX and MOX assemblies, for which standard reconstruction techniques are known to perform poorly. Furthermore, we show that our methods are highly and easily parallelizable. (authors)
NASA Astrophysics Data System (ADS)
Sinclair Yemini, Francis; Chenu, Claire; Monga, Olivier; Vieuble Gonond, Laure; Juarez, Sabrina; Pihneiro, Marc; otten, Wilfred; Garnier, Patricia
2014-05-01
Contaminant degradation by microorganisms is very variable in soils because of the very heterogeneous spatial relationship of contaminant/degraders. Repacked Soil columns were carried out to study the degradation of 2,4D pesticide labelled with C14 for different scenarios of microorganisms and pesticide initial location. Measurements of global C14-CO2 emission and C14 distribution in the soil column showed that the initial location play a crucial rule on the dissipation of the pollutant. Experiments were simulated using a 3D model able to model microbial degradation and substrate diffusion between aggregates by considering explicitly the 3D structure of soil from CT images. The initial version of the model (Monga et al., 2008) was improved in order to simulate diffusion in samples of large size. Partial differential equations were implemented using freefem++ solver. The model simulates properly the dynamics of 2,4D in the column for the different initial situations. CT images of the same soil but using undisturbed structure instead of repacked aggregates were also carried out. Significant differences of the simulated results were observed between the repacked and the undisturbed soil. The conclusion of our work is that the heterogeneity of the soil structure and location of pollutants and decomposers has a very strong influence on the dissipation of pollutants.
Fine-scale heterogeneity in the Earth's inner core
Vidale; Earle
2000-03-16
The seismological properties of the Earth's inner core have become of particular interest as we understand more about its composition and thermal state. Observations of anisotropy and velocity heterogeneity in the inner core are beginning to reveal how it has grown and whether it convects. The attenuation of seismic waves in the inner core is strong, and studies of seismic body waves have found that this high attenuation is consistent with either scattering or intrinsic attenuation. The outermost portion of the inner core has been inferred to possess layering and to be less anisotropic than at greater depths. Here we present observations of seismic waves scattered in the inner core which follow the expected arrival time of the body-wave reflection from the inner-core boundary. The amplitude of these scattered waves can be explained by stiffness variations of 1.2% with a scale length of 2 kilometres across the outermost 300 km of the inner core. These variations might be caused by variations in composition, by pods of partial melt in a mostly solid matrix or by variations in the orientation or strength of seismic anisotropy. PMID:10749207
Samluk, Jesse P.; Geiger, Cathleen A.; Weiss, Chester J.; Kolodzey, James
2015-10-01
In this article we explore simulated responses of electromagnetic (EM) signals relative to in situ field surveys and quantify the effects that different values of conductivity in sea ice have on the EM fields. We compute EM responses of ice types with a three-dimensional (3-D) finite-volume discretization of Maxwell's equations and present 2-D sliced visualizations of their associated EM fields at discrete frequencies. Several interesting observations result: First, since the simulator computes the fields everywhere, each gridcell acts as a receiver within the model volume, and captures the complete, coupled interactions between air, snow, sea ice and sea water as a function of their conductivity; second, visualizations demonstrate how 1-D approximations near deformed ice features are violated. But the most important new finding is that changes in conductivity affect EM field response by modifying the magnitude and spatial patterns (i.e. footprint size and shape) of current density and magnetic fields. These effects are demonstrated through a visual feature we define as 'null lines'. Null line shape is affected by changes in conductivity near material boundaries as well as transmitter location. Our results encourage the use of null lines as a planning tool for better ground-truth field measurements near deformed ice types.
Samluk, Jesse P.; Geiger, Cathleen A.; Weiss, Chester J.; Kolodzey, James
2015-10-01
In this article we explore simulated responses of electromagnetic (EM) signals relative to in situ field surveys and quantify the effects that different values of conductivity in sea ice have on the EM fields. We compute EM responses of ice types with a three-dimensional (3-D) finite-volume discretization of Maxwell's equations and present 2-D sliced visualizations of their associated EM fields at discrete frequencies. Several interesting observations result: First, since the simulator computes the fields everywhere, each gridcell acts as a receiver within the model volume, and captures the complete, coupled interactions between air, snow, sea ice and sea water asmore » a function of their conductivity; second, visualizations demonstrate how 1-D approximations near deformed ice features are violated. But the most important new finding is that changes in conductivity affect EM field response by modifying the magnitude and spatial patterns (i.e. footprint size and shape) of current density and magnetic fields. These effects are demonstrated through a visual feature we define as 'null lines'. Null line shape is affected by changes in conductivity near material boundaries as well as transmitter location. Our results encourage the use of null lines as a planning tool for better ground-truth field measurements near deformed ice types.« less
A Global 3D P-Velocity Model of the Earth's Crust and Mantle for Improved Event Location
NASA Astrophysics Data System (ADS)
Ballard, S.; Young, C. J.; Hipp, J. R.; Chang, M.; Lewis, J.; Begnaud, M. L.; Rowe, C. A.
2009-12-01
further refinement takes place around adjusted nodes to form a new model, and the process is repeated until no more improvement can be obtained. We thus produce a smooth, multi-resolution model with node density appropriate to both ray coverage and the velocity gradients required by the data. This scheme is computationally expensive, so we use a Java-based distributed computing framework developed by Sandia National Laboratories (SNL), providing us with 300+ processors having an efficiency of better than 90% for the calculations. We evaluate our model both in terms of travel time residual variance reduction and in location improvement for GT events. For the latter, we use a new multi-threaded version of the SNL-developed LocOO code modified to use 3D velocity models.
ERIC Educational Resources Information Center
Isik-Ercan, Zeynep; Zeynep Inan, Hatice; Nowak, Jeffrey A.; Kim, Beomjin
2014-01-01
This qualitative case study describes (a) the ways 3D visualization, coupled with other science and literacy experiences, supported young children's first exploration of the Earth-Sun-Moon system and (b) the perspectives of classroom teachers and children on using 3D visualization. We created three interactive 3D software modules that…
ERIC Educational Resources Information Center
Isik-Ercan, Zeynep; Zeynep Inan, Hatice; Nowak, Jeffrey A.; Kim, Beomjin
2014-01-01
This qualitative case study describes (a) the ways 3D visualization, coupled with other science and literacy experiences, supported young children's first exploration of the Earth-Sun-Moon system and (b) the perspectives of classroom teachers and children on using 3D visualization. We created three interactive 3D software modules that simulate day…
Re-Dimensional Thinking in Earth Science: From 3-D Virtual Reality Panoramas to 2-D Contour Maps
ERIC Educational Resources Information Center
Park, John; Carter, Glenda; Butler, Susan; Slykhuis, David; Reid-Griffin, Angelia
2008-01-01
This study examines the relationship of gender and spatial perception on student interactivity with contour maps and non-immersive virtual reality. Eighteen eighth-grade students elected to participate in a six-week activity-based course called "3-D GeoMapping." The course included nine days of activities related to topographic mapping. At the end…
NASA Astrophysics Data System (ADS)
Monge-Sanz, B. M.; Chipperfield, M. P.; Cariolle, D.; Feng, W.
2010-05-01
A detailed full-chemistry 3-D chemistry and transport model (CTM) is used to evaluate the current stratospheric O3 parameterisation in the European Centre for Medium-Range Weather Forecasts (ECMWF) model and to obtain an alternative version of the ozone scheme implicitly including heterogeneous chemistry. The approach avoids the inaccurate treatment currently given to heterogeneous ozone chemistry in the ECMWF model, as well as the uncertainties of a cold-tracer. The new O3 scheme (COPCAT) is evaluated within the same CTM used to calculate it. It is the first time such a comparison has been possible, providing direct information on the validity of the linear parameterisation approach for stratospheric ozone. Simulated total column and O3 profiles are compared against Total Ozone Mapping Spectrometer (TOMS) and Halogen Occultation Experiment (HALOE) observations. COPCAT successfully simulates polar loss and reproduces a realistic Antarctic O3 hole. The new scheme is comparable to the full-chemistry in many regions for multiannual runs. The parameterisation produces less ozone over the tropics around 10 hPa, compared to full-chemistry and observations, however, this problem can be ameliorated by choosing a different ozone climatology for the scheme. The new scheme is compared to the current ECMWF scheme in the same CTM runs. The Antarctic O3 hole with the current ECMWF scheme is weaker and disappears earlier than with the new COPCAT scheme. Differences between the current ECMWF scheme and COPCAT are difficult to explain due to the different approach used for heterogeneous chemistry and differences in the photochemical models used to calculate the scheme coefficients. Results with the new COPCAT scheme presented here show that heterogeneous and homogeneous ozone chemistry can be included in a consistent way in a linear ozone parameterisation, without any additional tunable parameters, providing a parameterisation scheme in much better agreement with the current
NASA Astrophysics Data System (ADS)
Monge-Sanz, B. M.; Chipperfield, M. P.; Cariolle, D.; Feng, W.
2011-02-01
A detailed full-chemistry 3-D chemistry and transport model (CTM) is used to evaluate the current stratospheric O3 parameterisation in the European Centre for Medium-Range Weather Forecasts (ECMWF) model and to obtain an alternative version of the ozone scheme implicitly including heterogeneous chemistry. The approach avoids the inaccurate treatment currently given to heterogeneous ozone chemistry in the ECMWF model, as well as the uncertainties of a cold-tracer. The new O3 scheme (COPCAT) is evaluated within the same CTM used to calculate it. It is the first time such a comparison has been possible, providing direct information on the validity of the linear parameterisation approach for stratospheric ozone. Simulated total column and O3 profiles are compared against Total Ozone Mapping Spectrometer (TOMS) and Halogen Occultation Experiment (HALOE) observations. COPCAT successfully simulates polar loss and reproduces a realistic Antarctic O3 hole. The new scheme is comparable to the full-chemistry in many regions for multiannual runs. The parameterisation produces less ozone over the tropics around 10 hPa, compared to full-chemistry and observations. However, this problem can be ameliorated by choosing a different ozone climatology for the scheme. The new scheme is compared to the current ECMWF scheme in the same CTM runs. The Antarctic O3 hole with the current ECMWF scheme is weaker and disappears earlier than with the new COPCAT scheme. Differences between the current ECMWF scheme and COPCAT are difficult to explain due to the different approach used for heterogeneous chemistry and differences in the photochemical models used to calculate the scheme coefficients. Results with the new COPCAT scheme presented here show that heterogeneous and homogeneous ozone chemistry can be included in a consistent way in a linear ozone parameterisation, without any additional tunable parameters, providing a parameterisation scheme in better agreement with the current knowledge
NASA Astrophysics Data System (ADS)
Qu, Wulin; Zhang, Bei; Huang, Luyuan; Cheng, Huihong; Shi, Yaolin
2016-04-01
The 26 December 2004 Sumatra-Andaman earthquake with moment magnitude (Mw) of 9.1 to 9.3 is the first great earthquake recorded by digital broadband, high-dynamic-range seismometers and global positioning system (GPS) equipment, which recorded many high-quality geophysical data sets. The spherical curvature is not negligible in far field especially for large event and the real Earth is laterally inhomogeneity and the analytical results still are difficult to explain the geodetic measurements. We use equivalent body force finite elements method Zhang et al. (2015) and mesh the whole earth, to compute global co-seismic displacements using four fault slip models of the 2004 Sumatra earthquake provided by different authors. Comparisons of calculated co-seismic displacements and GPS show that the confidences are well in near field for four models, and the confidences are according to different models. In the whole four models, the Chlieh model (Chlieh et al., 2007) is the best as this slip model not only accord well with near field data but also far field data. And then we use the best slip model, Chlieh model to explore influence of three dimensional lateral earth structure on both layered spherically symmetric (PREM) and real 3-D heterogeneous earth model (Crust 1.0 model and GyPSuM). Results show that the effects of 3-D heterogeneous earth model are not negligible and decrease concomitantly with increasing distance from the epicenter. The relative effects of 3-D crust model are 23% and 40% for horizontal and vertical displacements, respectively. The effects of the 3-D mantle model are much smaller than that of 3-D crust model but with wider impacting area.
NASA Astrophysics Data System (ADS)
Ballard, S.; Begnaud, M. L.; Hipp, J. R.; Chael, E. P.; Encarnacao, A.; Maceira, M.; Yang, X.; Young, C. J.; Phillips, W.
2013-12-01
SALSA3D is a global 3D P wave velocity model of the Earth's crust and mantle developed specifically to provide seismic event locations that are more accurate and more precise than are locations from 1D and 2.5D models. In this paper, we present the most recent version of our model, for the first time jointly derived from multiple types of data: body wave travel times, surface wave group velocities, and gravity. The latter two are added to provide information in areas with poor body wave coverage, and are down-weighted in areas where body wave coverage is good. To constrain the inversions, we invoked empirical relations among the density, S velocity, and P velocity. We demonstrate the ability of the new SALSA3D model to reduce mislocations and generate statistically robust uncertainty estimates for a large set of realizations derived from a carefully chosen set of globally-distributed ground truth events. We obtain path-dependent travel time prediction uncertainties for our model by computing the full 3D model covariance matrix of our tomographic system and integrating the model slowness variance and covariance along paths of interest. This approach yields very low travel time prediction uncertainties for well-sampled paths through the Earth and higher uncertainties for paths that are poorly represented in the data set used to develop the model. While the calculation of path-dependent prediction uncertainties with this approach is computationally expensive, uncertainties can be pre-computed for a network of stations and stored in 3D lookup tables that can be quickly and efficiently interrogated using GeoTess software.
NASA Astrophysics Data System (ADS)
Viscardy, S.; Errera, Q.; Pitts, M. C.; Daerden, F.
2012-04-01
A 3-D Chemical Transport Model (CTM), with full stratospheric chemistry and driven by the ECMWF temperature and wind fields, is coupled to the PSC microphysical model PSCBox. This interactively describes the formation and evolution of four types of PSC particles (STS, SAT, NAT, and ice) through relevant microphysical processes. The number density and composition of each type of particles are computed for a binned size distribution. As a result, the calculation of surface area densities is accurately performed, of which the computation of the heterogeneous reaction constants takes advantage. The explicit computation of the particle size distributions allows obtaining the same optical properties as those measured by CALIPSO. Hence, the evolution of PSC coverage and composition will be studied and compared to the CALIPSO observations during Antarctic winters. The relationship between the presence of PSCs and the heterogeneous chemistry will also be investigated. In particular, two issues will be considered: (i) how the ozone depletion is related to the PSC coverage, and (ii) how each PSC particle type contributes to the chlorine activation.
NASA Astrophysics Data System (ADS)
Jackson, B. V.; Hamilton, M. S.; Hick, P. P.; Buffington, A.; Bisi, M. M.; Clover, J. M.; Tokumaru, M.; Fujiki, K.
2011-07-01
SMEI and IPS remotely observe increased brightness and velocity enhancements behind interplanetary shocks that are also seen in situ. We use the UCSD time-dependent 3-D reconstruction technique to map these enhancements, and compare them with measurements at the SOHO, Wind, ACE, and STEREO spacecraft. The analyses of these shocks from hour-averaged in-situ data show that the enhanced density column associated with the shock response varies considerably between different instruments, even for in-situ instruments located at L1 near Earth. The relatively-low-resolution SMEI 3-D reconstructions generally show density enhancements, and within errors, the column excesses match those observed in situ. In these SMEI 3-D reconstructions from remotely-sensed data, the shock density enhancements appear not as continuous broad fronts, but as segmented structures. This may provide part of the explanation for the observed discrepancies between the various in-situ measurements at Earth and STEREO, but not between individual instruments near L1.
NASA Astrophysics Data System (ADS)
Kobayashi, H.; Suzuki, R.; Nagai, S.; Nakai, T.; Kim, Y.
2012-12-01
Over the last couple of decades, the three dimensional plant canopy radiative transfer models have been developed, improved and used for the retrievals of biophysical variables of vegetative surface. Fraction of absorbed photosynthetically active radiation (FAPAR) by plant canopy, a similar variable to heating rate in the atmosphere, is one of the important biophysical variables to infer the terrestrial plant canopy photosynthesis. FAPAR can be estimated by the radiative transfer model inversion or the empirical relationships between FAPAR and vegetation indices such as normalized difference vegetation index (NDVI). To date, some global FAPAR products are publicly available. These products are estimated from the moderate resolution satellites such as MODIS and SPOT-VEGETATION. One may apply the similar FAPAR algorithms to higher spatial resolution satellites if the ecosystem structures are horizontally homogeneous, which means that the adjacent satellite pixels have a similar spectral properties. If the vegetation surface is highly heterogeneous, "domain average FAPAR", which assumes no net horizontal radiation fluxes, can be unrealistically high (more than 1). In this presentation, we analyzed the characteristics of FAPAR in a heterogeneous landscape. As a case study, we selected our study site in a sparse black spruce forest in Alaska. We conducted the field campaigns to measure forest structural and optical properties that are used in the radiative transfer simulation. We used a 3D radiative transfer, FLiES (Kobayashi, H. and H. Iwabuchi (2008), A coupled 1-D atmosphere and 3-D canopy radiative transfer model for canopy reflectance, light environment, and photosynthesis simulation in a heterogeneous landscape, Remote Sensing of Environment, 112, 173-185) to create a high resolution simulated spectral reflectance and FAPAR images over the course of the growing season. From the analysis, we show (1) FAPAR with no net horizontal fluxes assumption can be higher than
NASA Astrophysics Data System (ADS)
Barrash, W.; Cardiff, M. A.; Kitanidis, P. K.
2012-12-01
The distribution of hydraulic conductivity (K) is a major control on groundwater flow and contaminant transport. Our limited ability to determine 3D heterogeneous distributions of K is a major reason for increased costs and uncertainties associated with virtually all aspects of groundwater contamination management (e.g., site investigations, risk assessments, remediation method selection/design/operation, monitoring system design/operation). Hydraulic tomography (HT) is an emerging method for directly estimating the spatially variable distribution of K - in a similar fashion to medical or geophysical imaging. Here we present results from 3D transient field-scale experiments (3DTHT) which capture the heterogeneous K distribution in a permeable, moderately heterogeneous, coarse fluvial unconfined aquifer at the Boise Hydrogeophysical Research Site (BHRS). The results are verified against high-resolution K profiles from multi-level slug tests at BHRS wells. The 3DTHT field system for well instrumentation and data acquisition/feedback is fully modular and portable, and the in-well packer-and-port system is easily assembled and disassembled without expensive support equipment or need for gas pressurization. Tests are run for 15-20 min and the aquifer is allowed to recover while the pumping equipment is repositioned between tests. The tomographic modeling software developed uses as input observations of temporal drawdown behavior from each of numerous zones isolated in numerous observation wells during a series of pumping tests conducted from numerous isolated intervals in one or more pumping wells. The software solves for distributed K (as well as storage parameters Ss and Sy, if desired) and estimates parameter uncertainties using: a transient 3D unconfined forward model in MODFLOW, the adjoint state method for calculating sensitivities (Clemo 2007), and the quasi-linear geostatistical inverse method (Kitanidis 1995) for the inversion. We solve for K at >100,000 sub-m3
3D Dynamics of Magnetic Flux Ropes Across Scales: Solar Eruptions and Sun-Earth Plasma Coupling
NASA Astrophysics Data System (ADS)
Chen, James
2012-10-01
Central to the understanding of the eruptive phenomena on the Sun and their impact on the terrestrial plasma environment is the dynamics of coronal mass ejections (CMEs)---a 3D magnetic flux rope configuration---and the evolution of their magnetic fields. I will discuss the basic physics of CME eruption and solar flare energy release in the context of the analytic erupting flux rope model of CMEs. In this ideal MHD model, a CME is treated as a 3D flux rope with its two stationary footpoints anchored in the Sun. The model structure is non-axisymmetric and embedded in a model corona/solar wind. The initial flux rope is driven out of equilibrium by ``injection'' of poloidal flux and propagates under the Lorentz hoop force from the Sun to 1 AU, across a wide range of spatial and temporal scales. Comparisons of the model results and recent STEREO observations show that the solutions that best fit the observed CME position-time data (to within 1-2% of data) also correctly replicate the temporal profiles of associated flare X-ray emissions (GOES data) and the in situ magnetic field and plasma data of the CME ejecta at 1 AU where such data are available (e.g., ACE and STEREO/IMPAXCT/PLASTIC data), providing a unified basis of understanding CME dynamics and flare energetics.
Time scales and heterogeneous structure in geodynamic earth models
Bunge; Richards; Lithgow-Bertelloni; Baumgardner; Grand; Romanowicz
1998-04-01
Computer models of mantle convection constrained by the history of Cenozoic and Mesozoic plate motions explain some deep-mantle structural heterogeneity imaged by seismic tomography, especially those related to subduction. They also reveal a 150-million-year time scale for generating thermal heterogeneity in the mantle, comparable to the record of plate motion reconstructions, so that the problem of unknown initial conditions can be overcome. The pattern of lowermost mantle structure at the core-mantle boundary is controlled by subduction history, although seismic tomography reveals intense large-scale hot (low-velocity) upwelling features not explicitly predicted by the models. PMID:9525864
Sharifi, Tiva; Gracia-Espino, Eduardo; Jia, Xueen; Sandström, Robin; Wågberg, Thomas
2015-12-30
We report efficient electrolysis of both water-splitting half reactions in the same medium by a bifunctional 3D electrode comprising Co3O4 nanospheres nucleated on the surface of nitrogen-doped carbon nanotubes (NCNTs) that in turn are grown on conductive carbon paper (CP). The resulting electrode exhibits high stability and large electrochemical activity for both oxygen and hydrogen evolution reactions (OER and HER). We obtain a current density of 10 mA/cm(2) in 0.1 M KOH solution at overpotentials of only 0.47 and 0.38 V for OER and HER, respectively. Additionally, the experimental observations are understood and supported by analyzing the Co3O4:NCNT and NCNT:CP interfaces by ab initio calculations. Both the experimental and the theoretical studies indicate that firm and well-established interfaces along the electrode play a crucial role on the stability and electrochemical activity for both OER and HER. PMID:26629887
Wei, Xiao-Hua; Yang, Lin-Yan; Liao, Sheng-Yun; Zhang, Ming; Tian, Jin-Lei; Du, Pei-Yao; Gu, Wen; Liu, Xin
2014-04-21
A series of metal-organic framework {Ln(BCPBA)(H2O)}n {Ln = Nd (1), Sm (2), Eu (3), Tb (4), Dy (5)}; {[Ln(BCPBA)(H2O)](H2O)}n {Ln = Pr (6), Gd (7)} have been synthesized through the hydrothermal synthesis method. These compounds possess non-interpenetrating 3D networks with 10.1438 Å× 17.9149 Å rhombic channels along the [001] direction. The results of temperature-dependent magnetic susceptibility measurements indicate that compounds 4 and 7 exhibit Ln(III)Ln(III) antiferromagnetic interactions, while compound 5 exhibits Ln(III)Ln(III) ferromagnetic interactions. Frequency dependent out-of-phase signals were observed in alternating current (ac) magnetic susceptibility measurements which indicate that they have slow magnetic relaxation characteristics. The luminescent properties of 1, 2, 3, 4, and 5 are also discussed. Due to the good match between the lowest triplet state of the ligand and the resonant energy level of the lanthanide ion, compound 4 has longer fluorescence lifetime (τ1 = 400.0000 ms, τ2 = 1143.469 ms) and higher quantum yield (Φ = 42%) compared with other compounds. PMID:24572766
NASA Astrophysics Data System (ADS)
Hazenberg, P.; Broxton, P. D.; Brunke, M.; Gochis, D.; Niu, G. Y.; Pelletier, J. D.; Troch, P. A. A.; Zeng, X.
2015-12-01
The terrestrial hydrological system, including surface and subsurface water, is an essential component of the Earth's climate system. Over the past few decades, land surface modelers have built one-dimensional (1D) models resolving the vertical flow of water through the soil column for use in Earth system models (ESMs). These models generally have a relatively coarse model grid size (~25-100 km) and only account for sub-grid lateral hydrological variations using simple parameterization schemes. At the same time, hydrologists have developed detailed high-resolution (~0.1-10 km grid size) three dimensional (3D) models and showed the importance of accounting for the vertical and lateral redistribution of surface and subsurface water on soil moisture, the surface energy balance and ecosystem dynamics on these smaller scales. However, computational constraints have limited the implementation of the high-resolution models for continental and global scale applications. The current work presents a hybrid-3D hydrological approach is presented, where the 1D vertical soil column model (available in many ESMs) is coupled with a high-resolution lateral flow model (h2D) to simulate subsurface flow and overland flow. H2D accounts for both local-scale hillslope and regional-scale unconfined aquifer responses (i.e. riparian zone and wetlands). This approach was shown to give comparable results as those obtained by an explicit 3D Richards model for the subsurface, but improves runtime efficiency considerably. The h3D approach is implemented for the Delaware river basin, where Noah-MP land surface model (LSM) is used to calculated vertical energy and water exchanges with the atmosphere using a 10km grid resolution. Noah-MP was coupled within the WRF-Hydro infrastructure with the lateral 1km grid resolution h2D model, for which the average depth-to-bedrock, hillslope width function and soil parameters were estimated from digital datasets. The ability of this h3D approach to simulate
NASA Astrophysics Data System (ADS)
Lucas, Yann; Younes, Anis; Ackerer, Philippe
2014-05-01
Knowing the distribution of fresh and saline groundwater is an important issue of the sustainable management of water resources. Here we deal with groundwater plume pollution coming from brines of salted waste heaps in the Upper Rhine aquifer, one of the most important groundwater resources for Western Europe. It has been contaminated in its upstream part by mining brine infiltration coming from potash extraction which started at the beginning of the twentieth century. The salty residues were stored in waste heaps near the mine shafts, and, due to rainwater percolation through the heaps, saline aqueous solutions mixed with groundwater and then progressed in the flow direction, forming high salinity plumes. In the present case, it is proved that density effects are all the more necessary to be considered in a modeling approach that brines density can reach high values just below waste heaps. For this purpose a 3D model is developed using the SEAWAT option of Modflow in a 30 layers model. Simulations are performed on about one hundred years, where chlorine concentration data on the last thirty years are used to constrain the model. In our case, conjugated to a heckled substratum topography, heterogeneous hydrodynamical conditions and the fact that pollution history alternates high pollution phases with relatively low pollution phases, density driven flows lead salted water into original directions downstream waste heaps. From a numerical point of view, it is also shown that discretization degree plays a major role in the density effect modeling.
NASA Astrophysics Data System (ADS)
Viscardy, Sébastien; Errera, Quentin; Pitts, Michael C.; Poole, Lamont R.; Chabrillat, Simon; Daerden, Frank
2013-04-01
A 3-D Chemical Transport Model (CTM), with full stratospheric chemistry and driven by the ECMWF temperature and wind fields, has been coupled to a detailed PSC microphysical model to simulate polar winters. The formation and evolution of four types of PSC particles (STS, SAT, NAT, and ice) are described through relevant microphysical processes which alter interactively the nitric acid and water vapor concentrations of the atmosphere. Each particle type is described by a binned size distribution for the number density and chemical composition. This set-up allows for detailed calculation of optical properties and surface area densities used to compute the heterogeneous reaction rates. After comparing the evolution of the stratospheric chemical structure to satellite observations, we will investigate how the model reproduces the PSC coverage detected by CALIPSO. A comparison between the model and CALIPSO optical properties will be used to discuss the PSC composition. Finally, we aim at estimating the contribution of each PSC particle type to the chlorine activation.
NASA Astrophysics Data System (ADS)
Lipatov, A. S.; Sibeck, D. G.
2016-09-01
We use a new hybrid kinetic model to simulate the response of ring current, outer radiation belt, and plasmaspheric particle populations to impulsive interplanetary shocks. Since particle distributions attending the interplanetary shock waves and in the ring current and radiation belts are non-Maxwellian, wave-particle interactions play a crucial role in energy transport within the inner magnetosphere. Finite gyroradius effects become important in mass loading the shock waves with the background plasma in the presence of higher energy ring current and radiation belt ions and electrons. Initial results show that shocks cause strong deformations in the global structure of the ring current, radiation belt, and plasmasphere. The ion velocity distribution functions at the shock front, in the ring current, and in the radiation belt help us determine energy transport through the Earth's inner magnetosphere.
NASA Astrophysics Data System (ADS)
Rollins, C.; Argus, D. F.; Avouac, J. P.; Landry, W.; Barbot, S.
2015-12-01
North-south compression across the Los Angeles basin is accommodated by slip on thrust faults beneath the basin that may present significant seismic hazard to Los Angeles. Previous geodesy-based efforts to constrain the distributions and rates of elastic strain accumulation on these faults [Argus et al 2005, 2012] have found that the elastic model used has a first-order impact on the inferred distribution of locking and creep, underlining the need to accurately incorporate the laterally heterogeneous elastic structure and complex fault geometries of the Los Angeles basin into this analysis. We are using Gamra [Landry and Barbot, in prep.], a newly developed adaptive-meshing finite-difference solver, to compute elastostatic Green's functions that incorporate the full 3D regional elastic structure provided by the SCEC Community Velocity Model. Among preliminary results from benchmarks, forward models and inversions, we find that: 1) for a modeled creep source on the edge dislocation geometry from Argus et al [2005], the use of the SCEC CVM material model produces surface velocities in the hanging wall that are up to ~50% faster than those predicted in an elastic halfspace model; 2) in sensitivity-modulated inversions of the Argus et al [2005] GPS velocity field for slip on the same dislocation source, the use of the CVM deepens the inferred locking depth by ≥3 km compared to an elastic halfspace model; 3) when using finite-difference or finite-element models with Dirichlet boundary conditions (except for the free surface) for problems of this scale, it is necessary to set the boundaries at least ~100 km away from any slip source or data point to guarantee convergence within 5% of analytical solutions (a result which may be applicable to other static dislocation modeling problems and which may scale with the size of the area of interest). Here we will present finalized results from inversions of an updated GPS velocity field [Argus et al, AGU 2015] for the inferred
NASA Astrophysics Data System (ADS)
De Donatis, Mauro; Susini, Sara
2014-05-01
-consistent in showing the structural features of the study area. The work was not so straightforward, but the result is more then satisfying, even if some limitations were not solved (i.e. visualisation of bedding attitudes). Geological maps are fundamental for knowledge transfer among experts but, if combined with the innovative digital methods from survey to 3D model, this knowledges could reach a much larger number of people, allowing a cultural growth and the establishment of a larger awareness of the Earth and Environment.
Tajbakhsh, Jian; Stefanovski, Darko; Tang, George; Wawrowsky, Kolja; Liu, Naiyou; Fair, Jeffrey H.
2015-03-15
Cell-surface markers and transcription factors are being used in the assessment of stem cell fate and therapeutic safety, but display significant variability in stem cell cultures. We assessed nuclear patterns of 5-hydroxymethylcytosine (5hmC, associated with pluripotency), a second important epigenetic mark, and its combination with 5-methylcytosine (5mC, associated with differentiation), also in comparison to more established markers of pluripotency (Oct-4) and endodermal differentiation (FoxA2, Sox17) in mouse embryonic stem cells (mESC) over a 10-day differentiation course in vitro: by means of confocal and super-resolution imaging together with 3D high-content analysis, an essential tool in single-cell screening. In summary: 1) We did not measure any significant correlation of putative markers with global 5mC or 5hmC. 2) While average Oct-4 levels stagnated on a cell-population base (0.015 lnIU/day), Sox17 and FoxA2 increased 22-fold and 3-fold faster, respectively (Sox17: 0.343 lnIU/day; FoxA2: 0.046 lnIU/day). In comparison, global DNA methylation levels increased 4-fold faster (0.068 lnIU/day), and global hydroxymethylation declined at 0.046 lnIU/day, both with a better explanation of the temporal profile. 3) This progression was concomitant with the occurrence of distinct nuclear codistribution patterns that represented a heterogeneous spectrum of states in differentiation; converging to three major coexisting 5mC/5hmC phenotypes by day 10: 5hmC{sup +}/5mC{sup −}, 5hmC{sup +}/5mC{sup +}, and 5hmC{sup −}/5mC{sup +} cells. 4) Using optical nanoscopy we could delineate the respective topologies of 5mC/5hmC colocalization in subregions of nuclear DNA: in the majority of 5hmC{sup +}/5mC{sup +} cells 5hmC and 5mC predominantly occupied mutually exclusive territories resembling euchromatic and heterochromatic regions, respectively. Simultaneously, in a smaller subset of cells we observed a tighter colocalization of the two cytosine variants, presumably
Waveform synthesis of surface waves in a laterally heterogeneous earth by the Gaussian beam method
NASA Technical Reports Server (NTRS)
Yomogida, K.; Aki, K.
1985-01-01
The present investigation is concerned with an application of the Gaussian beam method to surface waves in the laterally heterogeneous earth. The employed method has been developed for ray tracing and synthesizing seismograms of surface waves in cases involving the laterally heterogeneous earth. The procedure is based on formulations derived by Yomogida (1985). Vertical structure of the wave field is represented by the eigenfunctions of normal mode theory, while lateral variation is expressed by the parabolic equation as in two-dimensional acoustic waves or elastic body waves. It is demonstrated that a large-amplitude change can result from a slight perturbation in the phase velocity model.
Li, X. D.; Fang, Y. M.; Wu, S. Q. E-mail: wsq@xmu.edu.cn; Zhu, Z. Z. E-mail: wsq@xmu.edu.cn
2015-05-15
Single adsorption of different atoms on pristine two-dimensional monolayer MoS{sub 2} have been systematically investigated by using density functional calculations with van der Waals correction. The adatoms cover alkali metals, alkaline earth metals, main group metal, 3d-transition metals, coinage metal and nonmetal atoms. Depending on the adatom type, metallic, semimetallic or semiconducting behavior can be found in direct bandgap monolayer MoS{sub 2}. Additionally, local or long-range magnetic moments of two-dimensional MoS{sub 2} sheet can also attained through the adsorption. The detailed atomic-scale knowledge of single adsorption on MoS{sub 2} monolayer is important not only for the sake of a theoretical understanding, but also device level deposition technological application.
Wen, Zhong Quan; Li, Min; Li, Fei; Zhu, Shi Jin; Liu, Xiao Ying; Zhang, Yu Xin; Kumeria, Tushar; Losic, Dusan; Gao, Yang; Zhang, Wei; He, Shi Xuan
2016-01-21
3-Dimensional (3D) composites based on a unique combination of MnO2-nanostructures, graphene oxide nanosheets and porous Diatomaceous Earth (DE) microparticles (GO-DE@MnO2) were synthesized and explored for application in high-performance supercapacitors. To explore the influence of the structural properties of MnO2 nanostructures on supercapacitor performances, several MnO2 structures with nanosheet and nanowire morphologies were synthesized and characterized. The prepared GO-DE@MnO2 composites with MnO2 nanosheets due to their higher conductivity and higher surface area showed a larger specific capacitance of 152.5 F g(-1) and a relatively better cycle stability (83.3% capacitance retention after 2000 cycles at a scan rate of 2 A g(-1)), indicating great potential for application in supercapacitors. PMID:26645931
NASA Astrophysics Data System (ADS)
Vidal Royo, O.
2014-12-01
GIS and geological modeling software have radically changed the means by which geological mapping is produced, published and visualized. This type of software environment normally requires a spatially aware reference system to position data and interpretation, often referred as georeferenced data (i.e. geographic data referenced on the Earth). However, for this study we coin the term areoreferenced data (i.e. Mars-referenced "geographic" data). Thanks to the wealth of areoreferenced data made available by the NASA and the HiRise at University of Arizona it is now possible to carry out 3D areographic and areologic (i.e. related to the topography and geology of Mars, respectively) reconstructions in great detail. The present work benefits from the availability of software and areographic data, and presents the results of an areologic map and 3D model of the fault systems in the Meridiani Planum of Mars. The work has been carried out in Move™ (developed by Midland Valley Exploration), a geological modeling toolkit that allows for easy data loading in a wide range of formats as well as straightforward 2D/3D model building tools of geological bodies. Initial data consisted of Digital Terrain Model and orthoimages (NASA/JPL/University of Arizona/USGS). From these we have interpreted several structural domains: right-lateral strike-slip systems with associated releasing bends, which gave room to an extensional event causing a horizontal-axis rotation of the bedding. Bedding ranges from subhorizontal in the southern domain where strike-slip prevails to nearly 40º in the central and northern domains, where a more complex interaction between strike-slip and extensional faults is described. The stratigraphic sequence is mainly composed by moderately rounded well laminated basaltic sandstones (Squyres et al., 2004) in which a high component of sulfurs (e.g. sulfate anhydrate, hexahydrite, epsomite, gypsum) and salts (e.g. halite) has been described (Squyres et al., 2004
NASA Astrophysics Data System (ADS)
von Hebel, Christian; Rudolph, Sebastian; Huisman, Johan A.; van der Kruk, Jan; Vereecken, Harry
2013-04-01
three different coil offsets in HCP and VCP measurement modes. This resulted in six high spatial resolution data sets of approximately 60000 measurements with different sensing depths. A 5 m block-kriging was applied to all six data sets to re-grid the sampling points on the same regular grid. For each grid node, the six measured apparent conductivities were used in a three-layer inversion. The three-layer inversion results of electrical conductivity thus obtained were used to derive a three-dimensional (3D) model of subsurface heterogeneity, which clearly indicated lateral and vertical conductivity changes of the subsurface that are related to changes in soil texture and soil water content.
PKiKP Coda Observations Interpreted in Terms of the Earth's Inner Core Heterogeneities
NASA Astrophysics Data System (ADS)
Krasnoshchekov, D.; Kaazik, P.; Ovtchinnikov, V.
2007-12-01
Although resulting from gradual few billion years long crystallization process, Earth's inner core (IC) is rather heterogeneous than a single crystal of iron. IC fabric variations such as changes in crystallographic alignment or supposed increase in crystal size with depth often constitute physical grounds for hypothesizing IC heterogeneities of various scale lengths. To constrain upper IC heterogeneities, we analyze reflections from the Earth's inner core boundary and the coda waves following the reflections (PKiKP) on array records of underground nuclear explosions. In particular, we compare PKiKP codas observed after reflections with close bounce points on the surface of the IC and diverse ray paths in the Earth's crust, mantle and outer core. Such PKiKP coda doublets show similar shape, frequency content, intensity and duration, and feature powerful arrivals originating from IC heterogeneities. The proposed interpretation favors misaligned anisotropic iron crystals up to 10 km in size to cause the observed PKiKP codas rather than classical scattering on inclusions in the outermost IC. Also, the IC fabric's image observed as PKiKP coda appears to be time stationary in 20 years span, failing to exhibit any pronounced effect of IC differential rotation.
NASA Astrophysics Data System (ADS)
Rodrigo Rodríguez Cardozo, Félix; Hjörleifsdóttir, Vala
2015-04-01
One important ingredient in the study of the complex active tectonics in Mexico is the analysis of earthquake focal mechanisms, or the seismic moment tensor. They can be determined trough the calculation of Green functions and subsequent inversion for moment-tensor parameters. However, this calculation is gets progressively more difficult as the magnitude of the earthquakes decreases. Large earthquakes excite waves of longer periods that interact weakly with laterally heterogeneities in the crust. For these earthquakes, using 1D velocity models to compute the Greens fucntions works well. The opposite occurs for smaller and intermediate sized events, where the relatively shorter periods excited interact strongly with lateral heterogeneities in the crust and upper mantle and requires more specific or regional 3D models. In this study, we calculate Greens functions for earthquakes in Mexico using a laterally heterogeneous seismic wave speed model, comprised of mantle model S362ANI (Kustowski et al 2008) and crustal model CRUST 2.0 (Bassin et al 1990). Subsequently, we invert the observed seismograms for the seismic moment tensor using a method developed by Liu et al (2004) an implemented by Óscar de La Vega (2014) for earthquakes in Mexico. By following a brute force approach, in which we include all observed Rayleigh and Love waves of the Mexican National Seismic Network (Servicio Sismológico Naciona, SSN), we obtain reliable focal mechanisms for events that excite a considerable amount of low frequency waves (Mw > 4.8). However, we are not able to consistently estimate focal mechanisms for smaller events using this method, due to high noise levels in many of the records. Excluding the noisy records, or noisy parts of the records manually, requires interactive edition of the data, using an efficient tool for the editing. Therefore, we developed a graphical user interface (GUI), based on python and the python library ObsPy, that allows the edition of observed and
Caicedo-Carvajal, Carlos E.; Liu, Qing; Remache, Yvonne; Goy, Andre; Suh, K. Stephen
2011-01-01
Isolation and amplification of primary lymphoma cells in vitro setting is technically and biologically challenging task. To optimize culture environment and mimic in vivo conditions, lymphoma cell lines were used as a test case and were grown in 3-dimension (3D) using a novel 3D tissue culture polystyrene scaffold with neonatal stromal cells to represent a lymphoma microenvironment. In this model, the cell proliferation was enhanced more than 200-fold or 20,000% neoplastic surplus in 7 days when less than 1% lymphoma cells were cocultured with 100-fold excess of neonatal stroma cells, representing 3.2-fold higher proliferative rate than 2D coculture model. The lymphoma cells grew and aggregated to form clusters during 3D coculture and did not maintained the parental phenotype to grow in single-cell suspension. The cluster size was over 5-fold bigger in the 3D coculture by day 4 than 2D coculture system and contained less than 0.00001% of neonatal fibroblast trace. This preliminary data indicate that novel 3D scaffold geometry and coculturing environment can be customized to amplify primary cancer cells from blood or tissues related to hematological cancer and subsequently used for personalized drug screening procedures. PMID:22073378
NASA Astrophysics Data System (ADS)
Jauer, P. R.; Gonzalez, W. D.; de Souza Costa, C. L.; Souza, V. M.
2013-12-01
The interaction, transport and conversion of energy between the solar wind and Earth's magnetosphere have been studied for decades through in situ measurements and Magnetohydrodynamics simulation, (MHD). Nevertheless, due to the vast regions of space and nonlinearities of the physical processes there are many questions that still remain without conclusive answers. Currently, the MHD simulation is a powerful tool that helps other means of already existing research, even within its theoretical limitation; it provides information of the space regions where in situ measurements are rare or nonexistent. The aim of this work is the study of energy transfer from the solar wind through the calculation of the divergence of the Poynting vector for the inner regions of the Earth's magnetosphere, especially the magneto tail using 3D global MHD numerical code Space Weather Modelling Framework (SWMF) / (Block Adaptive Tree Solar wind Roe Upwind Scheme) (BATS-R-US), developed by the University of Michigan. We conducted a simulation study for the event that occurred on September 21-27, 1999, for which the peak value of the interplanetary magnetic field was -22 nT, and gave rise to an intense magnetic storm with peak Dst of -160 nT. Furthermore, we compare the results of the power estimated by the model - through the integration of the Poynting vector in rectangular region of the tail, with a domain -130
The relation of mantle heterogeneity to the bulk composition and origin of the earth
NASA Technical Reports Server (NTRS)
Smith, J. V.
1980-01-01
The heterogeneity of the mantle can be viewed in the context of models for accretion of the terrestrial planets from the solar nebula. Oxygen isotope ratios and mineralogy indicate the existence of hot planetesimals of diverse compositions. Assuming that nebular condensates range from a reduced state near the sun to an oxidized state near Jupiter, a new model is proposed for heterogeneous accretion of the earth beginning with hot, reduced condensates and ending with cool, oxidized condensates. The Ganapathy-Anders cosmochemical model for the bulk composition of the earth was tested by summing measured compositions for the three outer zones and likely compositions for the inner zones. Revisions are suggested, including reduction of the content of the early condensate from that suggested by taking the U concentration as 30 ng/g, as suggested by the naive interpretation of the heat flow.
How to Represent 100-meter Spatial Heterogeneity in Earth System Models
NASA Astrophysics Data System (ADS)
Chaney, Nathaniel; Shevliakova, Elena; Malyshev, Sergey
2016-04-01
Terrestrial ecosystems play a pivotal role in the Earth system; they have a profound impact on the global climate, food and energy production, freshwater resources, and biodiversity. One of the most fascinating yet challenging aspects of characterizing terrestrial ecosystems is their field-scale (~100 m) spatial heterogeneity. It has been observed repeatedly that the water, energy, and biogeochemical cycles at multiple temporal and spatial scales have deep ties to an ecosystem's spatial structure. Current Earth system models largely disregard this important relationship leading to an inadequate representation of ecosystem dynamics. In this presentation, we will show how existing hyperresolution environmental datasets can be harnessed to explicitly represent field-scale spatial heterogeneity in Earth system models. For each macroscale grid cell, these environmental data are clustered according to their field-scale soil and topographic attributes to define unique sub-grid tiles or hydrologic response units (HRUs). The novel Geophysical Fluid Dynamics Laboratory (GFDL) LM3-TiHy-PPA land model is then used to simulate these HRUs and their spatial interactions via the exchange of water, energy, and nutrients along explicit topographic gradients. Using historical simulations over the contiguous United States, we will show how a robust representation of field-scale spatial heterogeneity impacts modeled ecosystem dynamics including the water, energy, and biogeochemical cycles as well as vegetation composition and distribution.
NASA Astrophysics Data System (ADS)
Isik-Ercan, Zeynep; Zeynep Inan, Hatice; Nowak, Jeffrey A.; Kim, Beomjin
2014-01-01
This qualitative case study describes (a) the ways 3D visualization, coupled with other science and literacy experiences, supported young children's first exploration of the Earth-Sun-Moon system and (b) the perspectives of classroom teachers and children on using 3D visualization. We created three interactive 3D software modules that simulate day and night, Moon phases and seasons. These modules were used in a science and literacy unit for 35 second graders at an urban elementary school in Midwestern USA. Data included pre- and post-interviews, audio-taped lessons and classroom observations. Post-interviews demonstrated that children's knowledge of the shapes and the movements of the Earth and Moon, alternation of day and night, the occurrence of the seasons, and Moon's changing appearance increased. Second graders reported that they enjoyed expanding their knowledge through hands-on experiences; through its reality effect, 3D visualization enabled them to observe the space objects that move in the virtual space. The teachers noted that 3D visualization stimulated children's interest in space and that using 3D visualization in combination with other teaching methods-literacy experiences, videos and photos, simulations, discussions, and presentations-supported student learning. The teachers and the students still experienced challenges using 3D visualization due to technical problems with 3D vision and time constraints. We conclude that 3D visualization offers hands-on experiences for challenging science concepts and may support young children's ability to view phenomena that would typically be observed through direct, long-term observations in outer space. Results imply a reconsideration of assumed capabilities of young children to understand astronomical phenomena.
NASA Astrophysics Data System (ADS)
Moulik, P.; Ekstrom, G.
2012-12-01
We have developed a framework that can be used to investigate anisotropic velocity, density and anelastic heterogeneity in the Earth's mantle using a wide spectrum (0.3-50 mHz) of seismological observables. We start with the extensive dataset of surface-wave phase anomalies, long-period waveforms, and body-wave travel times collected by Kustowski et al. (2008) for the development of the global model S362ANI. The additional data included in our analysis are splitting functions of spheroidal and toroidal modes, which are analogous to phase velocity maps at low frequencies. We include in this set of observations a new dataset containing the splitting functions of 56 spheroidal fundamental modes and overtones, measured by Deuss et al. (2011, 2012) using data from large recent earthquakes. Apart from providing unique constraints on the long-wavelength elastic and density structure in the mantle, the overtone splitting data are especially sensitive to the velocity (and anisotropic) structure in the transition zone and in the deeper mantle. The detection of anisotropy, a marker of flow, in the transition zone has implications for our understanding of mantle convection. Our forward modeling of the splitting functions, like the other types of data, includes the effects of radial anisotropy (Mochizuki, 1986). We show that the upper-mantle shear-wave anisotropy of S362ANI generates a clear contribution to the splitting functions of the modes that are sensitive to the upper-mantle structure. We explore the tradeoffs between fitting the mode splitting functions and the travel-times of body waves that turn in the transition zone or in the lower mantle (e.g. SS), while observing that the waveforms and the surface wave phase-anomalies provide complementary information about the mantle. Our experiments suggest that the splitting data are sufficiently sensitive to the anisotropy in the mantle such that their inclusion may provide a better depth resolution of the anisotropic shear
Lay, T.; Wu, R.S.
1994-12-13
This document is the final report for this grant to develop new three-dimensional wave propagation techniques for high frequency waves in heterogeneous media. The report is divided into four sections, each being a published paper sponsored by this grant. In the first section we formulate a one-way wide-angle elastic wave propagation method for arbitrarily heterogeneous media in both the space and wavenumber domains using elastic Rayleigh integrals and local elastic Born scattering theory. In the second section this complex phase screen method is compared with fourth-order finite differences and exact eigenfunction expansion calculations for two-dimensional inhomogeneous media to assess the accuracy of the one-way propagation algorithm. In the third section, an observational study of continental margin structure influence on Lg propagation is presented, using data from the former Soviet stations for nuclear explosions at Novaya Zemlya. We find that bathymetric features can be correlated with energy levels of Lg, suggesting that waveguide structure influences regional phase energy partitioning. This idea is pursued in the fourth section, using Eurasian earthquake and nuclear explosion data along with information about the crustal structure in Eurasia. We develop empirical relations that reduce the scatter in the P/Lg discriminant at low frequency.
Heterogeneous growth of meteorites and planets, especially the earth and moon
NASA Technical Reports Server (NTRS)
Smith, J. V.
1982-01-01
It is pointed out that enough information remains in meteorites to show that the solar nebula was chemically inhomogeneous, and that it did not develop under physical and chemical equilibrium. It is found that the earth is not in chemical equilibrium, and that it still retains vestiges of its heterogeneous growth in spite of prolonged convection and volcanism. The present investigation has the aim to unify simple ideas into a comprehensive, complex scheme for heterogeneous growth of the planets and parent bodies of meteorites. A scheme is considered for the conversion of dust and gas in the solar nebula into the present population of bodies. Attention is given to possible processes for the growth and development of planetesimals from dust via clouds, the major episodes in the development of the solar system, the relation of meteorites to a solar nebula with monotonic variation of temperature and oxidation state, the significance of meteorite properties, the composition of the inner planets, and the development of the earth.
NASA Astrophysics Data System (ADS)
Audigane, Pascal; Chiaberge, Christophe; Mathurin, Frédéric; Lions, Julie; Picot-Colbeaux, Géraldine
2011-04-01
This paper is addressed to the TOUGH2 user community. It presents a new tool for handling simulations run with the TOUGH2 code with specific application to CO 2 geological storage. This tool is composed of separate FORTRAN subroutines (or modules) that can be run independently, using input and output files in ASCII format for TOUGH2. These modules have been developed specifically for modeling of carbon dioxide geological storage and their use with TOUGH2 and the Equation of State module ECO2N, dedicated to CO 2-water-salt mixture systems, with TOUGHREACT, which is an adaptation of TOUGH2 with ECO2N and geochemical fluid-rock interactions, and with TOUGH2 and the EOS7C module dedicated to CO 2-CH 4 gas mixture is described. The objective is to save time for the pre-processing, execution and visualization of complex geometry for geological system representation. The workflow is rapid and user-friendly and future implementation to other TOUGH2 EOS modules for other contexts (e.g. nuclear waste disposal, geothermal production) is straightforward. Three examples are shown for validation: (i) leakage of CO 2 up through an abandoned well; (ii) 3D reactive transport modeling of CO 2 in a sandy aquifer formation in the Sleipner gas Field, (North Sea, Norway); and (iii) an estimation of enhanced gas recovery technology using CO 2 as the injected and stored gas to produce methane in the K12B Gas Field (North Sea, Denmark).
Gigot, C.; de Vallavieille-Pope, C.; Huber, L.; Saint-Jean, S.
2014-01-01
Background and Aims Recent developments in plant disease management have led to a growing interest in alternative strategies, such as increasing host diversity and decreasing the use of pesticides. Use of cultivar mixtures is one option, allowing the spread of plant epidemics to be slowed down. As dispersal of fungal foliar pathogens over short distances by rain-splash droplets is a major contibutor to the spread of disease, this study focused on modelling the physical mechanisms involved in dispersal of a non-specialized pathogen within heterogeneous canopies of cultivar mixtures, with the aim of optimizing host diversification at the intra-field level. Methods Virtual 3-D wheat-like plants (Triticum aestivum) were used to consider interactions between plant architecture and disease progression in heterogeneous canopies. A combined mechanistic and stochastic model, taking into account splash droplet dispersal and host quantitative resistance within a 3-D heterogeneous canopy, was developed. It consists of four sub-models that describe the spatial patterns of two cultivars within a complex canopy, the pathway of rain-splash droplets within this canopy, the proportion of leaf surface area impacted by dispersal via the droplets and the progression of disease severity after each dispersal event. Key Results Different spatial organization, proportions and resistance levels of the cultivars of two-component mixtures were investigated. For the eight spatial patterns tested, the protective effect against disease was found to vary by almost 2-fold, with the greatest effect being obtained with the smallest genotype unit area, i.e. the ground area occupied by an independent unit of the host population that is genetically homogeneous. Increasing both the difference between resistance levels and the proportion of the most resistant cultivar often resulted in a greater protective effect; however, this was not observed for situations in which the most resistant of the two
Necessity of using heterogeneous ellipsoidal Earth model with terrain to calculate co-seismic effect
NASA Astrophysics Data System (ADS)
Cheng, Huihong; Zhang, Bei; Zhang, Huai; Huang, Luyuan; Qu, Wulin; Shi, Yaolin
2016-04-01
Co-seismic deformation and stress changes, which reflect the elasticity of the earth, are very important in the earthquake dynamics, and also to other issues, such as the evaluation of the seismic risk, fracture process and triggering of earthquake. Lots of scholars have researched the dislocation theory and co-seismic deformation and obtained the half-space homogeneous model, half-space stratified model, spherical stratified model, and so on. Especially, models of Okada (1992) and Wang (2003, 2006) are widely applied in the research of calculating co-seismic and post-seismic effects. However, since both semi-infinite space model and layered model do not take the role of the earth curvature or heterogeneity or topography into consideration, there are large errors in calculating the co-seismic displacement of a great earthquake in its impacted area. Meanwhile, the computational methods of calculating the co-seismic strain and stress are different between spherical model and plane model. Here, we adopted the finite element method which could well deal with the complex characteristics (such as anisotropy, discontinuities) of rock and different conditions. We use the mash adaptive technique to automatically encrypt the mesh at the fault and adopt the equivalent volume force replace the dislocation source, which can avoid the difficulty in handling discontinuity surface with conventional (Zhang et al., 2015). We constructed an earth model that included earth's layered structure and curvature, the upper boundary was set as a free surface and the core-mantle boundary was set under buoyancy forces. Firstly, based on the precision requirement, we take a testing model - - a strike-slip fault (the length of fault is 500km and the width is 50km, and the slippage is 10m) for example. Because of the curvature of the Earth, some errors certainly occur in plane coordinates just as previous studies (Dong et al., 2014; Sun et al., 2012). However, we also found that: 1) the co
NASA Astrophysics Data System (ADS)
Chang, Xin
This dissertation proposal is concerned with the use of fast and broadband full-wave electromagnetic methods for modeling high speed interconnects (e.g, vertical vias and horizontal traces) and passive components (e.g, decoupling capacitors) for structures of PCB and packages, in 3D IC, Die-level packaging and SIW based devices, to effectively modeling the designs signal integrity (SI) and power integrity (PI) aspects. The main contributions finished in this thesis is to create a novel methodology, which hybridizes the Foldy-Lax multiple scattering equations based fast full wave method, method of moment (MoM) based 1D technology, modes decoupling based geometry decomposition and cavity modes expansions, to model and simulate the electromagnetic scattering effects for the irregular power/ground planes, multiple vias and traces, for fast and accurate analysis of link level simulation on multilayer electronic structures. For the modeling details, the interior massively-coupled multiple vias problem is modeled most-analytically by using the Foldy-Lax multiple scattering equations. The dyadic Green's functions of the magnetic field are expressed in terms of waveguide modes in the vertical direction and vector cylindrical wave expansions or cavity modes expansions in the horizontal direction, combined with 2D MoM realized by 1D technology. For the incident field of the case of vias in the arbitrarily shaped antipad in finite large cavity/waveguide, the exciting and scattering field coefficients are calculated based on the transformation which converts surface integration of magnetic surface currents in antipad into 1D line integration of surface charges on the vias and on the ground plane. Geometry decomposition method is applied to model and integrate both the vertical and horizontal interconnects/traces in arbitrarily shaped power/ground planes. Moreover, a new form of multiple scattering equations is derived for solving coupling effects among mixed metallic
Fine scale heterogeneity in the Earth's upper mantle - observation and interpretation
NASA Astrophysics Data System (ADS)
Thybo, Hans
2014-05-01
High resolution seismic data has over the last decade provided significant evidence for pronounced fine scale heterogeneity in the Earth's mantle at an unprecedented detail. Seismic tomography developed tremendously during the last 20-30 years. The results show overall structure in the mantle which can be correlated to main plate tectonic features, such as oceanic spreading centres, continental rift zones and subducting slabs. Much seismological mantle research is now concentrated on imaging fine scale heterogeneity, which may be detected and imaged with high-resolution seismic data with dense station spacing and at high frequency, e.g. from the Russian Peaceful Nuclear Explosion (PNE) data set and array recordings of waves from natural seismic sources. Mantle body waves indicate pronounced heterogeneity at three depth levels whereas other depth intervals appear transparent, at least in the frequency band of 0.5-15 Hz: (1) The Mantle Low-Velocity Zone (LVZ) is a global feature which has been detected in more than 50 long-range seismic profiles (Thybo and Perchuc, Science, 1997). Since then numerous studies based on receiver functions, surface waves, and controlled source seismology have confirmed the presence of this zone. The data demonstrates that the top of the LVZ everywhere is at a depth of 100±20 km. A pronounced coda shows that the zone is highly heterogeneous at characteristic scale lengths of 5-15 by 2-6 km. We interpret that the rocks in the LVZ have a temperature close to the solidus or even may contain small fractions of partial melt. The solidus of mantle rocks is very low below a depth of ca. 90 km if volatiles are present due to a characteristic kink in the solidus which is much lower than for dry mantle rocks. We suggest that the rocks are in a totally solid state below the LVZ and that the depth to the interface to fully solid rocks is an indicator of the thermal state of the upper mantle. (2) Significant scattering from around the top of the
Fine scale heterogeneity in the Earth's upper mantle - observation and interpretation
NASA Astrophysics Data System (ADS)
Thybo, Hans
2013-04-01
High resolution seismic data has over the last decade provided significant evidence for pronounced fine scale heterogeneity in the Earth's mantle at an unprecedented detail. Seismic tomography developed tremendously during the last 20-30 years. The results show overall structure in the mantle which can be correlated to main plate tectonic features, such as oceanic spreading centres, continental rift zones and subducting slabs. Much seismological mantle research is now concentrated on imaging fine scale heterogeneity, which may be detected and imaged with high-resolution seismic data with dense station spacing and at high frequency, e.g. from the Russian Peaceful Nuclear Explosion (PNE) data set and array recordings of waves from natural seismic sources. Mantle body waves indicate pronounced heterogeneity at three depth levels whereas other depth intervals appear transparent, at least in the frequency band of 0.5-15 Hz: (1) The Mantle Low-Velocity Zone (LVZ) is a global feature which has been detected in more than 50 long-range seismic profiles (Thybo and Perchuc, Science, 1997). Since then numerous studies based on receiver functions, surface waves, and controlled source seismology have confirmed the presence of this zone. The data demonstrates that the top of the LVZ everywhere is at a depth of 100±20 km. A pronounced coda shows that the zone is highly heterogeneous at characteristic scale lengths of 5-15 by 2-6 km. We interpret that the rocks in the LVZ have a temperature close to the solidus or even may contain small fractions of partial melt. The solidus of mantle rocks is very low below a depth of ca. 90 km if volatiles are present due to a characteristic kink in the solidus which is much lower than for dry mantle rocks. We suggest that the rocks are in a totally solid state below the LVZ and that the depth to the interface to fully solid rocks is an indicator of the thermal state of the upper mantle. (2) Significant scattering from around the top of the
NASA Astrophysics Data System (ADS)
Pletinckx, D.
2011-09-01
The current 3D hype creates a lot of interest in 3D. People go to 3D movies, but are we ready to use 3D in our homes, in our offices, in our communication? Are we ready to deliver real 3D to a general public and use interactive 3D in a meaningful way to enjoy, learn, communicate? The CARARE project is realising this for the moment in the domain of monuments and archaeology, so that real 3D of archaeological sites and European monuments will be available to the general public by 2012. There are several aspects to this endeavour. First of all is the technical aspect of flawlessly delivering 3D content over all platforms and operating systems, without installing software. We have currently a working solution in PDF, but HTML5 will probably be the future. Secondly, there is still little knowledge on how to create 3D learning objects, 3D tourist information or 3D scholarly communication. We are still in a prototype phase when it comes to integrate 3D objects in physical or virtual museums. Nevertheless, Europeana has a tremendous potential as a multi-facetted virtual museum. Finally, 3D has a large potential to act as a hub of information, linking to related 2D imagery, texts, video, sound. We describe how to create such rich, explorable 3D objects that can be used intuitively by the generic Europeana user and what metadata is needed to support the semantic linking.
Tajbakhsh, Jian; Stefanovski, Darko; Tang, George; Wawrowsky, Kolja; Liu, Naiyou; Fair, Jeffrey H.
2015-01-01
Cell-surface markers and transcription factors are being used in the assessment of stem cell fate and therapeutic safety, but display significant variability in stem cell cultures. We assessed nuclear patterns of 5-hydroxymethylcytosine (5hmC, associated with pluripotency), a second important epigenetic mark, and its combination with 5-methylcytosine (5mC, associated with differentiation), also in comparison to more established markers of pluripotency (Oct-4) and endodermal differentiation (FoxA2, Sox17) in mouse embryonic stem cells (mESC) over a ten-day differentiation course in vitro: by means of confocal and super-resolution imaging together with high-content analysis, an essential tool in single-cell screening. In summary: 1) We did not measure any significant correlation of putative markers with global 5mC or 5hmC. 2) While average Oct-4 levels stagnated on a cell-population base (0.015 lnIU per day), Sox17 and FoxA2 increased 22-fold and 3-fold faster, respectively (Sox17:0.343 lnIU/day; FoxA2: 0.046 lnIU/day). In comparison, DNA global methylation levels increased 4-fold faster (0.068 lnIU/day), and global hydroxymethylation declined at 0.046 lnIU/day, both with a better explanation of the temporal profile. 3) This progression was concomitant with the occurrence of distinct nuclear codistribution patterns that represented a heterogeneous spectrum of states in differentiation; converging to three major coexisting 5mC/5hmC phenotypes by day 10: 5hmC+/5mC−, 5hmC+/5mC+, and 5hmC−/5mC+ cells. 4) Using optical nanoscopy we could delineate the respective topologies of 5mC/5hmC colocalization in subregions of nuclear DNA: in the majority of 5hmC+/5mC+ cells 5hmC and 5mC predominantly occupied mutually exclusive territories resembling euchromatic and heterochromatic regions, respectively. Simultaneously, in a smaller subset of cells we observed a tighter colocalization of the two cytosine variants, presumably delineating chromatin domains in remodeling. We
Where should fine-resolution spatial heterogeneity be captured within Earth System Models?
NASA Astrophysics Data System (ADS)
Adam, J. C.; Hull, R.; Tague, C.; Reyes, J. J.; Liu, M.
2015-12-01
Land-atmosphere interactions impact the environment in many ways, such as through partially driving our climate system, and in changing the availability and usability of our natural resources. Earth System Models (EaSMs) are being used increasingly to explore these coupled dynamics from watershed to global scales. However, many EaSMs do not adequately represent landscape-scale spatial heterogeneity that influences land surface response, as relatively coarse resolution simulations are necessitated by computational limitations. Research is needed to understand which types of spatial heterogeneity, over which biomes and climate types, should be represented such that an EaSM accurately captures the aggregate land surface response to a changing climate. Spatial heterogeneity in a landscape arises due to differences in model forcings; in underlying soil, vegetation, and topographic properties that control moisture, energy and nutrient fluxes; and in land surface responses that arise due to spatially-organized connections. While our long-term goal is to understand how each of these sources should be represented in an EaSM, in this study we focus first on parameter heterogeneity. We apply the Regional Hydro-Ecological Simulation System (RHESSys), a distributed process-based model that was originally developed for catchment-scale applications. We explore the functional form of the hydrologic response of a RHESSys "patch" (a 200-400 m element with homogenous landscape parameters) to an invoked change. According to scale transition theory, a linear response makes it is possible to upscale (or aggregate) the model resolution without biasing the model response. We perform RHESSys simulations for more than 500 individual catchments within the Willamette and Yakima River basins in the Pacific Northwest region of the U.S. Each catchment was imposed with incremental perturbations of temperature and precipitation. The response curves for hydrologic variables such as
3d-3d correspondence revisited
NASA Astrophysics Data System (ADS)
Chung, Hee-Joong; Dimofte, Tudor; Gukov, Sergei; Sułkowski, Piotr
2016-04-01
In fivebrane compactifications on 3-manifolds, we point out the importance of all flat connections in the proper definition of the effective 3d {N}=2 theory. The Lagrangians of some theories with the desired properties can be constructed with the help of homological knot invariants that categorify colored Jones polynomials. Higgsing the full 3d theories constructed this way recovers theories found previously by Dimofte-Gaiotto-Gukov. We also consider the cutting and gluing of 3-manifolds along smooth boundaries and the role played by all flat connections in this operation.
Ma, Xue; Tian, Jing; Yang, Hong-Y.; Zhao, Kai; Li, Xia
2013-05-01
The reaction of RE(NO)₃·6H₂O and formamide yielded the coordination polymers, [RE(HCOO)₄]⁻[NH₂CHNH₂]⁺ (RE=Y 1, Eu 2, Gd 3, Tb 4, Dy 5, Er 6, and Yb 7). They possess 3D porous frameworks with the 1D rhombic channels occupied by [NH₂CHNH₂]⁺ cations. Complexes 2 and 4 display the characteristic down-conversion emissions corresponding to ⁵D₀→⁷F_{J} (J=1–4) transitions of Eu(III) ion and ⁵D₄→⁷F_{J} (J=6–3) transitions of Tb(III) ion, respectively. Longer lifetime values of 2.128±0.002 ms (⁵D₀) for 2 and 2.132±0.002 ms (⁵D₄) for 4 have been observed. The up-conversion spectra of the Y:Yb,Er and Gd:Yb,Er codoped complexes exhibit three emission bands around 410 (⁴H_{9/2}→⁴I_{15/2}, blue), 518–570 (⁴S_{3/2}, ²H_{11/2}→⁴I_{15/2}, green), and 655 nm (⁴F_{9/2}→⁴I_{15/2}, red). - Graphical Abstract: The complexes [RE(HCOO)₄]⁻[NH₂CHNH₂]⁺ possess 3D porous frameworks. Eu(III) and Tb(III) complexes show characteristic emission of Ln(III) ions. The up-conversion emission of the Y:Yb,Er and Gd:Yb,Er codoped complexes was observed. Highlights: •The reaction of RE(NO)₃·6H₂O and formamide produced complexes [RE(HCOO)₄]⁻[NH₂CHNH₂]⁺. • The complexes possess 3D frameworks with the 1D channels occupied by [NH₂CHNH₂]^{+} cations. • Eu(III)/Tb(III) complexes display the characteristic down-conversion emission of Ln(III) ions. • The Y:Yb,Er and Gd:Yb,Er doped complexes exhibit the up-conversion emission.
NASA Technical Reports Server (NTRS)
Tanimoto, T.
1984-01-01
A simple modification of Gilbert's formula to account for slight lateral heterogeneity of the earth leads to a convenient formula to calculate synthetic long period seismograms. Partial derivatives are easily calculated, thus the formula is suitable for direct inversion of seismograms for lateral heterogeneity of the earth. Previously announced in STAR as N83-29893
NASA Technical Reports Server (NTRS)
Tanimoto, T.
1983-01-01
A simple modification of Gilbert's formula to account for slight lateral heterogeneity of the Earth leads to a convenient formula to calculate synthetic long period seismograms. Partial derivatives are easily calculated, thus the formula is suitable for direct inversion of seismograms for lateral heterogeneity of the Earth.
NASA Astrophysics Data System (ADS)
Meulien Ohlmann, Odile
2013-02-01
Today the industry offers a chain of 3D products. Learning to "read" and to "create in 3D" becomes an issue of education of primary importance. 25 years professional experience in France, the United States and Germany, Odile Meulien set up a personal method of initiation to 3D creation that entails the spatial/temporal experience of the holographic visual. She will present some different tools and techniques used for this learning, their advantages and disadvantages, programs and issues of educational policies, constraints and expectations related to the development of new techniques for 3D imaging. Although the creation of display holograms is very much reduced compared to the creation of the 90ies, the holographic concept is spreading in all scientific, social, and artistic activities of our present time. She will also raise many questions: What means 3D? Is it communication? Is it perception? How the seeing and none seeing is interferes? What else has to be taken in consideration to communicate in 3D? How to handle the non visible relations of moving objects with subjects? Does this transform our model of exchange with others? What kind of interaction this has with our everyday life? Then come more practical questions: How to learn creating 3D visualization, to learn 3D grammar, 3D language, 3D thinking? What for? At what level? In which matter? for whom?
Viscoelastic relaxation in a heterogeneous Earth following the 2004 Sumatra-Andaman earthquake
NASA Astrophysics Data System (ADS)
Wiseman, Kelly; Bürgmann, Roland; Freed, Andrew M.; Banerjee, Paramesh
2015-12-01
Consideration of the three-dimensional heterogeneity of mantle rheology allows models of viscoelastic relaxation following the 2004 Sumatra-Andaman earthquake to simultaneously fit both the observed far-field and near-field postseismic deformation. We use horizontal and vertical campaign and continuous GPS observations from the Andaman, Nicobar, and Sumatran forearc islands, mainland Sumatra, Thailand, the Malay Peninsula, the Indian Ocean, and southern India, spanning the first five years of postseismic deformation. The postseismic relaxation models consider contributions from the 2004 Mw 9.2 Sumatra-Andaman, the 2005 Mw 8.7 Nias, and 2007 Mw 8.4 Bengkulu earthquakes. Far-field motions to the east of the ruptures are equally well fit by homogeneous or laterally variable earth models. However, only models with contrasting rheology across the subducting slab, a ten-times higher mantle viscosity under the Indian Ocean lithosphere than the backarc mantle, can also produce the observed enduring postseismic uplift along the forearc and lack of far-field transient displacements in southern India. While postseismic uplift of forearc stations can also be produced by rapid and enduring down-dip afterslip, the inferred rheology structure is consistent with the distribution of mantle temperature inferred from seismic tomography.
NASA Astrophysics Data System (ADS)
Wolf, Eric T.; Kopparapu, Ravi; Haqq-Misra, Jacob; Toon, Owen Brian
2015-12-01
The host star imposes a primary control on terrestrial planet climate. Both the spectral energy distribution and the main sequence lifetime vary as a function of stellar type. Here we present recent results from three-dimensional climate system models describing the evolutionary sequence of Earth-analog planets throughout their habitable lifetimes. Climatological evolution is traced from snowball to moist greenhouse, representing the conventional end-member states of the habitable zone. For Earth the habitable period would have been tantalizingly short, if not for geological and biological regulation of greenhouse gases. Without active carbon cycling, an early snowball could not have been broken until late in Earth’s history. Abrupt solar driven deglaciation would soon be followed by the onset of the water vapor greenhouse feedback and a moist greenhouse climate, leaving little over 1 billion years of habitable surface conditions. Around bluer stars, the habitable period for terrestrial planets is constricted further due to their reduced main sequence lifetimes and thus more rapid brightening. Planets with long-lived habitable periods are most likely found around stars redder than the Sun due to their more gradual brightening.
NASA Astrophysics Data System (ADS)
Tung, Sui; Masterlark, Timothy
2016-05-01
We derive a coseismic slip model of the 2015 Mw7.8 Gorkha earthquake on the basis of GPS and line-of-sight displacements from ALOS-2 descending interferograms, using Green's functions calculated with a 3-D finite element model (FEM). The FEM simulates a nonuniform distribution of elastic material properties and a precise geometric configuration of the irregular topographical surface. The rupturing fault is modeled as a low-angle and north dipping surface within the Main Frontal Thrust along the convergent margin of the Himalayas. The optimal model that inherits heterogeneous material properties provides a significantly better solution than that in a homogenous domain at the 95% confidence interval. The best fit solution for the domain having a nonuniform distribution of material properties reveals a rhombus-shaped slip zone of three composite asperities. Slip is primarily concentrated at a depth of 15 km with both dip-slip (maximum 6.54 m) and strike-slip (maximum 2.0 m) components, giving rise to a geodetic-based moment of 1.09 × 1021 Nm in general agreement with the seismological estimate. The optimal relative weights among GPS and interferometric synthetic aperture radar (InSAR) are deduced from a new method, MC-HVCE which combines a Monte Carlo search and a Helmert Method of Variance Components Estimation. This method determines the relative weights in a systemic approach which preserves the intrinsic solution smoothness. The joint solution is significantly better than those inverted from each individual data set. This methodology allows us to integrate multiple data sets of geodetic observations with seismic tomography, in an effort to achieve a better understanding of seismic ruptures within crustal heterogeneity.
Zhao, Hongying; Qian, Lin; Guan, Xiaohong; Wu, Deli; Zhao, Guohua
2016-05-17
Novel iron-copper-carbon (FeCuC) aerogel was fabricated through a one-step process from metal-resin precursors and then activated with CO2 and N2 in environmentally friendly way. The activated FeCuC aerogel was applied in a heterogeneous electro-Fenton (EF) process and exhibited higher mineralization efficiency than homogeneous EF technology. High total organic carbon (TOC) removal of organic pollutants with activated FeCuC aerogel was achieved at a wide range of pH values (3-9). The chemical oxygen demand (COD) of real dyeing wastewater was below China's discharge standard after 30 min of treatment, and the specific energy consumption was low (9.2 kW·h·kg(-1)COD(-1)), corresponding to a power consumption of only ∼0.34 kW·h per ton of wastewater. The enhanced mineralization efficiency of FeCuC aerogel was mostly attributable to ultradispersed metallic Fe-Cu nanoparticles embedded in 3D carbon matrix and the CO2-N2 treatment. The CO2 activation enhanced the accessibility of the aerogel's pores, and the secondary N2 activation enlarged the porosity and regenerated the ultradispersed zerovalent iron (Fe(0)) with reductive carbon. Cu(0) acted as a reduction promoter for interfacial electron transfer. Moreover, activated FeCuC aerogel presented low iron leaching (<0.1 ppm) in acidic solution and can be molded into different sizes with high flexibility. Thus, this material could be used as a low-cost cathode and efficient heterogeneous EF technology for actual wastewater treatment. PMID:27082750
ERIC Educational Resources Information Center
Hastings, S. K.
2002-01-01
Discusses 3 D imaging as it relates to digital representations in virtual library collections. Highlights include X-ray computed tomography (X-ray CT); the National Science Foundation (NSF) Digital Library Initiatives; output peripherals; image retrieval systems, including metadata; and applications of 3 D imaging for libraries and museums. (LRW)
Constraints from Earth's heat budget on the origins of deep mantle heterogeneity
NASA Astrophysics Data System (ADS)
Kellogg, L. H.
2008-12-01
To understand the origin and evolution of mantle heterogeneity, a multidisciplinary approach is required. A successful model must satisfy the constraints provided by seismology, which reveals the current state of the mantle, and geochemistry, which reflects the time-integrated history of the planet. One of the most important observations is the prominence of the degree-2 structure in the lower mantle. The large structures in the lower mantle beneath Africa and the Pacific contribute to the degree-2 signal, but their origin is not yet understood. A number of mantle models have been proposed containing multiple reservoirs in various configurations to produce a depleted source of mid-ocean ridge basalts along with the range of long-lived sources inferred from the geochemistry of ocean islands. One important constraint is provided by the Earth's total heat loss. Heat is produced by radioactive decay in the depleted mantle and continental crust; an additional radioactive source is required to make up the heat production of the bulk silicate earth known from as revealed by cosmochemistry and geodynamical models. This heat source may be located in the lower mantle, possibly associated with the superplume structures. This requires a tradeoff between the radiogenic heat productivity of the reservoir and its size: a small volume must contain a higher concentration of heat producing elements to balance the global heat budget. A very small reservoir at the base of the mantle can likely be ruled out because it would become extremely hot (resulting in velocities that are incompatible with seismic models) or buoyant (and therefore short-lived.) The dynamical implications of a dense layer can be assessed using the criteria that the layer must be stable through time, shows topography at its interface, has an effective density profile compatible with seismic models, and has an appropriate heat flux across the core- mantle boundary. A reservoir that is enriched in radiogenic elements
NASA Astrophysics Data System (ADS)
McClusky, S.; Reilinger, R. E.; Uzel, T.; Eren, K.; Dindar, A. A.
2010-12-01
Over the last twenty-five years an increasing number of discrete GNSS networks have been installed across the globe for a variety of purposes. These networks have been built with specific scientific or civil goals in mind such as; monitoring atmospheric or ionospheric properties, definition of terrestrial or cadastral reference frames, monitoring or tracking of infrastructure, observation tectonic processes, etc… These networks were designed with specific site spacing and topologies, monument stability requirements, instrumentation types, data sampling and telemetry requirements. Currently there are hundreds or even thousands of these networks in operation - Can the data collected by this diverse heterogeneous set of networks be assimilated into operational crustal deformation network analyses and providing new information for constraining models of Earth processes? Are there new opportunities in Earth Science research presented by these combined super networks? As case study we will compare the performance characteristics of the recently installed Turkish CORS-TR real time network with the Mediterranean Crustal Deformation Monitoring network in operated by MIT since the mid 1980’s. Additionally we will examine how the integration of data from these two networks improves our ability to constrain crustal deformation models of the region. The CORS-TR network is a continuously operating ~150 station GNSS network installed at ~80-100 km spacing across Turkey and North Cyprus. The CORS-TR network was established by Istanbul Kultur University in corporation with General Directorate of Land Registry and Cadastre (GDLRC) and General Command of Mapping (GCM) between May 2006 and May 2009 and was designed to fulfill a range of civil and scientific goals including; real time access to the Turkish cadastral reference frame for surveyors, precise RTK tracking of infrastructure, weather, atmospheric and ionospheric studies and monitoring of earthquake and crustal deformation
NASA Astrophysics Data System (ADS)
Stanciu, A. C.; Russo, R. M.; Mocanu, V. I.; Foster, D. A.; Bremner, P. M.; Torpey, M. E.; Hongsresawat, S.; VanDecar, J. C.; Hole, J. A.; Tikoff, B.
2014-12-01
We present new images of lithospheric structure obtained from P-to-S receiver functions at 86 broadband seismic stations we deployed as part of the EarthScope IDOR experiment. The high-quality data of the IDOR survey allowed us to image in detail the crustal thickness, and the structure of the lithospheric mantle beneath the Blue Mountains accreted terranes and the Paleozoic margin of North America. At the surface, the ~110 Ma suture between the arc terranes and the North American Craton is marked by the dextral transpressional system of the Western Idaho Shear Zone (WISZ), which closely follows the Sr 0.706 isopleth. Long-lived magmatism (108 - 50 Ma) of the Idaho Batholith occurred during and after the WISZ deformation. We constructed more than 7000 P-to-S teleseismic receiver functions using iterative time domain deconvolution, and we used H-k grid search, inversion, and common conversion point stacking (CCP) to image the lithospheric structure. Moho depths vary from 35-40 km beneath the Blue Mountain terranes and the extended regions of central Idaho, and shallow to ~30 km in the central part of the IDOR network. We see a clear break in the continuity of the Moho across the WISZ, with depths increasing from 30 to 40 km east of the WISZ. We also see a strong mid-crustal converting interface at 15-20 km depth beneath the Idaho Batholith on the east side of the WISZ. To the east, beneath the Lost River Range and Pahsimeroi valley, we see an offset in Moho depth possibly related to Basin-and-Range extension. A strong shallow contrast in crustal velocities beneath the West Snake River Plain is probably due to Columbia River Basalts and sedimentary cover. We observe a horizontal, continuous positive amplitude phase at ~90 km across the entire network that we interpret as the Mid-Lithospheric Discontinuity. A second positive amplitude phase varies from ~120 km beneath the accreted terranes to ~140 km beneath the western margin of the North American Craton. The
NASA Astrophysics Data System (ADS)
Kershaw, Philip; Lawrence, Bryan; Lowe, Dominic; Norton, Peter; Pascoe, Stephen
2010-05-01
CEDA (Centre for Environmental Data Archival) based at STFC Rutherford Appleton Laboratory is host to the BADC (British Atmospheric Data Centre) and NEODC (NERC Earth Observation Data Centre) with data holdings of over half a Petabyte. In the coming months this figure is set to increase by over one Petabyte through the BADC's role as one of three data centres to host the CMIP5 (Coupled Model Intercomparison Project Phase 5) core archive of climate model data. Quite apart from the problem of managing the storage of such large volumes there is the challenge of collating the data together from the modelling centres around the world and enabling access to these data for the user community. An infrastructure to support this is being developed under the US Earth System Grid (ESG) and related projects bringing together participating organisations together in a federation. The ESG architecture defines Gateways, the web interfaces that enable users to access data and data serving applications organised into Data Nodes. The BADC has been working in collaboration with US Earth System Grid team and other partners to develop a security system to restrict access to data. This provides single sign-on via both OpenID and PKI based means and uses role based authorisation facilitated by SAML and OpenID based interfaces for attribute retrieval. This presentation will provide an overview of the access control architecture and look at how this has been implemented for CEDA. CEDA has developed an expertise in data access and information services over several years through a number of projects to develop and enhance these capabilities. Participation in CMIP5 comes at a time when a number of other software development activities are coming to fruition. New services are in the process of being deployed alongside services making up the system for ESG. The security system must apply access control across this heterogeneous environment of different data services and technologies. One strand
NASA Astrophysics Data System (ADS)
Erlström, M.; Niemi, A.; Lindström, S.; Gunnarsson, N.; Daher, S. Bou
2012-04-01
llviken Halfgraben. Relatively less sand in wells away from the main faults imply a distal position and/or a lack of accumulation space. In this study special emphasis is in building a 3D site model by using the simulation software Petrel, evaluating geostatistical data as well as stochastic simulations by using different geostatistical algorithms and evaluating the benefits in this. The primary aim has been to produce a 3D model of the distribution patterns of the different facies and the porosity. The results will be used for CO2 injection simulation purposes in the continuing work of CO2 Mustang (EU Fp 7 project).
Crandall, K.R.
1987-08-01
TRACE 3-D is an interactive beam-dynamics program that calculates the envelopes of a bunched beam, including linear space-charge forces, through a user-defined transport system. TRACE 3-D provides an immediate graphics display of the envelopes and the phase-space ellipses and allows nine types of beam-matching options. This report describes the beam-dynamics calculations and gives detailed instruction for using the code. Several examples are described in detail.
NASA Astrophysics Data System (ADS)
Matsumoto, Munehisa; Akai, Hisazumi; Harashima, Yosuke; Doi, Shotaro; Miyake, Takashi
2016-06-01
A classical spin model derived ab initio for rare-earth-based permanent magnet compounds is presented. Our target compound, NdFe12N, is a material that goes beyond today's champion magnet compound Nd2Fe14B in its intrinsic magnetic properties with a simpler crystal structure. Calculated temperature dependence of the magnetization and the anisotropy field agrees with the latest experimental results in the leading order. Having put the realistic observables under our numerical control, we propose that engineering 5d-electron-mediated indirect exchange coupling between 4f-electrons in Nd and 3d-electrons from Fe would most critically help enhance the material's utility over the operation-temperature range.
NASA Astrophysics Data System (ADS)
Withers, K.; Olsen, K. B.; Shi, Z.; Day, S. M.
2015-12-01
We model blind thrust scenario earthquakes matching the fault geometry of 1994 Mw 6.7 Northridge earthquake up to 8 Hz by first performing dynamic rupture propagation using a support operator method (SORD). We extend the ground motion by converting the slip-rate data to a kinematic source for the finite difference wave propagation code AWP-ODC, which incorporates an improved frequency-dependent attenuation approach. This technique has high accuracy for Q values down to 15. The desired Q function is fit to the 'effective' Q over the coarse grained-cell for low Q, and a simple interpolation formula is used to interpolate the weights for arbitrary Q. Here, we use a power-law model Q above a reference frequency in the form Q 0 f^n with exponents ranging from 0.0-0.9. We find envelope and phase misfits only slightly larger than that of the elastic case when compared with that of the frequency-wavenumber solution for both a homogenous and a layered model with a large-velocity contrast. We also include small-scale medium complexity in both a 1D layered model and a 3D medium extracted from SCEC CVM-S4 including a surface geotechnical layer (GTL). We model additional realizations of the scenario by varying the hypocenter location, and find that similar moment magnitudes are generated. We observe that while the ground motion pattern changes, the median ground motion is not affected significantly, when binned as a function of distance, and is within 1 interevent standard deviation from the median GMPEs. We find that intra-event variability for the layered model simulations is similar to observed values of single-station standard deviation. We show that small-scale heterogeneity can significantly affect the intra-event variability at frequencies greater than ~1 Hz, becoming increasingly important at larger distances from the source. We perform a parameter space study by varying statistical parameters and find that the variability is fairly independent of the correlation length
NASA Astrophysics Data System (ADS)
Oldham, Mark
2015-01-01
Radiochromic materials exhibit a colour change when exposed to ionising radiation. Radiochromic film has been used for clinical dosimetry for many years and increasingly so recently, as films of higher sensitivities have become available. The two principle advantages of radiochromic dosimetry include greater tissue equivalence (radiologically) and the lack of requirement for development of the colour change. In a radiochromic material, the colour change arises direct from ionising interactions affecting dye molecules, without requiring any latent chemical, optical or thermal development, with important implications for increased accuracy and convenience. It is only relatively recently however, that 3D radiochromic dosimetry has become possible. In this article we review recent developments and the current state-of-the-art of 3D radiochromic dosimetry, and the potential for a more comprehensive solution for the verification of complex radiation therapy treatments, and 3D dose measurement in general.
NASA Astrophysics Data System (ADS)
Iliesiu, Luca; Kos, Filip; Poland, David; Pufu, Silviu S.; Simmons-Duffin, David; Yacoby, Ran
2016-03-01
We study the conformal bootstrap for a 4-point function of fermions < ψψψψ> in 3D. We first introduce an embedding formalism for 3D spinors and compute the conformal blocks appearing in fermion 4-point functions. Using these results, we find general bounds on the dimensions of operators appearing in the ψ × ψ OPE, and also on the central charge C T . We observe features in our bounds that coincide with scaling dimensions in the GrossNeveu models at large N . We also speculate that other features could coincide with a fermionic CFT containing no relevant scalar operators.
The Digital Space Shuttle, 3D Graphics, and Knowledge Management
NASA Technical Reports Server (NTRS)
Gomez, Julian E.; Keller, Paul J.
2003-01-01
The Digital Shuttle is a knowledge management project that seeks to define symbiotic relationships between 3D graphics and formal knowledge representations (ontologies). 3D graphics provides geometric and visual content, in 2D and 3D CAD forms, and the capability to display systems knowledge. Because the data is so heterogeneous, and the interrelated data structures are complex, 3D graphics combined with ontologies provides mechanisms for navigating the data and visualizing relationships.
NASA Astrophysics Data System (ADS)
Greenwood, J.; Rucker, D.; Levitt, M.; Yang, X.; Lagmanson, M.
2007-12-01
High Resolution Resistivity data is currently used by hydroGEOPHYSICS, Inc to detect and characterize the distribution of suspected contaminant plumes beneath leaking tanks and disposal sites within the U.S. Department of Energy Hanford Site, in Eastern Washington State. The success of the characterization effort has led to resistivity data acquisition in extremely large survey areas exceeding 0.6 km2 and containing over 6,000 electrodes. Optimal data processing results are achieved by utilizing 105 data points within a single finite difference or finite element model domain. The large number of measurements and electrodes and high resolution of the modeling domain requires a model mesh of over 106 nodes. Existing commercially available resistivity inversion software could not support the domain size due to software and hardware limitations. hydroGEOPHYSICS, Inc teamed with Advanced Geosciences, Inc to advance the existing EarthImager3D inversion software to allow for parallel-processing and large memory support under a 64 bit operating system. The basis for the selection of EarthImager3D is demonstrated with a series of verification tests and benchmark comparisons using synthetic test models, field scale experiments and 6 months of intensive modeling using an array of multi-processor servers. The results of benchmark testing show equivalence to other industry standard inversion codes that perform the same function on significantly smaller domain models. hydroGEOPHYSICS, Inc included the use of 214 steel-cased monitoring wells as "long electrodes", 6000 surface electrodes and 8 buried point source electrodes. Advanced Geosciences, Inc. implemented a long electrode modeling function to support the Hanford Site well casing data. This utility is unique to commercial resistivity inversion software, and was evaluated through a series of laboratory and field scale tests using engineered subsurface plumes. The Hanford site is an ideal proving ground for these methods due
An aerial 3D printing test mission
NASA Astrophysics Data System (ADS)
Hirsch, Michael; McGuire, Thomas; Parsons, Michael; Leake, Skye; Straub, Jeremy
2016-05-01
This paper provides an overview of an aerial 3D printing technology, its development and its testing. This technology is potentially useful in its own right. In addition, this work advances the development of a related in-space 3D printing technology. A series of aerial 3D printing test missions, used to test the aerial printing technology, are discussed. Through completing these test missions, the design for an in-space 3D printer may be advanced. The current design for the in-space 3D printer involves focusing thermal energy to heat an extrusion head and allow for the extrusion of molten print material. Plastics can be used as well as composites including metal, allowing for the extrusion of conductive material. A variety of experiments will be used to test this initial 3D printer design. High altitude balloons will be used to test the effects of microgravity on 3D printing, as well as parabolic flight tests. Zero pressure balloons can be used to test the effect of long 3D printing missions subjected to low temperatures. Vacuum chambers will be used to test 3D printing in a vacuum environment. The results will be used to adapt a current prototype of an in-space 3D printer. Then, a small scale prototype can be sent into low-Earth orbit as a 3-U cube satellite. With the ability to 3D print in space demonstrated, future missions can launch production hardware through which the sustainability and durability of structures in space will be greatly improved.
NASA Astrophysics Data System (ADS)
Iizuka, Keigo
2008-02-01
In order to circumvent the fact that only one observer can view the image from a stereoscopic microscope, an attachment was devised for displaying the 3D microscopic image on a large LCD monitor for viewing by multiple observers in real time. The principle of operation, design, fabrication, and performance are presented, along with tolerance measurements relating to the properties of the cellophane half-wave plate used in the design.
Visualization of 3D Geological Data using COLLADA and KML
NASA Astrophysics Data System (ADS)
Choi, Yosoon; Um, Jeong-Gi; Park, Myong-Ho
2013-04-01
This study presents a method to visualize 3D geological data using COLLAborative Design Activity(COLLADA, an open standard XML schema for establishing interactive 3D applications) and Keyhole Markup Language(KML, the XML-based scripting language of Google Earth).We used COLLADA files to represent different 3D geological data such as borehole, fence section, surface-based 3D volume and 3D grid by triangle meshes(a set of triangles connected by their common edges or corners). The COLLADA files were imported into the 3D render window of Google Earth using KML codes. An application to the Grosmont formation in Alberta, Canada showed that the combination of COLLADA and KML enables Google Earth to visualize 3D geological structures and properties.
NASA Astrophysics Data System (ADS)
Al-Attar, D.; Woodhouse, J. H.
2011-12-01
Normal mode spectra provide a valuable data set for global seismic tomography, and, notably, are amongst the few geophysical observables that are sensitive to lateral variations in density structure within the Earth. Nonetheless, the effects of lateral density variations on mode spectra are rather subtle. In order, therefore, to reliably determine density variations with in the earth, it is necessary to make use of sufficiently accurate methods for calculating synthetic mode spectra. In particular, recent work has highlighted the need to perform 'full-coupling calculations' that take into account the interaction of large numbers of spherical earth multiplets. However, present methods for performing such full-coupling calculations require diagonalization of large coupling matrices, and so become computationally inefficient as the number of coupled modes is increased. In order to perform full-coupling calculations more efficiently, we describe a new implementation of the direct solution method for calculating synthetic spectra in laterally heterogeneous earth models. This approach is based on the solution of the inhomogeneous mode coupling equations in the frequency domain, and does not require the diagonalization of large matrices. Early implementations of the direct solution method used LU-decomposition to solve the mode coupling equations. However, as the number of coupled modes is increased, this method becomes impractically slow. To circumvent this problem, we solve the mode coupling equations iteratively using the preconditioned biconjugate gradient algorithm. We present a number of numerical tests to display the accuracy and efficiency of this method for performing large full-coupling calculations. In addition, we describe a frequency-domain formulation of the adjoint method for the calculation of Frechet kernels that show the sensitivity of normal mode observations to variations in earth structure. The calculation of such Frechet kernels involves one solution
NASA Astrophysics Data System (ADS)
Kostrzewski, Andrew A.; Aye, Tin M.; Kim, Dai Hyun; Esterkin, Vladimir; Savant, Gajendra D.
1998-09-01
Physical Optics Corporation has developed an advanced 3-D virtual reality system for use with simulation tools for training technical and military personnel. This system avoids such drawbacks of other virtual reality (VR) systems as eye fatigue, headaches, and alignment for each viewer, all of which are due to the need to wear special VR goggles. The new system is based on direct viewing of an interactive environment. This innovative holographic multiplexed screen technology makes it unnecessary for the viewer to wear special goggles.
NASA Technical Reports Server (NTRS)
1992-01-01
Ames Research Center research into virtual reality led to the development of the Convolvotron, a high speed digital audio processing system that delivers three-dimensional sound over headphones. It consists of a two-card set designed for use with a personal computer. The Convolvotron's primary application is presentation of 3D audio signals over headphones. Four independent sound sources are filtered with large time-varying filters that compensate for motion. The perceived location of the sound remains constant. Possible applications are in air traffic control towers or airplane cockpits, hearing and perception research and virtual reality development.
NASA Technical Reports Server (NTRS)
Richter, F. M.; Daly, S. F.; Nataf, H.-C.
1982-01-01
It is experimentally shown that, although steady convective flows are efficient means to heterogeneity within a single cell, they do not produce a dispersal of heterogeneous material over scales that are large by comparison to their depth, which requires that the flow be time-dependent on a time scale comparable to the overturn time. Convection in an internally heated layer does possess this property, and numerical solutions are presently used to study the way in which it disperses a set of neutrally bouyant particles initially confined to a small space. The derived concept of effective diffusivity is applied to the isotopic evolution of the Sm-Nd and Rb-Sr systems, with spatial variations generated by horizontal variations in degree of melting 1.8 billion years ago.
NASA Astrophysics Data System (ADS)
Belle, C. L.; Kress, M. E.; Iraci, L. T.
2009-12-01
On average, 10^8 g of solar system debris impinges on the Earth system each day. It is estimated that a few percent of this material is carbonaceous in nature, yet the fate of this organic material once it enters our atmosphere is unexplored. Much of this incoming material arrives in the form of micrometeoroids which are large enough to suffer drag heating and volatilize their organic material. Preliminary work shows that the organic material contained in particles with diameters on the order of 10-100 um is expected to be volatilized at altitudes of 100-120 km. Observed species include aromatic compounds such as alkybenzenes, phenol, benzonitrile, naphthalene, and styrene. Once liberated, these molecules may be transformed by processes at the boundary of space, or may survive to be mixed throughout the atmosphere. Sulfuric acid particles exist in Earth's upper atmosphere, and organic compounds often react strongly with this acid. We will report the results of laboratory and theoretical investigations of the interaction of aromatic compounds with surrogate matrices which mimic upper atmospheric particles. These studies will explore how exogenous organic compounds are altered after liberation at altitudes near 100 km and will determine if they survive to reach the surface of the Earth, where they may have provided the starting materials for the evolution of life on Earth or on other bodies.