We present global and regional correlations between whole-rock values of Sr/Y and La/Yb and crustal thickness for intermediate rocks from modern subduction-related magmatic arcs formed around the Pacific. These correlations bolster earlier ideas that various geochemical parameters can be used to track changes of crustal thickness through time in ancient subduction systems. Inferred crustal thicknesses using our proposed empirical fits are consistent with independent geologic constraints for the Cenozoic evolution of the central Andes, as well as various Mesozoic magmatic arc segments currently exposed in the Coast Mountains, British Columbia, and the Sierra Nevada and Mojave- Transverse Range regions of California. We propose that these geochemical parameters can be used, when averaged over the typical lifetimes and spatial footprints of composite volcanoes and their intrusive equivalents to infer crustal thickness changes over time in ancient orogens.

This dissertation investigates the crustal structure of the eastern Basin and Range Province in the western United States and its relationship with the present-day extensional regime governing this region. The use of combined results from different geophysical methods provide a better understanding of the subsurface crustal structure and the processes involved in this extensional deformation. Teleseismic receiver functions were used to create a uniformly sampled map of the crustal thickness variations and stacked images of the crust beneath the majority of the state of Utah, which provide additional constraints on the seismic characteristics of the crust and upper mantle. These results reveal crustal variations characterized by a distinct change in crustal thickness that closely follows the surface trace of the Wasatch fault, with differences in depth of up to 10 km across a distance of less than 55 km. Analysis of seismic reflection profiles, horizontal and vertical crustal velocities from continuous GPS, and surface geology provide new constraints on the relationships between interseismic strain accumulation, subsurface fault geometry, and geologic slip rates on seismogenic faults. Seismic reflection data show recent activity along high-angle normal faults that become listric with depth, sole into preexisting décollements, reactivating them, and appear to be connected at depth with a regionally extensive detachment horizon. GPS data reveal present-day crustal extension of ~3 mm/yr and no net vertical motion between the Colorado Plateau and eastern Basin and Range. Inverse modeling results of the crustal deformation data include a low-angle dislocation (~8-20°) at a locking depth of ~7-10 km, consistent with the interpreted seismic data, and slipping at 3.2±0.2 mm/yr, suggesting an active regionally extensive sub-horizontal surface beneath the eastern Basin and Range. A test of this hypothesis using seismic data interpretation as the basis for a forward strain accumulation model shows that displacement across a deep low-angle detachment imaged seismically is also consistent with geodetic velocities. Seismic and geodetic data support a model for eastern Basin and Range mechanics wherein diffuse permanent strain of the upper crust by multiple discrete faults is facilitated by displacement along a single low-angle detachment at midcrustal depth.

In this study, we present results from teleseismic receiver functions, in order to investigate the crustal thickness and Vp/Vs ratio beneath northern Peru. A total number of 981 receiver functions were analyzed, from data recorded by 28 broadband seismic stations from the Peruvian permanent seismic network, the regional temporary SisNort network and one CTBTO station. The Moho depth and average crustal Vp/Vs ratio were determined at each station using the H-k stacking technique to identify the arrival times of primary P to S conversion and crustal reverberations (PpPms, PpSs + PsPms). The results show that the Moho depth correlates well with the surface topography and varies significantly from west to east, showing a shallow depth of around 25 km near the coast, a maximum depth of 55-60 km beneath the Andean Cordillera, and a depth of 35-40 km further to the east in the Amazonian Basin. The bulk crustal Vp/Vs ratio ranges between 1.60 and 1.88 with the mean of 1.75. Higher values between 1.75 and 1.88 are found beneath the Eastern and Western Cordilleras, consistent with a mafic composition in the lower crust. In contrast values vary from 1.60 to 1.75 in the extreme flanks of the Eastern and Western Cordillera indicating a felsic composition. We find a positive relationship between crustal thickness, Vp/ Vs ratio, the Bouguer anomaly, and topography. These results are consistent with previous studies in other parts of Peru (central and southern regions) and provide the first crustal thickness estimates for the high cordillera in northern Peru.

