Assimilation of Geosat Altimeter Data for the Agulhas Current Using the Ensemble Kalman Filter with a Quasigeostrophic Model

The ring-shedding process in the Agulhas Current is studied using the ensemble Kalman filter to assimilate geosat altimeter data into a two-layer quasigeostrophic ocean model. The properties of the ensemble Kalman filter are further explored with focus on the analysis scheme and the use of gridded data. The Geosat data consist of 10 fields of gridded sea-surface height anomalies separated 10 days apart that are added to a climatic mean field. This corresponds to a huge number of data values, and a data reduction scheme must be applied to increase the efficiency of the analysis procedure. Further, it is illustrated how one can resolve the rank problem occurring when a too large dataset or a small ensemble is used.

Application of data assimilation to ocean and climate prediction

Ocean prediction systems are now able to analyse and predict temperature, salinity and velocity structures within the ocean by assimilating measurements of the ocean’s temperature and salinity into physically based ocean models. Data assimilation combines current estimates of state variables, such as temperature and salinity, from a computational model with measurements of the ocean and atmosphere in order to improve forecasts and reduce uncertainty in the forecast accuracy. Data assimilation generally works well with ocean models away from the equator but has been found to induce vigorous and unrealistic overturning circulations near the equator. A pressure correction method was developed at the University of Reading and the Met Office to control these circulations using ideas from control theory and an understanding of equatorial dynamics. The method has been used for the last 10 years in seasonal forecasting and ocean prediction systems at the Met Office and European Center for Medium-range Weather Forecasting (ECMWF). It has been an important element in recent re-analyses of the ocean heat uptake that mitigates climate change.

A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

A parallel data assimilation model for oceanographic observations

In this paper we describe the development of a program that aims at the optimal integration of observed data in an oceanographic model describ

Magnetotelluric deep soundings, gravity and geoid in the south Sao Francisco craton: Geophysical indicators of cratonic lithosphere rejuvenation and crustal underplating

In the south Sao Francisco craton a circular and 8-m amplitude geoid anomaly coincides with the outcropping terrain of an Archean-Paleoproterozoic basement. Broadband magnetotelluric (MT) data inversions of two radial profiles within the positive geoid and Bouguer gravity anomaly yield geo-electrical crustal sections, whereby the lower crust is locally more conductive (10 to 100 Omega m) in spatial coincidence with a denser lower crust modeled by the gravity data. This anomalous lower crust may have resulted from magmatic underplating, associated with Mesoarchean and Proterozoic episodes of tholeiitic dike intrusion. Long-period MT soundings reveal a low electrical resistivity mantle (20 to 200 Omega m) from depths beyond 120 km. Forward geoid modeling, using the scope of the low electrical resistivity region within the mantle as a constraint, entails a density increase (40 to 50 kg/m(3)) possibly due to Fe enrichment of mantle minerals. However, this factor alone does not explain the observed resistivity. A supplemented presence of small amounts of percolated carbonatite melting (similar to 0.005 vol.%), dissolved water and enhanced oxygen fugacity within the peridotitic mantle are viable agents that could explain the less resistive upper mantle. We propose that metasomatic processes confined in the sub-continental lithospheric mantle foster the conditions for a low degree melting with variable CO(2), H(2)O and Fe content. Even though the precise age of this metasomatism is unknown it might be older than the Early Cretaceous based on the evidence that a high-degree of melting in a lithospheric mantle impregnated with carbonatites originated the tholeiitic dike intrusions dispersed from the southeastern border of the Sao Francisco craton, during the onset of the lithosphere extension and break-up of the western Gondwana. The proxies are the NE Parana and Espinhaco (130 Ma, Ar/Ar ages) tholeiitic dikes, which contain (similar to 3%) carbonatites in their composition. The occurrence of a positive geoid anomaly (+ 10 m) and pre-tholeiites (age > 138 Ma), carbonatites and kimberlites along the west African continental margin (Angola and Namibia) reinforces the presumed age of the Sao Francisco-Congo craton rejuvenation to be prior to its fragmentation in the Lower Cretaceous. (C) 2010 Elsevier B.V. All rights reserved.

Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

Genetic algorithm inversion of the average 1D crustal structure using local and regional earthquakes

