The intraplate Porto dos Gauchos seismic zone in the Amazon craton - Brazil

The largest earthquake observed in the stable continental interior of the South American plate occurred in Serra do Tombador, Mato Grosso state - Brazil, on January 31,1955 with a magnitude of 6.2 m(b). Since then no other earthquake has been located near the 1955 epicentre. However, in Porto dos Gauchos, 100 km northeast of Serra do Tombador, a recurrent seismicity has been observed since 1959. Both Serra do Tombador and Porto dos Gauchos are located in the Phanerozoic Parecis basin. Two magnitude 5 earthquakes occurred in Porto dos Gauchos, in 1998 and 2005, with intensities up to VI and V, respectively. These two main shocks were followed by aftershock sequences lasting more than three years each. Local seismic stations have been deployed by the Seismological Observatory of the University of Brasilia since 1998 to study the ""Porto dos Gauchos"" seismic zone (PGSZ). A local seismic refraction survey was carried out with two explosions to help define the seismic velocity model. Both the 1998 and 2005 earthquake sequences occurred in the same WSW-ENE oriented fault zone with right-lateral strike-slip mechanisms. The epicentral zone is in the Parecis basin, near its northern border where there are buried grabens, generally trending WNW-ESE, such as the deep Mesoproterozoic Caiabis graben which lies partly beneath the Parecis basin. However, the epicentral distribution indicates that the 1998 and 2005 sequences are related to a N60 degrees E fault which probably crosses the entire Caiabis graben. The 1955 earthquake, despite the uncertainty in its epicentre, does not seem to be directly related to any buried graben either. The seismicity in the Porto dos Gauchos seismic zone, therefore, is not directly related to rifted crust. The probable direction of the maximum horizontal stress near Porto dos Gauchos is roughly E-W, consistent with other focal mechanisms further south in the Pantanal basin and Paraguay. but seems to be different from the NW-SE direction observed further north in the Amazon basin. The recurrent seismicity observed in Porto dos Gauchos, and the large 1955 earthquake nearby, make this area of the Parecis basin one of the most important seismic zones of Brazil. (C) 2009 Elsevier B.V. All rights reserved.

Basement depths in the Parecis Basin (Amazon) with receiver functions from small local earthquakes in the Porto dos Gauchos seismic zone

Receiver functions from small local earthquakes were used to determine sediment thicknesses in Porto dos Gauchos seismic zone (PGSZ), Parecis basin, Amazonian craton, Brazil. The high velocity contrast between basement and sediments (P-wave velocities of 6.1 and 3.0 km/s, respectively) favors the generation of clear P-to-S converted phases (Ps) seen in the radial component, and also S-to-P conversions (Sp) seen in the vertical component. A reference 10 velocity model determined with shallow refraction experiment in PGSZ helped to convert Ps P time differences to basement depths at 15 stations deployed for aftershocks studies. The results of receiver function integrated with the shallow refraction reveal that the basement depths in the PGSZ increases from the basin border in the north up to about 600 m depth in the south. The basement topography, however, does not vary smoothly and a basement high with a steep topography was detected near the epicentral area. A 400 m elevation difference within 1.7 km distance suggests a possible border fault of a buried graben. This feature seems to be oriented roughly WSW-ENE and could indicate basement structures related to the seismicity of the Porto dos Gauchos Seismic Zone. (C) 2011 Elsevier Ltd. All rights reserved.