Isostatic equilibrium is commonly defined as the state achieved when there are no lateral gradients in hydrostatic pressure, and thus no lateral flow, at depth within the lower viscosity mantle that underlies a planetary body's outer crust. In a constant-gravity Cartesian framework, this definition is equivalent to the requirement that columns of equal width contain equal masses. Here we show, however, that this equivalence breaks down when the spherical geometry of the problem is taken into account. Imposing the "equal masses" requirement in a spherical geometry, as is commonly done in the literature, leads to significant lateral pressure gradients along internal equipotential surfaces and thus corresponds to a state of disequilibrium. Compared with the "equal pressures" model we present here, the equal masses model always overestimates the compensation depth-by similar to 27% in the case of the lunar highlands and by nearly a factor of 2 in the case of Enceladus. Plain Language Summary "Isostasy" is the principle that, just as an iceberg floats on the water, crustal rocks effectively float on the underlying higher density mantle, which behaves essentially like a fluid on geologic timescales. Although there are subtle inconsistencies among the various ways isostasy can be defined, they have not been historically problematic for bodies like the Earth, where the crust is thin compared with the overall radius. When the thickness of the crust is a nonnegligible fraction of a planetary body's radius, however, it becomes important to take the spherical geometry into account. In this case, the inconsistencies in the definitions can lead to significant discrepancies. Here we argue that one of the most commonly used approaches, which requires equal width columns to contain equal masses, always results in overestimating the crustal thickness. In particular, we suggest that the lunar and Martian highlands crustal thickness may have been overestimated by similar to 27% and similar to 10%, respectively, and that the ice shell thickness for Saturn's small icy moon Enceladus may have been overestimated by nearly a factor of 2.

Para a realização dessa pesquisa foi desenvolvido um estudo da espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais (entre 10ºN e 9ºS). Foram usadas reflexões da Moho de seis sismos de profundidade média para achar a espessura crustal na área de estudo. As ondas pmP e smS (reflexões na Moho) são encontradas como precursoras das fases profundas pP e sS de sismos registrados em estações localizadas a distâncias telessísmicas. A metodologia utilizada exigia sismos com profundidade maiores que 100 km e magnitude maiores que 6. Entretanto, devido à baixa amplitude apresentada pelas ondas pmP e smS, foi necessário realizar um empilhamento de vários sismogramas para realçar as chegadas das ondas analisadas. A área estudada foi dividida em três partes: área norte entre 6ºN e 0º(na Colômbia), área central entre 2ºS e 5ºS (no Equador) e área sul entre 6ºS e 9ºS (na borda entre Peru e Brasil). Na área norte foram encontradas espessuras crustais entre 26 e 56 km em quatro pontos de reflexão de dois sismos; na área central foram encontradas espessuras da crosta entre 40 e 63 km para três pontos de reflexão de um sismo; finalmente, na região sul, foram encontradas espessuras crustais entre 35 e 40 km para 7 pontos de reflexão de 3 sismos. Foi obtido a relação Vp/Vs = 1.79 ± 0.16 para a região norte dos Andes, usando as diferenças de tempo sS-smS e pP-pmP. Os resultados obtidos vão aumentar o banco de dados de espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais, como também possibilitaram uma melhora do modelo de espessura crustal desenvolvido por Assumpção et al. (2013) para América do Sul, permitindo assim a analise e comparação entre os dados obtidos nesse estudo e o modelo crustal de Assumpção et al. (2013). / In this work we studied the crustal thickness in northern Andes and northern central Andes (between 10ºN and 9ºS). We analyzed the reflections from the underside of the Moho for six intermediate earthquakes occurred in the study area to estimate the crustal thickness at the bounce point. The pmP and smS phases (reflections at the Moho) are identied as precursors to the depth phase pP and sS, respectively, when recorded at teleseismic distances. This method require events of magnitude > 6 and depth > 100 km. In order to better identify those two reflections, it was necessary to stacking traces from different stations. The study area was divided in three sections: northern area between 6ºN - 0º (Colombia), in this area we analyzed two earthquakes and the crustal thicknesses in four bounce points were found; central area between 2ºS- 5ºS (Ecuador), in this area we studied one earthquake and the crustal thicknesses in three bounce points were found; southern area between 6ºS - 9ºS (Peru-Brazil border), where we analyzed three earthquakes and the crustal thicknesses in seven bounce points were found. Crustal thickness found in the southern area varied from 35 to 40 km, from 40 to 63 km in the central area and from 26 to 56 km in the northern area. We estimated a Vp/Vs ratio of 1.79±0.16 for the northern Andes using the time differences sS-smS and pP-pmP relation. Our analysis complements the data base of Moho depths in the Northern Andes and also, it allows a comparison with the model of crustal thickness estimated by Assumpção et al. (2013) confirming the Moho depth on both studies.

