Integração de métodos geofísicos na modelagem crustal da Bacia de Santos

A Bacia de Santos é uma bacia de margem passiva gerada durante o Neocomiano, associado com o evento de ruptura do Supercontinente de Gondwana. A espessura da seção sedimentar pode alcançar 15 km, incluindo a seqüência rifte, e foi depositada sobre uma crosta afinada cuja espessura original era de cerca de 35 km. Foi realizada uma modelagem crustal em 8 linhas, aproximadamente perpendiculares ao mergulho deposicional da bacia, objetivando o reconhecimento da geometria da fase rifte e a estimativa do estiramento crustal. Foram utilizados durante a modelagem quatro horizontes sísmicos em profundidade, mapas magnético e gravimétrico além de dados da profundidade do embasamento magnético. Foi reconhecido um pronunciado estiramento crustal em toda a bacia (fatores de estiramento crustal variando entre 1.2 e 3.1), evidenciando ampla acomodação da deformação. Os resultados da modelagem necessitaram, ainda, de espessas e contínuas camadas de rochas vulcânicas para ajuste do modelo. Estas cunhas de rochas vulcânicas, limitadas lateralmente por falhas normais, são feições comuns em bacias rifte. Os resultados da modelagem foram confrontados com dados da profundidade do embasamento magnético, obtidos anteriormente, e evidenciaram forte correlação, sobretudo nas áreas rasas da bacia.

Tensor3D: A computer graphics program to simulate 3D real-time deformation and visualization of geometric bodies

Tensor3D is a geometric modeling program with the capacity to simulate and visualize in real-time the deformation, specified through a tensor matrix and applied to triangulated models representing geological bodies. 3D visualization allows the study of deformational processes that are traditionally conducted in 2D, such as simple and pure shears. Besides geometric objects that are immediately available in the program window, the program can read other models from disk, thus being able to import objects created with different open-source or proprietary programs. A strain ellipsoid and a bounding box are simultaneously shown and instantly deformed with the main object. The principal axes of strain are visualized as well to provide graphical information about the orientation of the tensor's normal components. The deformed models can also be saved, retrieved later and deformed again, in order to study different steps of progressive strain, or to make this data available to other programs. The shape of stress ellipsoids and the corresponding Mohr circles defined by any stress tensor can also be represented. The application was written using the Visualization ToolKit, a powerful scientific visualization library in the public domain. This development choice, allied to the use of the Tcl/Tk programming language, which is independent on the host computational platform, makes the program a useful tool for the study of geometric deformations directly in three dimensions in teaching as well as research activities. (C) 2007 Elsevier Ltd. All rights reserved.

Comportamento fractal de tracos de fraturas na borda sul do macico granitico de itu, estado de Sao Paulo

Fractal geometry is relevant to understand and explain many natural complex geometries. Using the fractal set concept (fig. 1) many authors have shown that shorelines, landscapes and fractures follow a fractal behaviour. These authors have developed many methods, including the Cantor's Dust Method (CDM) (VELDE et al., 1992), a linear method of analysis adapted for the determination of two-dimensional phenomena. The Itu Granitic Complex (IGC) is a wide granitic body that that crops out at northwest of Cabreuva City, Sao Paulo State (fig. 2) and was affected in its south border by dextral Itu-Jundiuvira Shear Zone (IJSZ) that produced fractures and alignment of feldspars crystals. The different types of fractures (compression, distension and shear) was discriminated from the relationship between them and medium stress ellipsoid of IJSZ (fig. 3). A modified version of CDM was used to study a possible fractal behaviour of the fracture traces in the south border of IGC. The main modification was the use only one direction of analysis (NE/SW). Four parallel profiles were traced with lengths between 9.75km and 12.75km, each one them was divided into six classes of segments (x) with 375m, 500m, 750m, 1.000m, 1.250m and 1.500m. The parameter (N) is provided by he rate between profile length and choiced segment. For each x the number of intervals is counted with at least one event (fracture intersection) which supplied the parameter(n). The n/N rate provide the parameter (p) that represents the relationship between frequency of events and x. And finally the parameters p and x were plotted in a logarithmic graphics (fig. 4) that provide a line with such a declivity (1) which is related to effective dimension (De). In theory, granitics bodies are isotropics and they would have a same fractal dimension in all segments, but the logarithmic graphics (fig. 4) show that fracture traces of IGC has a fractal behaviour in a restrict interval. This fact probably occurs from the passage of a ductil-brittle deformation condition to a more brittle deformation condition of IGC.

