The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

成矿流体动力学与大地构造环境的关系(英文)

热液矿床的形成需要大量流体,而流体的运移需要包括地势差、岩石变形、热梯度及热异常等多种驱动力。这些流体驱动力与构造环境及过程有密切关系。流体压力状态和热场及地势差的组合决定了压性构造环境流体以向上运动为主,如造山型成矿系统,而张性构造环境流体以对流为主,如VMS和SEDEX成矿系统。造山作用造成的地势差及水平挤压作用所产生的超压可以驱动流体侧向迁移数百公里,如MVT成矿系统。但是,具体的成矿流体动力学过程比较复杂,且随构造演化而变化。在经历过多个大地构造演化阶段的地区,如地洼区,老的成矿流体动力系统不断被新的系统叠加或取代;新的构造单元需要来自地幔的流体源补充才有利于成矿。成矿学研究的是成矿作用与大地构造的关系,而流体动力学系统与大地构造环境密切相关,因此,成矿流体动力学应成为成矿学的一个重要组成部分。     

On core surface flows inferred from satellite magnetic data

Satellite-data allows the magnetic field produced by the dynamo within the Earth’s core to be imaged with much more accuracy than previously possible with only ground-based data. Changes in this magnetic field can in turn be used to make some inferences about the core surface flow responsible for them. In this paper, we investigate the improvement brought to core flow computation by new satellite-data based core magnetic field models. It is shown that the main limitation now encountered is no longer the (now high) accuracy of those models, but the “non-modelled secular variation” produced by interaction of the non-resolvable small scales of the core flow with the core field, and by interaction of the (partly) resolvable large scales of the core flow with the small scales of the core field unfortunately masked by the crustal field. We show how this non-modelled secular variation can be taken into account to recover the largest scales of the core flow in a consistent way. We also investigate the uncertainties this introduces in core flows computed with the help of the frozen-flux and tangentially geostrophic assumptions. It turns out that flows with much more medium and small scales than previously thought are needed to explain the satellite-data-based core magnetic field models. It also turns out that a significant fraction of this flow unfortunately happens to be non-recoverable (being either “non-resolvable” because too small-scale, or “invisible”, because in the kernel of the inverse method) even though it produces the detectable “non-modelled secular variation”. Applying this to the Magsat (1980) to Ørsted (2000) field changes leads us to conclude that a flow involving at least strong retrograde vortices below the Atlantic Hemisphere, some less-resolved prograde vortices below the Pacific Hemisphere, and some poorly resolved (and partly non-resolvable) polar vortices, is needed to explain the 1980-2000 satellite-era average secular variation. The characteristics of the fraction of the secular variation left unexplained by this flow are also discussed.     

In-line response of a cylinder in oscillatory flow

The most widely used mathematical model to represent flow-induced in-line forces on structures is based on the Morison¹ equation. The present paper investigates the validity of using an extension of Morison's equation for non-stationary structures, by comparing predictions with results from a simple laboratory experiment. An elastically-mounted circular cylinder is placed in the sinusoidal flow of a U-tube, and responds in-line with the flow. Cylinder forces and responses are recorded over a range of Keulegan Carpenter numbers up to 35. An equation of motion is solved simply by using relative coordinates and by employing equivalent linearisation. The linear results are compared over a wide variation of parameters with solutions using the full nonlinear equation. Thereafter experimental results are compared with linear predictions.     

Changes in overland flow and infiltration after a rangeland fire in a Mediterranean scrubland

Changes in overland flow and infiltration after a wildfire (summer 1989) in a typical Mediterranean scrubland were measured during the winters of 1990, 1991, 1992 and 1995 by means of simulated rainfall. Infiltration increases gradually from 1990 (sixth months after the forest fire) to 1995 (five and a half years after the forest fire). Overland flow decreases from 45% of rainfall after the forest fire to less than 6% five and a half years later. The reduction in overland flow was greatest in the first two years after the fire because of the quick recovery of vegetation. The steady-state infiltration capacity increased every year after the fire. Runoff and infiltration changes are mainly determined by the gradual recovery of vegetation. © 1998 John Wiley & Sons, Ltd.     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Changes in the abundance and composition of picophytoplankton in relation to the occurrence of a Kyucho and a bottom intrusion in the Bungo Channel,Japan

