Exploring the potential linkages between oil‐field brine reinjection,crystalline basement permeability,and triggered seismicity for the Dagger Draw Oil field,southeastern New Mexico,USA, using hydrologic modeling

We used hydrologic models to explore the potential linkages between oil‐field brine reinjection and increases in earthquake frequency (up to M_d 3.26) in southeastern New Mexico and to assess different injection management scenarios aimed at reducing the risk of triggered seismicity. Our analysis focuses on saline water reinjection into the basal Ellenburger Group beneath the Dagger Draw Oil field, Permian Basin. Increased seismic frequency (>M_d 2) began in 2001, 5 years after peak injection, at an average depth of 11 km within the basement 15 km to the west of the reinjection wells. We considered several scenarios including assigning an effective or bulk permeability value to the crystalline basement, including a conductive fault zone surrounded by tighter crystalline basement rocks, and allowing permeability to decay with depth. We initially adopted a 7 m (0.07 MPa) head increase as the threshold for triggered seismicity. Only two scenarios produced excess heads of 7m five years after peak injection. In the first, a hydraulic diffusivity of 0.1 m² s^?1 was assigned to the crystalline basement. In the second, a hydraulic diffusivity of 0.3 m² s^?1 was assigned to a conductive fault zone. If we had considered a wider range of threshold excess heads to be between 1 and 60 m, then the range of acceptable hydraulic diffusivities would have increased (between 0.1–0.01 m² s^?1 and 1–0.1 m² s^?1 for the bulk and fault zone scenarios, respectively). A permeability–depth decay model would have also satisfied the 5‐year time lag criterion. We also tested several injection management scenarios including redistributing injection volumes between various wells and lowering the total volume of injected fluids. Scenarios that reduced computed excess heads by over 50% within the crystalline basement resulted from reducing the total volume of reinjected fluids by a factor of 2 or more.     

Wetland assessment and impact mitigation decision support framework for linear development projects: The Louis Riel Trail,Highway 11 North project,Saskatchewan, Canada

Next to agriculture, road development is one of the most significant sources of stress to wetlands in Prairie Canada. However, there currently exists limited guidance for incorporating direct, indirect, and induced effects to wetlands in impact assessment and mitigation planning for small and often routine developments, including access roads or highway improvement initiatives. Based on the Louis Riel Trail, Highway 11 North twinning project in Saskatchewan, Canada, this article demonstrates a methodological approach and decision support framework for assessing and managing direct, indirect, and induced effects to wetlands from linear developments. No regulatory‐based environmental assessment was required for the highway project; effects were deemed to be insignificant under current wetland mitigation practices. However, our results show that 1115 ha of potentially affected wetlands are located within a 500 m impact zone on either side of the proposed highway. More than 50 percent of these wetlands are seasonal, less than 1 ha in size, and typically not included in mitigation planning. An expert‐based multi‐criteria evaluation of impact and mitigation options for wetlands in the study area indicated “no net loss” as a planning priority, and a preference for a spatially ambitious mitigation plan focused on direct, indirect, and potentially induced impacts. In practice, however, mitigation is often restrictive, focused on mitigating only direct impacts within the project's right‐of‐way, in this case less than 50 ha of wetlands, resulting in the potential for significant net loss of wetland habitat and function. If the risk to wetlands is to be given due consideration in project planning and development for roads and road improvement initiatives, then structured assessment methods and decision support frameworks should be sensitive to the time and resource constraints of small projects and screening‐type assessments. This requires also that wetland mitigation policies are developed and implementation plans formulated as part of project planning and assessment initiatives for linear developments.     

Snow-tracking and GIS: using multiple species-environment models to determine optimal wildlife crossing sites and evaluate highway mitigation plans on the Trans-Canada Highway

Snow-tracking data were collected for cougars (  Felis concolor ), lynx (  Lynx canadensis ), martens ( Martes americana ) and wolves ( Canis lupus ) and combined with remotely sensed imagery in a geographic information system (GIS) to identify wildlife crossing sites on the Trans-Canada Highway in Banff National Park, Alberta. We used logistic regression to assess the dependent (species presence/absence) relative to measures of topography and vegetation. The exponent form of each logistic regression equation was used to predict crossing sites in a GIS, which were then contrasted with mitigation sites proposed by Parks Canada. We found that: (1) cougars were influenced positively by normalized difference vegetation index (NDVI); negatively by northness and distance to ruggedness, (2) lynx were influenced positively by wetness, greenness, rugged terrain, eastness and distance to rugged terrain; negatively by slope, (3) martens were related positively to wetness, elevation, eastness, and distance to rugged terrain; negatively to northness, (4) wolves were influenced positively by distance to ruggedness; negatively by brightness, elevation, eastness and terrain ruggedness. There were few sympatric crossing sites for all species, supporting the use of species-specific mitigation or wide structures that capture multiple species needs. Inconsistencies were observed between the crossing sites predicted in this study and the Parks Canada proposal. The usefulness of GIS and track data to enhance mitigation projects is illustrated.     

