Roman Builders Facing the Risk of Disaster: Coupling Archaeological and Geochemical Analyses on a Section of the ‘Aqua Augusta’ (the Roman Serino Aqueduct,Southern Italy)

Four samples of carbonate sinters from two sites of the Roman Serino aqueduct were studied in order to better constrain the history of the aqueduct and to estimate to what extent natural hazards can be recorded in such deposits. Micromorphological observations, trace element, stable oxygen and carbon isotope analyses were performed on the samples. Together with new quantitative observations on the construction techniques of the aqueduct, our results highlight the imprint left by geological hazards on the Serino aqueduct. Damage, disruption, ground movement and a debris‐flow are evidenced.     

HIGH‐RESOLUTION (PIXE) ANALYSES OF CARBONATE DEPOSITS IN A ROMAN AQUEDUCT (FRÉJUS,SE FRANCE): IMPLICATIONS FOR THE STUDY OF PALAEOHYDROLOGICAL VARIABILITY AND WATER RESOURCES MANAGEMENT IN SOUTHERN GAUL DURING THE ROMAN PERIOD*

Urban water was supplied to the Roman city of Forum Julii (Fréjus, southeastern France) for at least 200 years by a 39.4 km long aqueduct, operating in ad 50. Two perennial springs were successively collected, the Foux and the Siagnole, located at the outlet of Triassic and Jurassic karstic reservoirs, respectively. In this study, we performed high‐resolution PIXE (Particle‐Induced X‐ray Emission) measurements of Ca, Si, Fe and Sr concentrations on selected cross‐sections of laminated carbonate sampled along the sidewalls, before and after the connection of the two collection channels. Seasonal variations of water composition, suspended sediment load and discharge are recorded by alternating clear and dark layers, and can be traced by their Sr and Fe contents. On the basis of an annual bimodal high‐discharge regime for the two karstic reservoirs, the concentration measurements allow the derivation of a chronological record (117.5 years) of water supply and maintenance activities for a part of the operational period of the aqueduct. The water level in the channel was mainly controlled by the extent of carbonate deposition and by dredging, repair and maintenance operations rather than by the past hydrological regime of the two springs.     

PALAEOCLIMATIC IDENTIFICATION BASED ON AN ISOTOPE STUDY OF TRAVERTINE FROM THE COPPER AGE SITE AT LOS MILLARES,SOUTH-EASTERN SPAIN*

Remains of an aqueduct from the Copper Age settlement of Los Millares have been studied. The isotope results (—7‰ < δ¹⁸O — 5.3‰) show that the aqueduct was used to transport fresh water at room temperature. The results are inconsistent with the aqueduct being used for hot water transport or as a container for water exposed to evaporation processes. Both isotopic and mineralogical results were useful in identifying the origin and use of building materials and also in determining climatic conditions during the time of the settlement occupation.     

Pore‐scale investigation on stress‐dependent characteristics of granular packs and the impact of pore deformation on fluid distribution

Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore‐scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore‐scale experiments of single‐ and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction‐dependent characteristics of granular packs and their impact on fluid flow. A series of drainage and imbibition cycles were conducted on a water‐wet, soda‐lime glass bead pack under varying confining stress conditions. Simultaneously, X‐ray micro‐CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro‐CT images were segmented using a gradient‐based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. However, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.     

Temporal Patterns and Effects of Surface‐Water Diversions on Daily Flows and Aquatic Habitats: Bega‐Bemboka River,New South Wales,Australia

An understanding of the nature and magnitude of hydrological, physical habitat and physico‐chemical effects resulting from surface‐water diversions in river systems is essential for effective management of water resources. For most coastal‐draining rivers in New South Wales, however, there are few data available on irrigation diversions, their hydrological impacts and environmental effects. This paper therefore presents an analysis of mean daily surface‐water diversions for pasture irrigation from approximately three years of metered data and the resultant effects on daily flows and aquatic habitats in the Bega‐Bemboka River. The period of analysis of the hydrological effects of irrigation diversions is extended to the full length of record (approximately five years) for gauging stations most affected by irrigation diversions, using Maintenance of Variance Extension Type 1 (MOVE.1) modelling techniques. The annual mean, median and peak daily rates of water diversion by metered surface‐water licences used to irrigate 925 ha of dairy pasture are 10.5 Ml d^?1, 7.3 Ml d^?1 and 41.5 Ml d^?1, respectively. Diversion effects on flow duration statistics are such that the measured 90th and 95th daily flow duration percentiles at the gauging station most affected by upstream irrigation diversions are equivalent to the 97th and 99th flow duration percentiles, respectively, under MOVE.1 modelled natural flow conditions. While diversions for irrigation over the three‐year data period account for only 6.6% of total flow volumes, diversions as a proportion of daily surface‐water inflows increase exponentially under decreasing flow rates. Median and maximum daily diversion rates attain 91% and 118%, respectively, of total surface‐water inflows to the diversion‐affected reach when upstream inflows range from 15 to 20 Ml d^?1. This exponentially‐increasing relationship between daily diversion rates and declining surface‐water inflows suggests that ‘rule of thumb’ guidelines on sustainable diversion limits based on mean or median annual percentage diversion volumes need to be applied cautiously to river systems with no or limited capacity to manipulate flows to meet downstream consumptive demands.     

