DEM 点位地形信息量化模型研究

针对DEM 点位, 首先应用微分几何法对其所负载的语法信息量进行测度, 其次根据地形特征点类型及地形结构特征确定其语义信息量, 然后基于信息学理论构建了DEM 点位地形信息综合量化模型。在此基础上, 以黄土丘陵沟壑区作为实验样区, 对DEM 点位地形信息量提取方法及其在地形简化中的初步实例应用进行了探讨和验证。实验结果显示, 所提出的DEM 点位地形信息量化方案可行;基于DEM 地形信息量指数的多尺度DEM 构建方案, 具有机理明确、易于实现的特点, 并通过优先保留地形骨架特征点, 可以有效减少地形失真, 从而满足不同层次的多尺度数字地形建模和表达要求。对DEM 点位地形信息进行有效量化, 为认识DEM 地形信息特征提供了一个新的切入点, 同时为多尺度数字地形建模提供理论依据与方法支持。     

基于反馈机制的城市扩张模拟研究进展

回顾了城市扩张模拟的研究进展,发现相关研究在表达城市系统的复杂性特征时存在较多问题,主要原因是未能较好地反映城市扩张与驱动因子间的相互作用与反馈机制。部分研究虽然隐含地考虑了这种"因-果"相互作用关系,但尚未形成一套成熟的框架体系,主要表现在:1)从研究深度看,反馈机制的实现尚无一套系统的方法;2)就研究广度而言,目前城市扩张模拟较少考虑多类别反馈机制的综合作用。新理论与方法(如复杂系统与网络)的引入及模型耦合等是未来基于反馈机制的城市扩张模拟研究的重要方向。     

Long-term ice sheet–climate interactions under anthropogenic greenhouse forcing simulated with a complex Earth System Model

Several multi-century and multi-millennia simulations have been performed with a complex Earth System Model (ESM) for different anthropogenic climate change scenarios in order to study the long-term evolution of sea level and the impact of ice sheet changes on the climate system. The core of the ESM is a coupled coarse-resolution Atmosphere–Ocean General Circulation Model (AOGCM). Ocean biogeochemistry, land vegetation and ice sheets are included as components of the ESM. The Greenland Ice Sheet (GrIS) decays in all simulations, while the Antarctic ice sheet contributes negatively to sea level rise, due to enhanced storage of water caused by larger snowfall rates. Freshwater flux increases from Greenland are one order of magnitude smaller than total freshwater flux increases into the North Atlantic basin (the sum of the contribution from changes in precipitation, evaporation, run-off and Greenland meltwater) and do not play an important role in changes in the strength of the North Atlantic Meridional Overturning Circulation (NAMOC). The regional climate change associated with weakening/collapse of the NAMOC drastically reduces the decay rate of the GrIS. The dynamical changes due to GrIS topography modification driven by mass balance changes act first as a negative feedback for the decay of the ice sheet, but accelerate the decay at a later stage. The increase of surface temperature due to reduced topographic heights causes a strong acceleration of the decay of the ice sheet in the long term. Other feedbacks between ice sheet and atmosphere are not important for the mass balance of the GrIS until it is reduced to 3/4 of the original size. From then, the reduction in the albedo of Greenland strongly accelerates the decay of the ice sheet.     

Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations

The Regional Atmospheric Modeling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES), developed and evaluated here, is used to explore the effects of three-dimensional canopy heterogeneity, at the individual tree scale, on the statistical properties of turbulence most pertinent to mass and momentum transfer. In RAFLES, the canopy interacts with air by exerting a drag force, by restricting the open volume and apertures available for flow (i.e. finite porosity), and by acting as a heterogeneous source of heat and moisture. The first and second statistical moments of the velocity and flux profiles computed by RAFLES are compared with turbulent velocity and scalar flux measurements collected during spring and winter days. The observations were made at a meteorological tower situated within a southern hardwood canopy at the Duke Forest site, near Durham, North Carolina, U.S.A. Each of the days analyzed is characterized by distinct regimes of atmospheric stability and canopy foliage distribution conditions. RAFLES results agreed with the 30-min averaged flow statistics profiles measured at this single tower. Following this intercomparison, two case studies are numerically considered representing end-members of foliage and midday atmospheric stability conditions: one representing the winter season with strong winds above a sparse canopy and a slightly unstable boundary layer; the other representing the spring season with a dense canopy, calm conditions, and a strongly convective boundary layer. In each case, results from the control canopy, simulating the observed heterogeneous canopy structure at the Duke Forest hardwood stand, are compared with a test case that also includes heterogeneity commensurate in scale to tree-fall gaps. The effects of such tree-scale canopy heterogeneity on the flow are explored at three levels pertinent to biosphere-atmosphere exchange. The first level (zero-dimensional) considers the effects of such heterogeneity on the common representation of the canopy via length scales such as the zero-plane displacement, the aerodynamic roughness length, the surface-layer depth, and the eddy-penetration depth. The second level (one-dimensional) considers the normalized horizontally-averaged profiles of the first and second moments of the flow to assess how tree-scale heterogeneities disturb the entire planar-averaged profiles from their canonical (and well-studied planar-homogeneous) values inside the canopy and in the surface layer. The third level (three-dimensional) considers the effects of such tree-scale heterogeneities on the spatial variability of the ejection-sweep cycle and its propagation to momentum and mass fluxes. From these comparisons, it is shown that such microscale heterogeneity leads to increased spatial correlations between attributes of the ejection-sweep cycle and measures of canopy heterogeneity, resulting in correlated spatial heterogeneity in fluxes. This heterogeneity persisted up to four times the mean height of the canopy (h _c ) for some variables. Interestingly, this estimate is in agreement with the working definition of the thickness of the canopy roughness sublayer (2h _c –5h _c ).     

