Potential Field,Geoelectrical and Reflection Seismic Investigations for Massive Sulphide Exploration in the Skellefte Mining District,Northern Sweden

Multi-scale geophysical studies were conducted in the central Skellefte district (CSD) in order to delineate the geometry of the upper crust (down to maximum ～ 4.5 km depth) for prospecting volcanic massive sulphide (VMS) mineralization. These geophysical investigations include potential field, resistivity/induced polarization (IP), reflection seismic and magnetotelluric (MT) data which were collected between 2009 and 2010. The interpretations were divided in two scales: (i) shallow (～ 1.5 km) and (ii) deep (～4.5 km). Physical properties of the rocks, including density, magnetic susceptibility, resistivity and chargeability, were also used to improve interpretations. The study result delineates the geometry of the upper crust in the CSD and new models were suggested based on new and joint geophysical interpretation which can benefit VMS prospecting in the area. The result also indicates that a strongly conductive zone detected by resistivity/IP data may have been missed using other geophysical data.     

High resolution sensors in space and time for determination saltation and creep intensity

The measurement of aeolian sand transport rates at high temporal and spatial resolution is crucial for further progress in testing and developing numerical models of sand movement by wind and in the modelling of sand dunes, ripples and so on. This paper reports the development and field testing of two sand transport sensors. The first one, a webcam commonly used with personal computers, is a new device in aeolian research. The webcam frame transfer is triggered by a sonic anemometer every 0·1 second. Consecutive frames are compared and analysed in real‐time by a computer program. Changes in pixel light intensity exceeding a threshold level are recorded and interpreted as grain movements. The second sensor is a small‐sized ‘Saltiphone’‐type device made of simple loudspeakers with a diameter of 15 mm as used in MP3 players. It can be deployed as a 2 × 3 array of six such devices distributed to enable horizontal and vertical spatial sampling of the sand flux. The devices are tested under field conditions. Both signals and the sum of microphone impacts over 15 minutes are compared to data gathered using a Guelph sand trap, and very good agreement is found. Measurements in a wind tunnel using sieved natural sand indicate that the webcam can be used to infer additional information about the grain size. As an application, the fluid and impact thresholds for aeolian sand transport are investigated in field measurements by analysing the onset and breakdown of saltation in gust and lull intervals of rising and falling wind speeds, respectively. In this way, constitutive equations for sand transport in terms of the wind speed can be tested. If viable, they can be employed to infer estimates for the thresholds by minimizing the root‐mean‐square error between measured and calculated transport data. Copyright © 2012 John Wiley & Sons, Ltd.     

ARIMA representation for daily solar irradiance and surface air temperature time series

Autoregressive integrated moving average (ARIMA) models are used to compare long-range temporal variability of the total solar irradiance (TSI) at the top of the atmosphere (TOA) and surface air temperature series. The comparison shows that one and the same type of the model is applicable to represent the TSI and air temperature series. In terms of the model type surface air temperature imitates closely that for the TSI. This may mean that currently no other forcing to the climate system is capable to change the random walk type variability established by the varying activity of the rotating Sun. The result should inspire more detailed examination of the dependence of various climate series on short-range fluctuations of TSI.     

Improved Recharge Estimation from Portable,Low‐Cost Weather Stations

