Stochastic ground water flow simulation with a fracture zone continuum model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.     

Multi-criteria modelling and clustering of spatial information

To improve the formal integration of verbally given spatial information in a geographical information system, a methodology for converting such information to a geographical representation was developed. This was applied to examples of a spatial-related analysis of natural language text messages in the domain of disaster management. This article presents an approach for such a conversion by using a suitable knowledge representation as well as formal modelling structures. The structures are provided by an ontology-supported knowledge base. With respect to spatial uncertainty, a formal representation based on possibility theory was defined. Moreover, a cluster algorithm was developed for handling information of different sources and building a spatial context.     

Post-depositional impacts on ‘Findlinge’ (erratic boulders) and their implications for surface-exposure dating

Understanding and interpretation of ‘numbers’ produced about the depositional age of an erratic boulder by cosmogenic nuclide surface-exposure dating is important in the construction of glacial chronology. We have sampled three ‘Findlinge’ (glacially transported boulders) located on the right-lateral margin of the Aare glacier at Möschberg, Grosshöchstetten, southeast of Bern, with the aim of shedding light on this topic. The boulders have the same depositional, but different post-depositional histories: simple exposure; exhumation; and human impact. This sampling is specially selected for this study, since the boulders showing exhumation and human impact would not have been sampled in a regular surface-exposure dating application. We measured cosmogenic ¹⁰Be concentrations and calculated apparent exposure ages that are 13.6 ± 0.5, 18.1 ± 0.8, and 7.5 ± 0.4 ka, respectively. The exposure age of the first boulder reflects exhumation. The apparent exposure age of 18.1 ± 0.8 ka (erosion-corrected exposure age 19.0 ± 0.9 ka) from the second boulder correlates well with the end of the Alpine and global last glacial maximum. The third boulder shows evidence of quarrying as it is surrounded by a rim of excavation material, which is also reflected by the 7.5 ± 0.4 ka apparent exposure age. We modeled the variation of ¹⁰Be concentrations with depth down into the sediment in which the first (exhumed) boulder was once buried in, and down into the third (quarried) boulder. According to our modeling, we determined that the exhumed ‘Findling’ was buried in sediment at a depth of around 0.5 m, and around 2 m of rock was quarried from the third ‘Findling’. Our results reveal the importance of sampling for surface-exposure dating within a well defined field context, as post-depositional impacts can easily hinder exposure-dating of surfaces.     

南极威德尔海海冰晶体类型分析

为了解南极威德尔海不同区域海冰形成过程,分析了2006年9—10月威德尔海海冰调查中的27支冰芯晶体并描述了冰芯反映的冰层形成过程;给出2支完整冰芯晶体照片和多边形晶体雪冰照片,以及对应照片的冰层描述和冰盐度、密度、粒径。根据冰芯晶体揭示的海冰形成过程,发现尽管现场考察均选择较大表面平整的浮冰作业,但冰芯晶体揭示出实际冰层可以是受动力和热力联合作用形成的重叠冰、固结冰脊、二年冰;而从表面到底面属于纯热力学生长形成的冰层所占比例较小。积雪经过密实和水分参与,再冻结形成多边形晶体雪冰也是南极海冰主要组成之一。粒状晶体、过渡区混合晶体和柱状晶体海冰分别占27支冰芯总长度的28.7%、14.4%和55.2%。海冰的区域分布表现出海冰边缘区以纯热力学生长的一年冰为主;在冰架前缘,受动力作用形成的重叠冰和固结冰脊的比例增加,并且存在二年冰和多边形晶体雪冰;冰间湖内生长有大面积纯热力学生长的海冰并向外输送。     

Climatology of ionospheric upper transition height derived from COSMIC satellites during the solar minimum of 2008

The electron density profiles retrieved from the Constellation Observing System for Meteorology, Ionosphere & Climate (COSMIC) satellite Radio Occultation (RO) observations during 2008 are used to derive ionospheric upper transition height, where the density of O⁺ is equal to that of light ions (mainly H⁺ and He⁺). It is found that the ionosphere upper transition height is very low, with significant local time, latitude and seasonal variations, during the extremely low solar minimum of 2008. The transition height is higher in the daytime than at night, except over middle latitude region of winter hemisphere, where the transition height has minimum in the morning. There is a pronounced peak over equator for all seasons. The transition height is higher in summer than in winter hemisphere. Our results have comparability with C/NOFS satellite observations around the equatorial region during June–August of 2008. However, the IRI model gives much higher transition height than those from COSMIC and cannot reproduce its latitude and season variations well during 2008.     

A hybrid approach to compute seismic travel times in three-dimensional tetrahedral meshes