Coda wave attenuation in the Parecis Basin, Amazon Craton, Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos GaA(0)chos, Amazon craton-Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q (c)) of coda waves modeled usingwhere Q (0) is the coda quality factor at frequency of 1 Hz and eta is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239-252, 2008), based in the temporal attenuation coefficient, chi(f) instead of Q(f), given by the equation for the calculation of the geometrical attenuation (gamma) and effective attenuation Q (c) values have been computed at central frequencies (and band) of 1.5 (1-2), 3.0 (2-4), 6.0 (4-8), 9.0 (6-12), 12 (8-16), and 18 (12-24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin for the surrounding shield and for the whole region of Porto dos GaA(0)chos Using the independent frequency model, we found: for the cratonic zone, gamma = 0.014 s (-aEuro parts per thousand 1), nu a parts per thousand 1.12; for the basin zone with sediments of similar to 500 m, gamma = 0.031 s (-aEuro parts per thousand 1), nu a parts per thousand 1.27; and for the Parecis basin with sediments of similar to 1,000 m, gamma = 0.047 s (-aEuro parts per thousand 1), nu a parts per thousand 1.42. Analysis of the attenuation factor (Q (c)) for different values of the geometrical spreading parameter (nu) indicated that an increase of nu generally causes an increase in Q (c), both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, (in the deepest part of the basin), than in the basement, (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.

Seismic Stability Analysis of a Himalayan Rock Slope

The seismic slope stability analysis of the right abutment of a railway bridge proposed at about 350 m above the ground level, crossing a river and connecting two huge hillocks in the Himalayas, India, is presented in this paper. The rock slopes are composed of highly jointed rock mass and the joint spacing and orientation are varying at different locations. Seismic slope stability analysis of the slope under consideration is carried out using both pseudo-static approach and time response approach as the site is located in seismic zone V as per the earth quake zonation maps of India. Stability of the slope is studied numerically using program FLAC. The results obtained from the pseudo-static analysis are presented in the form of Factor of Safety (FOS) and the results obtained from the time response analysis of the slope are presented in terms of horizontal and vertical displacements along the slope. The results obtained from both the analyses confirmed the global stability of the slope as the FOS in case of pseudo-static analysis is above 1.0 and the displacements observed in case of time response analysis are within the permissible limits. This paper also presents the results obtained from the parametric analysis performed in the case of time response analysis in order to understand the effect of individual parameters on the overall stability of the slope.

Pseudo-Static and Pseudo-Dynamic Stability Analysis of Tailings Dam Under Seismic Conditions

In this paper the seismic slope stability analyses are performed for a typical section of 44 m high water retention type tailings earthen dam located in the eastern part of India, using both the conventional pseudo-static and recent pseudo-dynamic methods. The tailings earthen dam is analyzed for different upstream conditions of reservoir like filled up with compacted and non-compacted dumped waste materials with different water levels of the pond tailings portion. Phreatic surface is generated using seepage analysis in geotechnical software SEEP/W and that same is used in the pseudo-static and pseudo-dynamic analyses to make the approach more realistic. The minimum values of factor of safety using pseudo-static and pseudo-dynamic method are obtained as 1.18 and 1.09 respectively for the chosen seismic zone in India. These values of factor of safety show clearly the demerits of conventional pseudo-static analysis compared to recent pseudo-dynamic analysis, where in addition to the seismic accelerations, duration, frequency of earthquake, body waves traveling during earthquake and amplification effects are considered.

Typological/Placement possibilities and limitations of ductile behavior of bridges in seismic areas

Consideraciones sobre la ductilidad en zonas sísmicas. This paper analyses the ductile behavior of a highway overpass located in a seismic zone. The paper presents the results of a pushover analysis that enables the design engineer to estimate the behavior of the bridge’s columns in two directions in an independent manner. The differences with the theoretical bilinear behavior are described and explained. Indications are given on the need and possibilities of taking advantage of ductility in different seismic events scenarios.

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test `N' Values for the City of Agartala, Tripura State, India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M (w) > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V (s)) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V (s) profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V (s)30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V (s)30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V (s) and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as `V (s)' is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile-quantile (Q-Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q-Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V (s) profiles of the study area for site response studies.

Probabilistic Models for Forecasting Earthquakes in the Northeast Region of India

Northeast India and its adjoining areas are characterized by very high seismic activity. According to the Indian seismic code, the region falls under seismic zone V, which represents the highest seismic-hazard level in the country. This region has experienced a number of great earthquakes, such as the Assam (1950) and Shillong (1897) earthquakes, that caused huge devastation in the entire northeast and adjacent areas by flooding, landslides, liquefaction, and damage to roads and buildings. In this study, an attempt has been made to find the probability of occurrence of a major earthquake (M-w > 6) in this region using an updated earthquake catalog collected from different sources. Thereafter, dividing the catalog into six different seismic regions based on different tectonic features and seismogenic factors, the probability of occurrences was estimated using three models: the lognormal, Weibull, and gamma distributions. We calculated the logarithmic probability of the likelihood function (ln L) for all six regions and the entire northeast for all three stochastic models. A higher value of ln L suggests a better model, and a lower value shows a worse model. The results show different model suits for different seismic zones, but the majority follows lognormal, which is better for forecasting magnitude size. According to the results, Weibull shows the highest conditional probabilities among the three models for small as well as large elapsed time T and time intervals t, whereas the lognormal model shows the lowest and the gamma model shows intermediate probabilities. Only for elapsed time T = 0, the lognormal model shows the highest conditional probabilities among the three models at a smaller time interval (t = 3-15 yrs). The opposite result is observed at larger time intervals (t = 15-25 yrs), which show the highest probabilities for the Weibull model. However, based on this study, the IndoBurma Range and Eastern Himalaya show a high probability of occurrence in the 5 yr period 2012-2017 with >90% probability.

Attenuation tomography. Modelling regional love waves: Imperial Valley to Pasadena

Abstract to Part I

The inverse problem of seismic wave attenuation is solved by an iterative back-projection method. The seismic wave quality factor, Q, can be estimated approximately by inverting the S-to-P amplitude ratios. Effects of various uncertain ties in the method are tested and the attenuation tomography is shown to be useful in solving for the spatial variations in attenuation structure and in estimating the effective seismic quality factor of attenuating anomalies.

Back-projection attenuation tomography is applied to two cases in southern California: Imperial Valley and the Coso-Indian Wells region. In the Coso-Indian Wells region, a highly attenuating body (S-wave quality factor (Q_β ≈ 30) coincides with a slow P-wave anomaly mapped by Walck and Clayton (1987). This coincidence suggests the presence of a magmatic or hydrothermal body 3 to 5 km deep in the Indian Wells region. In the Imperial Valley, slow P-wave travel-time anomalies and highly attenuating S-wave anomalies were found in the Brawley seismic zone at a depth of 8 to 12 km. The effective S-wave quality factor is very low (Q_β ≈ 20) and the P-wave velocity is 10% slower than the surrounding areas. These results suggest either magmatic or hydrothermal intrusions, or fractures at depth, possibly related to active shear in the Brawley seismic zone.

No-block inversion is a generalized tomographic method utilizing the continuous form of an inverse problem. The inverse problem of attenuation can be posed in a continuous form , and the no-block inversion technique is applied to the same data set used in the back-projection tomography. A relatively small data set with little redundancy enables us to apply both techniques to a similar degree of resolution. The results obtained by the two methods are very similar. By applying the two methods to the same data set, formal errors and resolution can be directly computed for the final model, and the objectivity of the final result can be enhanced.

Both methods of attenuation tomography are applied to a data set of local earthquakes in Kilauea, Hawaii, to solve for the attenuation structure under Kilauea and the East Rift Zone. The shallow Kilauea magma chamber, East Rift Zone and the Mauna Loa magma chamber are delineated as attenuating anomalies. Detailed inversion reveals shallow secondary magma reservoirs at Mauna Ulu and Puu Oo, the present sites of volcanic eruptions. The Hilina Fault zone is highly attenuating, dominating the attenuating anomalies at shallow depths. The magma conduit system along the summit and the East Rift Zone of Kilauea shows up as a continuous supply channel extending down to a depth of approximately 6 km. The Southwest Rift Zone, on the other hand, is not delineated by attenuating anomalies, except at a depth of 8-12 km, where an attenuating anomaly is imaged west of Puu Kou. The Ylauna Loa chamber is seated at a deeper level (about 6-10 km) than the Kilauea magma chamber. Resolution in the Mauna Loa area is not as good as in the Kilauea area, and there is a trade-off between the depth extent of the magma chamber imaged under Mauna Loa and the error that is due to poor ray coverage. Kilauea magma chamber, on the other hand, is well resolved, according to a resolution test done at the location of the magma chamber.

Abstract to Part II

Long period seismograms recorded at Pasadena of earthquakes occurring along a profile to Imperial Valley are studied in terms of source phenomena (e.g., source mechanisms and depths) versus path effects. Some of the events have known source parameters, determined by teleseismic or near-field studies, and are used as master events in a forward modeling exercise to derive the Green's functions (SH displacements at Pasadena that are due to a pure strike-slip or dip-slip mechanism) that describe the propagation effects along the profile. Both timing and waveforms of records are matched by synthetics calculated from 2-dimensional velocity models. The best 2-dimensional section begins at Imperial Valley with a thin crust containing the basin structure and thickens towards Pasadena. The detailed nature of the transition zone at the base of the crust controls the early arriving shorter periods (strong motions), while the edge of the basin controls the scattered longer period surface waves. From the waveform characteristics alone, shallow events in the basin are easily distinguished from deep events, and the amount of strike-slip versus dip-slip motion is also easily determined. Those events rupturing the sediments, such as the 1979 Imperial Valley earthquake, can be recognized easily by a late-arriving scattered Love wave that has been delayed by the very slow path across the shallow valley structure.

Observation and interpretation of tectonic strain release mechanisms

In four chapters various aspects of earthquake source are studied.

Chapter I

Surface displacements that followed the Parkfield, 1966, earthquakes were measured for two years with six small-scale geodetic networks straddling the fault trace. The logarithmic rate and the periodic nature of the creep displacement recorded on a strain meter made it possible to predict creep episodes on the San Andreas fault. Some individual earthquakes were related directly to surface displacement, while in general, slow creep and aftershock activity were found to occur independently. The Parkfield earthquake is interpreted as a buried dislocation.

Chapter II

The source parameters of earthquakes between magnitude 1 and 6 were studied using field observations, fault plane solutions, and surface wave and S-wave spectral analysis. The seismic moment, M_O, was found to be related to local magnitude, M_L, by log M_O = 1.7 M_L + 15.1. The source length vs magnitude relation for the San Andreas system found to be: M_L = 1.9 log L - 6.7. The surface wave envelope parameter AR gives the moment according to log M_O = log AR₃₀₀ + 30.1, and the stress drop, τ, was found to be related to the magnitude by τ = 0.54 M - 2.58. The relation between surface wave magnitude M_S and M_L is proposed to be M_S = 1.7 M_L - 4.1. It is proposed to estimate the relative stress level (and possibly the strength) of a source-region by the amplitude ratio of high-frequency to low-frequency waves. An apparent stress map for Southern California is presented.

Chapter III

Seismic triggering and seismic shaking are proposed as two closely related mechanisms of strain release which explain observations of the character of the P wave generated by the Alaskan earthquake of 1964, and distant fault slippage observed after the Borrego Mountain, California earthquake of 1968. The Alaska, 1964, earthquake is shown to be adequately described as a series of individual rupture events. The first of these events had a body wave magnitude of 6.6 and is considered to have initiated or triggered the whole sequence. The propagation velocity of the disturbance is estimated to be 3.5 km/sec. On the basis of circumstantial evidence it is proposed that the Borrego Mountain, 1968, earthquake caused release of tectonic strain along three active faults at distances of 45 to 75 km from the epicenter. It is suggested that this mechanism of strain release is best described as "seismic shaking."

Chapter IV

The changes of apparent stress with depth are studied in the South American deep seismic zone. For shallow earthquakes the apparent stress is 20 bars on the average, the same as for earthquakes in the Aleutians and on Oceanic Ridges. At depths between 50 and 150 km the apparent stresses are relatively high, approximately 380 bars, and around 600 km depth they are again near 20 bars. The seismic efficiency is estimated to be 0.1. This suggests that the true stress is obtained by multiplying the apparent stress by ten. The variation of apparent stress with depth is explained in terms of the hypothesis of ocean floor consumption.

Projeto de fundações e estrutura de ampliação de escola no Açores utilizando os eurocódigos estruturais

Trabalho de Projeto para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Estruturas

Projeto de execução de um viaduto rodoviário em betão armado pré-esforçado

Trabalho de Projeto para obtenção do grau de Mestre em Engenharia Civil

Design earthquake ground motion prediction for Perth metropolitan area with microtremor measurements for site characterization

Perth is the largest city in Western Australia and home to three-quarters of the state's residents. In recent decades, there have been a lot of earthquake activities just east of Perth in an area known as the South-West Seismic Zone. Previous numerical results of site response analyses based on limited available geology information for PMA indicated that Perth Basin might amplify the bedrock motion by more than 10 times at some frequencies and at some sites. Hence, more detailed studies on site characterization and amplification are necessary. The microtremor method using spatial autocorrelation (SPAC) processing is a useful tool for gaining thickness and shear wave velocity (SWV) of sediments and has been adopted in many previous studies. In this study, the response spectrum of rock site corresponding to the 475-year return period for PMA is defined according to the probabilistic seismic hazard analysis (PSHA) based on the latest ground motion attenuation model of Southwest Western Australia. Site characterization in PMA is performed using two microtremor measurements, namely SPAC technique and H/V method. The clonal selection algorithm (CSA) is introduced to perform direct inversion of SPAC curves to determine the soil profiles of representative PMA sites investigated in this study. Using the simulated bedrock motion as input, the responses of the soil sites are estimated using numerical method based on the shear-wave velocity vs. depth profiles determined from the SPAC technique. The response spectrum of the earthquake ground motion on surface of each site is derived from the numerical results of the site response analysis, and compared with the respective design spectrum defined in the Australian Earthquake Loading Code. The comparison shows that the code spectra are conservative in the short period range, but may slightly underestimate the response spectrum at some long period range.

Arcabouço geofísico, isostasia e causas do magmatismo cenozóico da província Borborema e de sua margem continental (Nordeste do Brasil)

The Borborema Province (BP) is a geologic domain located in Northeastern Brazil. The BP is limited at the south by the São Francisco craton, at the west by the Parnaíba basin, and both at the north and east by coastal sedimentary basins. Nonetheless the BP surface geology is well known, several key aspects of its evolution are still open, notably: i)its tectonic compartmentalization established after the Brasiliano orogenesis, ii) the architecture of its cretaceous continental margin, iii) the elastic properties of its lithosphere, and iv) the causes of magmatism and uplifting which occurred in the Cenozoic. In this thesis, a regional coverage of geophysical data (elevation, gravity, magnetic, geoid height, and surface wave global tomography) were integrated with surface geologic information aiming to attain a better understanding of the above questions. In the Riacho do Pontal belt and in the western sector of the Sergipano belt, the neoproterozoic suture of the collision of the Sul domain of the BP with the Sanfranciscana plate (SFP) is correlated with an expressive dipolar gravity anomaly. The positive lobule of this anomaly is due to the BP lower continental crust uplifting whilst the negative lobule is due to the supracrustal nappes overthrusting the SFP. In the eastern sector of the Sergipano belt, this dipolar gravity anomaly does not exist. However the suture still can be identified at the southern sector of the Marancó complex arc, alongside of the Porto da Folha shear zone, where the SFP N-S geophysical alignments are truncated. The boundary associated to the collision of the Ceará domain of the BP with the West African craton is also correlated with a dipolar gravity anomaly. The positive lobule of this anomaly coincides with the Sobral-Pedro II shear zone whilst the negative lobule is associated with the Santa Quitéria magmatic arc. Judging by their geophysical signatures, the major BP internal boundaries are: i)the western sector of the Pernambuco shear zone and the eastern continuation of this shear zone as the Congo shear zone, ii) the Patos shear zone, and iii) the Jaguaribe shear zone and its southwestern continuation as the Tatajuba shear zone. These boundaries divide the BP in five tectonic domains in the geophysical criteria: Sul, Transversal, Rio Grande do Norte, Ceará, and Médio Coreaú. The Sul domain is characterized by geophysical signatures associated with the BP and SFP collision. The fact that Congo shear zone is now proposed as part of the Transversal domain boundary implies an important change in the original definition of this domain. The Rio Grande do Norte domain presents a highly magnetized crust resulted from the superposition of precambrian and phanerozoic events. The Ceará domain is divided by the Senador Pompeu shear zone in two subdomains: the eastern one corresponds to the Orós-Jaguaribe belt and the western one to the Ceará-Central subdomain. The latter subdomain exhibits a positive ENE-W SW gravity anomaly which was associated to a crustal discontinuity. This discontinuity would have acted as a rampart against to the N-S Brasiliano orogenic nappes. The Médio Coreaú domain also presents a dipolar gravity anomaly. Its positive lobule is due to granulitic rocks whereas the negative one is caused by supracrustal rocks. The boundary between Médio Coreaú and Ceará domains can be traced below the Parnaíba basin sediments by its geophysical signature. The joint analysis of free air anomalies, free air admittances, and effective elastic thickness estimates (Te) revealed that the Brazilian East and Equatorial continental margins have quite different elastic properties. In the first one 10 km < Te < 20 km whereas in the second one Te ≤ 10 km. The weakness of the Equatorial margin lithosphere was caused by the cenozoic magmatism. The BP continental margin presents segmentations; some of them have inheritance from precambrian structures and domains. The segmentations conform markedly with some sedimentary basin features which are below described from south to north. The limit between Sergipe and Alagoas subbasins coincides with the suture between BP and SFP. Te estimates indicates concordantly that in Sergipe subbasin Te is around 20 km while Alagoas subbasin has Te around 10 km, thus revealing that the lithosphere in the Sergipe subbasin has a greater rigidity than the lithosphere in the Alagoas subbasin. Additionally inside the crust beneath Sergipe subbasin occurs a very dense body (underplating or crustal heritage?) which is not present in the crust beneath Alagoas subbasin. The continental margin of the Pernambuco basin (15 < Te < 25 km) presents a very distinct free air edge effect displaying two anomalies. This fact indicates the existence in the Pernambuco plateau of a relatively thick crust. In the Paraíba basin the free air edge effect is quite uniform, Te ≈ 15 km, and the lower crust is abnormally dense probably due to its alteration by a magmatic underplating in the Cenozoic. The Potiguar basin segmentation in three parts was corroborated by the Te estimates: in the Potiguar rift Te ≅ 5 km, in the Aracati platform Te ≅ 25 km, and in the Touros platform Te ≅ 10 km. The observed weakness of the lithosphere in the Potiguar rift segment is due to the high heat flux while the relatively high strength of the lithosphere in the Touros platform may be due to the existence of an archaean crust. The Ceará basin, in the region of Mundaú and Icaraí subbasins, presents a quite uniform free air edge effect and Te ranges from 10 to 15 km. The analysis of the Bouguer admittance revealed that isostasy in BP can be explained with an isostatic model where combined surface and buried loadings are present. The estimated ratio of the buried loading relative to the surface loading is equal to 15. In addition, the lower crust in BP is abnormally dense. These affirmations are particularly adequate to the northern portion of BP where adherence of the observed data to the isostatic model is quite good. Using the same above described isostatic model to calculate the coherence function, it was obtained that a single Te estimate for the entire BP must be lower than 60 km; in addition, the BP north portion has Te around 20 km. Using the conventional elastic flexural model to isostasy, an inversion of crust thickness was performed. It was identified two regions in BP where the crust is thickened: one below the Borborema plateau (associated to an uplifting in the Cenozoic) and the other one in the Ceará domain beneath the Santa Quitéria magmatic arc (a residue associated to the Brasiliano orogenesis). On the other hand, along the Cariri-Potiguar trend, the crust is thinned due to an aborted rifting in the Cretaceous. Based on the interpretation of free air anomalies, it was inferred the existence of a large magmatism in the oceanic crust surrounding the BP, in contrast with the incipient magmatism in the continent as shown by surface geology. In BP a quite important positive geoid anomaly exists. This anomaly is spatially correlated with the Borborema plateau and the Macaú-Queimadas volcanic lineament. The integrated interpretation of geoid height anomaly data, global shear velocity model, and geologic data allow to propose that and Edge Driven Convection (EDC) may have caused the Cenozoic magmatism. The EDC is an instability that presumably occurs at the boundary between thick stable lithosphere and oceanic thin lithosphere. In the BP lithosphere, the EDC mechanism would have dragged the cold lithospheric mantle into the hot asthenospheric mantle thus causing a positive density contrast that would have generated the main component of the geoid height anomaly. In addition, the compatibility of the gravity data with the isostatic model, where combined surface and buried loadings are present, together with the temporal correlation between the Cenozoic magmatism and the Borborema plateau uplifting allow to propose that this uplifting would have been caused by the buoyancy effect of a crustal root generated by a magmatic underplating in the Cenozoic