Construction and Interpretation of a 60‐year Flow Record for the river Tay at Kenmore

The backward extension of a river flow record for the River Tay, based on daily water level observations from 1937–76, has created the third‐longest sensibly continuous daily flow record for any Scottish river based on direct observations. The paper describes the derivation of the record, and assesses the variability of the Tay discharge over the 60‐year period covered. Unprecedented winter, spring and flood flows in the late 1980s/1990s are reported, in the context of improving the understanding of hydroclimatic variability.     

Equations of state for basin geofluids: algorithm review and intercomparison for brines

Physical properties of formation waters in sedimentary basins can vary by more than 25% for density and by one order of magnitude for viscosity. Density differences may enhance or retard flow driven by other mechanisms and can initiate buoyancy‐driven flow. For a given driving force, the flow rate and injectivity depend on viscosity and permeability. Thus, variations in the density and viscosity of formation waters may have or had a significant effect on the flow pattern in a sedimentary basin, with consequences for various basin processes. Therefore, it is critical to correctly estimate water properties at formation conditions for proper representation and interpretation of present flow systems, and for numerical simulations of basin evolution, hydrocarbon migration, ore genesis, and fate of injected fluids in sedimentary basins. Algorithms published over the years to calculate water density and viscosity as a function of temperature, pressure and salinity are based on empirical fitting of laboratory‐measured properties of predominantly NaCl solutions, but also field brines. A review and comparison of various algorithms are presented here, both in terms of applicability range and estimates of density and viscosity. The paucity of measured formation‐water properties at in situ conditions hinders a definitive conclusion regarding the validity of any of these algorithms. However, the comparison indicates the versatility of the various algorithms in various ranges of conditions found in sedimentary basins. The applicability of these algorithms to the density of formation waters in the Alberta Basin is also examined using a high‐quality database of 4854 water analyses. Consideration is also given to the percentage of cations that are heavier than Na in the waters.     

Transient groundwater impacts on the development of paleoclimatic lake records in semi-arid environments

Lacustrine sedimentary deposits are a primary archive of terrestrial paleoclimatic information, and groundwater fluxes under arid climates can dominate water budgets in lake and playa environments. In this paper, we investigate the effects of basin‐scale water table fluctuations on lacustrine paleoclimatic records in arid regions using the hydrogeologic model MWT3D. Boundary conditions for paleo‐recharge over the past 57 ka, land surface topography, the location of lakes, and subsurface geology were based on information available for a 600‐km cross section across the western Murray Basin, Australia. Results indicate the importance of the physical interaction between the water table and the land surface in controlling the location and the timing of formation of paleoclimatic lacustrine records. Lakes in different positions across the watershed preserved different records of climate (wet versus dry), although the same recharge history was applied. Changes in groundwater seepage to the surface lagged behind the climatic signal between 500 and 220 years across the basin. The size of the flow system, specifically the distance from the local discharge point to the adjacent groundwater divide, was the primary control on the time lag. With large fluctuations in water table elevation (>50 m) this distance changes considerably from wet to dry climatic conditions. Our results are largely consistent with observed paleoclimatic records derived from lakes across the Murray Basin.     

Numerical modeling of the origin of calcite mineralization in the Refugio‐Carneros fault,Santa Barbara Basin,California