Groundwater recharge estimation is a critical quantity for sustainable groundwater management. The feasibility and robustness of recharge estimation was evaluated using physical‐based modeling procedures, and data from a low‐cost weather station with remote sensor techniques in Southern Abbotsford, British Columbia, Canada. Recharge was determined using the Richards‐based vadose zone hydrological model, HYDRUS‐1D. The required meteorological data were recorded with a HOBO^TM weather station for a short observation period (about 1 year) and an existing weather station (Abbotsford A) for long‐term study purpose (27 years). Undisturbed soil cores were taken at two locations in the vicinity of the HOBO^TM weather station. The derived soil hydraulic parameters were used to characterize the soil in the numerical model. Model performance was evaluated using observed soil moisture and soil temperature data obtained from subsurface remote sensors. A rigorous sensitivity analysis was used to test the robustness of the model. Recharge during the short observation period was estimated at 863 and 816 mm. The mean annual recharge was estimated at 848 and 859 mm/year based on a time series of 27 years. The relative ratio of annual recharge‐precipitation varied from 43% to 69%. From a monthly recharge perspective, the majority (80%) of recharge due to precipitation occurred during the hydrologic winter period. The comparison of the recharge estimates with other studies indicates a good agreement. Furthermore, this method is able to predict transient recharge estimates, and can provide a reasonable tool for estimates on nutrient leaching that is often controlled by strong precipitation events and rapid infiltration of water and nitrate into the soil.     

Atmospheric and soil moisture controls on evapotranspiration from above and within a Western Boreal Plain aspen forest

The Western Boreal Plain of North Central Alberta comprises a mosaic of wetlands and aspen (Populus tremuloides) dominated uplands where precipitation (P) is normally exceeded by evapotranspiration (ET). As such these systems are highly susceptible to the climatic variability that may upset the balance between P and ET. Above canopy evapotranspiration (ET_C) and understory evapotranspiration (ET_B) were examined using the eddy covariance technique situated at 25.5 m (7.5 m above tree crown) and 4.0 m above the ground surface, respectively. During the peak period of the growing seasons (green periods), ET_C averaged 3.08 mm d^?1 and 3.45 mm d^?1 in 2005 and 2006, respectively, while ET_B averaged 1.56 mm d^?1 and 1.95 mm d^?1. Early in the growing season, ET_B was equal to or greater than ET_C once understory development had occurred. However, upon tree crown growth, ET_B was lessened due to a reduction in available energy. ET_B ranged from 42 to 56% of ET_C over the remainder of the snow‐free seasons. Vapour pressure deficit (VPD) and soil moisture (θ) displayed strong controls on both ET_C and ET_B. ET_C responded to precipitation events as the developed tree crown intercepted and held available water which contributed to peak ET_C following precipitation events >10 mm. While both ET_C and ET_B were shown to respond to VPD, soil moisture in the rooting zone is shown to be the strongest control regardless of atmospheric demand. Further, soil moisture and tension data suggest that rooting zone soil moisture is controlled by the redistribution of soil water by the aspen root system. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Complexity vs. Simplicity: Groundwater Model Ranking Using Information Criteria

A groundwater model characterized by a lack of field data about hydraulic model parameters and boundary conditions combined with many observation data sets for calibration purpose was investigated concerning model uncertainty. Seven different conceptual models with a stepwise increase from 0 to 30 adjustable parameters were calibrated using PEST. Residuals, sensitivities, the Akaike information criterion (AIC and AICc), Bayesian information criterion (BIC), and Kashyap's information criterion (KIC) were calculated for a set of seven inverse calibrated models with increasing complexity. Finally, the likelihood of each model was computed. Comparing only residuals of the different conceptual models leads to an overparameterization and certainty loss in the conceptual model approach. The model employing only uncalibrated hydraulic parameters, estimated from sedimentological information, obtained the worst AIC, BIC, and KIC values. Using only sedimentological data to derive hydraulic parameters introduces a systematic error into the simulation results and cannot be recommended for generating a valuable model. For numerical investigations with high numbers of calibration data the BIC and KIC select as optimal a simpler model than the AIC. The model with 15 adjusted parameters was evaluated by AIC as the best option and obtained a likelihood of 98%. The AIC disregards the potential model structure error and the selection of the KIC is, therefore, more appropriate. Sensitivities to piezometric heads were highest for the model with only five adjustable parameters and sensitivity coefficients were directly influenced by the changes in extracted groundwater volumes.     

Comparing quantitative precipitation forecast methods for prediction of sewer flows in a small urban area

Abstract