Knowing the best 1D model of the crustal and upper mantle structure is useful not only for routine hypocenter determination, but also for linearized joint inversions of hypocenters and 3D crustal structure, where a good choice of the initial model can be very important. Here, we tested the combination of a simple GA inversion with the widely used HYPO71 program to find the best three-layer model (upper crust, lower crust, and upper mantle) by minimizing the overall P- and S-arrival residuals, using local and regional earthquakes in two areas of the Brazilian shield. Results from the Tocantins Province (Central Brazil) and the southern border of the Sao Francisco craton (SE Brazil) indicated an average crustal thickness of 38 and 43 km, respectively, consistent with previous estimates from receiver functions and seismic refraction lines. The GA + HYPO71 inversion produced correct Vp/Vs ratios (1.73 and 1.71, respectively), as expected from Wadati diagrams. Tests with synthetic data showed that the method is robust for the crustal thickness, Pn velocity, and Vp/Vs ratio when using events with distance up to about 400 km, despite the small number of events available (7 and 22, respectively). The velocities of the upper and lower crusts, however, are less well constrained. Interestingly, in the Tocantins Province, the GA + HYPO71 inversion showed a secondary solution (local minimum) for the average crustal thickness, besides the global minimum solution, which was caused by the existence of two distinct domains in the Central Brazil with very different crustal thicknesses. (C) 2010 Elsevier Ltd. All rights reserved.

High T/P evolution and metamorphic ages of the migmatitic basement of northern Sierras Pampeanas, Argentina: Characterization of a mid-crustal segment of the Famatinian belt

New petrologic, thermobarometric and U-Pb monazite geochronologic information allowed to resolve the metamorphic evolution of a high temperature mid-crustal segment of an ancient subduction-related orogen. The EI Portezuelo Metamorphic-Igneous Complex, in the northern Sierras Pampeanas, is mainly composed of migmatites that evolved from amphibolite to granulite metamorphic facies, reaching thermal peak conditions of 670-820 degrees C and 4.5-5.3 kbar. The petrographic study combined with conventional and pseudosection thermobarometry led to deducing a short prograde metamorphic evolution within migmatite blocks. The garnet-absent migmatites represent amphibolite-facies rocks, whereas the cordierite-garnet-K-feldspar-sillimanite migmatites represent higher metamorphic grade rocks. U-Pb geochronology on monazite grains within leucosome record the time of migmatization between approximate to 477 and 470 Ma. Thus, the El Portezuelo Metamorphic-Igneous Complex is an example of exhumed Early Ordovician anatectic middle crust of the Famatinian mobile belt. Homogeneous exposure of similar paleo-depths throughout the Famatinian back-arc and isobaric cooling paths suggest slow exhumation and consequent longstanding crustal residence at high temperatures. High thermal gradients uniformly distributed in the Famatinian back-arc can be explained by shallow convection of a low-viscosity asthenosphere promoted by subducting-slab dehydration. (C) 2011 Elsevier Ltd. All rights reserved.

Cordierite-bearing lavas from Jaguarao, southern Brazil: Petrological evidence for crustal melts during early rifting of Gondwana

The Jaguarao stratoid dacites (Rio Grande do Sul, Brazil) are limited in areal extent, are comprised of about 3.2 km(3) of preserved erupted material, and outcrop only in areas of the region underlain by mylonitic and ultramylonitic rocks. They are S-type volcanic rocks containing cordierite, orthopyroxene, plagioclase, and ilmenite as liquidus phases, and partially melted granite, gneiss, and migmatite enclaves that are very similar to the Precambrian basement rocks. The Jaguarao lavas have distinct geochemical signatures and Sr-Nd isotopes with respect to other volcanic rocks of the region. Available geochronological data for Jaguarao dacites range between 157 +/- 5 Ma and 139.6 +/- 7.4 Ma. Considering the errors, the younger ages obtained for Jaguarao lavas overlap the 138-128 Ma age of rocks of the Serra Geral Group, and thus indicate that the dacites were erupted prior to the break-up of Gondwana in this region. Petrographic, mineralogical, and petrochemical data, as well as the tectonic context of the Jaguarao lavas, suggest that magma genesis was linked, at least in part, to friction melts. The dacitic magma was generated by partial melting reactions involving biotite breakdown in a dominantly quartz-feldspathic source terrane, leaving a granulite facies residue in subsurface. These melts were probably generated as a consequence of crustal thinning linked to simple shear extension just prior to Gondwana break-up and rifting of the southern Atlantic Ocean. (C) 2009 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

A new Home, far from Home : The assimilation process of women involved in intercultural marriages based on Internet meeting

The internet has revolutionized the way we socialize, and as a consequence the way to love. The new communication technologies have facilitated intercultural relationships. Nowadays family relations are one of the major factors in immigration to European countries. Family relations means persons who arrive as family dependents and in accordance with laws regulating family reunification. This thesis aims to apply the classical assimilation theory stated by Milton Gordon (1964), which formulates a series of assimilation stages through which an individual must pass in order to be completely assimilated. In accordance with this theory, marriage is the final phase for a newcomer to fully incorporate into the host society. Thus, based on this presumption and other contemporary theories, the present study has analysed how women who get involved in intercultural marriages based on internet meeting experience these assimilation stages and evaluated the resources used by respondents to incorporate themselves into Swedish society.The main goal of the study was to determine if jumping to the last stage of assimilation does assure the incorporation in the social or/and labour spheres and the findings demonstrate that even though husbands are a valuable resource for assimilation, several cultural issues in Swedish society make it difficult to assure success for the newcomers.On the other hand, Sweden is a country with a strong national sentiment and the assimilation of immigrants still is an important issue to deal with. The Swedish Integration Board has disappeared and major projects for integration have been left in the hands of the municipalities or the Migration Board, institutions that still do not know how to deal with this dilemma.