Micromechanical parameters in bonded particle method for modelling of brittle material failure

Bonded particle modelling (BPM) is nowadays being extensively used for simulating brittle material failure. In BPM, material is modelled as a dense assemblage of particles (grains) connected together by contacts (cement). This sort of modelling seriously depends on the mechanical properties of particle and contact, which are named here as micro‐parameters. However, a definite calibration methodology to obtain micro‐parameters has not been so far established; and many have reported some serious problems. In this research, a calibration procedure to find a unique set of micro‐parameters is established. To attain this purpose, discrete element code of UDEC is used to perform BPM. This code can be conveniently developed by the user. The proposed BPM is composed of rigid polygonal particles interacting at their contact points. These contacts can undergo a certain amount of tension, and their shear resistance is provided by cohesion and friction angle. The results demonstrate that each material macro‐property (i.e. Young's modulus, Poisson's ratio, internal friction angel, internal cohesion, and tensile strength) is directly originated from and distinctly related to the contact properties (i.e. normal and shear stiffness, friction angel, cohesion, and tensile strength). Copyright © 2010 John Wiley & Sons, Ltd.     

The integration of thermal infrared imaging,discharge measurements and numerical simulation to quantify the relative contributions of freshwater inflows to small estuaries in Atlantic Canada

Nutrient fluxes from developed catchments are often a significant factor in the declining water quality and ecological functioning in estuaries. Determining the relative contributions of surface water and groundwater discharge to nutrient‐sensitive estuaries is required because these two pathways may be characterized by different nutrient concentrations and temporal variability, and may thus require different remedial actions. Quantifying the volumetric discharge of groundwater, which may occur via diffuse seepage or springs, remains a significant challenge. In this contribution, the total discharge of freshwater, including groundwater, to two small nutrient‐sensitive estuaries in Prince Edward Island (Canada) is assessed using a unique combination of airborne thermal infrared imaging, direct discharge measurements in streams and shoreline springs, and numerical simulation of groundwater flow. The results of the thermal infrared surveys indicate that groundwater discharge occurs at discrete locations (springs) along the shoreline of both estuaries, which can be attributed to the fractured sandstone bedrock aquifer. The discharge measured at a sub‐set of the springs correlates well with the area of the thermal signal attributed to each discharge location and this information was used to determine the total spring discharge to each estuary. Stream discharge is shown to be the largest volumetric contribution of freshwater to both estuaries (83% for Trout River estuary and 78% for McIntyre Creek estuary); however, groundwater discharge is significant at between 13% and 18% of the total discharge. Comparison of the results from catchment‐scale groundwater flow models and the analysis of spring discharge suggest that diffuse seepage to both estuaries comprises only about 25% of the total groundwater discharge. The methods employed in this research provide a useful framework for determining the relative volumetric contributions of surface water and groundwater to small estuaries and the findings are expected to be relevant to other fractured sandstone coastal catchments in Atlantic Canada. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of the Takagi–Sugeno fuzzy system for combination forecasting of river flow in semiarid mountain regions

Because of high spatial heterogeneity and the degree of uncertainty about hydrological processes in large‐scale catchments of semiarid mountain areas, satisfactory forecasting of daily discharge is seldom available using a single model in many practical cases. In this paper the Takagi–Sugeno fuzzy system (TS) and the simple average method (SAM) are applied to combine forecasts of three individual models, namely, the simple linear model (SLM), the seasonally based linear perturbation model (LPM) and the nearest neighbour linear perturbation model (NNLPM) for modelling daily discharge, and the performance of modelling results is compared in five catchments of semiarid areas. It is found that the TS combination model gives good predictions. The results confirm that better prediction accuracy can be obtained by combining the forecasts of different models with the Takagi–Sugeno fuzzy system in semi‐arid mountain areas. Copyright © 2009 John Wiley & Sons, Ltd.     