Influence of visual feedback sampling on obstacle crossing behavior in people with Parkinson's disease

The purpose of the current study was to investigate the role of visual information on gait control in people with Parkinson's disease as they crossed over obstacles. Twelve healthy individuals, and 12 patients with mild to moderate Parkinson's disease, walked at their preferred speeds along a walkway and stepped over obstacles of varying heights (ankle height or half-knee height), under three visual sampling conditions: dynamic (normal lighting), static (static visual samples, similar to stroboscopic lighting), and voluntary visual sampling. Subjects wore liquid crystal glasses for visual manipulation. In the static visual sampling condition only, the patients with Parkinson's disease made contact with the obstacle more often than did the control subjects. In the successful trials, the patients increased their crossing step width in the static visual sampling condition as compared to the dynamic and voluntary visual sampling conditions; the control group maintained the same step width for all visual sampling conditions. The patients showed lower horizontal mean velocity values during obstacle crossing than did the controls. The patients with Parkinson's disease were more dependent on optic flow information for successful task and postural stability than were the control subjects. Bradykinesia influenced obstacle crossing in the patients with Parkinson's disease. © 2013 Elsevier B.V.

Análise de estruturas flexíveis com aplicação de materiais viscoelásticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Evolução tectônica dos altos estruturais de Pitanga, Artemis, Pau D'Alho e Jibóia - Centro do Estado de São Paulo

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo da geometria e cinemática das rochas sedimentares arqueanas da mina do Igarapé do Azul – Carajás-PA