Effect Of Different Hydrological Model Structures On The Assimilation Of Distributed Uncertain Observations

The reliable evaluation of the flood forecasting is a crucial problem for assessing flood risk and consequent damages. Different hydrological models (distributed, semi-distributed or lumped) have been proposed in order to deal with this issue. The choice of the proper model structure has been investigated by many authors and it is one of the main sources of uncertainty for a correct evaluation of the outflow hydrograph. In addition, the recent increasing of data availability makes possible to update hydrological models as response of real-time observations. For these reasons, the aim of this work it is to evaluate the effect of different structure of a semi-distributed hydrological model in the assimilation of distributed uncertain discharge observations. The study was applied to the Bacchiglione catchment, located in Italy. The first methodological step was to divide the basin in different sub-basins according to topographic characteristics. Secondly, two different structures of the semi-distributed hydrological model were implemented in order to estimate the outflow hydrograph. Then, synthetic observations of uncertain value of discharge were generated, as a function of the observed and simulated value of flow at the basin outlet, and assimilated in the semi-distributed models using a Kalman Filter. Finally, different spatial patterns of sensors location were assumed to update the model state as response of the uncertain discharge observations. The results of this work pointed out that, overall, the assimilation of uncertain observations can improve the hydrologic model performance. In particular, it was found that the model structure is an important factor, of difficult characterization, since can induce different forecasts in terms of outflow discharge. This study is partly supported by the FP7 EU Project WeSenseIt.

Parameterization Of Turbulence Models Using 3DVAR Data Assimilation

In this research the 3DVAR data assimilation scheme is implemented in the numerical model DIVAST in order to optimize the performance of the numerical model by selecting an appropriate turbulence scheme and tuning its parameters. Two turbulence closure schemes: the Prandtl mixing length model and the two-equation k-ε model were incorporated into DIVAST and examined with respect to their universality of application, complexity of solutions, computational efficiency and numerical stability. A square harbour with one symmetrical entrance subject to tide-induced flows was selected to investigate the structure of turbulent flows. The experimental part of the research was conducted in a tidal basin. A significant advantage of such laboratory experiment is a fully controlled environment where domain setup and forcing are user-defined. The research shows that the Prandtl mixing length model and the two-equation k-ε model, with default parameterization predefined according to literature recommendations, overestimate eddy viscosity which in turn results in a significant underestimation of velocity magnitudes in the harbour. The data assimilation of the model-predicted velocity and laboratory observations significantly improves model predictions for both turbulence models by adjusting modelled flows in the harbour to match de-errored observations. 3DVAR allows also to identify and quantify shortcomings of the numerical model. Such comprehensive analysis gives an optimal solution based on which numerical model parameters can be estimated. The process of turbulence model optimization by reparameterization and tuning towards optimal state led to new constants that may be potentially applied to complex turbulent flows, such as rapidly developing flows or recirculating flows.

Anisotropia sísmica crustal na região de Cascavel - CE

This dissertation presents a study on crustal seismic anisotropy in Cascavel - CE. The earthquake data employed here are from the Seismological Laboratory at Universidade Federal do Rio Grande do Norte (UFRN) and were colected from 29 September 1997 to 05 march 1998 using six three-component digital seismographic stations. In general, the cause of the observed seismic anisotropy in many regions of the world is interpreted in terms of fluid-filled stress aligned microcracks in the rockmass (EDA). In other words, the polarisation direction of the faster shear-wave splitting is parallel to SHmax. However, other researches on seismic anisotropy carried out in NE Brazil have shown a remarkable consistency of the faster shear-wave polarisation direction with the direction of the Precambrian fabric. The present work is another case study that is used to investigate this issue. In order to map the Precambrian fabric we used aeromagnetic data, since the study area is mostly covered with sediments (up to 50m thick) and in-situ field mapping would be very difficult to be carried out. According to the results from the present research, the observations of the faster shear-wave polarisation directions in two seismographic stations in Cascavel region are best explained in the framework of EDA. For the remaining two stations, the observed anisotropy may have two interpretions: (i) - 90_ flips of the direction of polarisation of the faster shear-wave, since that the event-to-station ray path would be through the fracture zone and hence would travel under a higher pore pressure and (ii) - the observed seismic anisotropy would agree with the direction on the ductile Precambrian fabric