Natural siliceous dust and human health： Pathways, toxicity, and impact

Fine atmospheric dust includes mineral particles and aggregates, fibrous minerals and fibrous organic material. Generation, dislodgement and transport （deflation） of natural dust with the finer （〈4 microns） components suspended as silt-size aggregates, is widespread in and adjacent to the world＇s drylands, as well as deriving from volcanic vents. Silica is a highly fibrogenic agent in lung tissue. Long-term inhaling of siliceous dusts can lead to a number of fibrotic lung diseases, including natural （non-occupational） pneumoconioses （notably silicosis, but including asbestosis and others）. Different polymorphs of silica show different levels of toxicity in interaction with lung tissue. Particles with highly active surfaces may release radicals, causing cell damage. Some types of inhaled particulates are degraded by macrophages, but many are highly resistant and persist in the lungs, some stimulating fibroblastic cells to deposit collagen. Silicosis is an inflammation of the lung commonly caused by silicate mineral particles, leading to fibrosis. Three types are recognized： nodular pulmonary fibrosis （simple or chronic silicosis）, acute silicosis, and accelerated silicosis. Generally, finer particulates have greater oxidative capacity than the coarser fractions. They contain more reactive oxygen species, their greater bioreactivity making them more toxic to pulmonary tissue. Nevertheless, inhalation of large dust particles （〉 10μm） may constitute a health risk if the mineralogy is toxic, regardless of where the grains lodge in the respiratory system. Dust may absorb harmful gases, disease-generating bacteria and carcinogenic hydrocarbon compounds. Silica-related respiratory disease may also an exacerbate cardiac problem, and epidemiology suggests a link with tuberculosis. Quantification of dust loading and exposure requires study of spatial and temporal patterns, supported by meteorological analysis, airflow modeling and satellite-borne imagery. Some acute, short-term health impacts have been assessed using atmospheric and health records both before and after a dust storm or by comparison of populations within and outside such events. Analysis of the size, shape, mineralogy and geochemistry of ambient dust particulates provides information on natural dust sources, dust concentrations, and potential particulate toxicity, as well as providing a datum for assessment of human exposure levels.     

Total gaseous mercury emissions from mercury-enriched soil in Guizhou, China

Guizhou is located in the Circum-Pacific Global Mercuriferous Belt, and mercury concentrations in soil in this area are enriched. In-situ total gaseous mercury （TGM） exchange fluxes between air and soil surface were intensively measured at four sampling sites in Guiyang from 21 May to 16 June, 2003, and five sites in the Lanmuchang mercury mining area in December 2002 and May 2003, respectively. The in-situ Hg flux measurement was conducted with a dynamic flux chamber （DFC） of quartz. Overall, net emissions were obtained from all sampling sites. Soil mercury concentration and solar radiation have been proved to be the two most important parameters to control mercury emissions from soil. Meanwhile, rain events can enhance mercury emission rate significantly.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Heavy metals in freshly deposited sediments of the river Subernarekha, India: an example of lithogenic and anthropogenic effects

Heavy metal distribution patterns in river sediments aid in understanding the exogenic cycling of elements as well as in assessing the effect of anthropogenic influences. In India, the Subernarekha river flows over the Precambrian terrain of the Singhbhum craton in eastern India. The rocks are of an iron ore series and the primary rock types are schist and quartzite. One main tributary, the Kharkhai, flows through granite rocks and subsequently flows through the schist and quartzite layers. The Subernarekha flows through the East Singhbhum district, which is one of India’s industrialised areas known for ore mining, steel production, power generation, cement production and other related activities. Freshly deposited river sediments were collected upstream and downstream the industrial zone. Samples were collected from four locations and analysed in <63-μm sediment fraction for heavy metals including Zn, Pb, Cd and Cu by anodic stripping voltammetry. Enrichment of these elements over and above the local natural concentration level has been calculated and reported. Sediments of the present study are classified by Muller’s geo-accumulation index (I _geo) and vary from element to element and with climatic seasons. During pre-monsoon period the maximum I _geo value for Zn is moderately to highly polluted and for Cu and Pb is moderately polluted, respectively, based on the Muller’s standard. Anthropogenic, lithogenic or cumulative effects of both components are the main reasons for such variations in I _geo values. The basic igneous rock layer through which the river flows or a seasonal rivulet that joins with the main river may be the primary source for lithogenic components.     

Analysis of a rigid footing lying on three-layered soil using the finite difference method

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.     

Performance of Mortar and Chemical Grout Injection into Surrounding Soil When Slurry Pipe-jacking Method is Used

Small-diameter shallow tunnels are often being built by using the slurry pipe-jacking method. This system involves the pushing or thrusting of a drivage machine and concrete pipes into the ground. Chemical grout injection into the surrounding soil around the tunnel is carried out after the drivage and pushing processes are finished. The purpose of the chemical grout injection is to maintain permanent stability of the surrounding soil. However, the behavior of the chemical grouting material in the surrounding soil around the tunnel and the amount of optimum injection is not clearly understood. From these points of view, this paper discusses the performance of the chemical grouting material, when it is injected into the surrounding soil around the tunnel, by means of 2-D Eulerian–Lagrangian seepage analysis. Moreover, the effectiveness of the chemical grout injection was evaluated by using the non-linear finite element method. This investigation show when the range of the grouted zone is designed; it is necessary that the relationship between Young’s modulus of the soil/grouted zone and the confining stress be taken into consideration in order to establish effective, economical and safe chemical grout injection system. Understanding the performance of the seepage/dispersion behavior of the chemical grout and the characteristics of soil/ grouted zone is also important.     

River bank erosion control by soil nailing

A study has been done for analysing soil nailed cuts with circular type wedge failure by friction circle method. Various parameters such as nail length, nail diameter, nail inclination, wall inclination and angle of internal friction of soil have been considered to determine the factor of safety of nailed open cuts. The study shows that for cohesionless soil nailed cut, factor of safety increase with increase of parameters like angle of internal friction of soil, length of nail (L) versus height of cut (H) ratio, cohesion of soil and nail inclination (upto 15°) with horizontal. The study revealed that nails grouted with cement perform better than driven nails. A case study further confirms the analytical findings. Received 7 October     

侵蚀红壤植被恢复后土壤呼吸日动态及日呼吸速率测定方法比较

采用密闭室红外气体分析仪法(IRGA法)观测了中亚热带红壤侵蚀裸地植被恢复后不同季节土壤呼吸速率的日动态变化,并比较了IRGA法与碱吸收法(AA法)测定的土壤呼吸速率.结果表明:侵蚀裸地植被恢复后土壤呼吸速率日动态呈单峰曲线,与土壤温度的昼夜变化基本一致,最高值一般出现在午后13∶00~17∶00,最低值出现在凌晨3∶00~7∶00;植被恢复显著提高了土壤日呼吸速率,但明显降低了土壤呼吸速率日变化幅度;马尾松林对土壤呼吸速率日变化幅度降低程度高于板栗园和百喜草地,且对夏季的降低程度影响最大.而IRGA法和AA法测定的土壤呼吸速率具有显著的幂函数关系,AA法测定的土壤呼吸速率为IRGA法的27.5%~218%,平均为76.2%.当土壤呼吸速率较低时,AA法比IRGA法高估了土壤呼吸速率;反之,AA法则低估了土壤呼吸速率.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.