Estudoda sismicidade na região de Sobral - CE em 2008

The town of Sobral is located at the northwest part of the Ceará State, 250 km away from its capital, Fortaleza. In January 2008, an intense seismic activity began near Sobral with one event with magnitude 4,2mb on May 21. Since the start of its seismic activity, all events were recorded by the SBBR station (located on EMBRAPA Caprinos Farm), which operates in the region since August 2007. After this event, monitoring the seismic activity was carried out with the deployment of a local three component digital seismographic network, from June 5, 2008 until September 24. Initially, this network was composed of six seismographic stations. Later additional five stations were deployed until August 26 2008. This local network detected approximately 2,800 earthquakes. In this study we analyzed 581 earthquakes recorded by at least three stations for hypocentral and focal mechanism determination, and to contribute to a better explanation of the seismicity which in this region. To determine the hypocenters, we used a half-space model, with vP = 6,00 km/s and vP/vS = 1,71. From the hypocentral determination, it was revealed an active seismic zone with depth ranging between 1 and 8 km, 6 km long in E - W direction. The determination of fault planes and focal mechanism was obtained using the programs FPFIT and PLAN, which allowed comparison between their respective results in order to obtain more accurate results. A set of 24 earthquakes were selected to determine fault using PLAN planes and focal mechanisms using FPFIT. With the aid of detailed map of hypocenters this set, it was possible to identify three structures. Therefore, the set of 24 earthquakes were divided into three subsets. The type of mechanism was predominantly strike-slip with a dextral direction. Although the region has two tectonic structures near the site of the study area: the Café- Ipueiras Fault (normal fault) and the Sobral-Pedro II Lineament (dextral strike-slip fault) it was not possible to correlate the seismicity founded with those structures