Community structure and dynamics in climate networks

We consider climate networks constructed from observed and model simulated fields of three climate variables and investigate their community structure. We find that for all fields the number of effective communities is rather small (four to five). We are able to trace the origin of these communities to certain dynamical properties of climate. Our results suggest that the complete complexity of the climate system condenses beyond the `weather` time scales into a small number of low-dimensional interacting components and provide clues as to the nature of the climate subsystems underlying these components.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

Dynamics of Space Particles and Spacecrafts Passing by the Atmosphere of the Earth

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Speleoclimate dynamics in Santana Cave (PETAR, S(a)over-tildeo Paulo State, Brazil): general characterization and implications for tourist management

Show caves provide tourists with the opportunity to have close contact with natural underground spaces. However, visitation to these places also creates a need for management measures, mainly the definition of tourist carrying capacity. The present work describes the results of climate monitoring and atmospheric profiling performed in Santana Cave (Alto Ribeira State and Tourist Park - PETAR, Brazil) between 2008 and 2011. Based on the results, distinct preliminary zones with different levels of thermal variation were identified, which classify Santana Cave as a warm trap. Two critical points along the tourist route (Cristo and Encontro Halls) were identified where the temperature of the locality increased by 1.3 degrees C when tourists were present. Air flow from the inner cave to the outside occurs during the austral summer, and the opposite flow occurs when the outside environment is colder than the air inside the cave during the austral winter. The temperature was used to establish thresholds to the tourist carrying capacity by computing the recovery time of the atmospheric conditions after the changes caused by the presence of tourists. This method suggests a maximum limit of approximately 350 visits per day to Santana Cave. The conclusion of the study is that Santana Cave has an atmosphere that is highly connected with the outside; daily variations in temperature and, to a lesser extent, in the relative humidity occur throughout the entire studied area of the cave. Therefore, the tourist carrying capacity in Santana Cave can be flexible and can be implemented based on the climate seasonality, the tourism demand and other management strategies.

Towards the construction of a validated numerical system to study the mesoscale dynamics of the North East Atlantic (2003-2006)

Máster en Oceanografía

New insights into POC dynamics in the subtropical northeast Atlantic Ocean

[EN]The capacity of the ocean to sequester atmospheric carbon (CO2) depends to a large extent on the dynamics of biogenic carbon in the water column. However, most current global and regional estimates of carbon balances are solely based on particles collected with drifting and moored sediment traps. As a consequence, construction of ocean carbon budgets has long been guided by the simplification introduced by sediment traps, which give a 1D vision of the whole picture. In this thesis we have assessed a quantitative analysis of the flux magnitude and the mechanisms of transport of the whole particle spectrum (suspended, slowly-sinking and sinking particles).

Lagrangian transport evaluation of atmospheric chemistry in the Mediterranean region

Global observations of the chemical composition of the atmosphere are essential for understanding and studying the present and future state of the earth's atmosphere. However, by analyzing field experiments the consideration of the atmospheric motion is indispensable, because transport enables different chemical species, with different local natural and anthropogenic sources, to interact chemically and so consequently influences the chemical composition of the atmosphere. The distance over which that transport occurs is highly dependent upon meteorological conditions (e.g., wind speed, precipitation) and the properties of chemical species itself (e.g., solubility, reactivity). This interaction between chemistry and dynamics makes the study of atmospheric chemistry both difficult and challenging, and also demonstrates the relevance of including the atmospheric motions in that context. In this doctoral thesis the large-scale transport of air over the eastern Mediterranean region during summer 2001, with a focus on August during the Mediterranean Intensive Oxidant Study (MINOS) measurement campaign, was investigated from a lagrangian perspective. Analysis of back trajectories demonstrated transport of polluted air masses from western and eastern Europe in the boundary layer, from the North Atlantic/North American area in the middle end upper troposphere and additionally from South Asia in the upper troposphere towards the eastern Mediterranean. Investigation of air mass transport near the tropopause indicated enhanced cross-tropopause transport relative to the surrounding area over the eastern Mediterranean region in summer. A large band of air mass transport across the dynamical tropopause develops in June, and is shifted toward higher latitudes in July and August. This shifting is associated with the development and the intensification of the Arabian and South Asian upper-level anticyclones and consequential with areas of maximum clear-air turbulence, hypothesizing quasi-permanent areas with turbulent mixing of tropospheric and stratospheric air during summer over the eastern Mediterranean as a result of large-scale synoptic circulation. In context with the latex knowledge about the transport of polluted air masses towards the Mediterranean and with increasing emissions, especially in developing countries like India, this likely gains in importance.

Influence of nitrogen and soil physical characteristics on belowground carbon flux dynamics of woody plants

At ecosystem level soil respiration (Rs) represents the largest carbon (C) flux after gross primary productivity, being mainly generated by root respiration (autotrophic respiration, Ra) and soil microbial respiration (heterotrophic respiration, Rh). In the case of terrestrial ecosystems, soils contain the largest C-pool, storing twice the amount of C contained in plant biomass. Soil organic matter (SOM), representing the main C storage in soil, is decomposed by soil microbial community. This process produces CO2 which is mainly released as Rh. It is thus relevant to understand how microbial activity is influenced by environmental factors like soil temperature, soil moisture and nutrient availability, since part of the CO2 produced by Rh, directly increases atmospheric CO2 concentration and therefore affects the phenomenon of climate change. Among terrestrial ecosystems, agricultural fields have traditionally been considered as sources of atmospheric CO2. In agricultural ecosystems, in particular apple orchards, I identified the role of root density, soil temperature, soil moisture and nitrogen (N) availability on Rs and on its two components, Ra and Rh. To do so I applied different techniques to separate Rs in its two components, the ”regression technique” and the “trenching technique”. I also studied the response of Ra to different levels of N availability, distributed either in a uniform or localized way, in the case of Populus tremuloides trees. The results showed that Rs is mainly driven by soil temperature, to which it is positively correlated, that high levels of soil moisture have inhibiting effects, and that N has a negligible influence on total Rs, as well as on Ra. Further I found a negative response of Rh to high N availability, suggesting that microbial decomposition processes in the soil are inhibited by the presence of N. The contribution of Ra to Rs was of 37% on average.