Due to the relatively small spatial scale, as well as rapid response, of urban drainage systems, the use of quantitative rainfall forecasts for providing quantitative flow and depth predictions is a challenging task. Such predictions are important when consideration is given to urban pluvial flooding and receiving water quality, and it is worthwhile to investigate the potential for improved forecasting. In this study, three quantitative precipitation forecast methods of increasing complexity were compared and used to create quantitative forecasts of sewer flows 0–3 h ahead in the centre of a small town in the north of England. The HyRaTrac radar nowcast model was employed, as well as two different versions of the more complex STEPS model. The STEPS model was used as a deterministic nowcasting system, and was also blended with the Numerical Weather Prediction (NWP) model MM5 to investigate the potential of increasing forecast lead-times (LTs) using high-resolution NWP. Predictive LTs between 15 and 90 min gave acceptable results, but were a function of the event type. It was concluded that higher resolution rainfall estimation as well as nowcasts are needed for prediction of both local pluvial flooding and combined sewer overflow spill events.

Editor D. Koutsoyiannis; Guest editor R.J. Moore     

A Comparative Study and Evaluation of Mixing-Height Estimation Based on Sodar-RASS, Ceilometer Data and Numerical Model Simulations

The unmanned aerial system SUMO (Small Unmanned Meteorological Observer) has been used for the observation of the structure and behaviour of the atmospheric boundary layer above the Advent Valley, Svalbard during a two-week period in early spring 2009. Temperature, humidity and wind profiles measured by the SUMO system have been compared with measurements of a small tethered balloon system that was operated simultaneously. It is shown that both systems complement each other. Above 200?m, the SUMO system outperforms the tethered balloon in terms of flexibility and the ability to penetrate strong inversion layers of the Arctic boundary layer. Below that level, the tethered balloon system provides atmospheric profiles with higher accuracy, mainly due to its ability to operate at very low vertical velocities. For the observational period, a numerical mesoscale model has been run at high resolution and evaluated with SUMO profiles reaching up to a height of 1500?m above the ground. The sensitivity to the choice of atmospheric boundary-layer schemes and horizontal resolution has been investigated. A new scheme especially suited for stable conditions slightly improves the temperature forecast in stable conditions, although all schemes show a warm bias close to the surface and a cold bias above the atmospheric boundary layer. During one cold and cloudless night, the SUMO system could be operated nearly continuously (every 30?C45?minutes). This allowed for a detailed case study addressing the structure and behaviour of the air column within and above Advent Valley and its interaction with the local topography. The SUMO measurements in conjunction with a 10-m meteorological mast enabled the identification of a very stable nocturnal surface layer adjacent to the valley bottom, a stable air column in the valley and a strong inversion layer above the summit height. The results indicate the presence of inertial-gravity waves during the night, a feature not captured by the model.     

Evaluation of the Interpretation of Ceilometer Data with RASS and Radiosonde Data

Since 2006 different remote monitoring methods for determining mixing-layer height have been operated in parallel in Augsburg (Germany). One method is based on the operation of eye-safe commercial mini-lidar systems (ceilometers). The optical backscatter intensities recorded with ceilometers provide information about the range-dependent aerosol concentration; gradient minima within this profile mark the tops of mixed layers. Special software for these ceilometers provides routine retrievals of lower atmospheric layering. A second method, based on sodar observations, detects the height of a turbulent layer characterized by high acoustic backscatter intensities due to thermal fluctuations and a high variance of the vertical velocity component. This information is extended by measurements with a radio-acoustic sounding system (RASS) that directly provides the vertical temperature profile from the detection of acoustic signal propagation and thus temperature inversions that mark atmospheric layers. Ceilometer backscatter information is evaluated by comparison with parallel measurements. Data are presented from 2 years of combined ceilometer and RASS measurements at the same site and from comparison with a nearby (60 km) radiosonde for larger-scale humidity information. This evaluation is designed to ensure mixing-layer height monitoring from ceilometer data more reliable.