Permeability of the continental crust: dynamic variations inferred from seismicity and metamorphism

The variation of permeability with depth can be probed indirectly by various means, including hydrologic models that use geothermal data as constraints and the progress of metamorphic reactions driven by fluid flow. Geothermal and metamorphic data combine to indicate that mean permeability ( k ) of tectonically active continental crust decreases with depth ( z ) according to log  k  ≈ −14–3.2 log  z , where k is in m² and z in km. Other independently derived, crustal-scale k – z relations are generally similar to this power-law curve. Yet there is also substantial evidence for local-to-regional-scale, transient, permeability-generation events that entail permeabilities much higher than these mean k – z relations would suggest. Compilation of such data yields a fit to these elevated, transient values of log  k  ≈ −11.5–3.2 log  z , suggesting a functional form similar to that of tectonically active crust, but shifted to higher permeability at a given depth. In addition, it seems possible that, in the absence of active prograde metamorphism, permeability in the deeper crust will decay toward values below the mean k – z curves. Several lines of evidence suggest geologically rapid (years to 10³ years) decay of high-permeability transients toward background values. Crustal-scale k – z curves may reflect a dynamic competition between permeability creation by processes such as fluid sourcing and rock failure, and permeability destruction by processes such as compaction, hydrothermal alteration, and retrograde metamorphism.     

The distribution and morphometry of lakes and reservoirs in British Columbia: a provincial inventory

An inventory of provincial lakes and reservoirs has been developed to characterise and assess the distribution and morphometry of standing water bodies in British Columbia. In the province, there are over 241,500 lakes and reservoirs greater than 1,000 m ² in size. These water bodies cover 2.37 percent of the province area and contain an estimated 521 km ³ of water (312 km ³ in natural lakes and 209 km ³ in reservoirs). A hypsometric relation suitable for order-of-magnitude estimates of lake volume from lake area is presented. Based on the distribution and morphometric attributes of lakes, several distinctive limnologic regions were identified, including the northeastern Alberta Plateau (highest proportion of circular lakes), the southwestern Alberta Plateau (lowest lake density/coverage and highest proportion of irregularly shaped lakes) and the Milbanke Strandflat (highest density of lakes). Observed regional and scale-related patterns in lake distribution/morphometry appear to be largely related to geomorphic controls, particularly tectonic and glacial history. Large-scale hydrologic implications of these standing water bodies and potential ecosystem/water resource management applications of the provincial inventory are also discussed.     

HILLSLOPE TERRACETTES

Several systems of lines of steps or terracettes, commonly known as “cattle tracks” or “sheep tracks”, were observed on treeless, grass-covered hillslopes in the Kislovodsk area of the North Caucasus. Three types are distinguished: (1) anthropogenic, or man-made, tracks, which suggest purposeful human activity and are distinguished by flat-bottomed or concave walkways; (2) animal tracks, connecting grazing areas and distinguished by deeper and narrower walkways produced by hooves; (3) delapsive terracettes (from the Latin delapsus = slip down), produced by a creeplike downslope movement of unconsolidated slope material on top of the bedrock surface. In contrast to the man-made and animal tracks, which have walkways bare of vegetation and grassy side walls, the naturally produced pseudotracks are distinguished by grassy top surfaces and bare side walls. Downslope rates of movement of these natural terracettes are estimated at several tens of millimeters to several centimeters a year.     

Modelling the Lost City hydrothermal field: influence of topography and permeability structure

