Modelling study of the impact of deep convection on the utls air composition - Part I: Analysis of ozone precursors

The aim of this work is to study the local impact on the upper troposphere/lower stratosphere air composition of an extreme deep convective system. For this purpose, we performed a simulation of a convective cluster composed of many individual deep convective cells that occurred near Bauru (Brazil). The simulation is performed using the 3-D mesoscale model RAMS coupled on-line with a chemistry model. The comparisons with meteorological measurements show that the model produces meteorological fields generally consistent with the observations. The present paper (part I) is devoted to the analysis of the ozone precursors (CO, NO x and non-methane volatile organic compounds) and HO x in the UTLS. The simulation results show that the distribution of CO with altitude is closely related to the upward convective motions and consecutive outflow at the top of the convective cells leading to a bulge of CO between 7 km altitude and the tropopause (around 17km altitude). The model results for CO are consistent with satellite-borne measurements at 700 hPa. The simulation also indicates enhanced amounts of NO x up to 2 ppbv in the 7-17 km altitude layer mainly produced by the lightning associated with the intense convective activity. For insoluble non-methane volatile organic compounds, the convective activity tends to significantly increase their amount in the 7-17km layer by dynamical effects. During daytime in the presence of lightning NO x, this bulge is largely reduced in the upper part of the layer for reactive species (e.g. isoprene, ethene) because of their reactions with OH that is increased on average during daytime. Lightning NO x also impacts on the oxydizing capacity of the upper troposphere by reducing on average HO x, HO 2, H 2O 2 and organic hydroperoxides. During the simulation time, the impact of convection on the air composition of the lower stratosphere is negligible for all ozone precursors although several of the simulated convective cells nearly reach the tropopause. There is no significant transport from the upper troposphere to the lower stratosphere, the isentropic barrier not being crossed by convection. The impact of the increase of ozone precursors and HO x in the upper troposphere on the ozone budget in the LS is discussed in part II of this series of papers.

Including the sub-grid scale plume rise of vegetation fires in low resolution atmospheric transport models

We describe and begin to evaluate a parameterization to include the vertical transport of hot gases and particles emitted from biomass burning in low resolution atmospheric-chemistry transport models. This sub-grid transport mechanism is simulated by embedding a 1-D cloud-resolving model with appropriate lower boundary conditions in each column of the 3-D host model. Through assimilation of remote sensing fire products, we recognize which columns have fires. Using a land use dataset appropriate fire properties are selected. The host model provides the environmental conditions, allowing the plume rise to be simulated explicitly. The derived height of the plume is then used in the source emission field of the host model to determine the effective injection height, releasing the material emitted during the flaming phase at this height. Model results are compared with CO aircraft profiles from an Amazon basin field campaign and with satellite data, showing the huge impact that this mechanism has on model performance. We also show the relative role of each main vertical transport mechanisms, shallow and deep moist convection and the pyro-convection (dry or moist) induced by vegetation fires, on the distribution of biomass burning CO emissions in the troposphere.

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

The influence of recent major crater impacts on the surrounding surfaces of (21) Lutetia

We present 3-D simulations of impacts into Asteroid 21 Lutetia, the subject of a fly-by by the European Space Agency's Rosetta mission to Comet 67P/Churyumov-Gerasimenko. Using a 3-D shape model of the asteroid, impacts of sizes sufficient to reproduce the observed craters in Lutetia's North Polar Crater Cluster (NPCC) as observed by the OSIRIS experiment have been simulated using the Smoothed Particle Hydrodynamics technique. The asteroid itself has been modelled both as a homogeneous body and as a body with an iron core.

Comparison of several models for calculating ozone stomatal fluxes on a Mediterranean wheat cultivar (Triticum durum Desf. cv. Camacho)

Ozone stomatal fluxes were modeled for a 3-year period following different approaches for a commercial variety of durum wheat (Triticum durum Desf. cv. Camacho) at the phenological stage of anthesis. All models performed in the same range, although not all of them afforded equally significant results. Nevertheless, all of them suggest that stomatal conductance would account for the main percentage of ozone deposition fluxes. A new modeling approach was tested, based on a 3-D architectural model of the wheat canopy, and fairly accurate results were obtained. Plant species-specific measurements, as well as measurements of stomatal conductance and environmental parameters, were required. The method proposed for calculating ozone stomatal fluxes (FO(3_3-D)) from experimental gs data and modeling them as a function of certain environmental parameters in conjunction with the use of the YPLANT model seems to be adequate, providing realistic estimates of the canopy FO(3_3-D), integrating and not neglecting the contribution of the lower leaves with respect to the flag leaf, although a further development of this model is needed.

