Investigation of submicron volcanic aerosol-particles by photoelectron emission

Volcanic aerosol emissions have been studied for the first time by in situ photoelectric charging. Explorative studies on Mt Etna reveal large concentrations of particles below 1000 nm with the spectrum peaking in the size range of 100–1000 nm diameter. Although a large fraction of the particles is already charged upon emission, the net electrical charge carried by the aerosol turned out to be close to zero.Particles with high photoelectric yield vary greatly in their relative abundance and seem to occur mainly at active points of the volcano.     

震后快速钻探:历史、现状与未来

大地震发生后立即在断层带上进行钻探可以帮助我们获取更多的地震信息,尤其是可获得决定断层动态破裂的摩擦水平和强度,观测断层的愈合过程及可能触发余震的应力变化,并可获取控制破裂过程的重要的物理和化学属性。在2008年11月由国际大陆科学钻探项目（ICDP）和南加州地震中心（SCEC）联合在日本东京举办的为期3天的“震后快速钻探：历史、现状与未来”研讨会上,     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

SEISMIC VIBRATOR MODELLING1

The wavefield in, and at the surface of, a homogeneous, isotropic, perfectly elastic half-space, excited by a traction distribution at the surface of the medium is investigated. The emitted wavefield is a spatial convolution of the surface tractions and the spatial impulse response. The properties of the wavefield in the far-field of the medium are derived and it is shown that the far-field particle velocity is essentially equal to a weighted sum of the time derivative of the integrated surface tractions, that is, of the components of the ‘ground force’. The theory is valid for an arbitrary geometry and orientation of the surface tractions, and is independent of the boundary conditions at the surface of the medium. The surface tractions are related to a source that consists of a mass distribution with an arbitrary force distribution imposed upon it. A boundary condition is introduced that accounts for the mass load and the forces applied to it but neglects vibrations within the mass. The boundary condition follows from the equation of motion of the surface mass load. The theory is applied to the Vibroseis configuration, using a P-wave vibrator model with a uniformly distributed force imposed on top of the baseplate, and assuming that horizontal surface traction components are absent. The distribution of displacement and stress directly underneath the baseplate of a single vibrator and an array of vibrators is investigated. Three different boundary conditions are used: (1) assuming uniform pressure, (2) assuming uniform displacement, (3) using the equation of motion of the baseplate as a boundary condition. The calculations of the distribution of stress and displacement over the plate for different elastic media and several frequencies of operation show that only the results obtained with the mixed boundary condition agree with measurements made in the field. The accuracy of three different phase-feedback signals is compared using synthetic data. Baseplate velocity phase-feedback leads to huge deviations in the determination of the far-field wavelet; reaction mass acceleration phase-feedback looks stable but neglects the differentiating earth filter; and phase-feedback to a weighted sum of baseplate and reaction mass accelerations becomes unstable with increasing frequency. The instability can be overcome using measurements over the whole baseplate. The model can be extended to a lossy layered earth.     

Rock-weathering rates as functions of time

The scarcity of documented numerical relations between rock weathering and time has led to a common assumption that rates of weathering are linear. This assumption has been strengthened by studies that have calculated long-term average rates. However, little theoretical or empirical evidence exists to support linear rates for most chemical-weathering processes, with the exception of congruent dissolution processes. The few previous studies of rock-weathering rates that contain quantitative documentation of the relation between chemical weathering and time suggest that the rates of most weathering processes decrease with time. Recent studies of weathering rinds on basaltic and andesitic stones in glacial deposits in the western United States also clearly demonstrate that rock-weathering processes slow with time. Some weathering processes appear to conform to exponential functions of time, such as the square-root time function for hydration of volcanic glass, which conforms to the theoretical predictions of diffusion kinetics. However, weathering of mineralogically heterogeneous rocks involves complex physical and chemical processes that generally can be expressed only empirically, commonly by way of logarithmic time functions. Incongruent dissolution and other weathering processes produce residues, which are commonly used as measures of weathering. These residues appear to slow movement of water to unaltered material and impede chemical transport away from it. If weathering residues impede weathering processes then rates of weathering and rates of residue production are inversely proportional to some function of the residue thickness. This results in simple mathematical analogs for weathering that imply nonlinear time functions. The rate of weathering becomes constant only when an equilibrium thickness of the residue is reached. Because weathering residues are relatively stable chemically, and because physical removal of residues below the ground surface is slight, many weathering features require considerable time to reach constant rates of change. For weathering rinds on volcanic stones in the western United States, this time is at least 0.5 my.     

Absorption properties of urban/suburban aerosols in China

The broadband diffuse radiation method is improved to retrieve the aerosol refractive index imaginary part （AIP） and broadband （400-1000 nm mean） single scattering albedo （SSA）. In this method, four sets of SSA selection criteria are proposed for quality control. The method is used to retrieve AIP, SSA and absorptive optical thickness （AbOT） from routine hourly-exposed pyrheliometer and paranometer measurements over 11 sites （meteorological observatories） in China during 1998-2003. Apart from one suburban site （Ejin Qi）, the other urban sites are all located around big or medium cities. As shown in the retrieval results, annual mean SSA during 1998-2003 changes from 0.941 （Wuhan） to 0.849 （Lanzhou）, and AIP from 0.0054 to 0.0203. The 11-site average annual mean SSA and AIP are 0.898 and 0.0119, respectively. SSA during winter is smaller for most sites. There is an evidently positive correlation between SSA and aerosol optical thickness （AOT） for all sites. There is also a positive correlation between SSA and relative humidity for most sites, but a negative correlation for a few sites, such as Kashi and ǚrǚmqi in Northwest China.     

