Aerosol microphysics simulations of the Mt.~Pinatubo eruption with the UM-UKCA composition-climate model

We use a stratosphere–troposphere composition–climate model with interactive sulfur chemistry and aerosol microphysics, to investigate the effect of the 1991 Mount Pinatubo eruption on stratospheric aerosol properties. Satellite measurements indicate that shortly after the eruption, between 14 and 23 Tg of SO2 (7 to 11.5 Tg of sulfur) was present in the tropical stratosphere. Best estimates of the peak global stratospheric aerosol burden are in the range 19 to 26 Tg, or 3.7 to 6.7 Tg of sulfur assuming a composition of between 59 and 77 % H2SO4. In light of this large uncertainty range, we performed two main simulations with 10 and 20 Tg of SO2 injected into the tropical lower stratosphere. Simulated stratospheric aerosol properties through the 1991 to 1995 period are compared against a range of available satellite and in situ measurements. Stratospheric aerosol optical depth (sAOD) and effective radius from both simulations show good qualitative agreement with the observations, with the timing of peak sAOD and decay timescale matching well with the observations in the tropics and mid-latitudes. However, injecting 20 Tg gives a factor of 2 too high stratospheric aerosol mass burden compared to the satellite data, with consequent strong high biases in simulated sAOD and surface area density, with the 10 Tg injection in much better agreement. Our model cannot explain the large fraction of the injected sulfur that the satellite-derived SO2 and aerosol burdens indicate was removed within the first few months after the eruption. We suggest that either there is an additional alternative loss pathway for the SO2 not included in our model (e.g. via accommodation into ash or ice in the volcanic cloud) or that a larger proportion of the injected sulfur was removed via cross-tropopause transport than in our simulations. We also critically evaluate the simulated evolution of the particle size distribution, comparing in detail to balloon-borne optical particle counter (OPC) measurements from Laramie, Wyoming, USA (41° N). Overall, the model captures remarkably well the complex variations in particle concentration profiles across the different OPC size channels. However, for the 19 to 27 km injection height-range used here, both runs have a modest high bias in the lowermost stratosphere for the finest particles (radii less than 250 nm), and the decay timescale is longer in the model for these particles, with a much later return to background conditions. Also, whereas the 10 Tg run compared best to the satellite measurements, a significant low bias is apparent in the coarser size channels in the volcanically perturbed lower stratosphere. Overall, our results suggest that, with appropriate calibration, aerosol microphysics models are capable of capturing the observed variation in particle size distribution in the stratosphere across both volcanically perturbed and quiescent conditions. Furthermore, additional sensitivity simulations suggest that predictions with the models are robust to uncertainties in sub-grid particle formation and nucleation rates in the stratosphere.

Validation of MIPAS IMK/IAA temperature, water vapor, and ozone profiles with MOHAVE-2009 campaign measurements

MIPAS observations of temperature, water vapor, and ozone in October 2009 as derived with the scientific level-2 processor run by Karlsruhe Institute of Technology (KIT), Institute for Meteorology and Climate Research (IMK) and CSIC, Instituto de Astrofísica de Andalucía (IAA) and retrieved from version 4.67 level-1b data have been compared to co-located field campaign observations obtained during the MOHAVE-2009 campaign at the Table Mountain Facility near Pasadena, California in October 2009. The MIPAS measurements were validated regarding any potential biases of the profiles, and with respect to their precision estimates. The MOHAVE-2009 measurement campaign provided measurements of atmospheric profiles of temperature, water vapor/relative humidity, and ozone from the ground to the mesosphere by a suite of instruments including radiosondes, ozonesondes, frost point hygrometers, lidars, microwave radiometers and Fourier transform infra-red (FTIR) spectrometers. For MIPAS temperatures (version V4O_T_204), no significant bias was detected in the middle stratosphere; between 22 km and the tropopause MIPAS temperatures were found to be biased low by up to 2 K, while below the tropopause, they were found to be too high by the same amount. These findings confirm earlier comparisons of MIPAS temperatures to ECMWF data which revealed similar differences. Above 12 km up to 45 km, MIPAS water vapor (version V4O_H2O_203) is well within 10% of the data of all correlative instruments. The well-known dry bias of MIPAS water vapor above 50 km due to neglect of non-LTE effects in the current retrievals has been confirmed. Some instruments indicate that MIPAS water vapor might be biased high by 20 to 40% around 10 km (or 5 km below the tropopause), but a consistent picture from all comparisons could not be derived. MIPAS ozone (version V4O_O3_202) has a high bias of up to +0.9 ppmv around 37 km which is due to a non-identified continuum like radiance contribution. No further significant biases have been detected. Cross-comparison to co-located observations of other satellite instruments (Aura/MLS, ACE-FTS, AIRS) is provided as well.

