Simulation of moisture profiles in a Latossol in Dourados region, in the State of Mato Grosso do Sul, Brazil

The understanding of unsaturated soil water flow at process-level is essential to develop proper management actions for environmental protection in agricultural systems. One important tool for simulation of soil water flow that has been used worldwide is the SWAP model. The aim of this work was to test and to calibrate the SWAP model by inverse modeling to describe moisture profiles in a Brazilian very clayey Latossol in Dourados, State of Mato Grosso do Sul, Brazil. The SWAP model was tested in an experimental field of 0.09 ha cultivated with soybean and soil profiles were sampled eight times between December 2006 and October 2007. The SWAP input values (i.e. soil water retention curves and meteorological data) were based on in-situ measurements. Simulations with uncalibrated soil water retention curves resulted in moisture profiles that were too wet for almost all sampling dates, in particular between 0-10 cm depth. After calibration of soil water retention curves, there was a good improvement in the simulated moisture profiles, which were within the range of measured values for almost all depths and sampling dates.

Sensibility analysis of the pearl model for pesticide leaching in the State of Mato Grosso do Sul, Brazil

For an accurate use of pesticide leaching models it is necessary to assess the sensitivity of input parameters. The aim of this work was to carry out sensitivity analysis of the pesticide leaching model PEARL for contrasting soil types of Dourados river watershed in the state of Mato Grosso do Sul, Brazil. Sensitivity analysis was done by carrying out many simulations with different input parameters and calculating their influence on the output values. The approach used was called one-at-a-time sensitivity analysis, which consists in varying independently input parameters one at a time and keeping all others constant with the standard scenario. Sensitivity analysis was automated using SESAN tool that was linked to the PEARL model. Results have shown that only soil characteristics influenced the simulated water flux resulting in none variation of this variable for scenarios with different pesticides and same soil. All input parameters that showed the greatest sensitivity with regard to leached pesticide are related to soil and pesticide properties. Sensitivity of all input parameters was scenario dependent, confirming the need of using more than one standard scenario for sensitivity analysis of pesticide leaching models.

On the observability of turbulent transport rates by Argo: supporting evidence from an inversion experiment

Although estimation of turbulent transport parameters using inverse methods is not new, there is little evaluation of the method in the literature. Here, it is shown that extended observation of the broad scale hydrography by Argo provides a path to improved estimates of regional turbulent transport rates. Results from a 20 year ocean state estimate produced with the ECCO v4 non-linear inverse modeling framework provide supporting evidence. Turbulent transport parameter maps are estimated under the constraints of fitting the extensive collection of Argo profiles collected through 2011. The adjusted parameters dramatically reduce misfits to in situ profiles as compared with earlier ECCO solutions. They also yield a clear reduction in the model drift away from observations over multi-century long simulations, both for assimilated variables (temperature and salinity) and independent variables (bio-geochemical tracers). Despite the minimal constraints imposed specifically on the estimated parameters, their geography is physically plausible and exhibits close connections with the upper ocean ocean stratification as observed by Argo. The estimated parameter adjustments furthermore have first order impacts on upper-ocean stratification and mixed layer depths over 20 years. These results identify the constraint of fitting Argo profiles as an effective observational basis for regional turbulent transport rates. Uncertainties and further improvements of the method are discussed.

The influence of petrogenetic processes on the composition of Ocean Island Basalts

