Wetland maps including open water extent dynamics based on ENVISAT ASAR WS for Siberia, 2007 and 2008, links to GeoTIFFs

Wetlands store large amounts of carbon, and depending on their status and type, they release specific amounts of methane gas to the atmosphere. The connection between wetland type and methane emission has been investigated in various studies and utilized in climate change monitoring and modelling. For improved estimation of methane emissions, land surface models require information such as the wetland fraction and its dynamics over large areas. Existing datasets of wetland dynamics present the total amount of wetland (fraction) for each model grid cell, but do not discriminate the different wetland types like permanent lakes, periodically inundated areas or peatlands. Wetland types differently influence methane fluxes and thus their contribution to the total wetland fraction should be quantified. Especially wetlands of permafrost regions are expected to have a strong impact on future climate due to soil thawing. In this study ENIVSAT ASAR Wide Swath data was tested for operational monitoring of the distribution of areas with a long-term SW near 1 (hSW) in northern Russia (SW = degree of saturation with water, 1 = saturated), which is a specific characteristic of peatlands. For the whole northern Russia, areas with hSW were delineated and discriminated from dynamic and open water bodies for the years 2007 and 2008. The area identified with this method amounts to approximately 300,000 km**2 in northern Siberia in 2007. It overlaps with zones of high carbon storage. Comparison with a range of related datasets (static and dynamic) showed that hSW represents not only peatlands but also temporary wetlands associated with post-forest fire conditions in permafrost regions. Annual long-term monitoring of change in boreal and tundra environments is possible with the presented approach. Sentinel-1, the successor of ENVISAT ASAR, will provide data that may allow continuous monitoring of these wetland dynamics in the future complementing global observations of wetland fraction.

(Table 1) Gas tie points used to correlate the ice cores NGRIP and EDML

Water isotope records from the EPICA Dronning Maud Land (EDML) and the NorthGRIP ice cores have revealed a one to one coupling between Antarctic Isotope Maxima (AIM) and Greenland Dansgaard-Oeschger (DO) events back to 50 kyr. In order to explore if this north-south coupling is persistent over Marine Isotopic Stage 5 (MIS 5), a common timescale must first be constructed. Here, we present new records of d18O of O2 (d18Oatm) and methane (CH4) measured in the air trapped in ice from the EDML (68-147 kyr) and NorthGRIP (70-123 kyr) ice cores. We demonstrate that, through the period of interest, CH4 records alone are not sufficient to construct a common gas timescale between the two cores. Millennial-scale variations of d18Oatm are evidenced over MIS 5 both on the Antarctic and Greenland ice cores and are coupled to CH4 profiles to synchronise the NorthGRIP and EDML records. They are shown to be a precious tool for ice core synchronisation. With this new dating strategy, we produce the first continuous and accurate sequence of the north-south climatic dynamics on a common ice timescale for the last glacial inception and the first DO events of MIS 5, reducing relative dating uncertainties to an accuracy of a few centuries at the onset of DO events 24 to 20. This EDML-NorthGRIP synchronisation provides new firm evidence that the bipolar seesaw is a pervasive pattern from the beginning of the glacial period. The relationship between Antarctic warming amplitudes and their concurrent Greenland stadial duration highlights the particularity of DO event 21 and its Antarctic counterpart. Our results suggest a smaller Southern Ocean warming rate for this long DO event compared to DO events of MIS 3.

Dynamics of halocarbons in coastal surface waters during short term mesocosm experiments

The aim of the present study was to evaluate the influence of different light quality, especially ultraviolet radiation (UVR), on the dynamics of volatile halogenated organic compounds (VHOCs) at the sea surface. Short term experiments were conducted with floating gas-tight mesocosms of different optical qualities. Six halocarbons (CH3I, CHCl3, CH2Br2, CH2ClI, CHBr3 and CH2I2), known to be produced by phytoplankton, together with a variety of biological and environmental variables were measured in the coastal southern Baltic Sea and in the Raunefjord (North Sea). These experiments showed that ambient levels of UVR have no significant influence on VHOC dynamics in the natural systems. We attribute it to the low radiation doses that phytoplankton cells receive in a normal turbulent surface mixed layer. The VHOC concentrations were influenced by their production and removal processes, but they were not correlated with biological or environmental parameters investigated. Diatoms were most likely the dominant biogenic source of VHOCs in the Baltic Sea experiment, whereas in the Raunefjord experiment macroalgae probably contributed strongly to the production of VHOCs. The variable stable carbon isotope signatures (d13C values) of bromoform (CHBr3) also indicate that different autotrophic organisms were responsible for CHBr3 production in the two coastal environments. In the Raunefjord, despite strong daily variations in CHBr3 concentration, the carbon isotopic ratio was fairly stable with a mean value of -26 per mil. During the declining spring phytoplankton bloom in the Baltic Sea, the d13C values of CHBr3 were enriched in 13C and showed noticeable diurnal changes (-12 per mil ± 4). These results show that isotope signature analysis is a useful tool to study both the origin and dynamics of VHOCs in natural systems.