A non-linear analysis of lying a floating pipeline on the seabed

In this article, a model for the determination of displacements, deformations and tensions of a submarine pipeline during the construction is presented. The process is carried out from an initial floating situation to the final laying position on the seabed. The existence of currents and small waves are also considered. Firstly, this technique, usually applied to polyethylene pipelines, is described in this paper as well as some real world examples, as well as the variables that can be modified to control the behavior of the structure. A detailed description of the actions in this process is considered, specially the ones related to marine environment, as Archimedes force, current and sea waves. The behavior of the pipeline is modeled with a non linear elasto dynamic model where geometric non linearities are taken into account. A 3-D beam model, without cross section deformation effects, is developed. Special care is taken in the numerical analysis, developed within an updated lagrangian formulation framework, with the sea bed contact, the follower forces due to the external water pressures and the dynamic actions. Finally, some subroutines are implemented into ANSYS to simulate the two dimensional case, where the whole construction process is achieved. With this software, a sensibility analysis of the bending moments, axial forces and stresses obtained with different values of the control variables in order to optimize the construction steps. These control variables are, the axial load in the pipe, the inundated inner length and the distance of the control barge from the coast.

Three-dimensional modeling of the human eye based on magnetic resonance imaging

PURPOSE. A methodology for noninvasively characterizing the three-dimensional (3-D) shape of the complete human eye is not currently available for research into ocular diseases that have a structural substrate, such as myopia. A novel application of a magnetic resonance imaging (MRI) acquisition and analysis technique is presented that, for the first time, allows the 3-D shape of the eye to be investigated fully. METHODS. The technique involves the acquisition of a T2-weighted MRI, which is optimized to reveal the fluid-filled chambers of the eye. Automatic segmentation and meshing algorithms generate a 3-D surface model, which can be shaded with morphologic parameters such as distance from the posterior corneal pole and deviation from sphericity. Full details of the method are illustrated with data from 14 eyes of seven individuals. The spatial accuracy of the calculated models is demonstrated by comparing the MRI-derived axial lengths with values measured in the same eyes using interferometry. RESULTS. The color-coded eye models showed substantial variation in the absolute size of the 14 eyes. Variations in the sphericity of the eyes were also evident, with some appearing approximately spherical whereas others were clearly oblate and one was slightly prolate. Nasal-temporal asymmetries were noted in some subjects. CONCLUSIONS. The MRI acquisition and analysis technique allows a novel way of examining 3-D ocular shape. The ability to stratify and analyze eye shape, ocular volume, and sphericity will further extend the understanding of which specific biometric parameters predispose emmetropic children subsequently to develop myopia. Copyright © Association for Research in Vision and Ophthalmology.

Shape of the posterior vitreous chamber in human emmetropia and myopia

PURPOSE - To compare posterior vitreous chamber shape in myopia to that in emmetropia. METHODS - Both eyes of 55 adult subjects were studied, 27 with emmetropia (MSE =-0.55; <+0.75D; mean +0.09 ±0.36D) and 28 with myopia (MSE -5.87 ±2.31D). Cycloplegic refraction was measured with a Shin Nippon autorefractor and anterior chamber depth and axial length with a Zeiss IOLMaster. Posterior vitreous chamber shapes were determined from T2-weighted MRI (3-Tesla) using procedures previously reported by our laboratory. 3-D surface model coordinates were assigned to nasal, temporal, superior and inferior quadrants and plotted in 2-D to illustrate the composite shape of respective quadrants posterior to the second nodal point. Spherical analogues of chamber shape were constructed to compare relative sphericity between refractive groups and quadrants. RESULTS - Differences in shape occurred in the region posterior to points of maximum globe width and were thus in general accord with an equatorial model of myopic expansion. Shape in emmetropia is categorised distinctly as that of an oblate ellipse and in myopia as an oblate ellipse of significantly less degree such that it approximates to a sphere. There was concordance between shape and retinotopic projection of respective quadrants into right, left, superior and inferior visual fields. CONCLUSIONS - The transition in shape from oblate ellipse to sphere with axial elongation supports the hypothesis that myopia may be a consequence of equatorial restriction associated with biomechanical anomalies of the ciliary apparatus. The synchronisation of quadrant shapes with retinotopic projection suggests that binocular growth is coordinated by processes that operate beyond the optic chiasm.