The Bungo Channel in southwestern Japan receives both warm, called Kyucho, and cold deep-water intrusions (bottom intrusion) from the Pacific Ocean. Abundances of Prochlorococcus, Synechococcus, and eukaryotic picophytoplankton were monitored from 18 July to 17 August 2001 to clarify whether advected picophytoplankton from the Pacific Ocean can grow in the channel or not. Synechococcus cells were further discriminated into low- and high-PUB types according to their fluorescence property in flow cytometry. From 18 to 25 July, the water temperature decreased by 3 °C at a 5-m depth at all stations, indicating the occurrence of a bottom intrusion. From 25 July to 4 August, a Kyucho occurred and the water temperature rapidly increased. From 4 to 17 August, a bottom intrusion and a Kyucho both occurred twice, although the intensities were smaller than those occurring until 4 August. From 18 to 30 July, the abundance of both Prochlorococcus and a high-PUB type of Synechococcus drastically decreased because of a bottom intrusion; however, the abundances rapidly increased due to the advection by a Kyucho. These advected cells increased from 4 to 17 August in the channel and Kitanada Bay. Changes in the abundance of low-PUB type of Synechococcus and eukaryotic picophytoplankton were less noticeable than those in the abundance of Prochlorococcus and high-PUB type. The present study demonstrated that oceanic picophytoplankton advected by the Kyucho could grow in the channel. However, abundances of low-PUB type and eukaryotic picophytoplankton increased higher than those of Prochlorococcus and high-PUB type did. Thus, these oceanic phytoplankters will be excluded when Kyucho does not occur for a long time. The co-occurrence of various types of picophytoplankton found in the channel is probably achieved by both Kyucho event and their growth capability in the channel.     

Evidence for deep groundwater flow and convective heat transport in mountainous terrain, Delta County, Colorado, USA

The Tongue Creek watershed lies on the south flank of Grand Mesa in western Colorado, USA and is a site with 1.5 km of topographic relief, heat flow of 100 mW/m², thermal conductivity of 3.3 W m^–1 °C^–1, hydraulic conductivity of 10^-8 m/s, a water table that closely follows surface topography, and groundwater temperatures 3–15°C above mean surface temperatures. These data suggest that convective heat transport by groundwater flow has modified the thermal regime of the site. Steady state three-dimensional numerical simulations of heat flow, groundwater flow, and convective transport were used to model these thermal and hydrological data. The simulations provided estimates for the scale of hydraulic conductivity and bedrock base flow discharge within the watershed. The numerical models show that (1) complex three-dimensional flow systems develop with a range of scales from tens of meters to tens of kilometers; (2) mapped springs are frequently found at locations where contours of hydraulic head indicate strong vertical flow at the water table, and; (3) the distribution of groundwater temperatures in water wells as a function of surface elevation is predicted by the model.     

新一轮二滩电站库区遥感调查

新一轮二滩电站库区遥感调查是在我国前20年滑坡、泥石流遥感调查技术方法基础上的一次革新。本文介绍了革新方法的数据源、工作平台、技术流程、及人机交互与自动分类相结合的遥感解译方法，并对获得的工作区滑坡、泥石流的数量、分布、特性及它们的环境背景作简要分析。     

Geodynamic modelling of the Century deposit,Mt Isa Province,Queensland

This paper is concerned with an example of quantitative modelling of orebody formation as a guide to reducing the risk for future mineral exploration. Specifically, the paper presents a detailed 3–D numerical model for the formation of the Century zinc deposit in northern Queensland. The model couples fluid flow with deformation, thermal transport and chemical reactions. The emphasis of the study is a systems approach where the holistic mineralising system is considered rather than concentrating solely on the mineral deposit. In so doing the complete plumbing system for mineralisation is considered with a view to specifying the critical conditions responsible for the ore deposit occurring where it does and having the size and metal grades that are observed. The numerical model is based on detailed geological, tectonic, isotopic and mineralogical data collected over the past 20 years. The conclusions are that the Century zinc deposit is located where it is because of the following factors: (i) a thermal anomaly is associated with the Termite Range Fault due to advection of heat from depth by fluid flow up the Termite Range Fault; (ii) bedding‐plane fissility in the shale rocks hosting the Century zinc deposit has controlled the wavelength and nature of D₁ folding in the vicinity of the deposit and has also controlled increases in permeability due to hydrofracture of the shales; such hydrofracture is also associated with the production of hydrocarbons as these shales passed through the ‘oil‐window’; (iii) Pb–Zn leached from crustal rocks in the stratigraphic column migrated up along faults normal to the Termite Range Fault driven by topographic relief associated with inversion at the end of the Isan Orogeny; these fluids mixed with H₂S derived at depth moving up the Termite Range Fault to mix with the crustal fluids to precipitate Pb–Zn in a plume downstream from the point of mixing. Critical factors to be used as exploration guides are high temperatures, carbonaceous fissile shales now folded into relatively tight D₁ folds, fault‐controlled plumbing systems that enable fluid mixing, depletion of metals upstream of the deposit and,in particular,a very wide Fe‐depletion halo upstream of the deposit.