Liquid lead and Lithium. A Comprehensive molecular dynamics study

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. One of the main issues in these programs is the problem of liquid metals breeder blanket behavior. Structural material of the blanket should meet high requirements because of extreme operating conditions. Therefore the knowledge of eutectic properties like optimal composition, physical and thermodynamic behavior or diffusion coefficients of Tritium are extremely necessary for current designs. In particular, the knowledge of the function linking the tritium concentration dissolved in liquid materials with the tritium partial pressure at a liquid/gas interface in equilibrium, CT=f(PT), is of basic importance because it directly impacts all functional properties of a blanket determining: tritium inventory, tritium permeation rate and tritium extraction efficiency. Nowadays, understanding the structure and behavior of this compound is a real goal in fusion engineering and materials science. Simulations of liquids can provide much information to the community; not only supplementing experimental data, but providing new tests of theories and ideas, making specific predictions that require experimental tests, and ultimately helping to lead to the deeper understanding and better predictive behavior.

Energy-Entropy-Momentum integration of discrete thermo-visco-elastic dynamics.

A novel time integration scheme is presented for the numerical solution of the dynamics of discrete systems consisting of point masses and thermo-visco-elastic springs. Even considering fully coupled constitutive laws for the elements, the obtained solutions strictly preserve the two laws of thermo dynamics and the symmetries of the continuum evolution equations. Moreover, the unconditional control over the energy and the entropy growth have the effect of stabilizing the numerical solution, allowing the use of larger time steps than those suitable for comparable implicit algorithms. Proofs for these claims are provided in the article as well as numerical examples that illustrate the performance of the method.

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos

Influence of the composition of natural gas-hydrogen blends (HCNG) on performance an CO2 emissions of a spark ignition engine

Addition of hydrogen to natural gas could be a short-term alternative to nowadays fossil fuels as the emissions of greenhouse gases may be reduced. The aim of this study is to evaluate the performance and emissions of a park ignition engine fuelled with pure natural gas, pure hydrogen and different blends of hydrogen and natural gas (HCNG). The increase of the hydrogen fraction leads to variations in the cylinder pressure and CO2 emissions. In this work, a combustion model based on thermodynamic equations is used considering separated zones for the burned and unburned gases. The results show that the maximum cylinder pressure gets higher as the fraction of hydrogen in the blend increases. The presence of hydrogen in the blend leads to a drecrease in the CO2 emissions. Due to hydrogen properties, leaner fuel-air mixtures can be used along with the appropiate spark timing, leading to an engine emissions improvement without a performance worsening.

Flujos de gases de efecto invernadero en suelos forestales y agroforestales del centro de la Península Ibérica=Greenhouse gas flows in forest and agroforestry soils in the center of the Iberian Peninsula

