Subsurface DC resistivity mapping: approximate 1-D interpretation

Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.     

Tar ball frequency data and analytical results from a long-term beach monitoring program

Following the spill of fuel oils from the New Carissa in February 1999, approximately 300 km of beaches on the Pacific coast of North America were surveyed. A long-term observation program focused on the documentation of stranded tar balls in the vicinity of the spill site. Systematic beach surveys which were conducted over the period March 1999 to April 2001 and semi-logarithmic scale, time-series plots proved the most useful format for identifying trends. Beach monitoring continued through to August 2001. by which time 212 tar balls had been analyzed by GC/MS for their chemical characteristics. The samples of tar balls collected between February 1999 and August 2001 were qualitatively compared with New Carissa source oils (NCSO) and 101 (48%) were not consistent with NSCO. The presence of tar balls that are not related to an incident can confound attempts to define cleanup or endpoint criteria and to assess possible injury to natural resources.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Physical partitioning of mercury in till: an example from central British Columbia, Canada

The physical partitioning of Hg into different grain size fractions of till is predominantly controlled by the primary bedrock mineralogy, the distance of glacial transport, and the relative stability of cinnabar (HgS) in the soil weathering environment. At sites located short distances down-ice from bedrock cinnabar mineralization, the highest Hg concentrations in unoxidized till were measured in the sand- and granule-sized fractions reflecting the abundance of cinnabar in those size ranges. Similar partitioning was measured in oxidized till as cinnabar was found to be relatively resistant to postglacial weathering. Discrete clay-sized cinnabar grains obtained from the unoxidized till were viewed under the scanning electron microscope and suggest that the terminal grade of cinnabar is in the clay-sized range. In till collected from areas barren of cinnabar mineralization, the highest Hg levels were found in the clay-sized fraction which is attributed to the high adsorption of Hg by clay minerals.     

Resolving climate change in the period 15–23 ka in Greenland ice cores: a new application of spectral trend analysis

Northern Hemisphere climate history through and following the Last Glacial Maximum is recorded in detail in ice cores from Greenland. However, the period between Greenland Interstadials 1 and 2 (15–23 ka), i.e. the period of deglaciation following the last major glaciation, has been difficult to resolve in great detail. We here offer a new subdivision of this in the NGRIP, GRIP and GISP2 ice cores, by newly introducing spectral trend analysis to the study of climate-related data series from ice cores. This analysis reveals patterns of change and discontinuity in the waveform properties of a data series, relating to the environmental (including climatic) history of accumulation of the rock or ice record. The application allows high-resolution correlation between the ice cores, and a greatly improved subdivision of the study interval. Nine climatic phases are recognized, within which more identifiable events can also be correlated between the three locations.     

Changes in stream water quality due to logging of the boreal forest in the Montmorency Forest,Québec

Summer stream water quality was monitored before and following the logging of 50% of the boreal forest within three small watersheds (<50 ha) nested in the ‘Ruisseau des Eaux‐Volées’ Experimental Watershed, Montmorency Forest (Québec, Canada). Logging was conducted in winter, on snow cover according to recommended best management practices (BMPs) to minimize soil disturbance and protect advance growth. A 20‐m forest buffer was maintained along perennial streams. In watershed 7·2, cut‐blocks were located near the stream network and logging was partially allowed within the riparian buffer zone. In watersheds 7·5 and 7·7, logging occurred farther away from the stream network. Observations were also made for watershed 7·3 that collected the runoff from watersheds 7·2 and 7·5, and watershed 7·6, the uproad portion of watershed 7·7. The control watershed 0·2 was contiguous to the impacted watersheds and remained undisturbed. Following clearcutting, changes in summer daily maximum and minimum stream temperatures remained within ± 1 °C while changes in diurnal variation did not decrease by more than 0·5 °C. Concentrations of NO₃^? greatly increased by up to 6000% and concentrations of K⁺ increased by up to 300% during the second summer after logging. Smaller increases were observed for Fe_total (up to 71%), specific conductance (up to 26%), and Mg²⁺ (up to 19%). Post‐logging pH decreased slightly by no more than 7% while PO₄^3? concentration remained relatively constant. Suspended sediment concentrations appeared to increase during post‐logging, but there was not enough pre‐logging data to statistically confirm this result. Logging of moderate intensity and respecting established BMPs may account for the limited changes of water quality parameters and the low exceedances of the criteria for the protection of aquatic life. The proximity of the cutover to the stream network and logging within the riparian zone did not appear to affect water quality. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrocarbon pollution of the Mediterranean coastline of Morocco

Organic geochemical study of hydrocarbon pollution of the Mediterranean coastline of Morocco (Tangier-Nador) was performed during March–April 2002. It corresponds to a preliminary work representing an initial assessment of oil contamination of Moroccan coast.Three coastal sites corresponding to Tangier, Tetouan and Nador, representing the main built-up area of this marine region, were selected for superficial sediments, water and organism sampling. Results showed high contamination levels in these three selected sites, known as the most exposed areas to urban and industrial discharges. The maximal contamination corresponding to the sum of non-aromatic (NAH) and polycyclic aromatic hydrocarbons (PAH) encountered in sediment samples was about 553 μg/g (dw). Nevertheless, even a small city agglomeration without any industrial activity showed high values reaching 370 μg/g (dw), which confirm the important hydrocarbons input coming from discharges of fishing ships and oil transporting tankers.The obtained results from this survey reflecting in part the important amounts drained via effluents toward the Moroccan Mediterranean coasts, presented to the concerned authorities drove to the realization of the project of treatment of the wastewater.They later expect to lower rates by 50%, which is still insufficient for a fragile ecosystem such as the Mediterranean and sensitive to anthropogenic effects.     

Limit analysis for saturated porous media without fluid flow calculation

Application of yield design to porous media usually requires a preliminary calculation of the fluid flow net. The stability analysis is then carried out with seepage forces associated with the flow. We assume here that the flow is steady and that the yield criterion is defined by a function of the effective stress tensor. The formulation that we propose here allows taking into account seepage force in the expression of the kinematic stability conditions by means of hydraulic boundary conditions without calculation of the fluid flow. One obtains a formulation of the kinematic condition similar to the case of classic, non-porous media. The method is illustrated by two examples: a cylinder subjected to fluid flow and a vertical cut. It can be adapted to various boundary conditions and to the case of a criterion defined by a function of a generalized effective stress tensor. We also give a method to derive rigorous lower bounds using approximate fluid pressure field. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical Simulation of a Solar Active Region. I: Bastille Day Flare

We present three-dimensional unsteady modeling and numerical simulations of a coronal active region, carried out within the compressible single-fluid MHD approximation. We focus on AR 9077 on 14 July 2000, and the triggering of the X5.7 GOES X-ray class “Bastille Day” flare. We simulate only the lower corona, although we include a virtual photosphere and chromosphere below. The boundary conditions at the base of this layer are set using temperature maps from line intensities and line-of-sight magnetograms (SOHO/MDI). From the latter, we generate vector magnetograms using the force-free approximation; these vector magnetograms are then used to produce the boundary condition on the velocity field using a minimum energy principle (Longcope, Astrophys. J. 612, 1181, 2004). The reconnection process is modeled through a dynamical hyper-resistivity which is activated when the current exceeds a critical value (Klimas et al., J. Geophys. Res. 109, 2218, 2004). Comparing the time series of X-ray fluxes recorded by GOES with modeled time series of various mean physical variables such as current density, Poynting energy flux, or radiative loss inside the active region, we can demonstrate that the model properly captures the evolution of an active region over a day and, in particular, is able to explain the initiation of the flare at the observed time.