Airborne laser altimetry survey of Glaciar Tyndall, Patagonia

The first airborne laser altimetry measurements of a glacier in South America are presented. Data were collected in November of 2001 over Glaciar Tyndall, Torres del Paine National Park, Chilean Patagonia, onboard a Twin Otter airplane of the Chilean Air Force. A laser scanner with a rotating polygon-mirror system together with an Inertial Navigation System (INS) were fixed to the floor of the aircraft, and used in combination with two dual-frequency GPS receivers. Together, the laser–INS–GPS system had a nominal accuracy of 30 cm after data processing. On November 23rd, a total of 235 km were flown over the ablation area of Glaciar Tyndall, with 5 longitudinal tracks with a mean swath width of 300 m, which results in a point spacing of approximately 2 m both along and across track. A digital elevation model (DEM) generated using the laser altimetry data was compared with a DEM produced from a 1975 map (1:50,000 scale — Instituto Geográfico Militar (IGM), Chile). A mean thinning of − 3.1 ± 1.0 m a^− 1 was calculated for the ablation area of Glaciar Tyndall, with a maximum value of − 7.7 ± 1.0 m a^− 1 at the calving front at 50 m a.s.l. and minimum values of between − 1.0 and − 2.0 ± 1.0 m a^− 1 at altitudes close to the equilibrium line altitude (900 m a.s.l.). The thinning rates derived from the airborne survey were similar to the results obtained by means of ground survey carried out at  600 m of altitude on Glaciar Tyndall between 1975 and 2002, yielding a mean thinning of − 3.2 m a^− 1 [Raymond, C., Neumann, T.A., Rignot, E., Echelmeyer, K.A., Rivera, A., Casassa, G., 2005. Retreat of Tyndall Glacier, Patagonia, over the last half century. Journal of Glaciology 173 (51), 239–247.]. A good agreement was also found between ice elevation changes measured with laser data and previous results obtained with Shuttle Radar Topography Mission (SRTM) data. We conclude that airborne laser altimetry is an effective means for accurately detecting glacier elevation changes in Patagonia, where an ice thinning acceleration trend has been observed during recent years, presumably in response to warming and possibly also drier conditions.     

Modeling Surface Water-Groundwater Interaction with MODFLOW: Some Considerations

The accuracy with which MODFLOW simulates surface water-groundwater interaction is examined for connected and disconnected losing streams. We compare the effect of different vertical and horizontal discretization within MODFLOW and also compare MODFLOW simulations with those produced by HydroGeoSphere. HydroGeoSphere is able to simulate both saturated and unsaturated flow, as well as surface water, groundwater and the full coupling between them in a physical way, and so is used as a reference code to quantify the influence of some of the simplifying assumptions of MODFLOW. In particular, we show that (1) the inability to simulate negative pressures beneath disconnected streams in MODFLOW results in an underestimation of the infiltration flux; (2) a river in MODFLOW is either fully connected or fully disconnected, while in reality transitional stages between the two flow regimes exist; (3) limitations in the horizontal discretization of the river can cause a mismatch between river width and cell width, resulting in an error in the water table position under the river; and (4) because coarse vertical discretization of the aquifer is often used to avoid the drying out of cells, this may result in an error in simulating the height of the groundwater mound. Conditions under which these errors are significant are investigated.     

Mineralogy of Erromango lavas (New Hebrides): Evidence of an early stage of fractionation in island arc basalts

The orogenic volcanic rocks of Erromango island are divided into a Quaternary theoleiitic group and an older Pliocene one showing calc-alkalic affinity. The microprobe mineralogy of these lavas agrees with this geochemical distinction and indicates a marked iron enrichment trend in the recent tholeiitic lavas, whereas titanomagnetite occurs as early formed phenocrysts (Al, Cr, Mg-rich) only in Pliocene basalts. Some Ni and Mg-rich olivine phenocrysts (Fo88) and unusual chemically zoned clinopyroxene megacrysts are described from Pliocene basic lavas and interpreted as deep-fractionated minerals. The large compositional range of these latter (diopside core to augite rim, through salite intermediate zone) leads us to discuss the physical conditions and the corresponding cation substitutions prevailing during their formation. Diopside cores crystallized at depths of about 20 km from a relatively ‘primitive’ melt. The formation of salite is considered to be the result of concomitant fractionation of Mg-rich phases and increase of f_O2. An early stage of fractionation of Cr-bearing diopside + Ni-bearing olivine would account for the typical Cr and Ni depletion of Erromango lavas.     