Studying an interaction between tropical cyclones and jet streams using the numerical simulation data

Discussed are the results of studying an evolution of tropical cyclones (TCs) in the Pacific Ocean using the data of computation of ETA and WRF NMM mesoscale numerical atmospheric models. Computed are the trajectories of TCs and the fields of meteorological variables in the typhoons, of the wind speed and kinetic energy in the subtropical jet stream during the development of Parma, Melor, and Lupit typhoons. Carried out are the analysis and comparison of computed fields of pressure, wind, kinetic energy, and trajectories of TCs obtained using these models and their comparison with the actual fields. It is demonstrated that both models computed rather well the complex trajectories and the fields of wind and kinetic energy varying in the course of the interaction. Proposed is an explanation of processes taking place during the interaction between the vortices and the subtropical jet stream and the polar front.     

Climatic change and Chinese population growth dynamics over the last millennium

The climate–population relationship has long been conceived. Although the topic has been repeatedly investigated, most of the related works are Eurocentric or qualitative. Consequently, the relationship between climate and population remains ambiguous. In this study, fine-grained temperature reconstructions and historical population data sets have been employed to statistically test a hypothesized relationship between temperature change and population growth (i.e., cooling associated with below average population growth) in China over the past millennium. The important results were: (1) Long-term temperature change significantly determined the population growth dynamics of China. However, spatial variation existed, whilst population growth in Central China was shown to be responsive to both long- and short-term temperature changes; in marginal areas, population growth was only sensitive to short-term temperature fluctuations. (2) Temporally, the temperature–population relationship was obscured in some periods, which was attributable to the factors of drought and social buffers. In summary, a temperature–population relationship was mediated by geographic factors, the aridity threshold, and social factors. Given the upcoming threat posed by climate change to human societies, this study seeks to improve our knowledge and understanding of the climate–society relationship.     

Use of atmospheric radiation measurement program data from Barrow, Alaska, for evaluation and development of snow-albedo parameterizations

Summary Snow albedo is determined from the ratio of out-going to incoming solar radiation using three years of broadband shortwave radiometer data obtained from the Barrow, Alaska, Atmospheric Radiation Measurement (ARM) site. These data are used for the evaluation of various types of snow-albedo parameterizations applied in numerical weather prediction or climate models. These snow-albedo parameterizations are based on environmental conditions (e.g., air or snow temperature), snow related characteristics (e.g., snow depth, snow age), or combinations of both. The ARM data proved to be well suited for snow-albedo evaluation purposes for a low-precipitation tundra environment. The evaluation confirms that snow-age dependent parameterizations of snow albedo work well during snowmelt, while parameterizations considering meteorological conditions often perform better during snow accumulation. Current difficulties in parameterizing snow albedo occur for long episodes of snow-event free conditions and episodes with a high frequency of snow events or strong snowfall. In a further step, the first two years of the ARM albedo dataset is used to develop a snow-albedo parameterization, and the third year’s data serves for its evaluation. This parameterization considers snow depth, wind speed, and air temperature which are found to be significant parameters for snow-albedo modeling under various conditions. Comparison of all evaluated snow-albedo parameterizations with this new parameterization shows improved snow-albedo prediction. Correspondence: Nicole M?lders, Geophysical Institute and College of Natural Science and Mathematics, University of Alaska Fairbanks, 903 Koyukuk Drive, P.O. Box 757320, Fairbanks, AK 99775-7320, USA     

Nonlinear interactions of spherical Rossby modes

Weakly nonlinear triad interactions between spherical Rossby harmonics are studied. Each of the harmonics has the form AP_n^m(sin θ)exp[i(mλ − σt)], where A is an amplitude and P_n^m is the associated Legendre function. Equations for the amplitudes are derived and resonance conditions are analysed. The resonance conditions differ substantially from the usual resonance conditions on a β-plane and include a Diophantine equation and a few inequalities for the integer wavenumbers n and m of the interacting modes. Particular analytical series of solutions to the resonance conditions are constructed, which show that modes with arbitrary large wavenumbers can participate in the interactions. A numerical study of the resonance conditions reveals that no more than 21% of the Rossby harmonics can participate in the triad interactions and that chains of the interacting triads soon break. Thus precise interactions (for which the resonance conditions hold exactly) do not result in any significant redistribution of energy over the spectrum. On the other hand, approximate interactions (for which the resonance conditions hold approximately) generate an intensive energy redistribution among short Rossby modes with typical scales much smaller than the Earth's radius. Thus the energy cascade is concentrated mainly in the short-wave part of the spectrum and is very weak in the large-scale domain. The efficiency of the triad interaction of Rossby modes with scales much smaller than the Earth's radius depends strongly on the existence of the so-called interaction latitude at which the local wave-vectors and frequencies of the interacting modes satisfy resonance conditions for plane Rossby waves on the β-plane approximating the neighbourhood of the latitude. If the interaction latitude exists, the interaction is intensive; in the opposite case, it is weak.