Many faults in active and exhumed hydrocarbon‐generating basins are characterized by thick deposits of carbonate fault cement of limited vertical and horizontal extent. Based on fluid inclusion and stable isotope characteristics, these deposits have been attributed to upward flow of formation water and hydrocarbons. The present study sought to test this hypothesis by using numerical reactive transport modeling to investigate the origin of calcite cements in the Refugio‐Carneros fault located on the northern flank of the Santa Barbara Basin of southern California. Previous research has shown this calcite to have low δ¹³C values of about ?40 to ?30‰PDB, suggesting that methane‐rich fluids ascended the fault and contributed carbon for the mineralization. Fluid inclusion homogenization temperatures of 80–125°C in the calcite indicate that the fluids also transported significant quantities of heat. Fluid inclusion salinities ranging from fresh water to seawater values and the proximity of the Refugio‐Carneros fault to a zone of groundwater recharge in the Santa Ynez Mountains suggest that calcite precipitation in the fault may have been induced by the oxidation of methane‐rich basinal fluids by infiltrating meteoric fluids descending steeply dipping sedimentary layers on the northern basin flank. This oxidation could have occurred via at least two different mixing scenarios. In the first, overpressures in the central part of the basin may have driven methane‐rich formation waters derived from the Monterey Formation northward toward the basin flanks where they mixed with meteoric water descending from the Santa Ynez Mountains and diverted upward through the Refugio‐Carneros fault. In the second scenario, methane‐rich fluids sourced from deeper Paleogene sediments would have been driven upward by overpressures generated in the fault zones because of deformation, pressure solution, and flow, and released during fault rupture, ultimately mixing with meteoric water at shallow depth. The models in the present study were designed to test this second scenario, and show that in order for the observed fluid inclusion temperatures to be reached within 200 m of the surface, moderate overpressures and high permeabilities were required in the fault zone. Sudden release of overpressure may have been triggered by earthquakes and led to transient pulses of accelerated fluid flow and heat transport along faults, most likely on the order of tens to hundreds of years in duration. While the models also showed that methane‐rich fluids ascending the Refugio‐Carneros fault could be oxidized by meteoric water traversing the Vaqueros Sandstone to form calcite, they raised doubts about whether the length of time and the number of fault pulses needed for mineralization by the fault overpressuring mechanism were too high given existing geologic constraints.     

The Conception of Hull Shape by Shell‐builders in the Ancient Mediterranean

Did ancient Mediterranean shipbuilders at the turn of the first millennium have methods to pre‐design the shape of their hulls prior to the construction of the ship? If so, can these methods be identified? This paper proposes that ancient shipwrights were in possession of such design methods, based on geometry, that are comparable to the medieval and post medieval hull‐design procedures called whole‐moulding. These methods are well suited to be transformed into rules‐of‐thumb, therefore allowing the storage and transmission of accumulated knowledge. The implications of this proposal on the present understanding of the transition from shell‐ to skeleton‐building are explored.     

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.     

Analysis of strain‐induced ground‐water fluctuations at Devils Hole,Nevada

The pool in Devils Hole is a sensitive indicator of crustal strain and fluctuates in response to changes in atmospheric pressure, earth tides, earthquakes, large‐scale tectonic activity and ground‐water development. Short‐term and cyclic water‐level fluctuations caused by atmospheric pressure and earth tides were found to be on the order of millimeters to centimeters. The 1992 Landers/Little Skull Mountain earthquake sequence and the 1999 Hector Mine earthquake induced water‐level offsets of greater than ?12 and ?3.6 cm, respectively. The results of a dislocation model used to compute volumetric strain for each earthquake indicates that the coseismic water‐level offsets are consistent in magnitude and sense with poroelastic responses to earthquake‐induced strain. Theoretical postseismic fluid‐flow modeling indicates that the diffusivity of the system is on the order of 0.03 m² sec^?1, and identified areas of anomalous water‐level fluctuations. Interpretation of model results suggests that while the persistent post‐Landers rise in water‐level can be attributed to deformation‐induced channeling of fluid to the Devils Hole fault zone, the cause of the pre‐Hector Mine water‐level rise may be related to postseismic excess fluid pressures or preseismic strain accumulation.     

flowAMOEBA: Identifying Regions of Anomalous Spatial Interactions

This study aims at developing a data‐driven and bottom‐up spatial statistic method for identifying regions of anomalous spatial interactions (clusters of extremely high‐ or low‐value spatial flows), based on which it creates a spatial flow weights matrix. The method, dubbed flowAMOEBA, upgrades a multidirectional optimum ecotope‐based algorithm (AMOEBA) from areal data to spatial flow data through a proper spatial flow neighborhood definition. The method has the potential to dramatically change the way we study spatial interactions. First, it breaks the convention that spatial interaction data are always collected and modeled between spatial entities of the same granularity, as it delineates the OD region of anomalous spatial interactions, regardless of the size, shape, scale, or administrative level. Second, the method creates an empirical spatial flow weights matrix that can handle network autocorrelation embedded in spatial interaction modeling, thus improving related policy‐making or problem‐solving strategies. flowAMOEBA is tested and demonstrated on a synthetic data set as well as a county‐to‐county migration data set.     

