Simulating the effects of the 1991 Mount Pinatubo volcanic eruption using the ARPEGE atmosphere general circulation model

The climate changes that occured following the volcanic eruption of Mount Pinatubo in the Phillippines on 15 June 1991 have been simulated using the ARPEGE atmosphere general circulation model （AGCM）. The model was forced by a reconstructed spatial-time distribution of stratospheric aerosols intended for use in long climate simulations. Four statistical ensembles of the AGCM simulations with and without volcanic aerosols over a period of 5 years following the eruption have been made, and the calculated fields have been compared to available observations. The model is able to reproduce some of the observed features after the eruption, such as the winter warming pattern that was observed over the Northern Hemisphere （NH） during the following winters. This pattern was caused by an enhanced Equator-to-pole temperature gradient in the stratosphere that developed due to aerosol heating of the tropics. This in turn led to a strengthening of the polar vortex, which tends to modulate the planetary wave field in such a way that an anomalously positive Arctic Oscillation pattern is produced in the troposphere and at the surface, favouring warm conditions over the NH. During the summer, the model produced a more uniform cooling over the NH.     

耦合模式长期积分中东亚夏季风与ENSO联系的不稳定性

通过对挪威卑尔根全球大气-海洋-海冰耦合模式300a控制积分结果进行交叉子波分析,揭示了东亚夏季风（EASM）与同期Nio3区（90°W～150°W,5°S～5°N）海洋表面温度异常的相关关系在长期变化中是不稳定的,呈现出明显的阶段性特征.气候要素场在二者联系的紧密（HCP）和微弱（LCP）时期差别显著,在HCP时期,西北太平洋对流层低层出现一对耦合的异常气旋和反气旋性环流系统;东南亚地区对流层低层表现为强东风异常,风速的年际变率加大;热带西太平洋对流层温度和位势高度场的年际变率普遍加强.此外,中国夏季降水与同期Nio3区海洋表面温度异常的相关关系在上述两种时期也存在较大差别.     

On the potential for CO2 mineral storage in continental flood basalts - PHREEQC batch- and 1D diffusion-reaction simulations

ABSTRACT: Continental flood basalts (CFB) are considered as potential CO2 storage sites because of their high reactivity and abundant divalent metal ions that can potentially trap carbon for geological timescales. Moreover, laterally extensive CFB are found in many place in the world within reasonable distances from major CO2 point emission sources.Based on the mineral and glass composition of the Columbia River Basalt (CRB) we estimated the potential of CFB to store CO2 in secondary carbonates. We simulated the system using kinetic dependent dissolution of primary basalt-minerals (pyroxene, feldspar and glass) and the local equilibrium assumption for secondary phases (weathering products). The simulations were divided into closed-system batch simulations at a constant CO2 pressure of 100?bar with sensitivity studies of temperature and reactive surface area, an evaluation of the reactivity of H2O in scCO2, and finally 1D reactive diffusion simulations giving reactivity at CO2 pressures varying from 0 to 100?bar.Although the uncertainty in reactive surface area and corresponding reaction rates are large, we have estimated the potential for CO2 mineral storage and identified factors that control the maximum extent of carbonation. The simulations showed that formation of carbonates from basalt at 40?C may be limited to the formation of siderite and possibly FeMg carbonates. Calcium was largely consumed by zeolite and oxide instead of forming carbonates. At higher temperatures (60 - 100?C), magnesite is suggested to form together with siderite and ankerite. The maximum potential of CO2 stored as solid carbonates, if CO2 is supplied to the reactions unlimited, is shown to depend on the availability of pore space as the hydration and carbonation reactions increase the solid volume and clog the pore space. For systems such as in the scCO2 phase with limited amount of water, the total carbonation potential is limited by the amount of water present for hydration of basalt.     

Variation in mercury content of Guillemot feathers over 150 years

The feathers of 277 Black, Common and Brünnich's Guillemot from the Baltic, Kattegat, Faroe Islands and Greenland were analysed for mercury. The levels were found to be higher in the Baltic and the Kattegat compared to the Faroe Islands and Greenland. While the levels were almost constant for the last two areas a substantial increase during this century was found for the Baltic and the Kattegat. In Common Guillemots from the Baltic a decrease was indicated after 1969. In general the mercury levels were higher for Black Guillemots living close to the coast compared to the Uria sp. living off-shore. The inhomogeniety of the mercury load near the shore was indicated by a large individual variation within a colony of Black Guillemots. This variation was reinforced by differences in migration and feeding habits. It has been shown by analysis of individuals recaptured during successive years that the year to year variation of one individual is considerably smaller.It is concluded that chronological series based on feathers from sea birds and museum collections may contribute to the elucidation of the long term trend of mercury pollution at sea.     

The response of the North Pacific Decadal Variability to strong tropical volcanic eruptions

