人流联系和经济联系视角下区域城市关联比较——基于手机信令数据和企业关联数据的研究

手机信令数据所代表的人流联系与国内外研究常用的企业关联数据在测度城市关联存在何种差异,尚缺乏明确的研究。本文以江西省北部地区为研究对象,分别以手机信令数据、企业关联数据测度城市关联,比较了两种城市关联与城市之间空间距离、人口规模的关系,比较了两种城市关联网络的层级、结构、腹地。主要研究发现：人流联系不仅集中在城市周边地区,而且集中于较高等级城市之间;经济联系则主要集中在较高等级城市之间。人流联系、经济联系受到城市人口规模影响均较小。人流联系显著地符合幂函数空间衰减规律。手机信令数据和企业关联数据是两种不同的城市关联数据,但可互为补充,较为全面地反映城市关联。     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

Hybrid finite element–finite volume discretization of complex geologic structures and a new simulation workflow demonstrated on fractured rocks

The generation of computational meshes of complex geological objects is a challenge: their shape needs to be retained, resolution has to adapt to local detail, and variations of material properties in the objects have to be represented. Also mesh refinement and adaptation must be sufficient to resolve variations in the computed variable(s). Here, we present an unstructured hybrid finite element, node‐centred finite‐volume discretization suitable for solving fluid flow, reactive transport, and mechanical partial differential equations on a complex geometry with inhomogeneous material domains. We show that resulting meshes accurately capture free‐form material interfaces as defined by non‐uniform rational B‐spline curves and surfaces. The mesh discretization error is analysed for the elliptic pressure equation and an error metric is introduced to guide mesh refinement. Finite elements and finite volumes are represented in parametric space and integrations are conducted numerically. Subsequently, integral properties are mapped to physical space using Jacobian transformations. This method even retains its validity when the mesh is deformed. The resulting generic formulation is demonstrated for a transport calculation performed on a complex discrete fracture model.     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

Heat transport at boiling,near‐critical conditions

Two‐phase flow and near‐critical phenomena are likely to enhance energy transport in high‐temperature hydrothermal systems. We present a series of two‐dimensional simulations of two‐phase flow of pure water at near‐critical conditions. The results show that at near‐critical conditions, two‐phase convection can be more efficient in transporting energy than single‐phase convection. The highest heat fluxes are attained when two‐phase heat‐pipes form near the bottom boundary, recharging the root of the upflow zone and thereby enabling the formation of broad upflow regions. When the system becomes more vapor‐dominated, it loses this ability, upflow zones become narrower and the energy efficiency drops to more moderate values.     

Modelling of hydrology and potential population levels at Bronze Age Jawa,Northern Jordan: a Monte Carlo approach to cope with uncertainty

An annual water balance model for Wadi Rajil, in Northern Jordan, is used to simulate the ancient water supply system for the Early Bronze Age site of Jawa. The model includes: water delivery from the catchment; local pond storage; and water demand for people, animals and irrigation. A Monte Carlo approach is used to incorporate the uncertainty associated with a range of factors including rainfall, evaporation, water losses and use. The stochastic simulation provides estimates of potential population levels sustainable by the water supply system. Historical precipitation estimates from a Global Circulation Model, with uncertainty bounds, are used to reconstruct the climate at Jawa in the Early Bronze Age (EBA). Model results indicate that the population levels in the predicted wetter conditions during parts of the EBA could have risen to ∼6000 and may have been higher in wet years. However, palaeoclimatic proxies also suggest prolonged droughts in the EBA; and during these periods the water management system was unable to provide adequate supply for a population of 6000. The utility of Monte Carlo based hydrological modelling as a tool within archaeological science is discussed.     

Resource abundance constraints on the early post-contact Cherokee population

We present the combination of an analysis of resource demand by the early post-contact (1721) Cherokee population with spatially explicit estimates of production for five key resources: architectural land, agricultural land, firewood, hard mast, and white-tailed deer (Odocoileus virginianus). We combine a recent synthesis of village location and population, a map of recognized Cherokee territory, digital terrain data, estimates of per capita demand, and productive capacity for each resource. Average, high, and low demands were estimated for each resource and assigned based on a weighted function of terrain and distance from each village. We conclude that Cherokee demands for architectural and agricultural land, hard mast, and fuelwood were easily met within a short proximity to each town under all combinations of production and demand. These resources were likely not limiting, and were satisfied for the entire Cherokee population by less than 1% of the entire recognized Cherokee territory in 1721. These resources likely exceeded demand even when sources were restricted to the convex hull of the Cherokee territory, or to near stream, flat regions. Deer resources were likely harvested over a much larger area and to a much greater extent. Our best estimate of deer resource demand was 32% of annual sustainable production in the Cherokee territory, with from 16 to 48% of estimated sustainable production harvested for low and high demand estimates, respectively. Our estimates vary in response to uncertainties in deer production, harvest proportion, deer density, and sustainable harvest rates. Deer demand was substantially higher under all combinations of conditions than that available within the convex hull of Cherokee villages, indicating significant travel was needed to furnish deer requirements. Spatially explicit models that consider terrain- and distance-related tradeoffs suggest that Cherokee demand for deer drove harvest over areas consisting of over half the recognized Cherokee territory.     

The significance of pockmarks to understanding fluid flow processes and geohazards

Underwater gas and liquid escape from the seafloor has long been treated as a mere curiosity. It was only after the advent of the side‐scan sonar and the subsequent discovery of pockmarks that the scale of fluid escape and the moon‐like terrain on parts of the ocean floor became generally known. Today, pockmarks ranging in size from the ‘unit pockmark’ (1–10 m wide, < 0.6 m deep) to the normal pockmark (10–700 m wide, up to 45 m deep) are known to occur in most seas, oceans, lakes and in many diverse geological settings. In addition to indicating areas of the seabed that are ‘hydraulically active’, pockmarks are known to occur on continental slopes with gas hydrates and in association with slides and slumps. However, possibly their potentially greatest significance is as an indicator of deep fluid pressure build‐up prior to earthquakes. Whereas only a few locations containing active (bubbling) pockmarks are known, those that become active a few days prior to major earthquakes may be important precursors that have been overlooked. Pockmark fields and individual pockmarks need to be instrumented with temperature and pressure sensors, and monitoring should continue over years. The scale of such research calls for a multinational project in several pockmark fields in various geological settings.     

Fluid flow and stability of the US continental slope offshore New Jersey from the Pleistocene to the present

We predict that portions of the New Jersey continental slope were unstable approximately 0.5 million years ago. This instability was caused by rapid sediment loading during a Pleistocene sea‐level lowstand and by flow focusing in underlying, permeable Miocene strata. The simulated instability is consistent with soft‐sediment deformation and small slumps in Pleistocene strata of the Hudson Apron. Stability of the New Jersey margin has increased since 0.3 Ma because sedimentation rate has decreased. Today, the modelled factor of safety (FS) for the upper slope is approximately 1.5 whereas in the lower slope it exceeds 3. We predict that sedimentation rate is a dominant factor on slope stability. When rapid and asymmetric loading of a highly permeable sedimentary layer occurs, the location of instability can shift seaward to regions where sedimentation rates are low. Stability calculations use pressures and effective stresses predicted by a coupled sedimentation‐fluid flow model. This hydrodynamic analysis demonstrates how the interplay of sedimentation and fluid migration affects the distribution, timing, and size of sedimentary failures.