Debido a la complejidad de los procesos que controlan el intercambio de gases de carbono (C) y nitrógeno (N) entre el suelo y la atmósfera, en los sistemas forestales y agroforestales, son comprensibles las incógnitas existentes respecto a la estimación de los flujos de los gases de efecto invernadero (GEI) y la capacidad como reservorios de carbono de los suelos, bajo diferentes formas de uso y regímenes de alteración a escala regional y global. Esta escasez de información justifica la necesidad de caracterizar la dinámica de intercambio de GEI en los ecosistemas Mediterráneos, en especial en el contexto actual de cambio climático, y el incremento asociado de temperatura y periodos de sequía, alteración de los patrones de precipitación, y el riesgo de incendios forestales; cuyas consecuencias afectarán tanto a los compartimentos de C y de N del suelo como a la capacidad de secuestro de C de estos ecosistemas. Dentro de este contexto se enmarca la presente tesis doctoral cuyo objetivo ha sido cuantificar y caracterizar los flujos de dióxido de carbono (CO2), de oxido nitroso (N2O) y de metano (CH4), junto con los stocks de C y N, en suelos forestales de Quercus ilex, Quercus pyrenaica y Pinus sylvestris afectados por incendios forestales; así como el estudiar el efecto de la gestión y la cubierta arbórea en la respiración del suelo y los stocks de C y N en una dehesa situada en el centro de la Península Ibérica. De manera que los flujos de CO2, N2O y CH4; y los parámetros físico-químicos y biológicos del suelo fueron estudiados en los diferentes tratamientos y ecosistemas a lo largo del trabajo que se presenta. Los resultados obtenidos muestran la existencia de variaciones temporales y espaciales de la respiración del suelo dentro de una escala geográfica pequeña, controladas principalmente por la temperatura y la humedad del suelo; y por los contenidos de C y N del suelo en un bosque de Pinus sylvestris en la vertiente norte de la Sierra de Guadarrama , en España. El análisis de los efectos de los incendios forestales a largo plazo (6-8 años) revela que las pérdidas anuales de C a través de la respiración del suelo en las zonas quemadas de Quercus ilex, Quercus pyrenaica y Pinus sylvestris fueron 450 gCm-2yr-1, 790 gCm-2yr-1 y 1220 gCm-2yr-1, respectivamente; lo que representa una reducción del 43%, 22% y 11% en comparación con las zonas no quemadas de dichas especies, debido a la destrucción de la masa arbórea. El efecto del fuego también alteró los flujos N2O y CH4 del suelo, de una forma diferente en los distintos ecosistemas y estacionalidades estudiadas. De tal modo, que los suelos quemados mostraron una mayor oxidación del CH4 en las masas de Q. ilex, y una menor oxidación en las de P. sylvestris; además de una disminución de los flujos de N2O en Q. pyrenaica. Los incendios también afectaron los parámetros microclimáticos de los suelos forestales, observándose un incremento de la temperatura del suelo y una disminución de la humedad en los emplazamientos quemados que en los no quemados. Los cationes intercambiables, el pH, el cociente C/N, el contenido en raicillas y la biomasa microbiana también disminuyeron en las zonas quemadas. Aunque el C orgánico del suelo no se alteró de manera significativa, si lo hizo la calidad de la materia orgánica, disminuyendo el carbono lábil y aumentando las formas recalcitrantes lo que se tradujo en menor sensibilidad de la respiración del suelo a la temperatura (valores de Q10) en las zonas quemadas. Los resultados del estudio realizado en la Dehesa muestran que las actividades silvopastorales estudiadas afectaron levemente y de forma no constante a la respiración del suelo y las condiciones microclimáticas del suelo. Se observó una reducción 12% de la respiración del suelo por efecto del pastoreo no intensivo. Sin embargo, se observaron incrementos de 3Mg/ha en los stocks de C y de 0.3 Mg/ha en los stocks de N en los suelos pastoreados en comparación con los no pastoreados. Aunque, no se observó un claro efecto de la labranza sobre la respiración del suelo en nuestro experimento, sin embargo si se observó una disminución de 3.5 Mg/ha en las reservas de C y de 0.3 Mg/ ha en las de N en los suelos labrados comparados con los no labrados. La copa del arbolado influyó de forma positiva tanto en la respiración del suelo, como en los stocks de C y N de los suelos. La humedad del suelo jugó un papel relevante en la sensibilidad de la respiración a la temperatura del suelo. Nuestros resultados ponen de manifiesto la sensibilidad de la respiración del suelo a cambios en la humedad y los parámetros edáficos, y sugieren que la aplicación de modelos estándar para estimar la respiración del suelo en áreas geográficas pequeñas puede no ser adecuada a menos que otros factores sean considerados en combinación con la temperatura del suelo. Además, las diferentes respuestas de los flujos de gases de efecto invernadero a los cambios, años después de la ocurrencia de incendios forestales, destaca la necesidad de incluir estos cambios en las futuras investigaciones de la dinámica del carbono en los ecosistemas mediterráneos. Por otra parte, las respuestas divergentes en los valores de respiración del suelo y en los contenidos de C y N del suelo observados en la dehesa, además de la contribución de la copa de los árboles en los nutrientes del suelo ilustran la importancia de mantener la gestión tradicional aplicada en beneficio de la capacidad de almacenar C en la dehesa estudiada. La información obtenida en este trabajo pretende contribuir a la mejora del conocimiento de la dinámica y el balance de C en los sistemas mediterráneos, además de ayudar a predecir el impacto del cambio climático en el intercambio de C entre los ecosistemas forestales y agroforestales y la atmósfera. ABSTRACT Due to the complexity of the processes that control the exchange of carbon (C) and nitrogen (N) gasses between soils and the atmosphere in forest and agroforestry ecosystems, understandable uncertainties exist as regards the estimation of greenhouse gas (GHG) fluxes and the soil sink capacity at regional and global scale under different forms of land use and disturbance regimes. These uncertainties justify the need to characterize the exchange dynamics of GHG between the atmosphere and soils in Mediterranean terrestrial ecosystems, particularly in the current context of climate change and the associated increase in temperature, drought periods, heavy rainfall events, and increased risk of wildfires, which affect not only the C and N pools but also the soil C sink capacity of these ecosystems. Within this context, the aims of the present thesis were, firstly, to quantify and characterize the fluxes of carbon dioxide (CO2), nitrous oxide (N2O) and methane (CH4) as well as the C and N stocks in Quercus ilex, Quercus pyrenaica and Pinus sylvestris stands affected by wildfires, and secondly, to study the effects of Quercus ilex canopy and management on both soil respiration and C and N pools in dehesa systems in the center of Iberian Peninsula. Soil CO2, N2O and CH4 fluxes, and soil physical-chemical and biological parameters were studied under the different treatments and ecosystems considered in this study. The results showed seasonal and spatial variations in soil respiration within small geographic areas, mainly controlled by soil temperature and moisture in addition to soil carbon and nitrogen stocks in mixed pine–oak forest ecosystems on the north facing slopes of the Sierra de Guadarrama in Spain. The analysis of long term effects of wildfires (6–8 years) revealed that annual carbon losses through soil respiration from burned sites in Quercus ilex, Quercus pyrenaica and Pinus sylvestris stands were 450 gCm-2yr-1, 790 gCm-2yr-1 and 1220 gCm-2yr-1, respectively; with burned sites emitting 43%, 22% and 11% less in burned as opposed to non-burned sites due the loss of trees. Fire may alter both N2O and CH4 fluxes although the magnitude of such variation depends on the site, soil characteristics and seasonal climatic conditions. The burned sites showed higher CH4 oxidation in Q.ilex stands, and lower oxidation rates in P. sylvestris stands. A reduction in N2O fluxes in Q. pyrenaica stands was detected at burned sites along with changes in soil microclimate; higher soil temperature and lower soil moisture content. Exchangeable cations, the C/N ratio, pH, fine root and microbial biomass were also found to decrease at burned sites. Although the soil organic carbon was not significantly altered, the quality of the organic matter changed, displaying a decrease in labile carbon and a relative increase in refractory forms, leading to lower sensitivity of soil respiration to temperature (Q10 values) at burned sites. The results from the dehesa study show that light grazing and superficial tilling practices used in the studied dehesa system in Spain had a slight but non-consistent impact on soil respiration and soil microclimate over the study period. The reduction in soil respiration in the dehesa system due to the effects of grazing was around 12 %. However, increments of 3Mg/ha in C stocks and 0.3 Mg/ha in N stocks in grazed soils were observed. Although no clear effect of tilling on soil respiration was found, a decrease of 3.5 Mg/ha in C stocks and 0.3 Mg/ha in N stocks was detected for tilled soils. The presence of a tree canopy induced increases in soil respiration, soil C and N stocks, while soil moisture was found to play an important role in soil respiration temperature response. Our results suggest that the use of standard models to estimate soil respiration in small geographical areas may not be adequate unless other factors are considered in addition to soil temperature. Furthermore, the different responses of GHG flux to climatic shifts, many years after the occurrence of wildfire, highlight the need to include these shifts in C dynamics in future research undertaken in Mediterranean ecosystems. Furthermore, divergent responses in soil respiration and soil C and N stocks to grazing or tilling practices in Dehesa systems, and the influence of tree canopy on soil respiration and soil nutrient content, illustrate the importance of maintaining beneficial management practices. Moreover, the carbon sequestration capacity of the Dehesa system studied may be enhanced through improvements in the management applied. It is hoped that the information obtained through this research will contribute towards improving our understanding of the dynamics and balance of C in Mediterranean systems, and help predict the impact of climate change on the exchange of C between forest and agroforestry ecosystems and the atmosphere.