Variable-density groundwater flow and solute transport in porous media containing nonuniform discrete fractures

Variations in fluid density can greatly affect fluid flow and solute transport in the subsurface. Heterogeneities such as fractures play a major role for the migration of variable-density fluids. Earlier modeling studies of density effects in fractured media were restricted to orthogonal fracture networks, consisting of only vertical and horizontal fractures. The present study addresses the phenomenon of 3D variable-density flow and transport in fractured porous media, where fractures of an arbitrary incline can occur. A general formulation of the body force vector is derived, which accounts for variable-density flow and transport in fractures of any orientation. Simulation results are presented that show the verification of the new model formulation, for the porous matrix and for inclined fractures. Simulations of variable-density flow and solute transport are then conducted for a single fracture, embedded in a porous matrix. The simulations show that density-driven flow in the fracture causes convective flow within the porous matrix and that the high-permeability fracture acts as a barrier for convection. Other simulations were run to investigate the influence of fracture incline on plume migration. Finally, tabular data of the tracer breakthrough curve in the inclined fracture is given to facilitate the verification of other codes.     

Channel Representation in Physically Based Models Coupling Groundwater and Surface Water: Pitfalls and How to Avoid Them

Recent models that couple three‐dimensional subsurface flow with two‐dimensional overland flow are valuable tools for quantifying complex groundwater/stream interactions and for evaluating their influence on watershed processes. For the modeler who is used to defining streams as a boundary condition, the representation of channels in integrated models raises a number of conceptual and technical issues. These models are far more sensitive to channel topography than conventional groundwater models. On all spatial scales, both the topography of a channel and its connection with the floodplain are important. For example, the geometry of river banks influences bank storage and overbank flooding; the slope of the river is a primary control on the behavior of a catchment; and at the finer scale bedform characteristics affect hyporheic exchange. Accurate data on streambed topography, however, are seldom available, and the spatial resolution of digital elevation models is typically too coarse in river environments, resulting in unrealistic or undulating streambeds. Modelers therefore perform some kind of manual yet often cumbersome correction to the available topography. In this context, the paper identifies some common pitfalls, and provides guidance to overcome these. Both aspects of topographic representation and mesh discretization are addressed. Additionally, two tutorials are provided to illustrate: (1) the interpolation of channel cross‐sectional data and (2) the refinement of a mesh along a stream in areas of high topographic variability.     

Temporal Geochemical Evolution of Neogene Magmatism in the Baguio Gold–Copper Mining District (Northern Luzon, Philippines)

Baguio, in the Central Cordillera of Northern Luzon, is a district that displays porphyry copper and epithermal gold mineralization, associated with Early Miocene–Pliocene–Quaternary calc‐alkaline and adakitic intrusions. Systematic sampling, K‐Ar dating, major and trace elements, and Sr, Nd, Pb isotopic analyses of fresh magmatic rocks indicate three magmatic pulses: an Early Miocene phase (21.2–18.7 Ma), a Middle–Late Miocene phase (15.3–8 Ma) and finally a Pliocene–Quaternary event (3–1 Ma). The first phase emplaced evolved calc‐alkaline magmas, essentially within the Agno Batholith complex, and is thought to be related to the westward‐dipping subduction of the West Philippine Basin. After a quiescence period during which the Kennon limestone was deposited, magmatic activity resumed at 15.3 Ma, in connection with the start of the subduction of the South China Sea along the Manila Trench. It emplaced first petrogenetically related and relatively unradiogenic low‐K calc‐alkaline lavas and intermediate adakites. Temporal geochemical patterns observed from 15.3 to 1 Ma include progressive enrichment in K and other large ion lithophile elements, increase in radiogenic Sr and Pb and corresponding decrease in radiogenic Nd. These features are thought to reflect the progressive addition to the Luzon arc mantle wedge of incompatible elements largely inherited from South China Sea sediments. The origin of the long quiescence period, from 8 to 3 Ma, remains problematic. It might represent a local consequence of the docking of the Zambales ophiolitic terrane to Northern Luzon. Then, magmatic activity resumed at 3 Ma, emplacing chemically diversified rocks ranging from low K to high K and including a large proportion of adakites, especially during the Quaternary (dacitic plugs). The authors tentatively relate this diversity to the development of a slab tear linked with the subduction of the fossil South China Sea ridge beneath the Baguio area.     