Para a realização dessa pesquisa foi desenvolvido um estudo da espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais (entre 10ºN e 9ºS). Foram usadas reflexões da Moho de seis sismos de profundidade média para achar a espessura crustal na área de estudo. As ondas pmP e smS (reflexões na Moho) são encontradas como precursoras das fases profundas pP e sS de sismos registrados em estações localizadas a distâncias telessísmicas. A metodologia utilizada exigia sismos com profundidade maiores que 100 km e magnitude maiores que 6. Entretanto, devido à baixa amplitude apresentada pelas ondas pmP e smS, foi necessário realizar um empilhamento de vários sismogramas para realçar as chegadas das ondas analisadas. A área estudada foi dividida em três partes: área norte entre 6ºN e 0º(na Colômbia), área central entre 2ºS e 5ºS (no Equador) e área sul entre 6ºS e 9ºS (na borda entre Peru e Brasil). Na área norte foram encontradas espessuras crustais entre 26 e 56 km em quatro pontos de reflexão de dois sismos; na área central foram encontradas espessuras da crosta entre 40 e 63 km para três pontos de reflexão de um sismo; finalmente, na região sul, foram encontradas espessuras crustais entre 35 e 40 km para 7 pontos de reflexão de 3 sismos. Foi obtido a relação Vp/Vs = 1.79 ± 0.16 para a região norte dos Andes, usando as diferenças de tempo sS-smS e pP-pmP. Os resultados obtidos vão aumentar o banco de dados de espessura crustal nos Andes setentrionais e na parte norte dos Andes centrais, como também possibilitaram uma melhora do modelo de espessura crustal desenvolvido por Assumpção et al. (2013) para América do Sul, permitindo assim a analise e comparação entre os dados obtidos nesse estudo e o modelo crustal de Assumpção et al. (2013). / In this work we studied the crustal thickness in northern Andes and northern central Andes (between 10ºN and 9ºS). We analyzed the reflections from the underside of the Moho for six intermediate earthquakes occurred in the study area to estimate the crustal thickness at the bounce point. The pmP and smS phases (reflections at the Moho) are identied as precursors to the depth phase pP and sS, respectively, when recorded at teleseismic distances. This method require events of magnitude > 6 and depth > 100 km. In order to better identify those two reflections, it was necessary to stacking traces from different stations. The study area was divided in three sections: northern area between 6ºN - 0º (Colombia), in this area we analyzed two earthquakes and the crustal thicknesses in four bounce points were found; central area between 2ºS- 5ºS (Ecuador), in this area we studied one earthquake and the crustal thicknesses in three bounce points were found; southern area between 6ºS - 9ºS (Peru-Brazil border), where we analyzed three earthquakes and the crustal thicknesses in seven bounce points were found. Crustal thickness found in the southern area varied from 35 to 40 km, from 40 to 63 km in the central area and from 26 to 56 km in the northern area. We estimated a Vp/Vs ratio of 1.79±0.16 for the northern Andes using the time differences sS-smS and pP-pmP relation. Our analysis complements the data base of Moho depths in the Northern Andes and also, it allows a comparison with the model of crustal thickness estimated by Assumpção et al. (2013) confirming the Moho depth on both studies.

