Thermal models, lithosphere thickness and heat flow in South Portugal. Some comments about the subject

Reinterpretation of old heat flow data or use of new data and new techniques of detection of the temperature under the surface have conducted to new heat flow density values in some regions of the globe. The problem of ice melting in Greenland and Antarctica caught the public's attention to the importance of knowledge on heat flow values and thermal structure of the globe. In the last years, several models were presented trying to obtain lithosphere and Moho thickness of the Iberia Peninsula. The work we intend to present is related with the SW part of the Iberia Peninsula ( south of the Ossa Morena zone, South Portuguese Zone and Algarve). The results obtained show a decrease in the thickness of the crust and the lithosphere in this region. Density anomalies in the crust are also referred. I intend to make the connection between the results of these models and the heat flow thermal conductivity, heat production and geological data available for the region, trying to explain the results of heat flow density data obtained.

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

Abstract In the last years, several models have been presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. Data obtained in the region was collected and deviations from average values used were detected. In this work, models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, were used to adjust the model. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. This work shows the importance of trying to know and understand the thermal structure of the region.

Lithosphere Thickness, Heat Flow and Moho Depth in the South of Portugal

ABSTRACT In the last years, several models were presented trying to obtain lithosphere and Moho thickness in the Iberian Peninsula, using data related to geoid elevation and topography, gravity, seismicity and thermal analysis. The results obtained show a decrease in the thickness of the crust and the lithosphere in the SW part of the Iberian Peninsula. Density anomalies in the crust are also referred. The work I intend to present is related with the south of the Ossa Morena Zone, the South Portuguese Zone and the Algarve, in the south of Portugal. Data obtained in the region was collected and deviations from average values used were detected. Models were made taking into account the specific characteristics of the region. Heat flow, thermal conductivity, heat production, topography, gravity, seismic and geological data available for the region, are used to adapt the models. A special attention will be given to the spatial variation of heat flow values and to Moho depth in the region. The results show that this region is different from other parts of the Iberian Peninsula and a special attention must be given to it. The different values obtained using seismic, gravity, and geoid height data, and the results obtained with models using thermal data shows the importance of trying to know and understand the thermal structure of the regions. Problems related with the use of average values will be focused.

Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

The Lithium, Boron and Beryllium content of serpentinized peridotites from ODP Leg 209 (Sites 1272A and 1274A): Implications for lithium and boron budgets of oceanic lithosphere

Despite the key importance of altered oceanic mantle as a repository and carrier of light elements (B, Li, and Be) to depth, its inventory of these elements has hardly been explored and quantified. In order to constrain the systematics and budget of these elements we have studied samples of highly serpentinized (>50%) spinel harzburgite drilled at the Mid-Atlantic Ridge (Fifteen-Twenty Fracture zone, ODP Leg 209, Sites 1272A and 1274A). In-situ analysis by secondary ion mass spectrometry reveals that the B, Li and Be contents of mantle minerals (olivine, orthopyroxene, and clinopyroxene) remain unchanged during serpentinization. B and Li abundances largely correspond to those of unaltered mantle minerals whereas Be is close to the detection limit. The Li contents of clinopyroxene are slightly higher (0.44-2.8 mu g g(-1)) compared to unaltered mantle clinopyroxene, and olivine and clinopyroxene show an inverse Li partitioning compared to literature data. These findings along with textural observations and major element composition obtained from microprobe analysis suggest reaction of the peridotites with a mafic silicate melt before serpentinization. Serpentine minerals are enriched in B (most values between 10 and 100 mu g g(-1)), depleted in Li (most values below I mu g g(-1)) compared to the primary phases, with considerable variation within and between samples. Be is at the detection limit. Analysis of whole rock samples by prompt gamma activation shows that serpentinization tends to increase B (10.4-65.0 mu g g(-1)), H2O and Cl contents and to lower Li contents (0.07-3.37 mu g g(-1)) of peridotites, implying that-contrary to alteration of oceanic crust-B is fractionated from Li and that the B and Li inventory should depend essentially on rock-water ratios. Based on our results and on literature data, we calculate the inventory of B and Li contained in the oceanic lithosphere, and its partitioning between crust and mantle as a function of plate characteristics. We model four cases, an ODP Leg 209-type lithosphere with almost no igneous crust, and a Semail-type lithosphere with a thick igneous crust, both at I and 75 Ma, respectively. The results show that the Li contents of the oceanic lithosphere are highly variable (17-307 kg in a column of I m x I m x thickness of the lithosphere (kg/col)). They are controlled by the primary mantle phases and by altered crust, whereas the B contents (25-904 kg/col) depend entirely on serpentinization. In all cases, large quantities of B reside in the uppermost part of the plate and could hence be easily liberated during slab dehydration. The most prominent input of Li into subduction zones is to be expected from Semail-type lithosphere because most of the Li is stored at shallow levels in the plate. Subducting an ODP Leg 209-type lithosphere would mean only very little Li contribution from the slab. Serpentinized mantle thus plays an important role in B recycling in subduction zones, but it is of lesser importance for Li. (C) 2008 Elsevier Ltd. All rights reserved.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

Lithospheric thickness of Pangea