Ocean Island Basalts (OIB) provide important information on the chemical and physical characteristics of their mantle sources. However, the geochemical composition of a generated magma is significantly affected by partial melting and/or subsequent fractional crystallization processes. In addition, the isotopic composition of an ascending magma may be modified during transport through the oceanic crust. The influence of these different processes on the chemical and isotopic composition of OIB from two different localities, Hawaii and Tubuai in the Pacific Ocean, are investigated here. In a first chapter, the Os-isotope variations in suites of lavas from Kohala Volcano, Hawaii, are examined to constrain the role of melt/crust interactions on the evolution of these lavas. As 187Os/188Os sensitivity to any radiogenic contaminant strongly depend on the Os content in the melt, Os and other PGE variations are investigated first. This study reveals that Os and other PGE behavior change during the Hawaiian magma differentiation. While PGE concentrations are relatively constant in lavas with relatively primitive compositions, all PGE contents strongly decrease in the melt as it evolved through ~ 8% MgO. This likely reflects the sulfur saturation of the Hawaiian magma and the onset of sulfide fractionation at around 8% MgO. Kohala tholeiites with more than 8% MgO and rich in Os have homogeneous 187Os/188Os values likely to represent the mantle signature of Kohala lavas. However, Os isotopic ratios become more radiogenic with decreasing MgO and Os contents in the lavas, which reflects assimilation of local crust material during fractional crystallization processes. Less than 8% upper oceanic crust assimilation could have produced the most radiogenic Os-isotope ratios recorded in the shield lavas. However, these small amounts of upper crust assimilation have only negligible effects on Sr and Nd isotopic ratios and therefore, are not responsible for the Sr and Nd isotopic heterogeneities observed in Kohala lavas. In a second chapter, fractional crystallization and partial melting processes are constrained using major and trace element variations in the same suites of lavas from Kohala Volcano, Hawaii. This inverse modeling approach allows the estimation of most of the trace element composition of the Hawaiian mantle source. The calculated initial trace element pattern shows slight depletion of the concentrations from LREE to the most incompatible elements, which indicates that the incompatible element enrichments described by the Hawaiian melt patterns are entirely produced by partial melting processes. The “Kea trend” signature of lavas from Kohala Volcano is also confirmed, with Kohala lavas having lower Sr/Nd and La/Th ratios than lavas from Mauna Loa Volcano. Finally, the magmatic evolution of Tubuai Island is investigated in a last chapter using the trace element and Sr, Nd, Hf isotopic variations in mafic lava suites. The Sr, Nd and Hf isotopic data are homogeneous and typical for the HIMU-type OIB and confirms the cogenetic nature of the different mafic lavas from Tubuai Island. The trace element patterns show progressive enrichment of incompatible trace elements with increasing alkali content in the lavas, which reflect progressive decrease in the degree of partial melting towards the later volcanic events. In addition, this enrichment of incompatible trace elements is associated with relative depletion of Rb, Ba, K, Nb, Ta and Ti in the lavas, which require the presence of small amount of residual phlogopite and of a Ti-bearing phase (ilmenite or rutile) during formation of the younger analcitic and nephelinitic magmas.

Real-time control of an Electrodynamic Shaker

Questa tesi è essenzialmente focalizzata sullo sviluppo di un sistema di controllo in tempo reale per uno Shaker Elettrodinamico usato per riprodurre profili di vibrazione ambientale registrati in contesti reali e di interesse per il recupero di energia. Grazie all'utilizzo di uno shaker elettrodinamico è quindi possibile riprodurre scenari di vibrazione reale in laboratorio e valutare più agevolmente le prestazioni dei trasduttori meccanici. Tuttavia, è richiesto un controllo dello Shaker non solo in termini di stabilità ma anche per garantire l'esatta riproduzione del segnale registrato nel contesto reale. In questa tesi, si è scelto di sviluppare un controllo adattivo nel dominio del tempo per garantire la corretta riproduzione del profilo di accelerazione desiderato. L'algoritmo è stato poi implementato sul sistema di prototipazione rapida dSPACE DS1104 basata su microprocessore PowerPC. La natura adattiva dell'algoritmo proposto permette di identificare cambiamenti nella risposta dinamica del sistema, e di regolare di conseguenza i parametri del controllore. Il controllo del sistema è stato ottenuto anteponendo al sistema un filtro adattivo la cui funzione di trasferimento viene continuamente adattata per rappresentare al meglio la funzione di trasferimento inversa del sistema da controllare. Esperimenti in laboratorio confermano l'efficacia del controllo nella riproduzione di segnali reali e in tipici test di sweep frequenziale.

Integrated Geophysical Investigation of the St. James Fault Complex: A Case Study