The Mohorovičič discontinuity (Moho), which is the surface separating the Earth’s crust from the mantle, is of great interest among geoscientists. The Moho depth can be determined by seismic and gravimetric methods. The seismic methods are expensive, time-consuming and suffer from lack of global coverage of data, while the gravimetric methods use inexpensive and mostly already available global and regional data based on an isostatic model. The main reasons for studying an isostatic model are on one hand the gaps and uncertainties of the seismic models, and, on the other hand, the generous availability of gravity data from global models for the gravimetric-isostatic model. In this study, we present a new gravimetric-isostatic Moho model, called the Vening Meinesz-Moritz (VMM) model. Also, a combined Moho model based on seismic and gravimetric models is presented. Classical isostatic hypotheses assume that the topographic potential is fully compensated at all wavelengths, while is not the case in reality. We found that the maximum degree of compensation for the topographic potential based on the new Moho model is 60, corresponding to the resolution of about 330 km. Other (dynamic) isostatic effects (such as temporal compensation, plate tectonics, post-glacial rebound, etc) should be considered as well, which are disregarded in this thesis. Numerical results imply that the dynamic phenomena affect mostly the long-wavelengths. The VMM model is applied for different purposes. The Moho density contrast is an important parameter for estimating the Moho depth, and we present a technique to simultaneously estimate Moho depth and density contrast by the VMM and seismic models. Another application is the recovery of gravity anomaly from Satellite Gravity Gradiometry (SGG) data by a smoothing technique, and we show that the VMM model performs better than the Airy-Heiskanen isostatic model. We achieved an rms difference of 4 mGal for the gravity anomaly estimated from simulated GOCE data in comparison with EGM08, and this result is better than direct downward continuation of the data without smoothing. We also present a direct method to recover Moho depth from the SGG mission, and we show that the recovered Moho is more or less of the same quality as that obtained from terrestrial gravimetric data (with an rms error of 2 km). Moreover, a strategy is developed for creating substitutes for missing GOCE data in Antarctica, where there is a polar gap of such data. The VMM model is further used for constructing a Synthetic Earth Gravity Model (SEGM). The topographic-isostatic potential is simple to apply for the SEGM, and the latter can be an excellent tool to fill data gaps, extending the EGMs to higher degrees and validating a recovery technique of the gravity field from a satellite mission. Regional and global tests of the SEGM yield a relative error of less than 3 % vs. EGM08 to degree 2160. / QC 20110405

Made available in DSpace on 2015-03-13T17:08:14Z (GMT). No. of bitstreams: 1 MariaFNB.pdf: 1315338 bytes, checksum: 0d0177ce82f07bd3aa4bb82fa98ec33d (MD5) Previous issue date: 2008-06-13 / Conselho Nacional de Desenvolvimento Cient?fico e Tecnol?gico / Crustal thickness and VP/VS estimates are essential to the studies of subsurface geological structures and also to the understanding of the regional tectonic evolution of a given area. In this dissertation, we use the Langston?s (1979) Receiver Function Method using teleseismic events reaching the seismographic station with angles close to the vertical. In this method, the information of the geologic structures close to the station is isolated so that effects related to the instrument response and source mechanics are not present. The resulting time series obtained after the deconvolution between horizontal components contains the larger amplitude referring to the P arrival, followed by smaller arrival caused by the reverberation and conversion of the P-wave at the base of the crust. We also used the HK-Stacking after Zhu & Kanamori (2000) to obtain crustal thickness and Vp/VS estimates. This method works stacking receiver functions so that the best estimates of crustal thickness and Vp/VS are found when the direct P, the Ps wave and the first multiple are coherently stacked. We used five broadband seismographic stations distributed over the Borborema Province, NE Brazil. Crustal thickness and Vp/VS estimates are consistent with the crust-mantle interface obtained using gravity data. We also identified crutal thickening in the NW portion of the province, close to Sobral/CE. Towards the center-north portion of the province, there is an evident crustal thinning which coincides with a geological feature consisting of an alignment of sedimentary basins known as the Cariris-Potiguar trend. Towards the NE portion of the province, in Sol?nea/PB and Agrestina/PE regions, occurs a crustal thickening and a systematic increase in the VP/VS values which suggest the presence of mafic rocks in the lower crust also consistent with the hypothesis of underplating in the region / A estimativa da espessura da crosta terrestre e raz?o VP/VS s?o essenciais para o detalhamento de estruturas e fei??es geol?gicas, al?m de corroborarem para o entendimento da evolu??o tect?nica regional. Neste trabalho, utiliza-se o m?todo da Fun??o do Receptor de Langston (1979) com eventos teless?smicos que incidem sob uma esta??o sismogr?fica com um ?ngulo pr?ximo da vertical. ? necess?rio o isolamento de informa??es das estruturas pr?ximas ao receptor, eliminando informa??es relacionadas ? resposta do instrumento, e ao mecanismo da fonte. O sismograma sint?tico obtido ap?s uma deconvolu??o entre as componentes horizontais possui um pico maior referente ? onda P, seguido por picos menores da onda Ps e m?ltiplas. A Onda Ps ? considerada como onda P convertida em onda S refratada na descontinuidade Moho. Para o c?lculo das estimativas de espessura crustal e raz?es Vp/VS utilizou-se o procedimento HK-Stacking de Zhu & Kanamori (2000). Esse m?todo trabalha com o empilhamento das fun??es do receptor. As melhores estimativas de espessura crustal e raz?o VP/VS s?o encontradas quando as tr?s fases P, Ps e primeira m?ltipla s?o empilhadas coerentemente. Foram utilizadas cinco esta??es sismogr?ficas banda-larga distribu?das estrategicamente pela Prov?ncia Borborema, nordeste do Brasil. As estimativas de espessura crustal e raz?o VP/VS s?o consistentes com o modelo atual da interface crosta-manto usados pela gravimetria. Foi identificado um espessamento crustal na por??o NW da prov?ncia, pr?ximo de Sobral/CE, na borda leste da Bacia do Parna?ba. Em dire??o a por??o centro-norte da Prov?ncia Borborema ? evidente um afinamento crustal, coincidente com a fei??o geol?gica que consiste de um alinhamento de bacias sedimentares conhecido como o trend Cariris-Potiguar. Na por??o NE, nas regi?es de Sol?nea/PB e Agrestina/PE, ocorre um espessamento crustal e aumento da raz?o VP/VS sugerindo presen?a de rochas mais m?ficas na crosta inferior, consistente com a hip?tese de underplating na regi?o