Computing geochemical mass transfer and water/rock ratios in submarine hydrothermal systems: implications for estimating the vigour of convection

A method of calculating chemical water/rock ratios is presented that enables the estimation of fluid velocities in open, flow‐through hydrologic systems. The approach is based on relating the gain/loss of a chemical species per kilogram of solid phase to the loss/gain of that species in the fluid phase, integrated across a specified length of the flowpath. After examining the underlying approximations of the approach using a one‐dimensional model of seawater moving through a basalt under nonisothermal conditions, the method is applied to representative zones within a two‐dimensional hydrothermal convective system. The method requires that regions within the flow system can be identified in which the direction of flow is steady for an extended period of time. Estimates of fluid velocity are spatial and temporal averages for the length of the flowpath used in the calculation. The location within the flow system and the nature of the alteration reactions determine which species can provide reliable values of the chemical water/rock ratio and useful estimates of fluid velocities. Over the length of the flowpath considered, the calculation of water/rock ratios works best when a species is controlled by a single reaction. Accurate estimates are favoured if the concentration profile of a species along the flowpath increases or decreases monotonically. If the length of the flowpath extends over more than one reaction zone, then erroneous estimates of the water/rock ratio and fluid velocity are more likely. Model calculations suggest that the quartz/silica system should provide reliable estimates for fluid velocity under a wide range of temperature and flow conditions, in particular in those regions of a system at or near quartz equilibrium, so that the aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.     

Inferring intermediate‐scale fluid flow in a heterogeneous metasedimentary multilayer sequence during progressive deformation: evidence from the Monts d’Arrée slate belt (Brittany,France)

Quartz veins in the early Variscan Monts d’Arrée slate belt (Central Armorican Terrane, Western France), have been used to determine fluid‐flow characteristics. A combination of a detailed structural analysis, fluid inclusion microthermometry and stable isotope analyses provides insights in the scale of fluid flow and the water–rock interactions. This research suggests that fluids were expelled during progressive deformation and underwent an evolution in fluid chemistry because of changing redox conditions. Seven quartz‐vein generations were identified in the metasedimentary multilayer sequence of the Upper Silurian to Lower Devonian Plougastel Formation, and placed within the time frame of the deformation history. Fluid inclusion data of primary inclusions in syn‐ to post‐tectonic vein generations indicate a gradual increase in methane content of the aqueous–gaseous H₂O–CO₂–NaCl–CH₄–N₂ fluid during similar P–T conditions (350–400°C and 2–3.5 kbar). The heterogeneous centimetre‐ to metre‐scale multilayer sequence of quartzites and phyllites has a range of oxygen‐isotope values (8.0–14.1‰ Vienna Standard Mean Ocean Water), which is comparable with the range in the crosscutting quartz veins (10.5–14.7‰ V‐SMOW). Significant differences between oxygen‐isotope values of veins and adjacent host rock (Δ = ?2.8‰ to +4.9‰ V‐SMOW) suggest an absence of host‐rock buffering on a centimetre scale, but based on the similar range of isotope values in the Plougastel Formation, an intraformational buffering and an intermediate‐scale fluid‐flow system could be inferred. The abundance of veins, their well‐distributed and isolated occurrence, and their direct relationship with the progressive deformation suggests that the intermediate‐scale fluid‐flow system primarily occurred in a dynamically generated network of temporarily open fractures.     

Engaged Universals and Community Economies: The (Human) Right to Water in Colombia

While the United Nations' sanctioning of the human right to water was widely celebrated, many debate the adequacy and political potency of the rights discourse to frame water justice. Drawing on multi‐sited, ethnographic‐based fieldwork in Colombia in 2010 and 2011, and prioritizing activists' reflexivity, the paper explores how water activists in the 2007–2011 referendum campaign engaged the universal human right while making user‐run community aqueducts more visible as place‐based, not‐for‐profit, culturally attuned, and valid alternatives to the corporate model of water supply. This case study suggests that the human right to water cannot be separated from water commons, and that communal users and activists engage the universal under their own terms. It also suggests we think of these water models as “economic communities” in Gibson‐Graham's sense: ethical spaces to make explicit our social relations with water, and to cultivate selves and practices that enact alternative socio‐natural relations through water's circulations.     

Salinity structure of the central North Slope foreland basin, Alaska, USA: implications for pathways of past and present topographically driven regional fluid flow