We noninvasively detected the characteristics and location of a regional fault in an area of poor bedrock exposure complicated by karst weathering features in the subsurface. Because this regional fault is associated with sinkhole formation, its location is important for hazard avoidance. The bedrock lithologies on either side of the fault trace are similar; hence, we chose an approach that capitalized on the complementary strengths of very low frequency (VLF) electromagnetic, resistivity, and gravity methods. VLF proved most useful as a first-order reconnaissance tool, allowing us to define a narrow target area for further geophysical exploration. Fault-related epikarst was delineated using resistivity. Ultimately, a high-resolution gravity survey and subsequent inverse modeling using the results of the resistivity survey helped to further constrain the location and approximate orientation of the fault. The combined results indicated that the location of the fault trace needed to be adjusted 53 m south of the current published location and was consistent with a north-dipping thrust fault. Additionally, a gravity low south of the fault trace agreed with the location of conductive material from the resistivity and VLF surveys. We interpreted these anomalies to represent enhanced epikarst in the fault footwall. We clearly found that a staged approach involving a progression of methods beginning with a reconnaissance VLF survey, followed by high-resolution gravity and electrical resistivity surveys, can be used to characterize a fault and fault-related karst in an area of poor bedrock surface exposure.

ATMOSPHERIC RETRIEVAL ANALYSIS OF THE DIRECTLY IMAGED EXOPLANET HR 8799b

Directly imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b, for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity, and effective temperature of HR 8799b. We show that cloud-free models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediate cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity, and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.

(Table 1) Bulk and grain density, porosity, and p- and s-wave velocity characteristics of ODP Hole 163-990A sediments

Seismic velocities in rocks are influenced by the properties of the solid, the pore fluid, and the pore space. Cracks dramatically affect seismic velocities in rocks; their influence on the effective elastic moduli of rocks depends on their shape and concentration. Thin cracks (or fractures) substantially lower the moduli of a rock relative to the effect of spherical voids (or vesicles), and lower moduli are reflected by lower P- and S-wave velocities. The objective of this research is to determine the types and concentrations of cracks and their influence on the seismic properties of subaerially erupted basalts drilled from Hole 990A on the Southeast Greenland margin during Ocean Drilling Program Leg 163. Ellipsoidal cracks are used to model the voids in the rocks. The elastic moduli of the solid (grains) are also free parameters in the inverse modeling procedure. The apparent grain moduli reflect a weighted average of the moduli of the constituent minerals (e.g., plagioclase, augite, and clay minerals). The results indicate that (1) there is a strong relationship between P-wave velocity and porosity, suggesting a similarity of pore shape distributions, (2) the distribution of crack types within the massive, central region of aa flows from Hole 990A is independent of total porosity, (3) thin cracks are the first to be effectively sealed by alteration products, and (4) grain densities (an alteration index) and apparent grain moduli of the basalt samples are directly related.

Grain size and XRD analysis on Site 162-983

The grain size of deep-sea sediments provides an apparently simple proxy for current speed. However, grain size-based proxies may be ambiguous when the size distribution reflects a combination of processes, with current sorting only one of them. In particular, such sediment mixing hinders reconstruction of deep circulation changes associated with ice-rafting events in the glacial North Atlantic because variable ice-rafted detritus (IRD) input may falsely suggest current speed changes. Inverse modeling has been suggested as a way to overcome this problem. However, this approach requires high-precision size measurements that register small changes in the size distribution. Here we show that such data can be obtained using electrosensing and laser diffraction techniques, despite issues previously raised on the low precision of electrosensing methods and potential grain shape effects on laser diffraction. Down-core size patterns obtained from a sediment core from the North Atlantic are similar for both techniques, reinforcing the conclusion that both techniques yield comparable results. However, IRD input leads to a coarsening that spuriously suggests faster current speed. We show that this IRD influence can be accounted for using inverse modeling as long as wide size spectra are taken into account. This yields current speed variations that are in agreement with other proxies. Our experiments thus show that for current speed reconstruction, the choice of instrument is subordinate to a proper recognition of the various processes that determine the size distribution and that by using inverse modeling meaningful current speed reconstructions can be obtained from mixed sediments.