Developmental vitamin D deficiency alters adult behaviour: An informative animal model of schizophrenia

Background: There is growing evidence that vitamin D is active in the brain but until recently there was a lack of evidence about its role during brain development. Guided by certain features of the epidemiology of schizophrenia, we have explored the role of vitamin D in the developing brain and behaviour using whole animal models. Methods: Sprague-Dawley rats were fed a vitamin D deficient diet (DVD) or control diet 6 weeks prior to mating and housed under UVB-free lighting conditions. On the day of birth all rats were fed a control diet for the remainder of the study. We observed behaviour at two timepoints; on the day of birth to study maternal behaviour, and at 10 weeks of age to study offspring behaviour in adulthood, under baseline and drug induced conditions (MK-801, haloperidol, amphetamine). Results: Prenatal vitamin D deficiency results in subtle alterations in maternal behaviour as well as long lasting effects on the adult offspring, despite a return to normal vitamin D levels during postnatal life. These affects were specific to transient prenatal vitamin D depletion as adult vitamin D depletion, combined prenatal and chronic postnatal vitamin D depletion, or ablation of the vitamin D receptor in mice led to markedly different outcomes. Conclusions: The developmental vitamin D (DVD) model now draws strength from epidemiological evidence of schizophrenia and animal experiments. Although the DVD model does not replicate every aspect of schizophrenia, it has several attractive features: (1) the exposure is based on clues from epidemiology; (2) it reproduces the increase in lateral ventricles; (3) it reproduces well-regarded behavioural phenotypes associated with schizophrenia (e.g. MK- 801 induced hyperlocomotion); and (4) it implicates a disturbance in dopamine signaling. In summary, low prenatal levels of vitamin D can influence critical components of orderly brain development and that this has a long lasting effect on behaviour.

Signatures of doubly charged Higgs in the SU(3)(L) circle times U(1)(N) model at the CERN LHC

The scalar sector of the simplest version of the 3-3-1 electroweak model is constructed with three Higgs triplets only. We show that a relation involving two of the constants of the model, two vacuum expectation values of the neutral scalars, and the mass of the doubly charged Higgs boson leads to important information concerning the signals of this scalar particle.

Calculable lepton masses, seesaw relations, and four neutrino mixings in a 3-3-1 model with an extra U(1) symmetry

We propose a scheme in which the masses of the heavier leptons obey seesaw type relations. The light lepton masses, except the electron and the electron neutrino ones, are generated by one loop level radiative corrections. We work in a version of the 3-3-1 electroweak model that predicts singlets (charged and neutral) of heavy leptons beyond the known ones. An extra U(1)(Omega) symmetry is introduced in order to avoid the light leptons getting masses at the tree level. The electron mass induces an explicit symmetry breaking at U(1)(Omega). We discuss also the mixing matrix among four neutrinos. The new energy scale required is not higher than a few TeV.

Supersymmetric 3-3-1 model

We build a complete supersymmetric version of a 3-3-1 gauge model using the superfield formalism. We point out that a discrete symmetry, similar to R symmetry in the minimal supersymmetric standard model, is possible to be defined in this model. Hence we have both R-conservina and R-violating possibilities. Analysis of the mass spectrum of the neutral real scalar fields show that in this model the lightest scalar Higgs boson has a mass upper limit, and at the tree level it is 124.5 GeV for a given illustrative set of parameters.

Lepton masses in a supersymmetric 3-3-1 model

We consider the mass generation for both charginos and neutralinos in a 3-3-1 supersymmetric model. We show that R-parity breaking interactions leave the electron and one of the neutrinos massless at the tree level. However, the same interactions induce masses for these particles at the 1-loop level. Unlike the similar situation in the minimal supersymmetric standard model, the masses of the neutralinos are related to the masses of the charginos.

Yang-Lee zeros of the two- and three-state Potts model defined on π3 Feynman diagrams