Temporal dynamics of benthic assemblages along a gradient of ocean acidification at a CO2 vent’s system of the Ischia Island

Global warming and ocean acidification, due to rising atmospheric levels of CO2, represent an actual threat to terrestrial and marine environments. Since Industrial Revolution, in less of 250 years, pH of surface seawater decreased on average of 0.1 unit, and is expected to further decreases of approximately 0.3-0.4 units by the end of this century. Naturally acidified marine areas, such as CO2 vent systems at the Ischia Island, allow to study acclimatation and adaptation of individual species as well as the structure of communities, and ecosystems to OA. The main aim of this thesis was to study how hard bottom sublittoral benthic assemblages changed trough time along a pH gradient. For this purpose, the temporal dynamics of mature assemblages established on artificial substrates (volcanic tiles) over a 3 year- period were analysed. Our results revealed how composition and dynamics of the community were altered and highly simplified at different level of seawater acidification. In fact, extreme low values of pH (approximately 6.9), affected strongly the assemblages, reducing diversity both in terms of taxa and functional groups, respect to lower acidification levels (mean pH 7.8) and ambient conditions (8.1 unit). Temporal variation was observed in terms of species composition but not in functional groups. Variability was related to species belonging to the same functional group, suggesting the occurrence of functional redundancy. Therefore, the analysis of functional groups kept information on the structure, but lost information on species diversity and dynamics. Decreasing in ocean pH is only one of many future global changes that will occur at the end of this century (increase of ocean temperature, sea level rise, eutrophication etc.). The interaction between these factors and OA could exacerbate the community and ecosystem effects showed by this thesis.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1 clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Circulation dynamics and its influence on European and Mediterranean January-April climate over the past half millennium: results and insights from instrumental data, documentary evidence and coupled climate models

We use long instrumental temperature series together with available field reconstructions of sea-level pressure (SLP) and three-dimensional climate model simulations to analyze relations between temperature anomalies and atmospheric circulation patterns over much of Europe and the Mediterranean for the late winter/early spring (January–April, JFMA) season. A Canonical Correlation Analysis (CCA) investigates interannual to interdecadal covariability between a new gridded SLP field reconstruction and seven long instrumental temperature series covering the past 250 years. We then present and discuss prominent atmospheric circulation patterns related to anomalous warm and cold JFMA conditions within different European areas spanning the period 1760–2007. Next, using a data assimilation technique, we link gridded SLP data with a climate model (EC-Bilt-Clio) for a better dynamical understanding of the relationship between large scale circulation and European climate. We thus present an alternative approach to reconstruct climate for the pre-instrumental period based on the assimilated model simulations. Furthermore, we present an independent method to extend the dynamic circulation analysis for anomalously cold European JFMA conditions back to the sixteenth century. To this end, we use documentary records that are spatially representative for the long instrumental records and derive, through modern analogs, large-scale SLP, surface temperature and precipitation fields. The skill of the analog method is tested in the virtual world of two three-dimensional climate simulations (ECHO-G and HadCM3). This endeavor offers new possibilities to both constrain climate model into a reconstruction mode (through the assimilation approach) and to better asses documentary data in a quantitative way.

Hydration of the Bisulfate Ion: Atmospheric Implications

Using molecular dynamics configurational sampling combined with ab initio energy calculations, we determined the low energy isomers of the bisulfate hydrates. We calculated the CCSD(T) complete basis set (CBS) binding electronic and Gibbs free energies for 53 low energy isomers of HSO4–(H2O)n=1–6 and derived the thermodynamics of adding waters sequentially to the bisulfate ion and its hydrates. Comparing the HSO4–/H2O system to the neutral H2SO4/H2O cluster, water binds more strongly to the anion than it does to the neutral molecules. The difference in the binding thermodynamics of HSO4–/H2O and H2SO4/H2O systems decreases with increasing number of waters. The thermodynamics for the formation of HSO4–(H2O)n=1–5 is favorable at 298.15 K, and that of HSO4–(H2O)n=1–6 is favorable for T < 273.15 K. The HSO4– ion is almost always hydrated at temperatures and relative humidity values encountered in the troposphere. Because the bisulfate ion binds more strongly to sulfuric acid than it does to water, it is expected to play a role in ion-induced nucleation by forming a strong complex with sulfuric acid and water, thus facilitating the formation of a critical nucleus.

Atmospheric Implications for Formation of Clusters of Ammonium and 1−10 Water Molecules

A mixed molecular dynamics/quantum mechanics model has been applied to the ammonium/water clustering system. The use of the high level MP2 calculation method and correlated basis sets, such as aug-cc-pVDZ and aug-cc-pVTZ, lends confidence in the accuracy of the extrapolated energies. These calculations provide electronic and free energies for the formation of clusters of ammonium and 1−10 water molecules at two different temperatures. Structures and thermodynamic values are in good agreement with previous experimental and theoretical results. The estimated concentration of these clusters in the troposphere was calculated using atmospheric amounts of ammonium and water. Results show the favorability of forming these clusters and implications for ion-induced nucleation in the atmosphere.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.