MAXIPOL: a balloon-borne experiment for measuring the polarization anisotropy of the cosmic microwave background radiation

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E-mode polarization anisotropy of the cosmic microwave background radiation (CMB) on angular scales of 10^′ to 2°. MAXIPOL is the first CMB experiment to collect data with a polarimeter that utilizes a rotating half-wave plate and fixed wire-grid polarizer. We present the instrument design, elaborate on the polarimeter strategy and show the instrument performance during flight with some time domain data. Our primary dataset was collected during a 26 h turnaround flight that was launched from the National Scientific Ballooning Facility in Ft. Sumner, New Mexico in May 2003. During this flight five regions of the sky were mapped. Data analysis is in progress.     

Study on the kinetics of iron oxide leaching by oxalic acid

The presence of iron oxides in clay or silica raw materials is detrimental to the manufacturing of high quality ceramics. Although iron has been traditionally removed by physical mineral processing, acid washing has been tested as it is more effective, especially for extremely low iron (of less than 0.1% w/w). However, inorganic acids such as sulphuric or hydrochloric acids easily contaminate the clay products with SO₄²⁻ and Cl⁻, and therefore should be avoided as much as possible. On the other hand, if oxalic acid is used, any acid left behind will be destroyed during the firing of the ceramic products. The characteristics of dissolution of iron oxides were therefore investigated in this study.The dissolution of iron oxides in oxalic acid was found to be very slow at temperatures within the range 25–60 °C, but its rate increases rapidly above 90 °C. The dissolution rate also increases with increasing oxalate concentration at the constant pH values set within the optimum range of pH2.5–3.0. At this optimum pH, the dissolution of fine pure hematite (Fe₂O₃) (105–140 μm) follows a diffusion-controlled shrinking core model. The rate expression expressed as 1 − (2 / 3)x − (1 − x)^2 / 3 where x is a fraction of iron dissolution was found to be proportional to [oxalate]^1.5.The addition of magnetite to the leach liquor at 10% w/w hematite was found to enhance the dissolution rate dramatically. Such addition of magnetite allows coarser hematite in the range 0.5–1.4 mm to be leached at a reasonable rate.     

Drought in Bulgaria and atmospheric synoptic conditions over Europe

Drought in Bulgaria is analyzed from the multiple viewpoints of statistical occurrence, spatial patterns, and synoptic conditions. A new index of drought, the SD (spatial-dryness) index, characterizes drought by both intensity and spatial extent. The occurrence of the SD index is analyzed using global gridded data sets. Examination of transitional probabilities of multiple months and years with drought occurrence suggests persistence is sufficiently frequent to be important for climate-related environmental planning. Finally, it is shown that specific seasonal synoptic patterns are associated with wet and dry conditions in Bulgaria. This revised version was published online in July 2006 with corrections to the Cover Date.     

Flood-prone area mapping using machine learning techniques: a case study of Quang Binh province,Vietnam

Natural Hazards - Vietnam’s central coastal region is the most vulnerable and always at flood risk, severely affecting people’s livelihoods and socio-economic development. In...     

Predicting rock size distribution in mine blasting using various novel soft computing models based on meta-heuristics and machine learning algorithms

Blasting is well-known as an effective method for fragmenting or moving rock in open-pit mines. To evaluate the quality of blasting, the size of rock distribution is used as a critical criterion in blasting operations. A high percentage of oversized rocks generated by blasting operations can lead to economic and environmental damage. Therefore, this study proposed four novel intelligent models to predict the size of rock distribution in mine blasting in order to optimize blasting parameters, as well as the efficiency of blasting operation in open mines. Accordingly, a nature-inspired algorithm (i.e., firefly algorithm – FFA) and different machine learning algorithms (i.e., gradient boosting machine (GBM), support vector machine (SVM), Gaussian process (GP), and artificial neural network (ANN)) were combined for this aim, abbreviated as FFA-GBM, FFA-SVM, FFA-GP, and FFA-ANN, respectively. Subsequently, predicted results from the abovementioned models were compared with each other using three statistical indicators (e.g., mean absolute error, root-mean-squared error, and correlation coefficient) and color intensity method. For developing and simulating the size of rock in blasting operations, 136 blasting events with their images were collected and analyzed by the Split-Desktop software. In which, 111 events were randomly selected for the development and optimization of the models. Subsequently, the remaining 25 blasting events were applied to confirm the accuracy of the proposed models. Herein, blast design parameters were regarded as input variables to predict the size of rock in blasting operations. Finally, the obtained results revealed that the FFA is a robust optimization algorithm for estimating rock fragmentation in bench blasting. Among the models developed in this study, FFA-GBM provided the highest accuracy in predicting the size of fragmented rocks. The other techniques (i.e., FFA-SVM, FFA-GP, and FFA-ANN) yielded lower computational stability and efficiency. Hence, the FFA-GBM model can be used as a powerful and precise soft computing tool that can be applied to practical engineering cases aiming to improve the quality of blasting and rock fragmentation.     

An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Construction and analysis of higher order finite difference schemes for the 1D wave equation

In this article, we propose to: 1. Establish most of the properties conjectured in [2] about the higher order finite difference approximation of the 1D Laplace operator. 2. Generalize to any order the fourth-order accurate scheme in space and time of Shubin and Bell [20] and Cohen [6]. For this new family of 2m–2m schemes, we establish, via elementary mathematics, various stability and dispersion results that are helpful to compare these schemes to the 2–2m schemes of Anné et al. [2]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Assessment with satellite data of the urban heat island effects in Asian mega cities

This study focuses on using remote sensing for comparative assessment of surface urban heat island (UHI) in 18 mega cities in both temperate and tropical climate regions. Least-clouded day- and night-scenes of TERRA/MODIS acquired between 2001 and 2003 were selected to generate land-surface temperature (LST) maps. Spatial patterns of UHIs for each city were examined over its diurnal cycle and seasonal variations. A Gaussian approximation was applied in order to quantify spatial extents and magnitude of individual UHIs for inter-city comparison. To reveal relationship of UHIs with surface properties, UHI patterns were analyzed in association with urban vegetation covers and surface energy fluxes derived from high-resolution Landsat ETM+ data. This study provides a generalized picture on the UHI phenomena in the Asian region and the findings can be used to guide further study integrating satellite high-resolution thermal data with land-surface modeling and meso-scale climatic modeling in order to understand impacts of urbanization on local climate in Asia.