Interference of parallel cylindrical Langmuir probes

Current to a cylindrical probe of arbitrary cross section is discussed. Previous results for circular cylinders at the high bias and moderate radius R of interest for electrodynamic bare tethers, for which space charge may be ignored over a large neighborhood of the probe, depend in separate ways on both R and perimeter p. These results are extended to a general convex cross section by introducing certain equivalent radius Req. For any concave cross section, results use a proper equivalent perimeter peq , in addition to Req. Finally, for the joint cross section of separate parallel probes, certain effective perimeter peff replaces peq. Rules to determine Req. peq. and peff are used to discuss collection interference among two or more parallel cylinders when brought from far away to contact

The orbital-motion-limited regime of cylindrical Langmuir probes

An asymptotic analysis of electron collection at high bias Fp serves to determine the domain of validity of the orbital-motion-limited regime of cylindrical Langmuir probes, which is basic for the workings of conductive bare tethers. The radius of a wire collecting OML current in an unmagnetized plasma at rest cannot exceed a value, Rmax , which is found to exhibit a minimum as a function of Fp ; atFp values of interest, Rmax is already increasing and is larger than the electron Debye length lDe . The breakdown of the regime relates to conditions far fromthe probe, at electron energies comparable to the ion thermal energy, kTi ; Rmax is found to increase with Ti . It is also found that ~1! the maximumwidth of a thin tape, if used instead of a wire, is 4Rmax ; ~2! the electron thermal gyroradius must be larger than both R and lDe for magnetic effects to be negligible; and ~3! conditions applying to the tether case are such that trapped-orbit effects are negligible.

Bare-tether sheath and current: comparison of asymptotic theory and kinetic simulations in stationary plasma

Analytical expressions for current to a cylindrical Langmuir probe at rest in unmagnetized plasma are compared with results from both steady-state Vlasov and particle-in-cell simulations. Probe bias potentials that are much greater than plasma temperature (assumed equal for ions and electrons), as of interest for bare conductive tethers, are considered. At a very high bias, both the electric potential and the attracted-species density exhibit complex radial profiles; in particular, the density exhibits a minimum well within the plasma sheath and a maximum closer to the probe. Excellent agreement is found between analytical and numerical results for values of the probe radiusR close to the maximum radius Rmax for orbital-motion-limited (OML) collection at a particular bias in the following number of profile features: the values and positions of density minimum and maximum, position of sheath boundary, and value of a radius characterizing the no-space-charge behavior of a potential near the high-bias probe. Good agreement between the theory and simulations is also found for parametric laws jointly covering the following three characteristic R ranges: sheath radius versus probe radius and bias for Rmax; density minimum versus probe bias for Rmax; and (weakly bias-dependent) current drop below the OML value versus the probe radius for R > Rmax.

A theory of electron collection by a grid sphere

Use of a spherical grid as electron collector at the anodic end of a tether, as recently proposed, is considered. The standard analysis of space-charge limited current to a solid sphere (with neither magnetic nor plasma-motion effects), which has been shown to best fit TSS1R in-orbit results at very high bias, is used to determine effects from grid transparency on current collected; the analysis is first reformulated in the formalism recently introduced in the two-dimensional analysis of bare-tethers. A discussion of the electric potential created by a spherical grid in vacuum is then carried out; it is shown that each grid-wire collects current well below its maximum OML current, the effective grid transparency being close to its optical value. Formulae for the current to a spherical grid, showing the effects of grid transparency, is determined. A fully consistent analysis of electric potential and electron density, outside and inside the grid, is completed.

Low Work-Function Thermionic Emission and Orbital-Motion-Limited Ion Collection at Bare-Tether Cathodic Contact

