Outlier identification and correction for GRACE aggregated data

Gravity measurements within the Gravity Recovery and Climate Experiment (GRACE) provide a direct measure of monthly changes in mass over the Earth’s land masses. As such changes in mass mainly correspond to water storage changes, these measurements allow to close the continental water balance on large spatial scales and on a monthly time scale within the respective error bounds. When quantifying uncertainties, positive and negative peaks are detected in GRACE aggregated monthly time series (from different data providers) that do not correspond to hydrological or hydro-meteorological signals. These peaks must be interpreted as outliers, which carry the danger of signal degradation. In this paper an algorithm is developed to identify outliers and replace them with hydrologically plausible values. The algorithm is based on a statistical approach in which hydrological and hydro-meteorological signals are used to control the algorithm. The procedure of outlier detection is verified by evaluating catchment based aggregated GRACE monthly signals with ground truth from hydrology and hydro-meteorological signals. The results show improvement in the correlation of GRACE versus hydrometeorological and hydrological signals in most catchments. Also, the noise level is significantly reduced over 255 largest catchments.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

A spatio-temporal three-dimensional conceptualization and simulation of Dera Ismail Khan alluvial aquifer in visual MODFLOW: a case study from Pakistan

Dera Ismail Khan (DIK) is situated in the Lower Indus Basin of Pakistan. The land use has been changed in the canal command area due to irrigation activities near the Indus River. To check the current status and predict the groundwater levels in the area, the unconfined aquifer has been simulated in Visual MODFLOW for a period of 35 years, i.e., from 1985 to 2020. The 2900-km² area has been modeled with a grid of 500 by 500 m and the depth set to 100 m. The aquifer in the study area has been divided vertically and laterally into three and ten zones, respectively, for the characterization. Water wells and streams were used as the sinks and hydrologic boundaries, respectively. The model was successfully calibrated in steady and the non-steady state. The simulation revealed that the whole simulation can be divided into two phases, i.e., before and after the construction of the Chashma Right Bank Canal (CRBC), whereas the results were summarized in the form of water table depth maps and groundwater budget calculations. To determine the groundwater sustainability, a conjunctive use scenario has been employed to simulate the aquifer dynamics till 2020. The simulation revealed incremental drawdowns till the end.     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

Ground viscosity and stiffness measurements for near surface seismic velocity

In land surveys, the weathering layer can often distort the seismic signal due to it passing through rapid velocity and density changes, dispersion, scattering and inelastic absorption. In a simple spring‐dashpot model for the earth response, an equivalent medium groups these complex phenomena into two parameters only; these are called ground viscosity and ground stiffness. The most recent controllers for vibrators can estimate both parameters. To validate these measurements, Saudi Aramco conducted an experiment measuring ground viscosity and stiffness from two different vibrator control systems over an area of varying terrain conditions, including unconsolidated sand and limestone outcrop. The two systems measured different values, but detected similar trends that correlated well with weathering conditions and surface geology, e.g. lower viscosity values on the outcrop than on the sand. The ratio of ground viscosity to ground stiffness can approximate the shallow S‐wave velocity, which we converted into P‐wave velocity through calibration with sparse uphole data. Static corrections incorporating this velocity information somewhat improved the focusing of seismic time sections. This new approach does not require additional acquisition efforts, and can model shallow complex formations in arid areas where classical static methods often fail.     

Air pollution forecasting from sky images with shallow and deep classifiers

Air pollution is one of the most important problems in the new era. Detecting the level of air pollution from an image taken by a camera can be informative for the people who are not aware of exact air pollution level be declared daily by some organizations like municipalities. In this paper, we propose a method to predict the level of the air pollution of a location by taking an image by a camera of a smart phone then processing it. We collected an image dataset from city of Tehran. Afterward, we proposed two methods for estimation of level of air pollution. In the first method, the images are preprocessed and then Gabor transform is used to extract features from the images. At the end, two shallow classification methods are employed to model and predict the level of air pollution. In the second proposed method, a Convolutional Neural Network(CNN) is designed to receive a sky image as an input and result a level of air pollution. Some experiments have been done to evaluate the proposed method. The results show that the proposed 9 method has an acceptable accuracy in detection of the air pollution level. Our deep classifier achieved accuracy about 59.38% which is 10 about 6% higher than traditional combination of feature extraction and classification methods.     

On transfer equation in a semi-infinite atmosphere with albedo ω>1

In this note we derive an exact solution of transfer equation in a plane-parallel semiinfinite atmosphere with albedo >1, by the method of Laplace transform and Wiener-Hopf technique. The emergent intensityI(0, ) is obtained in terms of theH ⁰-functionH ⁰() (Das Gupta, 1978) for which some good approximations are given. Intensity at any depth is also obtained.I(0, )/I(0, 0) is plotted in graphs against [0,1], and shows a maximum which drops and shifts towards the origin as increases.     

Role of hydrological regime and floodplain sediments in channel instability of the Bhagirathi River,Ganga-Brahmaputra Delta,India

Abstract

This research deals with the surface dynamics and key factors – hydrological regime, sediment load, and erodibility of floodplain facies – of frequent channel shifting, intensive meandering, and lateral instability of the Bhagirathi River in the western part of the Ganga-Brahmaputra Delta (GBD). At present, the floodplain of the Bhagirathi is categorized as a medium energy (specific stream power of 10–300 W m^?2), non-cohesive floodplain, which exhibits a mixed-load and a meandering channel, an entrenchment ratio >2.2, width–depth ratio >12, sinuosity >1.4, and channel slope <0.02. In the study area, since 1975, four meander cutoffs have been shaped at an average rate of one in every 9–10 years. In the active meander belt and sand-silt dominated floodplains of GBD, frequent shifting of the channel and meander migration escalate severe bank erosion (e.g. 2.5 × 10⁶ m³ of land lost between 1999 and 2004) throughout the year. Remote sensing based spatio-temporal analysis and stratigraphic analysis reveal that the impact of the Farakka barrage, completed in 1975, is not the sole factor of downstream channel oscillation; rather, hydrogeomorphic instability induced by the Ajay–Mayurakshi fluvial system and the erodibility of floodplain sediments control the channel dynamics of the study area.     

Modulating the Impacts of Stochastic Uncertainties Linked to Deposit Locations in Data-Driven Predictive Mapping of Mineral Prospectivity

Natural Resources Research - The operation of large-scale ore-forming processes triggers the development of neighboring mineral deposits of the same or related types in a metallogenic province....