Crystal-rich mass flow deposits related to the eruption of a sublacustrine silicic cryptodome (Early Permian Collio Basin, Italian Alps)

Sedimentological, granulometric and petrographic data are presented from a detailed study on a crystal-rich mass flow deposit, which is presumably related to the eruption of a sublacustrine cryptodome. The deposit forms a prominent intercalation in the Lower Permian Collio Formation in the Italian Alps north of Brescia. Outcrops of the 10–20-m-thick volcaniclastic deposit (Dasdana I Beds, DB) can be traced over 12 km from east to west. The DB consists of a thick, crystal-rich, sandy–gravelly lower subunit representing a sequence of amalgamated Bouma-a(b) divisions overlain by a thin, well-bedded, sandy–muddy subunit that is rich in outsize porphyritic silicic fragments. Modal and computer-aided image analyses reveal that the crystal-rich lower subunit contains up to 80% of volcanogenic crystals. Some samples contain up to 60% of porphyritic fragments, which have a phenocryst content of about 20%. The wide textural range from cryptocrystalline, poikilomosaic, to rarer medium-grained granophyric groundmass, the irregular to lensoid shapes of the porphyritic fragments, and the presence of basement and sedimentary clasts suggest that the DB originated from a sublacustrine eruption of a partially extrusive cryptodome (ca. 1.6 km³). Two other porphyritic felsic cryptodomes (Dosso dei Lupi, Dosso del Bue), described briefly here, emplaced into the Collio Formation sometime after the DB event, and expose flat bases and tilted sediments at their sides. Textures observed in these domes are comparable to those found in the DB porphyritic fragments.     

Modelling of meteorological conditions at an urban scale for the PUMA winter campaign

For the winter 2000 campaign of the Pollution of Urban Midlands Atmosphere project, observation and numerical modelling of meteorological conditions over the West Midlands conurbation, UK, was undertaken. Modelling was performed using the regional atmospheric meteorological system (RAMS). This paper presents a comparison of modelled and observed wind and temperature for 25 and 26 January 2000. The RAMS model uses two nested grids with a mesh size of 2 km for the inner grid which is embedded in the outer grid with a mesh size of 8 km. Statistical evaluation of the model results against the observational data of wind speed, direction and temperature at 10 m was conducted. In general, the modelling results are in a reasonable agreement with observation. The statistical evaluation suggests that model performance is poorer for the inner grid than the outer grid as the model uncertainties (mainly mean bias) transfer from the outer to inner one. The low indices of agreement of temperature and wind are mainly associated with the systematic root-mean-square-difference values. For temperature, the systematic bias may also be affected by representation of cloud amount by the model. For wind, the model tends to have a poor performance for calm conditions, as under a stable anti-cyclonic situation local wind patterns associated with topography may develop, although the topography of the region is relatively flat. The results for the inner grid reveal some subtle spatial patterns at a scale smaller than 10 km near hills and valleys with differences in elevation of a few hundred metres.     

Development and testing of numerical methods for two-way nested air pollution modelling

A new 3-D model REGINA (REGIonal high resolutioN Air pollution model) is under development at the National Environmental Research Institute (NERI). The model is based on pieces from several models developed over the last decades at NERI. The aim of the work is to develop a nested model which can operate with very high resolution in both space and time. To fulfill this aim the choice and implementation of accurate numerical methods is crucial. The model will be applied for studying air pollution phenomena (both monitoring, forecasting and scenarios) over Denmark. The present paper is focussed on the modification, implementation and testing of a numerical method for treating the horizontal advection in the model as well as the implementation of two-way nesting techniques. The horizontal transport in the model is solved using an accurate space derivative algorithm. This method is traditionally implemented with periodic boundary conditions, however, this is not an option for nested modelling. A new method for calculating non-periodic boundary conditions has been developed in order to overcome this problem. Extensive testing of the numerical solution of the advection and the coupling of the solution of advection and chemistry in the model using Molenkamp–Crowley rotation tests have been carried out. The results show that the model with the current implementation of numerical methods is suitable for calculating air pollution levels with high resolution.     

A model of equilibrium bed topography for meander bends with erodible banks

Channel curvature produces secondary currents and a transverse sloping channel bed, along which the depth increases towards the outer bank. As a result deep pools tend to form adjacent to the outer bank, promoting bank collapse. The interaction of sediment grains with the primary and secondary flow and the transverse sloping bed also causes meanders to move different grain sizes in different proportions and directions, resulting in a consistent sorting pattern. Several models have been developed to describe this process, but they all have the potential to over‐predict pool depth because they cannot account for the influence of erodible banks. In reality, bank collapse might lead to the development of a wider, shallower cross‐section and any resulting flow depth discrepancy can bias associated predictions of flow, sediment transport, and grain‐size sorting. While bed topography, sediment transport and grain sorting in bends will partly be controlled by the sedimentary characteristics of the bank materials, the magnitude of this effect has not previously been explored. This paper reports the development of a model of flow, sediment transport, grain‐size sorting, and bed topography for river bends with erodible banks. The model is tested via intercomparison of predicted and observed bed topography in one low‐energy (5·3 W m^?2 specific stream power) and one high‐energy (43·4 W m^?2) study reach, namely the River South Esk in Scotland and Goodwin Creek in Mississippi, respectively. Model predictions of bed topography are found to be satisfactory, at least close to the apices of bends. Finally, the model is used in sensitivity analyses that provide insight into the influence of bank erodibility on equilibrium meander morphology and associated patterns of grain‐size sorting. The sensitivity of meander response to bank cohesion is found to increase as a function of the available stream power within the two study bends. Copyright © 2002 John Wiley & Sons, Ltd.     