Origin of the Grande Ronde Formation flows, Columbia River flood basalt group

Lavas belonging to the Grande Ronde Formation (GRB) constitute about 63% of the Columbia River Basalt Group (CRBG), a flood basalt province in the NW United States. A puzzling feature is the lack of phenocrysts (< 5%) in these chemically evolved lavas. Based mainly on this observation it has been hypothesized that GRB lavas were nearly primary melts generated by large-scale melting of eclogite. Another recent hypothesis holds that GRB magmas were extremely hydrous and rose rapidly from the mantle such that the dissolved water kept the magmas close to their liquidi. I present new textural and chemical evidence to show that GRB lavas were neither primary nor hydrous melts but were derived from other melts via efficient fractional crystallization and mixing in shallow intrusive systems. Texture and chemical features further suggest that the melt mixing process may have been exothermic, which forced variable melting of some of the existing phenocrysts. ^ Finally, reported here are the results of efforts to simulate the higher pressure histories of GRB using COMAGMAT and MELTS softwares. The intent was to evaluate (1) whether such melts could be derived from primary melts formed by partial melting of a peridotite source as an alternative to the eclogite model, or if bulk melting of eclogite is required; and (2) at what pressure such primary melts could have been in equilibrium with the mantle. I carried out both forward and inverse modeling. The best fit forward model indicates that most primitive parent melts related to GRB could have been multiply saturated at ∼1.5--2.0 GPa. I interpret this result to indicate that the parental melts last equilibrated with a peridotitic mantle at 1.5--2.0 GPa and such partial melts rose to ∼0.2 GPa where they underwent efficient mixing and fractionation before erupting. These models suggest that the source rock was not eclogitic but a fertile spinel lherzolite, and that the melts had ∼0.5% water. ^

Results of Closed-Flow Column Experiments (zip-archive 2.3 MB)

The identification of transport parameters by inverse modeling often suffers from equifinality or parameter correlation when models are fitted to observations of the solute breakthrough in column outflow experiments. This parameters uncertainty can be approached by the application of multiple experimental designs such as column experiments in open-flow mode and the recently proposed closed-flow mode. Latter are characterized by the recirculation of the column effluent into the solution supply vessel that feeds the inflow. Depending on the experimental conditions, the solute concentration in the solution supply vessel and the effluent follows a damped sinusoidal oscillation. As a result, the closed-flow experiment provides additional observables in the breakthrough curve. The evaluation of these emergent features allows intrinsic control over boundary conditions and impacts the uncertainty of parameters in inverse modeling. We present a comprehensive sensitivity analysis to illustrate the potential application of closed-flow experiments. We show that the sensitivity with respect to the apparent dispersion can be controlled by the experimenter leading to a decrease in parameter uncertainty as compared to classical experiments by an order of magnitude for optimal settings. With these finding we are also able to reduce the equifinality found for situations, where rate-limited interactions impede a proper determination of the apparent dispersion and rate coefficients. Furthermore, we show the expected breakthrough curve for equilibrium and kinetic sorption, the latter showing strong similarities to the behavior found for completely mixed batch reactor experiments. This renders the closed-flow mode a useful complementary approach to classical column experiments.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

Inverse problems in dynamic cognitive modeling

Inverse problems for dynamical system models of cognitive processes comprise the determination of synaptic weight matrices or kernel functions for neural networks or neural/dynamic field models, respectively. We introduce dynamic cognitive modeling as a three tier top-down approach where cognitive processes are first described as algorithms that operate on complex symbolic data structures. Second, symbolic expressions and operations are represented by states and transformations in abstract vector spaces. Third, prescribed trajectories through representation space are implemented in neurodynamical systems. We discuss the Amari equation for a neural/dynamic field theory as a special case and show that the kernel construction problem is particularly ill-posed. We suggest a Tikhonov-Hebbian learning method as regularization technique and demonstrate its validity and robustness for basic examples of cognitive computations.