The Lost City hydrothermal field (LCHF) is hosted in serpentinite at the crest of the Atlantis Massif, an oceanic core complex close to the mid‐Atlantic Ridge. It is remarkable for its longevity and for venting low‐temperature (40–91°C) alkaline fluids rich in hydrogen and methane. IODP Hole U1309D, 5 km north of the LCHF, penetrated 1415 m of gabbroic rocks and contains a near‐conductive thermal gradient close to 100°C km^?1. This is remarkable so close to an active hydrothermal field. We present hydrothermal modelling using a topographic profile through the vent field and IODP site U1309. Long‐lived circulation with vent temperatures similar to the LCHF can be sustained at moderate permeabilities of 10^?14 to 10^?15 m² with a basal heatflow of 0.22 W m^?2. Seafloor topography is an important control, with vents tending to form and remain in higher topography. Models with a uniform permeability throughout the Massif cannot simultaneously maintain circulation at the LCHF and the near‐conductive gradient in the borehole, where permeabilities <10^?16 m² are required. A steeply dipping permeability discontinuity between the LCHF and the drill hole is required to stabilize venting at the summit of the massif by creating a lateral conductive boundary layer. The discontinuity needs to be close to the vent site, supporting previous inferences that high permeability is most likely produced by faulting related to the transform fault. Rapid increases in modelled fluid temperatures with depth beneath the vent agree with previous estimates of reaction temperature based on geochemical modelling.     

Effective gas permeability of Tight Gas Sandstones as a function of capillary pressure – a non‐steady‐state approach

Single‐ and two‐phase (gas/water) fluid transport in tight sandstones has been studied in a series of permeability tests on core plugs of nine tight sandstones of the southern North Sea. Absolute (Klinkenberg‐corrected) gas permeability coefficients (k_gas__inf) ranged between 3.8 × 10^?16 and 6.2 × 10^?19 m² and decreased with increasing confining pressure (10–30 MPa) by a factor 3–5. Klinkenberg‐corrected (intrinsic) gas permeability coefficients were consistently higher by factors from 1.4 to 10 than permeability coefficients determined with water. Non‐steady‐state two‐phase (He/water) flow experiments conducted up to differential pressures of 10 MPa document the dynamically changing conductivity for the gas phase, which is primarily capillary‐controlled (drainage and imbibition). Effective gas permeability coefficients in the two‐phase flow tests ranged between 1.1 × 10^?17 and 2.5 × 10^?22 m², corresponding to relative gas permeabilities of 0.03% and 10%. In the early phase of the nonstationary flow regime (before establishment of steady‐state conditions), they may be substantially (>50%) lower. Effective gas permeability measurements are affected by the following factors: (i) Capillary‐controlled drainage/imbibition, (ii) viscous–dynamic effects (iii) and slip flow.     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Effects of episodic fluid flow on hydrocarbon migration in the Newport‐Inglewood Fault Zone,Southern California

Fault permeability may vary through time due to tectonic deformations, transients in pore pressure and effective stress, and mineralization associated with water‐rock reactions. Time‐varying permeability will affect subsurface fluid migration rates and patterns of petroleum accumulation in densely faulted sedimentary basins such as those associated with the borderland basins of Southern California. This study explores the petroleum fluid dynamics of this migration. As a multiphase flow and petroleum migration case study on the role of faults, computational models for both episodic and continuous hydrocarbon migration are constructed to investigate large‐scale fluid flow and petroleum accumulation along a northern section of the Newport‐Inglewood fault zone in the Los Angeles basin, Southern California. The numerical code solves the governing equations for oil, water, and heat transport in heterogeneous and anisotropic geologic cross sections but neglects flow in the third dimension for practical applications. Our numerical results suggest that fault permeability and fluid pressure fluctuations are crucial factors for distributing hydrocarbon accumulations associated with fault zones, and they also play important roles in controlling the geologic timing for reservoir filling. Episodic flow appears to enhance hydrocarbon accumulation more strongly by enabling stepwise build‐up in oil saturation in adjacent sedimentary formations due to temporally high pore pressure and high permeability caused by periodic fault rupture. Under assumptions that fault permeability fluctuate within the range of 1–1000 millidarcys (10^?15–10^?12 m²) and fault pressures fluctuate within 10–80% of overpressure ratio, the estimated oil volume in the Inglewood oil field (approximately 450 million barrels oil equivalent) can be accumulated in about 24 000 years, assuming a seismically induced fluid flow event occurs every 2000 years. This episodic petroleum migration model could be more geologically important than a continuous‐flow model, when considering the observed patterns of hydrocarbons and seismically active tectonic setting of the Los Angeles basin.     

Spatial variations in the salinity of pore waters in northern deep water Gulf of Mexico sediments: implications for pathways and mechanisms of solute transport