Surface wave tomography, using the fundamental Rayleigh wave velocities and those of higher modes between 1 and 4 and periods between 50 and 160 s, is used to image structures with a horizontal resolution of ~250 km and a vertical resolution of ~50 km to depths of ~300 km in the mantle. A new model, PM_v2_2012, obtained from 3×10**6 seismograms, agrees well with earlier lower resolution models. It is combined with temperature estimates from oceanic plate models and with pressure and temperature estimates from the mineral compositions of garnet peridotite nodules to generate a number of estimates of SV(P,T) based on geophysical and petrological observations alone. These are then used to estimate the unrelaxed shear modulus and its derivatives with respect to pressure and temperature, which agree reasonably with values from laboratory experiments. At high temperatures relaxation occurs, causing the shear wave velocity to depend on frequency. This behaviour is parameterised using a viscosity to obtain a Maxwell relaxation time. The relaxation behaviour is described using a dimensionless frequency, which depends on an activation energy E and volume Va. The values of E and Va obtained from the geophysical models agree with those from laboratory experiments on high temperature creep. The resulting expressions are then used to determine the lithospheric thickness from the shear wave velocity variations. The resolution is improved by about a factor of two with respect to earlier models, and clearly resolves the thick lithosphere beneath active intracontinental belts that are now being shortened. The same expressions allow the three dimensional variations of the shear wave attenuation and viscosity to be estimated.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the center of Iberia fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. This heterogeneity has been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Total sediment thickness grid of the Southern Pacific Ocean off West Antarctica, links to NetCDF files

Paleotopographic models of the West Antarctic margin, which are essential for robust simulations of paleoclimate scenarios, lack information on sediment thickness and geodynamic conditions, resulting in large uncertainties. A new total sediment thickness grid spanning the Ross Sea-Amundsen Sea-Bellingshausen Sea basins is presented and is based on all the available seismic reflection, borehole, and gravity modeling data offshore West Antarctica. This grid was combined with NGDC's global 5 arc minute grid of ocean sediment thickness (Whittaker et al., 2013, doi:10.1002/ggge.20181) and extends the NGDC grid further to the south. Sediment thickness along the West Antarctic margin tends to be 3-4 km larger than previously assumed. The sediment volume in the Bellingshausen, Amundsen, and Ross Sea basins amounts to 3.61, 3.58, and 2.78 million km³, respectively. The residual basement topography of the South Pacific has been revised and the new data show an asymmetric trend over the Pacific-Antarctic Ridge. Values are anomalously high south of the spreading ridge and in the Ross Sea area, where the topography seems to be affected by persistent mantle processes. In contrast, the basement topography offshore Marie Byrd Land cannot be attributed to dynamic topography, but rather to crustal thickening due to intraplate volcanism. Present-day dynamic topography models disagree with the presented revised basement topography of the South Pacific, rendering paleotopographic reconstructions with such a limited dataset still fairly uncertain.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.

Effects Of Misalignments In The Retinal Nerve Fiber Layer Thickness Measurements With Spectral Domain Optical Coherence Tomography.

Purpose. To investigate misalignments (MAs) on retinal nerve fiber layer thickness (RNFLT) measurements obtained with Cirrus(©) SD-OCT. Methods. This was a retrospective, observational, cross-sectional study. Twenty-seven healthy and 29 glaucomatous eyes of 56 individuals with one normal exam and another showing MA were included. MAs were defined as an improper alignment of vertical vessels in the en face image. MAs were classified in complete MA (CMA) and partial MA (PMA), according to their site: 1 (superior, outside the measurement ring (MR)), 2 (superior, within MR), 3 (inferior, within MR), and 4 (inferior, outside MR). We compared RNFLT measurements of aligned versus misaligned exams in all 4 sectors, in the superior area (sectors 1 + 2), inferior area (sectors 3 + 4), and within the measurement ring (sectors 2 + 3). Results. RNFLT measurements at 12 clock-hour of eyes with MAs in the superior area (sectors 1 + 2) were significantly lower than those obtained in the same eyes without MAs (P = 0.043). No significant difference was found in other areas (sectors 1 + 2 + 3 + 4, sectors 3 + 4, and sectors 2 + 3). Conclusion. SD-OCT scans with superior MAs may present lower superior RNFLT measurements compared to aligned exams.

The Effect of Prostaglandin Analogs and Prostamide on Central Corneal Thickness

Objective: The aim of this study was to verify the influence of prostaglandin analogs and prostamide on central corneal thickness (CCT). Methods: A prospective analysis was done of CCT in glautomatous patients submitted to monotherapy with prostaglandin analogs (latanoprost 0.005% or travoprost 0.004%) or prostamide (bimatoprost 0.03%) during an 8-week period. A control group of patients without any ocular medication was also evaluated. CCT measurements were performed with a commercially available ultrasound pachymeter. A total of 73 patients were included in this study. Mean age was 68.5 +/- 9.2 (range, 48-85) years old. Results: A statistically significant reduction in CCT was observed in all groups, except the control group (n = 21): Bimatoprost 0.03% group (n = 21): 544.41 +/- 35.4 vs. 540.35 +/- 35.9 mu m (P = 0.039); travoprost 0.004% group (n = 17): 538.47 +/- 32.0 vs. 532.25 +/- 30.4 mu m (P = 0.009); latanoprost 0.005% group (n = 14): 548.57 +/- 32.4 vs. 543.88 +/- 35.6 mu m (P = 0.036). Conclusion: Topical therapy with prostaglandin analogs and bimatoprost is associated with CCT reduction over a period of at least 8 weeks.