Una amarra electrodinámica (electrodynamic tether) opera sobre principios electromagnéticos intercambiando momento con la magnetosfera planetaria e interactuando con su ionosfera. Es un subsistema pasivo fiable para desorbitar etapas de cohetes agotadas y satélites al final de su misión, mitigando el crecimiento de la basura espacial. Una amarra sin aislamiento captura electrones del plasma ambiente a lo largo de su segmento polarizado positivamente, el cual puede alcanzar varios kilómetros de longitud, mientras que emite electrones de vuelta al plasma mediante un contactor de plasma activo de baja impedancia en su extremo catódico, tal como un cátodo hueco (hollow cathode). En ausencia de un contactor catódico activo, la corriente que circula por una amarra desnuda en órbita es nula en ambos extremos de la amarra y se dice que ésta está flotando eléctricamente. Para emisión termoiónica despreciable y captura de corriente en condiciones limitadas por movimiento orbital (orbital-motion-limited, OML), el cociente entre las longitudes de los segmentos anódico y catódico es muy pequeño debido a la disparidad de masas entre iones y electrones. Tal modo de operación resulta en una corriente media y fuerza de Lorentz bajas en la amarra, la cual es poco eficiente como dispositivo para desorbitar. El electride C12A7 : e−, que podría presentar una función de trabajo (work function) tan baja como W = 0.6 eV y un comportamiento estable a temperaturas relativamente altas, ha sido propuesto como recubrimiento para amarras desnudas. La emisión termoiónica a lo largo de un segmento así recubierto y bajo el calentamiento de la operación espacial, puede ser más eficiente que la captura iónica. En el modo más simple de fuerza de frenado, podría eliminar la necesidad de un contactor catódico activo y su correspondientes requisitos de alimentación de gas y subsistema de potencia, lo que resultaría en un sistema real de amarra “sin combustible”. Con este recubrimiento de bajo W, cada segmento elemental del segmento catódico de una amarra desnuda de kilómetros de longitud emitiría corriente como si fuese parte de una sonda cilíndrica, caliente y uniformemente polarizada al potencial local de la amarra. La operación es similar a la de una sonda de Langmuir 2D tanto en los segmentos catódico como anódico. Sin embargo, en presencia de emisión, los electrones emitidos resultan en carga espacial (space charge) negativa, la cual reduce el campo eléctrico que los acelera hacia fuera, o incluso puede desacelerarlos y hacerlos volver a la sonda. Se forma una doble vainas (double sheath) estable con electrones emitidos desde la sonda e iones provenientes del plasma ambiente. La densidad de corriente termoiónica, variando a lo largo del segmento catódico, podría seguir dos leyes distintas bajo diferentes condiciones: (i) la ley de corriente limitada por la carga espacial (space-charge-limited, SCL) o (ii) la ley de Richardson-Dushman (RDS). Se presenta un estudio preliminar sobre la corriente SCL frente a una sonda emisora usando la teoría de vainas (sheath) formada por la captura iónica en condiciones OML, y la corriente electrónica SCL entre los electrodos cilíndricos según Langmuir. El modelo, que incluye efectos óhmicos y el efecto de transición de emisión SCL a emisión RDS, proporciona los perfiles de corriente y potencial a lo largo de la longitud completa de la amarra. El análisis muestra que en el modo más simple de fuerza de frenado, bajo condiciones orbitales y de amarras típicas, la emisión termoiónica proporciona un contacto catódico eficiente y resulta en una sección catódica pequeña. En el análisis anterior, tanto la transición de emisión SCL a RD como la propia ley de emisión SCL consiste en un modelo muy simplificado. Por ello, a continuación se ha estudiado con detalle la solución de vaina estacionaria de una sonda con emisión termoiónica polarizada negativamente respecto a un plasma isotrópico, no colisional y sin campo magnético. La existencia de posibles partículas atrapadas ha sido ignorada y el estudio incluye tanto un estudio semi-analítico mediante técnica asintóticas como soluciones numéricas completas del problema. Bajo las tres condiciones (i) alto potencial, (ii) R = Rmax para la validez de la captura iónica OML, y (iii) potencial monotónico, se desarrolla un análisis asintótico auto-consistente para la estructura de plasma compleja que contiene las tres especies de cargas (electrones e iones del plasma, electrones emitidos), y cuatro regiones espaciales distintas, utilizando teorías de movimiento orbital y modelos cinéticos de las especies. Aunque los electrones emitidos presentan carga espacial despreciable muy lejos de la sonda, su efecto no se puede despreciar en el análisis global de la estructura de la vaina y de dos capas finas entre la vaina y la región cuasi-neutra. El análisis proporciona las condiciones paramétricas para que la corriente sea SCL. También muestra que la emisión termoiónica aumenta el radio máximo de la sonda para operar dentro del régimen OML y que la emisión de electrones es mucho más eficiente que la captura iónica para el segmento catódico de la amarra. En el código numérico, los movimientos orbitales de las tres especies son modelados para potenciales tanto monotónico como no-monotónico, y sonda de radio R arbitrario (dentro o más allá del régimen de OML para la captura iónica). Aprovechando la existencia de dos invariante, el sistema de ecuaciones Poisson-Vlasov se escribe como una ecuación integro-diferencial, la cual se discretiza mediante un método de diferencias finitas. El sistema de ecuaciones algebraicas no lineal resultante se ha resuelto de con un método Newton-Raphson paralelizado. Los resultados, comparados satisfactoriamente con el análisis analítico, proporcionan la emisión de corriente y la estructura del plasma y del potencial electrostático. ABSTRACT An electrodynamic tether operates on electromagnetic principles and exchanges momentum through the planetary magnetosphere, by continuously interacting with the ionosphere. It is a reliable passive subsystem to deorbit spent rocket stages and satellites at its end of mission, mitigating the growth of orbital debris. A tether left bare of insulation collects electrons by its own uninsulated and positively biased segment with kilometer range, while electrons are emitted by a low-impedance active device at the cathodic end, such as a hollow cathode, to emit the full electron current. In the absence of an active cathodic device, the current flowing along an orbiting bare tether vanishes at both ends and the tether is said to be electrically floating. For negligible thermionic emission and orbital-motion-limited (OML) collection throughout the entire tether (electron/ion collection at anodic/cathodic segment, respectively), the anodic-to-cathodic length ratio is very small due to ions being much heavier, which results in low average current and Lorentz drag. The electride C12A7 : e−, which might present a possible work function as low as W = 0.6 eV and moderately high temperature stability, has been proposed as coating for floating bare tethers. Thermionic emission along a thus coated cathodic segment, under heating in space operation, can be more efficient than ion collection and, in the simplest drag mode, may eliminate the need for an active cathodic device and its corresponding gas-feed requirements and power subsystem, which would result in a truly “propellant-less” tether system. With this low-W coating, each elemental segment on the cathodic segment of a kilometers-long floating bare-tether would emit current as if it were part of a hot cylindrical probe uniformly polarized at the local tether bias, under 2D probe conditions that are also applied to the anodic-segment analysis. In the presence of emission, emitted electrons result in negative space charge, which decreases the electric field that accelerates them outwards, or even reverses it, decelerating electrons near the emitting probe. A double sheath would be established with electrons being emitted from the probe and ions coming from the ambient plasma. The thermionic current density, varying along the cathodic segment, might follow two distinct laws under different con ditions: i) space-charge-limited (SCL) emission or ii) full Richardson-Dushman (RDS) emission. A preliminary study on the SCL current in front of an emissive probe is presented using the orbital-motion-limited (OML) ion-collection sheath and Langmuir’s SCL electron current between cylindrical electrodes. A detailed calculation of current and bias profiles along the entire tether length is carried out with ohmic effects considered and the transition from SCL to full RDS emission is included. Analysis shows that in the simplest drag mode, under typical orbital and tether conditions, thermionic emission provides efficient cathodic contact and leads to a short cathodic section. In the previous analysis, both the transition between SCL and RDS emission and the current law for SCL condition have used a very simple model. To continue, considering an isotropic, unmagnetized, colissionless plasma and a stationary sheath, the probe-plasma contact is studied in detail for a negatively biased probe with thermionic emission. The possible trapped particles are ignored and this study includes both semianalytical solutions using asymptotic analysis and complete numerical solutions. Under conditions of i) high bias, ii) R = Rmax for ion OML collection validity, and iii) monotonic potential, a self-consistent asymptotic analysis is carried out for the complex plasma structure involving all three charge species (plasma electrons and ions, and emitted electrons) and four distinct spatial regions using orbital motion theories and kinetic modeling of the species. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between the sheath and the quasineutral region. The parametric conditions for the current to be space-chargelimited are obtained. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers. In the numerical code, the orbital motions of all three species are modeled for both monotonic and non-monotonic potential, and for any probe radius R (within or beyond OML regime for ion collection). Taking advantage of two constants of motion (energy and angular momentum), the Poisson-Vlasov equation is described by an integro differential equation, which is discretized using finite difference method. The non-linear algebraic equations are solved using a parallel implementation of the Newton-Raphson method. The results, which show good agreement with the analytical results, provide the results for thermionic current, the sheath structure, and the electrostatic potential.