Mina de Manganês do Igarapé Azul posiciona-se geologicamente no interior do feixe da Falha Carajás, na porção central do Sistema Transcorrente de Carajás. O depósito do Manganês do Azul relaciona-se a rochas sedimentares pelíticas do Membro Azul, na base da Formação Águas Claras (Arqueano), em contato discordante, acima do Grupo Grão Pará (Nogueira et al, 1995). Três frentes de lavra a céu aberto estão atualmente em andamento na área: (1) Mina Principal (Mina 1), (2) Mina 2 e (3) Mina 3. Nestes locais encontram-se excelentes afloramentos de siltitos intercalados com argilitos e arenitos finos, intercalados com níveis manganesíferos. Essas rochas estão organizadas em conjuntos de dobras e falhas normais e inversas sob deformação heterogênea, particionada em diferentes escalas. As seções geológicas realizadas nas frentes de lavra mostram a predominância de siltitos intercalados com argilitos em contato com rochas pelíticas manganesíferas e minério (bióxido de Mn). Nessas rochas são comuns estruturas primárias tipo hummocky, estratificações cruzadas, e laminações plano-paralelas. O acamamento centimétrico a métrico (em média de 30 a 50 cm ) representa a principal estrutura primária, usada como marcador de deformação, observada nas rochas. A Mina do Igarapé Azul encontra-se dividida em dois blocos, separados por falha normal com rejeito de até dezenas de metros, com o bloco norte alto em relação ao bloco sul. O bloco sul encontra-se pouco deformado, apresentado uma regularidade no acamamento que mergulha com ângulos suaves para sul, colocando a camada de minério sucessivamente em níveis mais profundos na direção S. No bloco norte o acamamento apresenta um comportamento heterogêneo. A deformação é mais expressiva nessa região, estando o nível de minério deformado por dobras e falhas inversas. Além da cinemática vertical, as falhas apresentam deslocamento conjugado destral dando a essas feições um caráter oblíquo. Essa região pode ser definida como um corredor de deformação. O corredor observado no bloco norte, de acordo com os domínios principais separados pelas falhas anteriormente descritas, possui orientação NW-SE, com aproximadamente um quilômetro de extensão, sendo caracterizado por dobras assimétricas curvilineares com eixos de mergulhos suaves (10° a 25°) para NW e SE. Essas dobras são seccionadas por falhas normais sinuosas NW-SE e/ou E-W, com baixo ângulo de mergulho (em torno de 10° a 30°), subordinadas a transcorrências destrais, gerando em escala de detalhes feições como drag folds. Observam-se ainda falhas inversas retas e/ou sinuosas NW-SE e zonas de falhas sub-verticalizadas WNW -ESE. As dobras individuais nesta área são estruturas do tipo reversas, flexurais e com geometria en echelon com orientação semelhante às dobras curvilineares: eixos com baixo ângulo de mergulho (10° a 25°) e caimento para SE. O conjunto de feições anteriormente descrito desenha, em escala quilométrica, um antiforme aberto, provavelmente resultante da acomodação do acamamento em resposta a deformação dessas falhas. O paralelismo entre feições observadas na área da Mina do Igarapé Azul e os lineamentos maiores que desenham a Falha Carajás em planta sugere uma relação com dois importantes episódios deformacionais ocorridos durante a história tectônica da Falha Carajás. As falhas normais associadas a componente direcional destral, de maior expressividade da área da mina, estariam relacionadas ao episódio de transtensão destral responsável pela instalação da Falha Carajás anterior a 2.6 Ga (Pinheiro, 1997). As dobras, as falhas de cavalgamento e as zonas de falhas sub-verticalizadas estariam relacionados a deformações sob regime de transpressão sinistral, um segundo evento atuante na região, responsável pela reativação e inversão da maioria das estruturas próximas à zona da Falha Carajás (Pinheiro, 1997; Pinheiro & Holdsworth, 2000; Lima, 2002).

Deformação das rochas siliciclásticas paleoproterozoicas do Grupo Araí como exemplo das reativações de falhas do embasamento, Serra do Tepequém, Roraima, norte do Brasil

As rochas siliciclásticas da Serra do Tepequém são correlacionadas aos Grupos Araí e Suapi, pertencentes ao Supergrupo Roraima de idade paleoproterozoica. A análise estrutural destas rochas revelou que o acamamento exibe arranjos com mergulhos preferencialmente para SE e NW, individualizados em domínios limitados por zonas de falhas oblíquas sinistrais com rejeitos normais e inversos, com direção NE-SW. Essa estruturação é formada por dobras forçadas quilométricas do tipo kink bands e chevrons. O novo arcabouço geométrico observado é característico de um ambiente de deformação de nível crustal raso a médio. A história tectônica é controlada por reativações dos planos de fraqueza preexistentes nas tramas antigas, dúcteis, do embasamento. Esse modelo diverge dos regionais prévios para a região, os quais consideram as dobras existentes como produtos de ambiente dúctil sob tectônica colisional. Os resultados evidenciam a importância da presença de estruturas antigas do embasamento, relacionado ao Escudo das Guianas, como controladoras da geometria das rochas da Serra do Tepequém, em ambiente rúptil.

Propriedades genéricas de sistemas hamiltonianos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

An integrable deformation of the AdS(5) x S-5 superstring

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Myocardial Deformation by Speckle Tracking in Severe Dilated Cardiomyopathy