Spatial variations in the salinity of pore waters in sedimentary basins can provide important insight into basin-scale hydrogeologic processes. Although there have been numerous studies of brine seeps in the deep water Gulf of Mexico, much less is known about porewater salinities in the vast areas between seeps. A study has been made of spatial variation in pore water salinities in sediments in an approximately 500 km by 200 km area of the northern deep water (water depth >500 m) Gulf of Mexico sedimentary basin (GOM) to provide insight into pathways and mechanisms of solute transport in this portion of the basin. A second objective was to document salinities in the upper 500 m of the sedimentary section, the approximate depth to which methane hydrates, a potential future energy resource, may be stable. Elevated salinities would reduce the P – T stability range of hydrates. Salinities were calculated from borehole logs using a dual electrical conductivity model. Even though much of the northern GOM is underlain by allochthonous salt most of the undisturbed shallow sedimentary section has not been permeated by hypersaline waters. These waters are limited to areas near brine seeps. Hypersaline waters having salinities in excess of 100 g l⁻¹ become more common at subseafloor depths of 2 km and greater. A field study at Green Canyon 65 and published numerical simulations of fluid flow above tabular salt bodies suggest that brines produced by salt dissolution migrate laterally and pond above salt and/or within minibasins and that the dominant mechanism of vertical solute transport is a combination of compaction-driven advection and diffusion, not large-scale thermohaline overturn. Superimposed on this diffuse upward flux of dissolved salt is the more focused and localized expulsion of saline fluids up faults.     

Hiding to eat: the role of carnivores in the early Middle Pleistocene from the TD8 level of Gran Dolina (Sierra de Atapuerca,Burgos, Spain)

Damage generated by large and small carnivores is common in many Middle Pleistocene sites. However, identifying the predator that produces the faunal accumulations is often a difficult task. In order to recognize the main type of carnivore that acts on a faunal assemblage, a combination of several characteristics should be taken into account: taxonomic and skeletal element representation, age profiles, carnivore damage (location, frequencies and dimensions of tooth-marks, bone breakage and digested bones), degree of fragmentation and frequencies of coprolites. But, adding environmental characteristics and the ethology of non-human predators/scavengers is also important. All these aspects are applied to the faunal assemblage from the TD8 level of the Gran Dolina site (Sierra de Atapuerca, Burgos, Spain). Paleomagnetic data combined with ESR and U-series place the TD8 level at the beginning of the Middle Pleistocene, specifically circa 700 kyr ago. The TD8 level contains a large faunal accumulation primarily composed of ungulate skeletal elements, and to a lesser extent carnivore remains. This assemblage is characterized by an overrepresentation of fallow deer (Dama vallonetensis), a skeletal profile biased towards cranial remains and limb bones, diversity of ages at death, a high proportion of carnivore damage and tooth-marks of large size, and an absence of human activity. According to these data, the accumulation seems to have been produced primarily by large carnivores, possibly hyenas. This observation does not rule out the possible occasional activity by other carnivores. Nevertheless, the characteristics of the TD8 assemblage do not correlate entirely with those traditionally used to define carnivore dens. In TD8, there are (1) no immature carnivore remains (remains of just one young Mosbach wolf); (2) scarce traces related to the end stages of consumption and some anatomical connexions; (3) few coprolites; (4) high proportion of adult ungulates and; (5) high quantities of whole bones and epiphyses. From this perspective, the TD8 faunal assemblage seems to correspond to a succession of carnivore occupations that allows the development of a suite of features to identify the activities of several species of predators that may have used the cave in different ways and durations. This study aims to emphasize the importance of these analyzes in order to know the behaviour of different non-human predators/scavengers in the European Middle Pleistocene sites.     

Reconsideration of criteria for differentiating faunal assemblages accumulated by hyenas and hominids

Zooarchaeologists have established several criteria for differentiating hominid‐ and hyena‐derived faunal assemblages. In some cases, however, the patterns of skeletal part representation and bone surface modification on which these criteria are ultimately based have been observed in fossil bone assemblages of unknown origin, rather than in modern assemblages of known origin. When the proposed criteria are evaluated within an actualistic framework, only three are able to differentiate between hominid‐ and hyena‐created faunal assemblages. I suggest that only these three criteria—proportions of carnivores to ungulates in the assemblage, the preserved condition of long bone specimens (either as whole cylinders or as splintered shaft fragments), and the types of bone surface modifications—should be retained as important factors in a diagnosis of the ancient bone‐collector. The remaining four criteria—the relative proportion of horn pieces in the assemblage, the relative representation of podial bones, the relative representation of small and large bovid skeletal parts, and bovid mortality profiles—are not relevant or applicable to the problem of differentiating hominid‐ from hyena‐derived faunal assemblages. Copyright © 2002 John Wiley & Sons, Ltd.     