Low work-function thermionic emission and orbital-motion-limited ion collection at bare-tether cathodic contact

With a thin coating of low-work-function material, thermionic emission in the cathodic segment of bare tethers might be much greater than orbital-motion-limited (OML) ion collection current. The space charge of the emitted electrons decreases the electric field that accelerates them outwards, and could even reverse it for high enough emission, producing a potential hollow. In this work, at the conditions of high bias and relatively low emission that make the potential monotonic, an asymptotic analysis is carried out, extending the OML ion-collection analysis to investigate the probe response due to electrons emitted by the negatively biased cylindrical probe. At given emission, the space charge effect from emitted electrons increases with decreasing magnitude of negative probe bias. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between sheath and the quasineutral region. The space-charge-limited condition is located. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers.

Electrical characterization of the rectifying contact between aluminium and electrodeposited poly(3-methylthiophene)

A detailed investigation both of the DC and of the AC electrical properties of the Schottky barrier formed between aluminium and electrodeposited poly(3-methylthiophene) is reported. The devices show rectification ratios up to 2 x 10(4) which can be increased further after post-metal annealing. The reverse characteristics of the devices follow predictions based on the image-force lowering of the Schottky barrier, from which the doping density can be estimated, As the forward voltage increases, the device current is limited by the bulk resistance of the polymer with some evidence for injection limitation at the gold counter-electrode at high bias. In the bulk-limited regime, the device current is thermally activated near room temperature with an activation energy in the range 0.2-0.3 eV. Below about 150 K the device current is almost independent of temperature. Capacitance-voltage plots obtained at frequencies well below the device relaxation frequency indicate the presence of two distinct acceptor states, A set of shallow acceptor states are active in forward bias and are believed to determine the bulk conductivity of the polymer. A set of deeper accepters are active only for very small forward voltages and for all reverse voltages, namely when band banding causes the Fermi energy to cross these states. The density of these deeper states is approximately an order of magnitude greater than that of the shallow states. Evidence is presented also for the influence of fabrication conditions on the formation of an insulating interfacial layer at the rectifying interface. The presence of such a layer leads to inversion at the polymer surface and a modification of the I-V characteristics.