We present both analytical and numerical results on the position of partition function zeros on the complex magnetic field plane of the q=2 state (Ising) and the q=3 state Potts model defined on phi(3) Feynman diagrams (thin random graphs). Our analytic results are based on the ideas of destructive interference of coexisting phases and low temperature expansions. For the case of the Ising model, an argument based on a symmetry of the saddle point equations leads us to a nonperturbative proof that the Yang-Lee zeros are located on the unit circle, although no circle theorem is known in this case of random graphs. For the q=3 state Potts model, our perturbative results indicate that the Yang-Lee zeros lie outside the unit circle. Both analytic results are confirmed by finite lattice numerical calculations.

Development of a 3-Dimensional Mathematical Collision Risk Model Based on Recorded Aircraft Trajectories to Estimate the Safety Level in High Density En-Route Airspaces

Como en todos los medios de transporte, la seguridad en los viajes en avión es de primordial importancia. Con los aumentos de tráfico aéreo previstos en Europa para la próxima década, es evidente que el riesgo de accidentes necesita ser evaluado y monitorizado cuidadosamente de forma continúa. La Tesis presente tiene como objetivo el desarrollo de un modelo de riesgo de colisión exhaustivo como método para evaluar el nivel de seguridad en ruta del espacio aéreo europeo, considerando todos los factores de influencia. La mayor limitación en el desarrollo de metodologías y herramientas de monitorización adecuadas para evaluar el nivel de seguridad en espacios de ruta europeos, donde los controladores aéreos monitorizan el tráfico aéreo mediante la vigilancia radar y proporcionan instrucciones tácticas a las aeronaves, reside en la estimación del riesgo operacional. Hoy en día, la estimación del riesgo operacional está basada normalmente en reportes de incidentes proporcionados por el proveedor de servicios de navegación aérea (ANSP). Esta Tesis propone un nuevo e innovador enfoque para evaluar el nivel de seguridad basado exclusivamente en el procesamiento y análisis trazas radar. La metodología propuesta ha sido diseñada para complementar la información recogida en las bases de datos de accidentes e incidentes, mediante la provisión de información robusta de los factores de tráfico aéreo y métricas de seguridad inferidas del análisis automático en profundidad de todos los eventos de proximidad. La metodología 3-D CRM se ha implementado en un prototipo desarrollado en MATLAB © para analizar automáticamente las trazas radar y planes de vuelo registrados por los Sistemas de Procesamiento de Datos Radar (RDP) e identificar y analizar todos los eventos de proximidad (conflictos, conflictos potenciales y colisiones potenciales) en un periodo de tiempo y volumen del espacio aéreo. Actualmente, el prototipo 3-D CRM está siendo adaptado e integrado en la herramienta de monitorización de prestaciones de Aena (PERSEO) para complementar las bases de accidentes e incidentes ATM y mejorar la monitorización y proporcionar evidencias de los niveles de seguridad. ABSTRACT As with all forms of transport, the safety of air travel is of paramount importance. With the projected increases in European air traffic in the next decade and beyond, it is clear that the risk of accidents needs to be assessed and carefully monitored on a continuing basis. The present thesis is aimed at the development of a comprehensive collision risk model as a method of assessing the European en-route risk, due to all causes and across all dimensions within the airspace. The major constraint in developing appropriate monitoring methodologies and tools to assess the level of safety in en-route airspaces where controllers monitor air traffic by means of radar surveillance and provide aircraft with tactical instructions lies in the estimation of the operational risk. The operational risk estimate normally relies on incident reports provided by the air navigation service providers (ANSPs). This thesis proposes a new and innovative approach to assessing aircraft safety level based exclusively upon the process and analysis of radar tracks. The proposed methodology has been designed to complement the information collected in the accident and incident databases, thereby providing robust information on air traffic factors and safety metrics inferred from the in depth assessment of proximate events. The 3-D CRM methodology is implemented in a prototype tool in MATLAB © in order to automatically analyze recorded aircraft tracks and flight plan data from the Radar Data Processing systems (RDP) and identify and analyze all proximate events (conflicts, potential conflicts and potential collisions) within a time span and a given volume of airspace. Currently, the 3D-CRM prototype is been adapted and integrated in AENA’S Performance Monitoring Tool (PERSEO) to complement the information provided by the ATM accident and incident databases and to enhance monitoring and providing evidence of levels of safety.