Asphalt Dreams: Road Construction and Environmental Citizenship in Peru

Peru today is one of the main staging grounds for a continent‐wide integration effort. Launched in 2000, the Initiative for the Integration of Regional Infrastructure in South America (IIRSA) calls for an enormous expansion of the continent's transport and energy networks and an effort to increase the region's economic competitiveness. Among its most controversial projects is the Interoceanic Highway linking western Brazil with the Pacific coast of Peru. The highway has attracted fierce criticism from NGOs who point to major environmental impacts, an inadequate mitigation process, and a lack of transparency in funding flows and decision making. In an effort to voice their concerns, these groups engage the idea of ‘environmental governance’ to increase public participation in the development process and promote ecological sustainability. This alternative framework in turn opens up space for ‘environmental citizenship’. This article takes a closer look at how Peruvian NGOs employ this idea and suggests that while the group's advocacy of governance has had success, the building of environmental citizenship will require a move beyond urban Peruvian NGOs as technical experts.     

Rates of retrograde metamorphism and their implications for crustal rheology

The rate of reaction of a natural hornblende garnet granulite with water under a range of mid- to lower crustal conditions has been investigated experimentally. In runs of between 7 and 84 days small but measurable amounts of water were consumed, and sheet silicates (300°C, 300 MPa and 400°C, 400 MPa) and/or secondary actinolite (400°C, 400 MPa and 500°C, 500 MPa) were observed to have grown. When normalized to the surface area of the starting materials, hydration rates were in the range of 2–5 × 10⁻⁸ g m⁻² sec⁻¹. These reaction rates imply that a film of water that infiltrated a planar crack with a half width of 100 μm would be completely consumed within c. 100 years. These results imply that where water infiltrates the crust along faults or underlying shear zones in response to a deformation, it will remain as a free phase for only a finite period of time, which in some cases will be less than the repeat time for major earthquakes in the fault system. Thus, the rheology of fault zones and shear zones is likely to be cyclical, with the zone becoming stronger with time as water is consumed, and then weakened by infiltration of water after each rupture.     

Permeability changes in simulated granite faults during and after frictional sliding

Thermal–hydrological–mechanical coupling processes suggest that fault permeability should undergo dynamic change as a result of seismic slip. In igneous rocks, a fault's slip surface may have much higher permeability than the surrounding rock matrix and therefore operate as a conduit for fluids. We conducted laboratory experiments to investigate changes in fracture permeability (or transmissivity) of a fault in granite due to shear slip and cyclic heating and cooling. Our experiments showed that high initial fracture transmissivity (>10^?18 m³) was associated with a high friction coefficient and that transmissivity decreased during slip. We propose that this reduction in transmissivity reflects the presence of gouge in fracture voids, increasing the area of contact in the fault plane and reducing the hydraulic aperture. In contrast, when initial fracture transmissivity was low (<10^?18 m³), we observed that friction was lower and transmissivity increased during slip. The high transmissivity and high friction may be explained by large areas of bare rock being in contact on the slip surface. Slip velocity had little influence on the evolution of permeability, probably because gouge produced at different slip velocities had similar grain size distributions, or because gouge leaked from the slip surface. Transmissivity decreased with increasing temperature in heating tests, probably due to thermal expansion increasing normal stress on the fracture. Frictional heating did not influence transmissivity during the shearing tests.     

The effect of deformation on permeability development in anhydrite and implications for caprock integrity during geological storage of CO2