Utilizamos neste trabalho duas metodologias distintas, a função do receptor com ondas P e a função do receptor com ondas S, para mapear variações da crosta e interfaces do manto (litosfera-astenosfera, 410 km e 660 km) em diferentes estações sismográficas na placa Sul-Americana. No estudo da interface litosfera-astenosfera, por ser o primeiro realizado nesta região, utilizamos as estações temporárias do IAG/USP em conjunto com as estações permanentes da rede mundial cobrindo toda a placa Sul-Americana. O estudo para as outras interfaces (Crosta-Manto, 410 km e 660 km) foi feito com caráter regional, buscando detalhar características da crosta e manto na região estável da placa. Para ambos os métodos os traços (sismogramas) foram rotacionados para o sistema LQT, deconvolvidos, agrupados por pontos de perfuração e por estações, e finalmente empilhados. Nos traços empilhados as fases convertidas de interesse (Ps, Ppps, Ppss+Psps e Sp) foram identificadas e interpretadas. Para a parte estável da placa obtivemos um valor médio de espessura da crosta de 39.4±0.6 km, variando desde 31.0±0.5 km para a província Borborema, até 41.3±1.0 km para a bacia do Paraná, onde aplicamos uma correção para descontar o efeito do sedimento. A razão de velocidade para a crosta, Vp/Vs, apresentou valores mais altos para a bacia do Paraná (~1.75±0.08) e região litorânea oriental (>1.74), enquanto que as regiões cratônicas (cráton São Francisco e Amazônico) apresentaram valores de Vp/Vs baixos (<1.72), chegando até 1.68. O valor médio de Vp/Vs para todas as estações analisadas foi de 1.73±0.02. As variações dos tempos para as interfaces do manto mostraram boa correlação com resultados de tomografia sísmica de outros trabalhos, indicando alterações de até 5% na velocidade das ondas sísmicas para o manto superior sob os crátons, uma deflexão de até 15 km na interface de 660 km para a região Sul da bacia do Paraná e se mostraram bem correlacionadas com as médias globais para as outras região estudadas. Por fim, a espessura da litosfera apresentou valores desde ~40 km, sob as regiões de ilhas oceânicas, até ~160 km, sob as regiões mais estáveis. Para as regiões oceânicas a espessura da litosfera se mostra correlacionada com a idade da placa. À medida que adentramos a parte continental, o limite litosfera-astenosfera se torna menos proeminente, atingindo profundidades maiores no interior dos continentes e menores para as regiões marginais. Para a zona de subducção, observamos duas possíveis litosferas, uma oceânica, subduzindo junto com a placa de Nazca, e outra pertencente à parte continental. / Two distinct methodologies, the P- and S-wave receiver functions, are used to map variations in the crustal parameters (thickness and Vp/Vs) and mantle interfaces (lithosphere-asthenosphere, 410 km and 660 km) on a number of different seismograph stations located in the South American plate. The results of the S receiver function for the lithosphere-asthenosphere boundary are the first of this kind ever performed in South American continent and showed the large scale variations of this interface. To perform this study we analyze data from various global permanent stations together with all available data from temporary stations operated by the IAG/USP during the last15 years. For both methods the traces (seismograms) were rotated to the LQT system, deconvolved, grouped by piercing points and stations, and finally stacked. In the stacked traces, the converted phases (Ps, Ppps, Ppss+Psps and Sp) were identified and interpreted. Inside the stable part of the plate we found a mean crustal thickness of 39.4±0.6 km, ranging from 31.0±0.5 km in Borborema Province up to 41.3±1.0 km in the Paraná Basin, where we applied a correction to remove the sediment effects on the crustal estimates. The crustal velocity ratios, Vp/Vs, showed higher values for the Paraná Basin (~1.75±0.08) and Ribeira belt (>1.74), while the cratonic regions (São Francisco and Amazon cratons) showed low values of Vp/Vs (<1.72), down to 1.68. The average Vp/Vs obtained for all stations was equal to 1.73±0.02. The observed times of the converted mantle phases presented a good correlation with other tomographic studies, indicating that the upper mantle for the cratonic roots may be characterized by a variation up to 5% in seismic velocities, a 15 km deflection in the South Paraná 660 km discontinuity (probably due to a decreased temperature caused by the subducted slab); for other regions the converted times were close to the global average. As a final result, the lithospheric thickness presented values ranging from ~40 km under oceanic islands, to ~160 km under the stable continental regions. We found that for the oceanic islands the thickness of the lithosphere is correlated with the age of the plate. When we go further inside the continents, the lithosphere-asthenosphere boundary becomes less sharp, reaching larger depths inside the continents and shallower depths near the continental margin. In the Andean subduction area, we observed two possibles lithospheres, one oceanic, subducting together with the Nazca plate, and another belonging to the Continent, parallel to the crust interface.