Background: The high and increasing prevalence of Dilated Cardiomyopathy (DCM) represents a serious public health issue. Novel technologies have been used aiming to improve diagnosis and the therapeutic approach. In this context, speckle tracking echocardiography (STE) uses natural myocardial markers to analyze the systolic deformation of the left ventricle (LV). Objective: Measure the longitudinal transmural global strain (GS) of the LV through STE in patients with severe DCM, comparing the results with normal individuals and with echocardiographic parameters established for the analysis of LV systolic function, in order to validate the method in this population. Methods: Seventy-one patients with severe DCM (53 +/- 12 years, 72% men) and 20 controls (30 +/- 8 years, 45% men) were studied. The following variables were studied: LV volumes and ejection fraction calculated by two and three-dimensional echocardiography, Doppler parameters, Tissue Doppler Imaging systolic and diastolic LV velocities and GS obtained by STE. Results: Compared with controls, LV volumes were higher in the DCM group; however, LVEF and peak E-wave velocity were lower in the latter. The myocardial performance index was higher in the patient group. Tissue Doppler myocardial velocities (S', e', a') were significantly lower and E/e' ratio was higher in the DCM group. GS was decreased in the DCM group (-5.5% +/- 2.3%) when compared with controls (-14.0% +/- 1.8%). Conclusion: In this study, GS was significantly lower in patients with severe DCM, bringing new perspectives for therapeutic approaches in this specific population. (Arq Bras Cardiol 2012;99(3):834-842)

Influence of normal and radial contributions of local current density on local electrochemical impedance spectroscopy

A new tri-electrode probe is presented and applied to local electrochemical impedance spectroscopy (LEIS) measurements. As opposed to two-probe systems, the three-probe one allows measurement not only of normal, but also of radial contributions of local current densities to the local impedance values. The results concerning the cases of the blocking electrode and the electrode with faradaic reaction are discussed from the theoretical point of view for a disk electrode. Numerical simulations and experimental results are compared for the case of the ferri/ferrocyanide electrode reaction at the Pt working electrode disk. At the centre of the disk, the impedance taking into account both normal and radial contributions was in good agreement with the local impedance measured in terms of only the normal contribution. At the periphery of the electrode, the impedance taking into account both normal and radial contributions differed significantly from the local impedance measured in terms of only the normal contribution. The radial impedance results at the periphery of the electrode are in good agreement with the usual explanation that the associated larger current density is attributed to the geometry of the electrode, which exhibits a greater accessibility at the electrode edge. (C) 2011 Elsevier Ltd. All rights reserved.

A variational approach to 3D geometry reconstruction from two or multiple views

[EN] In the last years we have developed some methods for 3D reconstruction. First we began with the problem of reconstructing a 3D scene from a stereoscopic pair of images. We developed some methods based on energy functionals which produce dense disparity maps by preserving discontinuities from image boundaries. Then we passed to the problem of reconstructing a 3D scene from multiple views (more than 2). The method for multiple view reconstruction relies on the method for stereoscopic reconstruction. For every pair of consecutive images we estimate a disparity map and then we apply a robust method that searches for good correspondences through the sequence of images. Recently we have proposed several methods for 3D surface regularization. This is a postprocessing step necessary for smoothing the final surface, which could be afected by noise or mismatch correspondences. These regularization methods are interesting because they use the information from the reconstructing process and not only from the 3D surface. We have tackled all these problems from an energy minimization approach. We investigate the associated Euler-Lagrange equation of the energy functional, and we approach the solution of the underlying partial differential equation (PDE) using a gradient descent method.

Estimating crustal deformation from geodetic data

By the end of the 19th century, geodesy has contributed greatly to the knowledge of regional tectonics and fault movement through its ability to measure, at sub-centimetre precision, the relative positions of points on the Earth’s surface. Nowadays the systematic analysis of geodetic measurements in active deformation regions represents therefore one of the most important tool in the study of crustal deformation over different temporal scales [e.g., Dixon, 1991]. This dissertation focuses on motion that can be observed geodetically with classical terrestrial position measurements, particularly triangulation and leveling observations. The work is divided into two sections: an overview of the principal methods for estimating longterm accumulation of elastic strain from terrestrial observations, and an overview of the principal methods for rigorously inverting surface coseismic deformation fields for source geometry with tests on synthetic deformation data sets and applications in two different tectonically active regions of the Italian peninsula. For the long-term accumulation of elastic strain analysis, triangulation data were available from a geodetic network across the Messina Straits area (southern Italy) for the period 1971 – 2004. From resulting angle changes, the shear strain rates as well as the orientation of the principal axes of the strain rate tensor were estimated. The computed average annual shear strain rates for the time period between 1971 and 2004 are γ˙1 = 113.89 ± 54.96 nanostrain/yr and γ˙2 = -23.38 ± 48.71 nanostrain/yr, with the orientation of the most extensional strain (θ) at N140.80° ± 19.55°E. These results suggests that the first-order strain field of the area is dominated by extension in the direction perpendicular to the trend of the Straits, sustaining the hypothesis that the Messina Straits could represents an area of active concentrated deformation. The orientation of θ agree well with GPS deformation estimates, calculated over shorter time interval, and is consistent with previous preliminary GPS estimates [D’Agostino and Selvaggi, 2004; Serpelloni et al., 2005] and is also similar to the direction of the 1908 (MW 7.1) earthquake slip vector [e.g., Boschi et al., 1989; Valensise and Pantosti, 1992; Pino et al., 2000; Amoruso et al., 2002]. Thus, the measured strain rate can be attributed to an active extension across the Messina Straits, corresponding to a relative extension rate ranges between < 1mm/yr and up to ~ 2 mm/yr, within the portion of the Straits covered by the triangulation network. These results are consistent with the hypothesis that the Messina Straits is an important active geological boundary between the Sicilian and the Calabrian domains and support previous preliminary GPS-based estimates of strain rates across the Straits, which show that the active deformation is distributed along a greater area. Finally, the preliminary dislocation modelling has shown that, although the current geodetic measurements do not resolve the geometry of the dislocation models, they solve well the rate of interseismic strain accumulation across the Messina Straits and give useful information about the locking the depth of the shear zone. Geodetic data, triangulation and leveling measurements of the 1976 Friuli (NE Italy) earthquake, were available for the inversion of coseismic source parameters. From observed angle and elevation changes, the source parameters of the seismic sequence were estimated in a join inversion using an algorithm called “simulated annealing”. The computed optimal uniform–slip elastic dislocation model consists of a 30° north-dipping shallow (depth 1.30 ± 0.75 km) fault plane with azimuth of 273° and accommodating reverse dextral slip of about 1.8 m. The hypocentral location and inferred fault plane of the main event are then consistent with the activation of Periadriatic overthrusts or other related thrust faults as the Gemona- Kobarid thrust. Then, the geodetic data set exclude the source solution of Aoudia et al. [2000], Peruzza et al. [2002] and Poli et al. [2002] that considers the Susans-Tricesimo thrust as the May 6 event. The best-fit source model is then more consistent with the solution of Pondrelli et al. [2001], which proposed the activation of other thrusts located more to the North of the Susans-Tricesimo thrust, probably on Periadriatic related thrust faults. The main characteristics of the leveling and triangulation data are then fit by the optimal single fault model, that is, these results are consistent with a first-order rupture process characterized by a progressive rupture of a single fault system. A single uniform-slip fault model seems to not reproduce some minor complexities of the observations, and some residual signals that are not modelled by the optimal single-fault plane solution, were observed. In fact, the single fault plane model does not reproduce some minor features of the leveling deformation field along the route 36 south of the main uplift peak, that is, a second fault seems to be necessary to reproduce these residual signals. By assuming movements along some mapped thrust located southward of the inferred optimal single-plane solution, the residual signal has been successfully modelled. In summary, the inversion results presented in this Thesis, are consistent with the activation of some Periadriatic related thrust for the main events of the sequence, and with a minor importance of the southward thrust systems of the middle Tagliamento plain.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.