Post-solidification cooling and the age of Io's lava flows

The modeling of thermal emission from active lava flows must account for the cooling of the lava after solidification. Models of lava cooling applied to data collected by the Galileo spacecraft have, until now, not taken this into consideration. This is a flaw as lava flows on Io are thought to be relatively thin with a range in thickness from ∼1 to 13 m. Once a flow is completely solidified (a rapid process on a geological time scale), the surface cools faster than the surface of a partially molten flow. Cooling via the base of the lava flow is also important and accelerates the solidification of the flow compared to the rate for the ‘semi-infinite’ approximation (which is only valid for very deep lava bodies). We introduce a new model which incorporates the solidification and basal cooling features. This model gives a superior reproduction of the cooling of the 1997 Pillan lava flows on Io observed by the Galileo spacecraft. We also use the new model to determine what observations are necessary to constrain lava emplacement style at Loki Patera. Flows exhibit different cooling profiles from that expected from a lava lake. We model cooling with a finite-element code and make quantitative predictions for the behavior of lava flows and other lava bodies that can be tested against observations both on Io and Earth. For example, a 10-m-thick ultramafic flow, like those emplaced at Pillan Patera in 1997, solidifies in ∼450 days (at which point the surface temperature has cooled to ∼280 K) and takes another 390 days to cool to 249 K. Observations over a sufficient period of time reveal divergent cooling trends for different lava bodies [examples: lava flows and lava lakes have different cooling trends after the flow has solidified (flows cool faster)]. Thin flows solidify and cool faster than flows of greater thickness. The model can therefore be used as a diagnostic tool for constraining possible emplacement mechanisms and compositions of bodies of lava in remote-sensing data.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Applications of particle-tracking techniques to bank infiltration: a case study from El Paso,Texas, USA

This paper presents results of a small scale study that utilized particle-tracking techniques to evaluate transport of river water through an alluvial aquifer in a bank infiltration testing site in El Paso, Texas, USA. The particle-tracking survey was used to better define filtration parameters. Several simulations were generated to allow visualization of the effects of well placement and pumping rate on flow paths, travel time, the size of the pumping influence zone, and proportion of river-derived water and groundwater mixing in the pumping well. Simulations indicate that migration of river water into the aquifer is generally slow. Most water does not arrive at the well by the end of an 18-day pumping period at 0.54 m³/min pumping rate for a well located 18 m from the river. Forty-four percent of the water pumped from the well was river water. The models provided important information needed to design appropriate sampling schedules for bank filtration practices and ensured meeting adequate soil-retention times. The pumping rate has more effect on river water travel time than the location of the pumping well from the river. The examples presented in this paper indicate that operating the pumping well at a doubled distance from the river increased the time required for the water to travel to the well, but did not greatly change the capture zone.     

Integrating diverse methods to understand climate-land interactions in East Africa

The questions of how land use change affects climate, and how climate change affects land use, require examination of societal and environmental systems across space at multiple scales, from the global climate to regional vegetative dynamics to local decision making by farmers and herders. It also requires an analysis of causal linkages and feedback loops between systems. These questions and the conceptual approach of the research design of the Climate-Land Interaction Project (CLIP) are rooted in the classical human-environment research tradition in Geography.This paper discusses a methodological framework to quantify the two-way interactions between land use and regional climate systems, using ongoing work by a team of multi-disciplinary scientists examining climate-land dynamics at multiple scales in East Africa. East Africa is a region that is undergoing rapid land use change, where changes in climate would have serious consequences for people’s livelihoods, and requiring new coping and land use strategies. The research involves exploration of linkages between two important foci of global change research, namely, land use/land cover (LULC) and climate change. These linkages are examined through modeling agricultural systems, land use driving forces and patterns, the physical properties of land cover, and the regional climate. Both qualitative and quantitative methods are being used to illustrate a diverse pluralism in scientific discovery.     

Geoacoustic Inversion via Genetic Algorithm and Its Application to Manganese Sediment Identification

An acoustic inversion method using a wide-band signal and two near field receivers is proposed and applied to multiple layered seabed models including a manganese sediment. The inversion problem can be formulated into a probabilistic model comprised of signals, a forward model, and additive noise. The forward model simulates wide-band signals, such as chirp signals, and is chosen to be the source-waveletconvolution plane wave modeling method. The wavelet matching technique, using weighted least-squares fitting, estimates the sediment sound-speed and thickness on which determination of the possible numerical ranges for a priori uniform distribution is based. The genetic algorithm is applied to a global optimization problem to find a maximum a posteriori solution for determined a priori search space. Here the object function is defined by an L 2 norm of the difference between measured and modeled signals. Not only the marginal pdf but also its statistics are calculated by numerical evaluation of integrals using the samples selected during importance sampling process of the genetic algorithm.     

Integration of a Multi-Stage Bulk Materials Acquisition System from the Viewpoint of Dock Operations

A multi-stage bulk materials acquisition system is examined from the viewpoint of dock operations. The whole process can be divided into four stages and treated as four subsystems: ships scheduling planning, material yard planning, dock arrangement, and material discharging. Since the structure and the complexity of each subsystem is different, the solution approaches applied are also different with respect to each. To increase operations efficiency and data consistency, an integration of the subsystems by coordinating the operations sequence and data communications for the four subsystems is discussed. The four subsystems are illustrated by using the material docks of China Steel Corporation, and the models developed in this research are also validated. The results indicate that the models developed are capable of deriving solutions better than the existing ones. Although the discussions are confined to China Steel Corporation, the models can be applied to other companies with similar operations.     

Comparison of two three-dimensional hydrodynamic modeling systems for coastal tidal motion

A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.     

渤海、黄海强温跃层数值预报模式

文中对前人关于水温垂直结构计算、模拟和预报研究成果进行了评价，并在此基础上建立了考虑海面吸收辐射和透射辐射、地形、风、海流、界面摩擦及其引起的混合作用的强温跃层三维数值预报模式。该模式运行的驱动量仅为风和气温场资料，因而具有较好的实用性；此外，从试报的结果看，效果是令人满意的。     

复杂网络河的数值模拟

依据复杂网络河的特点及其内在联系，本文建立了其特征矩阵及由Newton-Raphson格式表述的隐式圣－维南差分方程组的系数矩阵，经输入边界条件及解方程组，可得到不同各支段上，下游端点上的水位和流量值。     

Numerical simulation of migration–accumulation of oil resources

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.