Glacial isostatic adjustment (GIA) is the process by which the solid Earth responds to past and present-day changes in glaciers, ice caps, and ice sheets. This thesis focuses on vertical crustal motion of the Earth caused by GIA, which is influenced by several factors including lithosphere thickness, mantle viscosity profile, and changes to the thickness and extent of surface ice. The viscosity of the mantle beneath Antarctica is a poorly constrained quantity due to the rarity of relative sea-level and heat flow observations. Other methods for obtaining a better-constrained mantle viscosity model must be investigated to obtain more accurate GIA model predictions. The first section of this study uses seismic wave tomography to determine mantle viscosity. By calculating the deviation of the P- and S-wave velocities relative to a reference Earth model (PREM), the viscosity can be determined. For Antarctica mantle viscosities obtained from S20A (Ekstrom and Dziewonski, 1998) seismic tomography in the asthenosphere range from 1016 Pa∙s to 1023 Pa∙s, with smaller viscosities beneath West Antarctica and higher viscosities beneath East Antarctica. This agrees with viscosity expectations based on findings from the Basin and Range area of North America, which is an analogue to the West Antarctic Rift System. Section two compares bedrock elevations in Antarctica to crustal thicknesses, to infer mantle temperatures and draw conclusions about mantle viscosity. Data from CRUST 2.0 (Bassin et al., 2000), BEDMAP (Lythe and Vaughan, 2001) and specific studies of crustal thickness in Antarctica were examined. It was found that the regions of Antarctica that are expected to have low viscosities agree with the hot mantle trend found by Hyndman (2010) while the regions expected to have high viscosity are in better agreement with the trend for cold mantle. Bevis et al. (2009) described new GPS observations of crustal uplift in Antarctica and compared the results to GIA model predictions, including IJ05 (Ivins and James, 2005). Here, we have generated IJ05 predictions for a three layered mantle (viscosities ranging over more than four orders of magnitude) and compared them to the GPS observations using a χ2 measure of goodness-of-fit. The IJ05 predictions that agree best with the Bevis et al. observations have a χ2 of 16, less than the null hypothesis value of 42. These large values for the best-fit model indicate the need for model revisions and/or that uncertainties are too optimistic. Equally important, the mantle viscosities of the best-fit models are much higher than expected for West Antarctica. The smallest χ2 values are found for an asthenosphere viscosity of 1021 Pa•s, transition zone viscosity of 1023 Pa∙s and lower mantle viscosity of 2 x 1023 Pa∙s, whereas the expected viscosity of the asthenosphere beneath West Antarctica is probably less than 1020 Pa∙s. This suggests that revisions to the IJ05 ice sheet history are required. Simulated annealing was performed on the ice sheet history and it was found that changes to the recent ice load history have the strongest effect on GIA predictions. / Graduate