Hydro-ecological modelers often use spatial variation of soil information derived from conventional soil surveys in simulation of hydro-ecological processes over watersheds at mesoscale (10–100 km2). Conventional soil surveys are not designed to provide the same level of spatial detail as terrain and vegetation inputs derived from digital terrain analysis and remote sensing techniques. Soil property layers derived from conventional soil surveys are often incompatible with detailed terrain and remotely sensed data due to their difference in scales. The objective of this research is to examine the effect of scale incompatibility between soil information and the detailed digital terrain data and remotely sensed information by comparing simulations of watershed processes based on the conventional soil map and those simulations based on detailed soil information across different simulation scales. The detailed soil spatial information was derived using a GIS (geographical information system), expert knowledge, and fuzzy logic based predictive mapping approach (Soil Land Inference Model, SoLIM). The Regional Hydro-Ecological Simulation System (RHESSys) is used to simulate two watershed processes: net photosynthesis and stream flow. The difference between simulation based on the conventional soil map and that based on the detailed predictive soil map at a given simulation scale is perceived to be the effect of scale incompatibility between conventional soil data and the rest of the (more detailed) data layers at that scale. Two modeling approaches were taken in this study: the lumped parameter approach and the distributed parameter approach. The results over two small watersheds indicate that the effect does not necessarily always increase or decrease as the simulation scale becomes finer or coarser. For a given watershed there seems to be a fixed scale at which the effect is consistently low for the simulated processes with both the lumped parameter approach and the distributed parameter approach. 相似文献
The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens
floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so
the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer
modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled
using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following
hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of
sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward
the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The
model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head
values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered
despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge,
and the well-developed karst system of the region has a very high hydraulic conductivity.
Electronic Publication 相似文献
The Xisha Trough, located in the northwest of the South China Sea (SCS) mainly rifted 30 Ma ago, has been a failed rift since the cessation of the seafloor spreading of the NW subbasin. Based on the velocity–depth model along Profile OBH-4 across the Xisha Trough, a seven-layer density–depth model is used to estimate density structure for the profile. The relationship between seismic velocity and radiogenic heat production is used to estimate the vertical distribution of heat sources in the lower crust. The 2-D temperature field is calculated by applying a 2-D numerical solution of the heat conduction equation and the thermal lithosphere thickness is obtained from the basalt dry solidus (BDS). The rheology of the profile is estimated on the basis of frictional failure in the brittle regime and power-law steady-state creep in the ductile regime. Rheological model is constructed for a three-layer model involving a granitic upper crust, a quartz diorite lower crust and an olivine upper mantle. Gravity modeling supports basically the velocity–depth model. The Moho along Profile OBH-4 is of relatively high heat flow ranging from 46 to 60 mW/m2 and the Moho heat flow is higher in the trough than on the flanks. The depth of the “thermal” lithospheric lower boundary is about 54 km in the center, deepens toward two sides, and is about 75 km at the northern slope area and about 70 km at the southern Xisha–Zhongsha Block. Rheological calculation indicates that the two thinnest ductile layers in the crust and the thickest brittle layer in the uppermost mantle lie in the central region, showing that the Xisha Trough has been rheologically strengthened, which are mainly due to later thermal relaxation. In addition, the strengthening in rheology during rifting was not the main factor in hampering the breakup of the Xisha Trough. 相似文献
We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ1 towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ3 rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures. 相似文献
Magma mixing structures from the lava flow of Lesbos (Greece) are analyzed in three dimensions using a technique that, starting from the serial sections of rock cubes, allows the reconstruction of the spatial distribution of magmas inside rocks. Two main kinds of coexisting structures are observed: (i) “active regions” (AR) in which magmas mix intimately generating wide contact surfaces and (ii) “coherent regions” (CR) of more mafic magma that have a globular shape and do not show large deformations. The intensity of mingling is quantified by calculating both the interfacial area (IA) between interacting magmas and the fractal dimension of the reconstructed structures. Results show that the fractal dimension is linearly correlated with the logarithm of interfacial area allowing discrimination among different intensities of mingling.
The process of mingling of magmas is simulated using a three-dimensional chaotic dynamical system consisting of stretching and folding processes. The intensity of mingling is measured by calculating the interfacial area between interacting magmas and the fractal dimension, as for natural magma mixing structures. Results suggest that, as in the natural case, the fractal dimension is linearly correlated with the logarithm of the interfacial area allowing to conclude that magma mixing can be regarded as a chaotic process.
Since chemical exchange and physical dispersion of one magma inside another by stretching and folding are closely related, we performed coupled numerical simulations of chaotic advection and chemical diffusion in three dimensions. Our analysis reveals the occurrence in the same system of “active mixing regions” and “coherent regions” analogous to those observed in nature. We will show that the dynamic processes are able to generate magmas with wide spatial heterogeneity related to the occurrence of magmatic enclaves inside host rocks in both plutonic and volcanic environments. 相似文献