Electrical characterization of the rectifying contact between aluminium and electrodeposited poly(3-methylthiophene)

A detailed investigation both of the DC and of the AC electrical properties of the Schottky barrier formed between aluminium and electrodeposited poly(3-methylthiophene) is reported. The devices show rectification ratios up to 2 x 10(4) which can be increased further after post-metal annealing. The reverse characteristics of the devices follow predictions based on the image-force lowering of the Schottky barrier, from which the doping density can be estimated, As the forward voltage increases, the device current is limited by the bulk resistance of the polymer with some evidence for injection limitation at the gold counter-electrode at high bias. In the bulk-limited regime, the device current is thermally activated near room temperature with an activation energy in the range 0.2-0.3 eV. Below about 150 K the device current is almost independent of temperature. Capacitance-voltage plots obtained at frequencies well below the device relaxation frequency indicate the presence of two distinct acceptor states, A set of shallow acceptor states are active in forward bias and are believed to determine the bulk conductivity of the polymer. A set of deeper accepters are active only for very small forward voltages and for all reverse voltages, namely when band banding causes the Fermi energy to cross these states. The density of these deeper states is approximately an order of magnitude greater than that of the shallow states. Evidence is presented also for the influence of fabrication conditions on the formation of an insulating interfacial layer at the rectifying interface. The presence of such a layer leads to inversion at the polymer surface and a modification of the I-V characteristics.

Isospin symmetry breaking at high spins in the mirror pair Se-67 and As-67

Recent experimental data have revealed large mirror energy differences (MED) between high-spin states in the mirror nuclei Se-67 and As-67, the heaviest pair where MED have been determined so far. The MED are generally attributed to the isospin symmetry breaking caused by the Coulomb force and by the isospin-nonconserving part of the nucleon-nucleon residual interaction. The different contributions of the various terms have been extensively studied in the fp shell. By employing large-scale shell-model calculations, we show that the inclusion of the g(9/2) orbit causes interference between the electromagnetic spin-orbit and the Coulomb monopole radial terms at high spin. The large MED are attributed to the aligned proton pair excitations from the p(3/2) and f(5/2) orbits to the g(9/2) orbit. The relation of the MED to deformation is discussed.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Experimental and numerical studies on rock breaking with TBM tools under high stress confinement

The understanding of rock breaking and chipping due to the TBM cutter disks mechanism in deep tunnels is considered in this paper. The interest stems from the use of TBMs for the excavation of long Trans-Alpine tunnels. Some tests that simulate the disk cutter action at the tunnel face by means of an indenter, acting on a rock specimen are proposed. The rock specimen is confined through a flat-jack and a confinement-free area on one side of the specimen simulates the formation of a groove near the indenter, like it occurs in TBM excavation conditions. Results show a limited influence of the confinement stress versus the thrust increment required for breaking the rock between the indenter and the free side of the specimen. Numerical modelling of the cutter disk action on confined material has also been carried out in order to investigate further aspects of the fracture initiation. Also in this case the importance of the relative position between disk cutter and groove is pointed out. © 2006 Springer-Verlag.

Safe operation of capacitance meters using high applied-bias voltage /

Includes bibliographical references (p. 49).