Pumping dry: an increasing groundwater budget deficit induced by urbanization,industrialization, and climate change in an over-exploited volcanic aquifer

A methodology is proposed to improve the groundwater budget model by determining the past, present, and future recharge and discharge rates. The model is applied to an increasingly urbanized and industrialized region with drying tendencies: the Toluca Valley, Mexico. This study includes spatially variable recharge determined from the historical climate data, the climate change predictions, and the multiple parameters used in the Hydrologic Evaluation of Landfill Performance (HELP3) model. Using HELP3 a spatial discretization for the average recharge is obtained and estimated at 376 million cubic meters per year (Mm³/year). When considering climate change predictions, by 2050 the average scenario projects recharge to decrease by 15 Mm³/year (from 376 to 361 Mm³/year), and in a worst case scenario up to a maximum decrease of 88 Mm³/year (from 376 to 288 Mm³/year). Groundwater pumping has increased steadily since 1970 and is estimated at 495 Mm³/year for 2010. The current average deficit estimated for 2010 is 172 Mm³/year with average projections increasing to over 292 Mm³/year by 2050. This study of two of the most important components of the water cycle (recharge and discharge) clearly shows that the decreasing water availability in the Toluca basin is due mainly to groundwater pumping and that the current pumping rates are not sustainable. The current deficit can be considered problematic and projections based on expected water consumption and climate change reinforce the need for management of the water resources to be addressed.     

Decision Analysis for Leachate Control at a Fractured Rock Landfill

The municipal landfill at the Complexe Environnemental de Saint-Michel (CESM) in Montreal, which is the third largest in North America, is located in a former quarry in fractured limestone. Impressive measures are taken to monitor and control biogas and leachate generated at the site. Leachate containment is presently performed with a pumping well completed within the waste. The efficiency of the well in controlling off-site leachate migration is questioned because field observations strongly suggest that the nearby former Francon quarry is diverting local ground water flow. To address this issue, four additional hydraulic control options are considered: (1) increased pumping at the existing waste well; (2) new pumping wells in the rock on the eastern limit of the site; (3) new injection wells in the rock on the eastern limit; and (4) combination of new injection wells at the same location and new water supply wells upgradient of the landfill. We evaluated the four hydraulic control options at the CESM using two coupled models: (1) a decision model based on an objective function weighting the risk, costs, and benefits of each option translated into dollar units; and (2) a numerical ground water flow model to represent the effect of operational conditions and ascertain success. Decision analysis offers a quantitative unbiased tool to evaluate the potential and relative cost of each option, but qualitative considerations and judgment still must be used for a complete evaluation. Our analysis confirms that scenario 4, which was the intuitively favored option, represents the best containment strategy.     

Airborne Gravimetry Survey for the Marine Area of the United Arab Emirates

The Military Survey Department (MSD) of the United Arab Emirates (UAE) undertook an airborne gravity survey project for the marine area of the country in 2009, especially to strengthen the marine and coastal geoid in the near-shore regions. For the airborne gravity survey, 5 km spacing coast-parallel flight lines were planned and surveyed. These lines were supplemented by cross-lines in order to assess the quality of the airborne gravity surveys. The flight lines were extended 10 km, spacing lines further offshore. A Beech King Air 350 aircraft was used for the surveys, collecting data at a typical flight speed of 170 knots and a typical flight elevation of 900–1500 m, depending on weather conditions and topography. Gravity was measured with a ZLS-modified LaCoste and Romberg gravimeter (S-99), augmented with a Honeywell strap-down inertial navigation system unit. The estimated accuracy for the airborne gravity data is better than 2.0 mGal r.m.s., as judged from the airborne track crossovers. The new airborne gravimetry data changed the UAE coastal geoid by up to 30 cm in some regions, highlighting the importance of airborne gravity coastal surveys.