Geological storage of CO₂ in depleted oil and gas reservoirs is one of the most promising options to reduce atmospheric CO₂ concentrations. Of great importance to CO₂ mitigation strategies is maintaining caprock integrity. Worldwide many current injection sites and potential storage sites are overlain by anhydrite‐bearing seal formations. However, little is known about the magnitude of the permeability change accompanying dilatation and failure of anhydrite under reservoir conditions. To this extent, we have performed triaxial compression experiments together with argon gas permeability measurements on Zechstein anhydrite, which caps many potential CO₂ storage sites in the Netherlands. Our experiments were performed at room temperature at confining pressures of 3.5–25 MPa. We observed a transition from brittle to semi‐brittle behaviour over the experimental range, and peak strength could be described by a Mogi‐type failure envelope. Dynamic permeability measurements showed a change from ‘impermeable’ (<10^?21 m²) to permeable (10^?16 to 10^?19 m²) as a result of mechanical damage. The onset of measurable permeability was associated with an increase in the rate of dilatation at low pressures (3.5–5 MPa), and with the turning point from compaction to dilatation in the volumetric versus axial strain curve at higher pressures (10–25 MPa). Sample permeability was largely controlled by the permeability of the shear faults developed. Static, postfailure permeability decreased with increasing effective mean stress. Our results demonstrated that caprock integrity will not be compromised by mechanical damage and permeability development. Geofluids (2010) 10 , 369–387     

Silver phiale said to be from al‐Juba (al‐Wusṭa Governorate) — an archaeological puzzle

In March 2014, while recording finds in the Ministry of Heritage and Culture in al‐Khuwair, artefacts unexpectedly came to light reportedly from al‐Juba in Oman's Bar al‐?ikmān, in al‐Wus?a Governorate, until recently an archaeologically little‐researched part of the Sultanate. Some of the pieces could be attributed to the Samad LIA or perhaps the PIR, both from the centuries at the turning point of the ages from BCE to CE. Such finds have never before come to light in this part of Oman. Samad LIA sites are generally located some 220 km to the north on the southern flank of the al‐?ajjar mountains in a zone c.160 x 105 km in area. Diagnostic pottery finds spread from Wadi Bānī Ruwāhah (UTM 40Q 620570 m E, 2561848 m N) eastwards to the coast, a smaller area than previously believed.     

Analytical and numerical models of hydrothermal fluid flow at fault intersections

Fault intersections are the locus of hot spring activity and Carlin‐type gold mineralization within the Basin and Range, USA. Analytical and numerical solutions to Stokes equation suggest that peak fluid velocities at fault intersections increase between 20% and 47% when fracture apertures have identical widths but increase by only about 1% and 8% when aperture widths vary by a factor of 2. This suggests that fault zone intersections must have enlarged apertures. Three‐dimensional finite element models that consider intersecting 10‐ to 20‐m wide fault planes resulted in hot spring activity being preferentially located at fault zone intersections when fault zones were assigned identical permeabilities. We found that the onset of convection at the intersections of the fault zones occurred in our hydrothermal model over a narrow permeability range between 5 × 10^?13 and 7 × 10^?13 m². Relatively high vertical fluid velocities (0.3–3 m year^?1) extended away from the fault intersections for about 0.5–1.5 km. For the boundary conditions and fault plane dimensions used, peak discharge temperatures of 112°C at the water table occurred with an intermediate fault zone permeability of 5 × 10^?13 m². When fault plane permeability differed by a factor of 2 or more, the locus of hot spring activity shifted away from the intersections. However, increasing the permeability at the core of the fault plane intersection by 40% shifted the discharge back to the intersections. When aquifer units were assigned a permeability value equal to those of the fault planes, convective rolls developed that extend about 3 km laterally along the fault plane and into the adjacent aquifer.     

Dietary ecology of extant guanaco (Lama guanicoe) from Southern Patagonia: seasonal leaf browsing and its archaeological implications

The analysis of dietary traits of ungulates through tooth microwear and mesowear has been applied to archaeological sites to investigate seasonal changes in settlements by hunter–gatherers. In this paper we propose to test the hypothesis that tooth microwear (combined to mesowear) is able to indicate seasonality in the diet of extant ungulates in arid habitats (semi-deserts or steppe). The material analyzed comes from six faunal monospecific assemblages of guanaco (Lama guanicoe) resulting from a mass mortality event in winter 2000 near the Cardiel Lake in Southern Patagonia (Province of Santa Cruz, Argentina). Mesowear results indicate that the guanacos from the Cardiel Lake area are mixed feeders, and thus, have a diet that shifts seasonally. Moreover, microwear analysis supports the hypothesis that tooth microwear is able to indicate seasonality in the diet of extant guanaco in arid habitats. The pattern is clear for the winter sample and needs to be confirmed for a summer sample. Consequently, tooth microwear is proposed as a new potential proxy for detecting seasonal occupation in archaeological sites in Patagonia and other arid environments.