Dynamical models for the Earth's geoid

The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.

Correlations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.

Temperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.

Measurements of mantle velocities of P waves with a large array

A large array has been used to investigate the P-wave velocity structure of the lower mantle. Linear array processing methods are reviewed and a method of nonlinear processing is presented. Phase velocities, travel times, and relative amplitudes of P waves have been measured with the large array at the Tonto Forest Seismological Observatory in Arizona for 125 earthquakes in the distance range of 30 to 100 degrees. Various models are assumed for the upper 771 km of the mantle and the Wiechert-Herglotz method applied to the phase velocity data to obtain a velocity depth structure for the lower mantle. The phase velocity data indicates the presence of a second-order discontinuity at a depth of 840 km, another at 1150 km, and less pronounced discontinuities at 1320, 1700 and 1950 km. Phase velocities beyond 85 degrees are interpreted in terms of a triplication of the phase velocity curve, and this results in a zone of almost constant velocity between depths of 2670 and 2800 km. Because of the uncertainty in the upper mantle assumptions, a final model cannot be proposed, but it appears that the lower mantle is more complicated than the standard models and there is good evidence for second-order discontinuities below a depth of 1000 km. A tentative lower bound of 2881 km can be placed on the depth to the core. The importance of checking the calculated velocity structure against independently measured travel times is pointed out. Comparisons are also made with observed PcP times and the agreement is good. The method of using measured values of the rate of change of amplitude with distances shows promising results.

Using Earth Deformation Caused by Surface Mass Loading to Constrain the Elastic Structure of the Crust and Mantle

Surface mass loads come in many different varieties, including the oceans, atmosphere, rivers, lakes, glaciers, ice caps, and snow fields. The loads migrate over Earth's surface on time scales that range from less than a day to many thousand years. The weights of the shifting loads exert normal forces on Earth's surface. Since the Earth is not perfectly rigid, the applied pressure deforms the shape of the solid Earth in a manner controlled by the material properties of Earth's interior. One of the most prominent types of surface mass loading, ocean tidal loading (OTL), comes from the periodic rise and fall in sea-surface height due to the gravitational influence of celestial objects, such as the moon and sun. Depending on geographic location, the surface displacements induced by OTL typically range from millimeters to several centimeters in amplitude, which may be inferred from Global Navigation and Satellite System (GNSS) measurements with sub-millimeter precision. Spatiotemporal characteristics of observed OTL-induced surface displacements may therefore be exploited to probe Earth structure. In this thesis, I present descriptions of contemporary observational and modeling techniques used to explore Earth's deformation response to OTL and other varieties of surface mass loading. With the aim to extract information about Earth's density and elastic structure from observations of the response to OTL, I investigate the sensitivity of OTL-induced surface displacements to perturbations in the material structure. As a case study, I compute and compare the observed and predicted OTL-induced surface displacements for a network of GNSS receivers across South America. The residuals in three distinct and dominant tidal bands are sub-millimeter in amplitude, indicating that modern ocean-tide and elastic-Earth models well predict the observed displacement response in that region. Nevertheless, the sub-millimeter residuals exhibit regional spatial coherency that cannot be explained entirely by random observational uncertainties and that suggests deficiencies in the forward-model assumptions. In particular, the discrepancies may reveal sensitivities to deviations from spherically symmetric, non-rotating, elastic, and isotropic (SNREI) Earth structure due to the presence of the South American craton.

Ichthyoplankton community in the Columbia River plume off Oregon: effects of fluctuating oceanographic conditions

Ichthyoplankton samples were collected at approximately 2-week intervals, primarily during spring and summer 1999−2004, from two stations located 20 and 30 km from shore near the Columbia River, Oregon. Northern anchovy (Engraulis mordax) was the most abundant species collected, and was the primary species associated with summer upwelling conditions, but it showed significant interannual and seasonal fluctuations in abundance and occurrence. Other abundant taxa included sanddabs (Citharichthys spp.), English sole (Parophrys vetulus), and blacksmelts (Bathylagidae). Two-way cluster analysis revealed strong species associations based primarily on season (before or after the spring transition date). Ichthyoplankton abundances were compared to biological and environmental data, and egg and larvae abundances were found to be most correlated with sea surface temperature. The Pacific Decadal Oscillation changed sign (from negative to positive) in late 2002 and indicated overall warmer conditions in the North Pacific Ocean. Climate change is expected to alter ocean upwelling, temperatures, and Columbia River flows, and consequently fish eggs and larvae distributions and survival. Long-term research is needed to identify how ichthyoplankton and fish recruitment are affected by regional and largescale oceanographic proces

Geocronologia (40Ar/39Ar e U-Pb), petrografia e litogeoquímica da intrusão alcalina do Marapicu RJ

A intrusão alcalina do Marapicu é uma intrusão localizada no maciço Marapicu-Gericinó-Mendanha situado na região metropolitana do Rio de Janeiro. Este maciço é formado por dois corpos alcalinos: Marapicu e Mendanha que fazem parte do lineamento magmático Poços de Caldas-Cabo Frio. Este lineamento inclui dezenas de corpos ígneos alcalinos de idade Cretácea com uma direção preferencial WNW-ESE. Os litotipos mais abundantes do Maciço Marapicu são representados por nefelina sienitos e sienitos de caráter plutônico, além de, fonolitos caracterizados por intrusões rasas geralmente em forma de diques. Além desses litotipos foram amostradas duas rochas com características químicas de magma parental (lamprófiro e fonolito tefrítico), porém, essas duas amostras não apresentam relação genética com as demais. Também foi amostrado um nefelina sienito que possui sodalita azul como feldspatóide, sendo assim, chamado de nefelina sodalita sienito. Entre os fonolitos coletados para esse trabalho, uma amostra apresenta granada melanita em sua assembleia mineralógica, e esta foi então denominada melanita fonolito. Quimicamente as rochas do Marapicu formam uma série alcalina predominantemente insaturada em sílica, miaskítica e metaluminosa. Dentro desta série se observam duas suítes sendo uma potássica (predominante) e outra sódica. A evolução química do corpo se deu por processo de cristalização fracionada com ou sem assimilação de crosta continental provavelmente dentro de uma fonte mantélica enriquecida. Duas idades de cristalização foram obtidas para o Maciço do Marapicu sendo uma idade 40Ar/39Ar de 80,46 0,58 Ma em hornblenda, e uma idade U-Pb em zircão bastante concordante de 78,0 2,1 Ma. Os dados apresentados aqui em conjunto com dados da literatura apontam para dois modelos geodinâmicos de geração dos corpos alcalinos do sudeste brasileiro, um considera a existência de uma pluma mantélica gerada na astenosfera, o outro tem por base a hipótese de flexura crustal e considera que a carga de sedimentos depositados na plataforma continental exerceria esforços que provocariam fraturas profundas permitindo a ascenção desses magmas. O presente trabalho vem para contribuir no entendimento do alojamento dos corpos alcalinos do sudeste brasileiro através do estudo especifico do Maciço Marapicu em conjunto com dados da literatura