2:00pm  2:20pmKeynote: With GOCE and Swarm Towards the 3D Lithosphere
Jörg Ebbing
Kiel University, Germany
Satellite measurements of the Earth gravity and magnetic field as from the GOCE and Swarm satellite missions have an increasing resolution. Their global coverage and specific characteristics (e.g. gradient measurements for GOCE) make them an ideal data sets for studying the lithosphere in combination with seismological models. Such a combined lithosphere is one of the goals of the ESA STSE “3D Earth”.
The prime objective of our project is to integrate seismological models and satellite observation towards a consistent image of the crust and upper mantle in 3D. Satellite gravity and (electro) magnetic data help to transfer velocity images towards composition and temperature that reflect the tectonic state and evolution of the Earth and offer a novel understanding of the processes that shape our planet.
We will analyse the limitations and sensitivities of the different geophysical methods in the context of their imaging capability and in combination for forward and inverse modelling of the Earth’s internal structure. Such analysis will for example help to assess the role of isostatic (lithospheric) and dynamic (deep Earth) effects in shaping the surface of the Earth.
2:20pm  2:35pmLargeScale Constraints on Tectonics and Macroscopic Magnetic Properties of the Earth's Lithosphere from the Swarm Constellation and CHAMP
Michael E Purucker^{1}, Suzanne McEnroe^{2}
^{1}NASA, Goddard Space Flight Center; ^{2}Norwegian University of Science and Technology
We constrain largescale macroscopic magnetic properties (susceptibility, remanence and their geometry) of the earth's lithosphere in three regimes (continent, ocean, and subducted slab) using the latest high degree CHAOS and CM models, and compare those constraints to measured compilations of magnetic properties. At the highest degrees, these models still reflect a significant contribution from CHAMP, but the difference with the latest Swarm models is narrowing. In the oceanic realm, the magnetic symmetry across spreading centers, and presence of transform faults, are clearly evident only in the fastest spreading ridges (e.g. East Pacific Rise). Both the CHAOS and CM models can now extract alongtrack features in the lithosphere. We discuss the implications for largescale oceanic and continental tectonics in these latest models, and anticipate some of the results when the Swarm satellites are at lower altitudes.
2:35pm  2:50pmA New Lithospheric Field Model based on CHAMP and Swarm Magnetic Satellite Data
Nils Olsen^{1}, Dhananjay Ravat^{2}, Christopher C Finlay^{1}
^{1}Technical University of Denmark, Denmark; ^{2}University of Kentucky
We used magnetic field observations from the last two years of the CHAMP satellite mission (at altitudes between 280 and 350 km), augmented with Swarm satellite gradient data, to determine a model of the lithospheric field. We first subtracted predictions of the core and largescale magnetospheric field as given by the CHAOS6 model from the observations. No further data treatment was applied (e.g., no orbitbyorbit filtering or “line leveling” of the individual satellite tracks was done). The lithospheric field is described by 35,000 point sources at 100 km depth below Earth’s surface. We estimate the amplitudes of these point sources from the observations using an Iteratively Reweighted Least Squares approach with robust weighting, to account for nonGaussian data errors. The model is regularized by minimization of the L1 norm of the vertical magnetic field at the surface (WGS84 ellipsoid). In a final step we expand our pointsource model in series of spherical harmonics up to degree and order N=185, accounting for div B = 0.
We will present our modeling approach and discuss the obtained lithospheric model, in particular how it compares with other lithospheric field models and with independent data sets.
2:50pm  3:05pmAn Updated Global lithospheric Model by Implementing 3D Joint Inversion of Gravity, Geoid, Topography and the Gravity Gradients in Spherical Coordinates
Farshad Salajegheh, Juan Carlos Afonso
CCFS, Department of Earth and Planetary Sciences,Macquarie University, Sydney, New South Wales, Australia
The everincreasing interest in the generation of highresolution models of the whole of the lithosphere by industry and academia has stimulated the development of probabilistic inversions of multiple geophysical datasets. These types of inversions are computationally expensive and rely on efficient algorithms to solve the forward problems. Potential field data (e.g. gravity, geoid and GOCE gravity gradients) is of particular interest given their complementary sensitivities to both crustal and deep lithospheric structure, particularly when jointly inverted with seismic data. When the inversion is performed in 3D Cartesian geometry, very efficient solvers exist, and this why current implementations of multiobservable probabilistic inversions rely on such geometry. This, however, precludes the application of the inversion method to largescale domains (e.g. continental or global studies).
For the forward calculation, we have implemented the prism approximation for calculating the tesseroid’s gravity effects. Thus, the lack of modelling singularities at the prism’s surface allows us to compute the gravity effects on the surface. If a prism instead of a tesseroid is used in the spherical system, the effect of the potential or its first and second derivatives is obtained in the local edge system of the prism. Therefore, these effects have to be transformed to the local reference frame of the observation point.
We have applied a smoothnessconstrained leastsquares algorithm based on the joint interpretation of freeair gravity, geoid, topography and gravity gradients data. The used inversion method is a linearized iterative inversion procedure in order to obtain variations in MOHO depth, average crustal density, average asthenosphere density and LAB depth. This method is able to provide in relatively efficient time a simplified global lithospheric model, which may utilize as the lateral extension of the model for the large scale study in a more complicated modelling using probabilistic inversion algorithm (e.g. LitMod3D_inv) to directly image the thermal and compositional structure of the lithosphere and sublithospheric upper mantle in the spherical coordinates. The prior information for this study has been extracted from the LITHO1.0 model.
In the meantime, the ScaLAPCK package has been used for parallelizing the inversion calculations. ScaLAPACK is a collection of routines written in Fortran77 for solving dense and banded linear systems, least squares problems, eigenvalue problems, and singular value problems in the distributed memory concurrent supercomputers.
3:05pm  3:20pmImaging Smallscale Seafloor and Subseafloor Tectonic Fabric Using Satellite Altimetry
David T. Sandwell^{1}, R. Dietmar Müller^{2}, Kara J. Matthews^{3}
^{1}UCSD, United States of America; ^{2}Univ. of Sydney, Australia; ^{3}Univ. of Oxford, United Kingdom
Marine gravity anomalies derived from satellite radar altimetry now provide an unprecedented resolution for mapping smallscale seafloor and subseafloor tectonic fabric. Most of the new information comes from the CryoSat2 satellite, which has routinely collected altimetry data over ice, land, and ocean since July 2010. The satellite has a long 369day repeat cycle resulting in an average ground track spacing of 3.5 km at the equator. To date it has completed more than 6 geodetic mappings of the ocean surface. These data are augmented by a complete 14month geodetic mapping of the ocean surface by Jason1 from its lower inclination orbit of 66 ̊ that compliments the higher inclination orbit CryoSat2 (88 ̊). The most recent global marine gravity anomaly map based on all of these geodetic mission data with 2pass retracking for optimal range precision has an accuracy that is 2 times better than the maps derived from Geosat and ERS1. The new data reveal the detailed fabric of fracture zones, previously unmapped, now extinct oceanic microplates in the central Pacific, and fault networks buried beneath thick sediments along continental margins. By combining satellite altimetry with marine magnetic anomalies and seafloor age dates from rock samples we are able to pinpoint the geometry and age of major plate reorganisations, particularly the enigmatic 100 Ma event, which occurred during the Cretaceous Magnetic Superchron.