流动成像测井研究进展

流动成像测井属于对非均匀介质动态的非线性测量，通过对油气井内多相流体信息的投影扫描和反演处理，实时成像显示流体分布及流动状况，国外研制的电导法，电容法仪器采用电导探针或电容元件构成阵列测量探头，分别利用油气与水的导电特性和介电特性差异辨识井内流体，对流动截面的测量局限于个别点上，已在水平井中见到应用效果，但对物场信息投影测量的数据量和分辨率未能满足成像要求。我国提出的电磁法采用环状阵列电极构成测量探头，综合利用油气与水的导电和介电特性差异辨识井内流体，可以测量获得多相流体流动截面的清晰图像，为研究开发全新的流动成像测井系统奠定了良好基础。     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Power of the Mann–Whitney test for detecting a shift in median or mean of hydro-meteorological data

 The non-parametric Mann–Whitney (MW) statistic test has been popularly used to assess the significance of a shift in median or mean of hydro-meteorological time series. It has been considered that the test is more suitable for non-normally distributed data and it may be not sensitive to the distribution type of sample data. However, no evidence has been provided to demonstrate these. This study investigates the power of the test in various circumstances by means of Monte Carlo simulation. Simulation results demonstrate that the power of the test is very sensitive to various properties of sample data. The power depends on the pre-assigned significance level, magnitude of a shift, sample size, and its occurrence position within a time series; and it is also strongly affected by the variation, skewness, and distribution type of a time series. The bigger the magnitude of a shift, the more powerful the test is; the larger the sample size, the more powerful the test is; and the bigger the variation within a time series, the less power the test has. The test has the highest power if a shift occurs at the midpoint of a time series. For the samples with different distribution types, the power of the test is dramatically different. The test has the highest power for time series with the extreme value type III (EV3) distribution while it indicates the lowest power for time series with the lognormal distribution.     

不同损失条件下的泥石流预报模型

泥石流预报是泥石流减灾的重要手段之一，然而泥石流形成的复杂性使泥石流预报准确度低，误报和漏报率较高。泥石流误报和漏报都会造成损失，但二者造成的损失有很大的差别。为了减少泥石流误报或漏报造成的损失，应当考虑两种错报造成损失的不同。根据使总平均损失达到最小的原则，建立了不同损失条件下的泥石流预报模型，并将该模型应用到云南东川蒋家沟。     

泥石流多发干旱河谷区植被恢复研究

云南东北部东川小江流域的干旱河谷土地退化十分严重，土地表层的砾石含量已超过60％，有机质含量、全氮、有效氮、有效磷、有效钾的含量相对于残存燥红土明显下降。草坡是当地的主要植被类型，并未起到防治土地退化的应有作用。以合欢（Leucaena glauca)、马桑（Coriaria sinica)种植为主的恢复方式对土地的恢复有明显的改善作用。高密度的乔木树种合欢的种植使林下的草本物种的数量和组成发生了极大变化，较耐荫的植物占据了林分下层，且物种数量较草坡减少，群落多样性指数降低。三种马桑种植方式下，不但原有的草坡群落的物种仍有大量的生存空间，而且为一些新物种创造了生存环境，从而使物种的种类较原有草坡地增加，提高了群落的多样性指数。合欢种植的生物量最高，对解决当地薪柴严重短缺意义重大，马桑种植在增加灌木层生物量的同时，并未改变草坡的物种组成，还可兼顾解决当地的饲草和放牧，应大力推广。     

改进的灰色模型在中国能源消费及其CO2排放预测中的研究

中国高速的经济发展导致了快速的能源消费和地球温暖化气体的排放量。本文应用灰色理论建立了中国的终端能源消费量的动态预测模型，同时应用等维递补及带有马尔科夫链符号估计的残差修正，对原始模型进行了改进。后验检验的结果证明了改进的预测模型具有较高的计算精度。预测结果表明，两种经济发展情况（高速及低速）下，从2000年到2030年期间的终端能源消费量的平均增长率将分别达到3．06％和2．18％，2030年CO2 的排放量将分别达到2000年的 2． 15及 1．60倍。其中工业及居民生活的能源消费量的增加速度快于其他部门。这表明了中国将要加快其工业化及都市化的进程．终端消费中，煤炭的消费量所占的比例将要逐年减少，而电力的比例则将逐年增大。