Relative Dynamics and Control of an Ion Beam Shepherd Satellite

The ion beam shepherd (IBS) is a recently proposed concept for modifying the orbit and/or attitude of a generic orbiting body in a contactless manner, which makes it a candidate technology for active space debris removal. In this paper we deal with the problem of controlling the relative position of a shepherd satellite coorbiting at small separation distance with a target debris. After deriving the orbit relative motion equations including the effect of the ion beam perturbation we study the system stability and propose different control strategies.

Web based interactive educational software introducing semiconductor laser dynamics: Sound Of Lasers (SOL)

In this work, educational software for intuitive understanding of the basic dynamic processes of semiconductor lasers is presented. The proposed tool is addressed to the students of optical communication courses, encouraging self consolidation of the subjects learned in lectures. The semiconductor laser model is based on the well known rate equations for the carrier density, photon density and optical phase. The direct modulation of the laser is considered with input parameters which can be selected by the user. Different options for the waveform, amplitude and frequency of thpoint. Simulation results are plotted for carrier density and output power versus time. Instantaneous frequency variations of the laser output are numerically shifted to the audible frequency range and sent to the computer loudspeakers. This results in an intuitive description of the “chirp” phenomenon due to amplitude-phase coupling, typical of directly modulated semiconductor lasers. In this way, the student can actually listen to the time resolved spectral content of the laser output. By changing the laser parameters and/or the modulation parameters,consequent variation of the laser output can be appreciated in intuitive manner. The proposed educational tool has been previously implemented by the same authors with locally executable software. In the present manuscript, we extend our previous work to a web based platform, offering improved distribution and allowing its use to the wide audience of the web.