Metamorphic evolution,partial melting and rapid exhumation above an ancient flat slab: insights from the San Emigdio Schist,southern California

The San Emigdio and related Pelona, Orocopia, Rand and Sierra de Salinas schists of southern California were underplated beneath the southern Sierra Nevada batholith and adjacent southern California batholith along a shallow segment of the subducting Farallon plate in Late Cretaceous to early Tertiary time. These subduction accretion assemblages represent a regional, deeply exhumed, shallowly dipping domain from an ancient slab segmentation system and record the complete life cycle of the segmentation process from initial flattening and compression to final extensional collapse. An important unresolved question regarding shallow subduction zones concerns how the thermal structure evolves during the slab flattening process. New field relationships, thermobarometry, thermodynamic modelling and garnet diffusion modelling are presented that speak to this issue and elucidate the tectonics of underplating and exhumation of the San Emigdio Schist. We document an upsection increase in peak temperature (i.e. inverted metamorphism), from 590 to 700 °C, peak pressures ranging from 8.5 to 11.1 kbar, limited partial melting, microstructural evidence for large seismic events, rapid cooling (825–380 °C Myr^?1) from peak conditions and an ‘out and back’P–T path. While inverted metamorphism is a characteristic feature of southern California schists, the presence of partial melt and high temperatures (>650 °C) are restricted to exposures with maximum depositional ages between 80 and 90 Ma. Progressive cooling and tectonic underplating beneath an initially hot upper plate following the onset of shallow subduction provide a working hypothesis explaining high temperatures and partial melting in San Emigdio and Sierra de Salinas schists, inverted metamorphism in the schist as a whole, and the observed P–T trajectory calculated from the San Emigdio body. Lower temperatures in Pelona, Orocopia and Rand schists are likewise explained in the context of this overarching model. These results are consistent with an inferred tectonic evolution from shallow subduction beneath the then recently active Late Cretaceous arc to exhumation by rapid trench‐directed channelized extrusion in the subducted schist.     

Tectono‐metamorphic history recorded in garnet porphyroblasts: insights from thermodynamic modelling and electron backscatter diffraction analysis of inclusion trails

In a Barrovian metamorphic sequence, garnetiferous mica schists document a heterogeneously developed superposition of sub‐orthogonal fabrics and multiple garnet growth episodes. In the variably deformed domains, four types of garnet porphyroblasts have been defined based on inclusion trail patterns. Modelled garnet zoning in the MnNCKFMASHTO system indicates a prograde evolution from 4–4.5 kbar and 490–510 °C to 5–6 kbar and 520–550 °C in the earliest subhorizontal fabric progressing towards 6.5–7.5 kbar and 560–590 °C in the subsequent subvertical foliation. This fabric is heterogeneously deformed into a shallow‐dipping retrograde foliation associated with garnet resorption. In situ electron backscatter diffraction measurements of ilmenite inclusions in individual garnet grains yield precise data on included planar and linear elements. Consistent orientations of internal foliations, lineations and foliation intersection axis sets indicate a superposition of three sub‐orthogonal foliation systems. Weak variations of internal records with increasing intensity of deformation suggest that a moderate buckling stage occurred, but apparent lack of porphyroblast rotation is interpreted as a result of dominant passive flow. Coupling the orientation of internal fabric sets with P–T estimates is used to complement the tectono‐metamorphic evolution of the thickened crust. We demonstrate that garnet porphyroblasts preserve features which reflect large‐scale tectonic processes in orogens.     

Modelling coastal ground‐ and surface‐water interactions using an integrated approach

Beach water table fluctuations have an impact on the transport of beach sediments and the exchange of solute and mass between coastal aquifer and nearby water bodies. Details are given of the refinement of a dynamically integrated ground‐ and surface‐water model, and its application to study ground‐ and surface‐water interactions in coastal regions. The depth‐integrated shallow‐water equations are used to represent the surface‐water flow, and the extended Darcy's equation is used to represent the groundwater flow, with a hydrostatic pressure distribution being assumed to apply for both these two types of flows. At the intertidal region, the model has two layers, with the surface‐water layer being located on the top of the groundwater layer. The governing equations for these two types of flows are discretized in a similar manner and they are combined to give one set of linear algebraic equations that can be solved efficiently. The model is used to predict water level distributions across sloping beaches, where the water table in the aquifer may or may not decouple from the free water surface. Five cases are used to test the model for simulating beach water table fluctuations induced by tides, with the model predictions being compared with existing analytical solutions and laboratory and field data published in the literature. The numerical model results show that the integrated model is capable of simulating the combined ground‐ and surface‐water flows in coastal areas. Detailed analysis is undertaken to investigate the capability of the model. Copyright © 2009 John Wiley & Sons, Ltd.