We propose an optimized method to compute travel times for seismic inversion problems. It is a hybrid method combining several approaches to deal with travel time computation accuracy in unstructured meshes based on tetrahedral elementary cells. As in the linear travel time interpolation method, the proposed approach computes travel times using seismic ray paths. The method operates in two sequential steps: At a first stage, travel times are computed for all nodes of the mesh using a modified version of the shortest path method. The difference with the standard version is that additional secondary nodes (called tertiary nodes) are added temporarily around seismic sources in order to improve accuracy with a reasonable increase in computational cost. During the second step, the steepest travel time gradient method is used to trace back ray paths for each source–receiver pair. Travel times at each receiver are then recomputed using slowness values at the intersection points between the ray path and the traversed cells. A number of numerical tests with an array of different velocity models, mesh resolutions and mesh topologies have been carried out. These tests showed that an average relative error in the order of 0.1% can be achieved at a computational cost that is suitable for travel time inversion.     

Urban tourism and small tourism enterprise development in Johannesburg: The case of township tourism

Urban tourism is a focus of vibrant research in developed countries. In this paper the focus is upon an element of urban tourism in a developing world city, viz., Johannesburg. Specific attention centers on an analysis of the role of small enterprise development in urban tourism. The major issues surrounding the development of small tourism businesses in two of Johannesburg's former black townships, namely Soweto and Alexandra, which are emerging as centers for so-termed ‘township tourism’ in South Africa, are investigated. This revised version was published online in August 2006 with corrections to the Cover Date.     

Shocked rocks and impact glasses from the El'gygytgyn impact structure,Russia

Abstract— The El'gygytgyn impact structure is about 18 km in diameter and is located in the central part of Chukotka, arctic Russia. The crater was formed in volcanic rock strata of Cretaceous age, which include lava and tuffs of rhyolites, dacites, and andesites. A mid‐Pliocene age of the crater was previously determined by fission track (3.45 ± 0.15 Ma) and ⁴⁰Ar/³⁹Ar dating (3.58 ± 0.04 Ma). The ejecta layer around the crater is completely eroded. Shock‐metamorphosed volcanic rocks, impact melt rocks, and bomb‐shaped impact glasses occur in lacustrine terraces but have been redeposited after the impact event. Clasts of volcanic rocks, which range in composition from rhyolite to dacite, represent all stages of shock metamorphism, including selective melting and formation of homogeneous impact melt. Four stages of shocked volcanic rocks were identified: stage I (≤35 GPa; lava and tuff contain weakly to strongly shocked quartz and feldspar clasts with abundant PFs and PDFs; coesite and stishovite occur as well), stage II (35–45 GPa; quartz and feldspar are converted to diaplectic glass; coesite but no stishovite), stage III (45–55 GPa; partly melted volcanic rocks; common diaplectic quartz glass; feldspar is melted), and stage IV (>55 GPa; melt rocks and glasses). Two main types of impact melt rocks occur in the crater: 1) impact melt rocks and impact melt breccias (containing abundant fragments of shocked volcanic rocks) that were probably derived from (now eroded) impact melt flows on the crater walls, and 2) aerodynamically shaped impact melt glass “bombs” composed of homogeneous glass. The composition of the glasses is almost identical to that of rhyolites from the uppermost part of the target. Cobalt, Ni, and Ir abundances in the impact glasses and melt rocks are not or only slightly enriched compared to the volcanic target rocks; only the Cr abundances show a distinct enrichment, which points toward an achondritic projectile. However, the present data do not allow one to unambiguously identify a meteoritic component in the El'gygytgyn impact melt rocks.     

Texture and anisotropy analysis of Qusaiba shales

Scanning and transmission electron microscopy, synchrotron X‐ray diffraction, microtomography and ultrasonic velocity measurements were used to characterize microstructures and anisotropy of three deeply buried Qusaiba shales from the Rub’al‐Khali basin, Saudi Arabia. Kaolinite, illite‐smectite, illite‐mica and chlorite show strong preferred orientation with (001) pole figure maxima perpendicular to the bedding plane ranging from 2.4–6.8 multiples of a random distribution (m.r.d.). Quartz, feldspars and pyrite crystals have a random orientation distribution. Elastic properties of the polyphase aggregate are calculated by averaging the single crystal elastic properties over the orientation distribution, assuming a nonporous material. The average calculated bulk P‐wave velocities are 6.2 km/s (maximum) and 5.5 km/s (minimum), resulting in a P‐wave anisotropy of 12%. The calculated velocities are compared with those determined from ultrasonic velocity measurements on a similar sample. In the ultrasonic experiment, which measures the effects of the shale matrix as well as the effects of porosity, velocities are smaller (P‐wave maximum 5.3 km/s and minimum 4.1 km/s). The difference between calculated and measured velocities is attributed to the effects of anisotropic pore structure and to microfractures present in the sample, which have not been taken into account in the matrix averaging.