Chaos in non lineal Alfven waves using the DNLS equations

The electro-dynamical tethers emit waves in structured denominated Alfven wings. The Derivative Nonlineal Schrödinger Equation (DNLS) possesses the capacity to describe the propagation of circularly polarized Alfven waves of finite amplitude in cold plasmas. The DNLS equation is truncated to explore the coherent, weakly nonlinear, cubic coupling of three waves near resonance, one wave being linearly unstable and the other waves damped. In this article is presented a theoretical and numerical analysis when the growth rate of the unstable wave is next to zero considering two damping models: Landau and resistive. The DNLS equation presents a chaotic dynamics when is consider only three wave truncation. The evolution to chaos possesses three routes: hard transition, period-doubling and intermittence of type I.

Numerical simulation of group combustion of pulverized coal

A mathematical model for the group combustion of pulverized coal particles was developed in a previous work. It includes the Lagrangian description of the dehumidification, devolatilization and char gasification reactions of the coal particles in the homogenized gaseous environment resulting from the three fuels, CO, H2 and volatiles, supplied by the gasification of the particles and their simultaneous group combustion by the gas phase oxidation reactions, which are considered to be very fast. This model is complemented here with an analysis of the particle dynamics, determined principally by the effects of aerodynamic drag and gravity, and its dispersion based on a stochastic model. It is also extended to include two other simpler models for the gasification of the particles: the first one for particles small enough to extinguish the surrounding diffusion flames, and a second one for particles with small ash content when the porous shell of ashes remaining after gasification of the char, non structurally stable, is disrupted. As an example of the applicability of the models, they are used in the numerical simulation of an experiment of a non-swirling pulverized coal jet with a nearly stagnant air at ambient temperature, with an initial region of interaction with a small annular methane flame. Computational algorithms for solving the different stages undergone by a coal particle during its combustion are proposed. For the partial differential equations modeling the gas phase, a second order finite element method combined with a semi-Lagrangian characteristics method are used. The results obtained with the three versions of the model are compared among them and show how the first of the simpler models fits better the experimental results.