Previous studies of the areal variation in heat flow in the National Petroleum Reserve Alaska (NPRA) support the existence of an active topographically driven regional fluid flow regime in this central part of the North Slope foreland basin. Drilling records and wireline logs for over 30 wells drilled in the NPRA provide additional field information, which can be used to further constrain interpretation of the pattern of regional flow of basinal waters within the NPRA. Hydraulic heads estimated from drilling mud weights show that ground water flow occurs generally from south to north, but with divergence to the north‐east and north‐west away from the central part of the NPRA towards coastal areas of elevated shallow heat flow. Salinities calculated from SP logs range from less than 1 g L^?1, to marine values of 35 g L^?1, to hypersaline values of over 150 g L^?1. The entire upper sedimentary section to a depth of 2 km or more in the eastern part of the NPRA has been preferentially flushed with meteoric water through an area corresponding to the sandiest portion of the Nanushuk group. Deeper areas of low salinity occur within the Sadlerochit and Lisburne sections. The pattern of regional flow in the east is complicated, however, by the presence of a large mass of hypersaline water at depth. It is not known whether these brines are being displaced laterally and upward towards the discharge end of the basin or whether fresher waters are simply riding up over the top. Deep, hypersaline waters also occur in fault slices in the Brooks Range and have survived meteoric flushing. The brines were probably formed at the time of deposition of the Lisburne carbonates. The fluid flow regime to the west is different. Low‐salinity waters may be flowing northward underneath this section through the Ellesmerian section and discharging upward nearer the coast. However, sparse well log control severely limits what can be deduced about the details of flow paths in the central and western parts of the NPRA.     

Evaluation of soluble benzene migration in the Uinta Basin

Field sampling and mathematical modeling are used to study the long‐distance transport and attenuation of petroleum‐derived benzene in the Uinta Basin, Utah. Benzene concentration was measured from oil and oil field formation waters of the Altamont‐Bluebell and Pariette Bench oil fields in the basin. It was also measured from springs located in the regional groundwater discharge areas, hydraulically down‐gradient from the oil fields sampled. The average benzene concentration in oils and co‐produced waters is 1946 and 4.9 ppm at the Altamont‐Bluebell field and 1533 and 0.6 ppm at the Pariette Bench field, respectively. Benzene concentration is below the detection limit in all springs sampled. Mathematical models are constructed along a north–south trending transect across the basin through both fields. The models represent groundwater flow, heat transfer and advective/dispersive benzene transport in the basin, as well as benzene diffusion within the oil reservoirs. The coupled groundwater flow and heat transfer model is calibrated using available thermal and hydrologic data. We were able to reproduce the observed excess fluid pressure within the lower Green River Formation and the observed convective temperature anomalies across the northern basin. Using the computed best‐fit flow and temperature, the coupled transport model simulates water washing of benzene from the oil reservoirs. Without the effect of benzene attenuation, dissolved benzene reaches the regional groundwater discharge areas in measurable concentration (>0.01 ppm); with attenuation, benzene concentration diminishes to below the detection limit within 1–4 km from the reservoirs. Attenuation also controls the amount of water washing over time. In general, models that represent benzene attenuation in the basin produce results more consistent with field observations.     

Factors controlling free thermal convection in faults in sedimentary basins: implications for the formation of zinc–lead mineral deposits

Stratiform sediment‐hosted Zn–Pb–Ag mineral deposits constitute about 40% of the Earth's zinc resources ( Allen 2001 ), and in most cases their genesis involves the discharge of basinal brines near or on the seafloor through syndepositional faults ( Sangster 2002 ). From the point of view of base metal exploration, it is therefore essential to identify all possible faults that formerly carried the upwelling ore‐forming solutions during mineralising events. This paper presents a numerical investigation of the relative importance of various physical parameters in controlling fluid discharge, recharge and heat transport in faults. A two‐dimensional, free convection of pure water, hydrogeological model is developed for the McArthur basin in northern Australia based on the surface geology, known stratigraphic and structural relationships and regional geophysical interpretations. Numerical experiments and sensitivity analyses reveal that faults with strong initial heat input, due to depth of penetration or magmatic activity, are the most likely candidates to carry discharge fluids to the sites of metal precipitation. Deeper, wider and more permeable faults are more likely to behave as the fluid discharge pathways, whereas shallow, narrow or less permeable faults act as marine water recharge pathways. Compared with these fault‐related factors, aquifer physical properties are less important in determining fluid flow patterns and the geothermal regime. These results are an important step in understanding hydrothermal fluid flow in sedimentary basins in order to develop effective exploration criteria for the location of stratiform Zn–Pb–Ag deposits.