In this study, the effects of volcanic forcing on North Pacific climate variability, on interannual to decadal time scales, are examined using climate model simulations covering the last 600?years. The model used is the Bergen Climate Model, a fully coupled atmosphere–ocean general circulation model. It is found that natural external forcings, such as tropical strong volcanic eruptions (SVEs) and variations in total solar irradiance, play an important role in regulating North Pacific Decadal Variability (NPDV). In response to tropical SVEs the lower stratospheric pole–to–equator temperature gradient is enhanced. The North polar vortex is strengthened, which forces a significant positive Arctic Oscillation. At the same time, dipole zonal wind anomalies associated with strong polar vortex propagate downward from the lower stratosphere. Through positive feedbacks in the troposphere, the surface westerly winds across the central North Pacific are significantly weakened, and positive sea level pressure anomalies are formed in the North Pacific. This anomalous surface circulation results in changes in the net heat fluxes and the oceanic advection across the North Pacific. As a result of this, warm water converges in the subtropical western North Pacific, where the surface waters in addition are heated by significantly reduced latent and sensible heat fluxes from the ocean. In the eastern and high–latitude North Pacific the ocean loses more heat, and large–scale decreases in sea surface temperatures are found. The overall response of this chain of events is that the North Pacific enters a negative phase of the Pacific decadal oscillation (PDO), and this negative phase of the PDO is maintained for several years. It is thus concluded that the volcanic forcing plays a key role in the phasing of the PDO. The model results furthermore highlight the important role of troposphere–stratosphere coupling, tropical–extratropical teleconnections and extratropical ocean–atmosphere interactions for describing NPDV.     

Discontinuous Galerkin methods for advective transport in single-continuum models of fractured media

Accurate simulation of flow and transport processes in fractured rocks requires that flow in fractures and shear zones to be coupled with flow in the porous rock matrix. To this end, we will herein consider a single-continuum approach in which both fractures and the porous rock are represented as volumetric objects, i.e., as cells in an unstructured triangular grid with a permeability and a porosity value associated with each cell. Hence, from a numerical point of view, there is no distinction between flow in the fractures and the rock matrix. This enables modelling of realistic cases with very complex structures. To compute single-phase advective transport in such a model, we propose to use a family of higher-order discontinuous Galerkin methods. Single-phase transport equations are hyperbolic and have an inherent causality in the sense that information propagates along streamlines. This causality is preserved in our discontinuous Galerkin discretization. We can therefore use a simple topological sort of the graph of discrete fluxes to reorder the degrees-of-freedom such that the discretized linear system gets a lower block-triangular form, from which the solution can be computed very efficiently using a single-pass forward block substitution. The accuracy and utility of the resulting transport solver is illustrated through several numerical experiments.     

1000 m long gas blow-out pipes

This study presents seismic observation of pipe anomalies from offshore Nigeria, outcrops of blow-out pipes from Rhodes, Greece, and geophysical modelling of an acoustic pipe. The studies give insight into how pipes form, their internal structure, the seismic image and geophysical artefacts related to the pipes. Over one hundred seafloor craters, 100 m-700 m wide and up to 30 m deep, have been observed on the seafloor offshore Nigeria. They are underlain by interpreted cones and seismic pipe anomalies that can be traced down to reservoir zones at 1000 m-1300 m below the seafloor. The seismic pipe anomalies are 50 m-150 m wide and almost vertical. They are interpreted as up-scaled pipes found in outcrops on Rhodes, Greece. The outcrops show pipe-related structures at three levels. Lowest, the reservoir rock contains metre-sized cavities which are filled with a mixture of clay derived from the overlying cap rock. In the middle, several circular to oval structures in plane view of pipes are observed in the lowest part of the cap rock. Highest, 15 m into the clay cap rock, strongly sheared country rock forms circular structures with a core of structureless clay. Based on outcrop observation on Rhodes we constructed an acoustic model of a 50 m wide and 1000 m long pipe. Seismic modelling proves that such pipes would be expressed in seismic data, that they are similar to the seismic pipe anomalies offshore Nigeria but this study also revealed that prominent intra-pipe reflections are artefacts. A formation model for the pipes is suggested: High fluid overpressure in the reservoir generated hydro fractures from the reservoir to seafloor where a mixture of gas and fluid flowed at high speed to form pipes, cones and seafloor craters. After hours to weeks of gas and fluid flow through the pipe the pore pressure in the reservoir dropped and the blow-out terminated. Muddy slurry fell back and plugged the cavity in the reservoir and the pipe.     

A2111: a z=0.23 Butcher–Oemler cluster with a non-isothermal atmosphere and normal metallicity

We report results from an X-ray study of the Abell 2111 galaxy cluster using the Advanced Satellite for Astrophysics and Cosmology ( ASCA ) and the ROSAT Position Sensitive Proportional Counter (PSPC). By correcting for the energy-dependent point-spread function of the ASCA instruments, we have examined the temperature structure of the cluster. The core of the cluster within 3 arcmin is found to have a temperature of 6.46±0.87 keV, significantly higher than the value of 3.10±1.19 keV in the surrounding region of r =3–6 arcmin. This radially decreasing temperature structure can be parametrized by a polytropic index of γ ≃ 1.45. The X-ray morphology of the cluster appears elongated and clumpy on scales ≤1 arcmin. These results, together with earlier ROSAT and optical studies which revealed that the X-ray centroid and ellipticity of A2111 shift with spatial scale, are consistent with the hypothesis that the cluster is a dynamically young system. Most likely, the cluster has recently undergone a merger, which may also be responsible for the high fraction of blue galaxies observed in the cluster. Alternatively, the temperature structure may also be caused by the gravitational potential of the cluster. We have also measured the emission-weighted abundance of the X-ray-emitting intracluster medium as 0.25±0.14 solar. This value is similar to those of nearby clusters that do not show a large blue galaxy fraction, indicating that star formation in disc galaxies and subsequent loss to the medium do not drastically alter the average abundance of a cluster. This is consistent with recent results, which indicate that cluster abundances have remained constant since at least z ∼0.3.