Rapid computation of magnetic anomalies with demagnetization included, for arbitrarily shaped magnetic bodies

Summary. The potential function ø for a magnetic body of susceptibility μ in a medium of susceptibility μ* satisfies the integral equation
Here Φ* is the potential function for the region without the heterogeneity and R is the distance from the point of observation to the point on the surface, s , of the body. δΦ /δn is the normal derivative, in the direction of the outward normal. The equation allows for the effects of demagnetization. For numerical purposes the surfaces can be divided into N facets over which δΦ/δ n is a constant. The unknown quantities δΦ/δn_j can be found from the system of equations defined by:
The prime on the summation sign denotes that the summation does not include the i th element. The magnetic field in the direction of the unit vector P( P ₁, P ₂, P₃ ) is given by:      

Modeling the effects of riparian planting strategies on stream temperature: Increasing suitable habitat for endangered Formosan Landlocked Salmon in Shei‐Pa National Park,Taiwan

Reducing or stabilizing the stream temperature of ChiChiaWan Creek is a crucial work for Formosan Landlocked Salmon because ChiChiaWan Creek is the only one habitat for this endangered species. Planting trees in the riparian zone would be one of the alternatives. The purpose of this study was to evaluate the effects of several planting strategies on daily maximum stream temperature along the river. The results showed the effective vegetative shading angles should be more than 50° along ChiChiaWan Creek to reduce the direct solar radiation heating effectively. Upstream planting with 70° vegetative shading angle could be the most effective way among all the scenarios. However, this planting strategy could not improve the worst situations in summer because of the large solar elevation angles. The upstream planting in ChiChiaWan Creek was strongly recommended because the canopies could be easier to extend to totally cover the narrow width of river producing the most effective shades. Practicing the upstream planting with 90° vegetative shading angle can increase more than 1 km long suitable habitats for the endangered Salmon in summer. Alternatively, the west‐side planting scenario was the second effective way for temperature reduction. Our result provided a useful suggestion for the authorities in charge of saving the Formosan Landlocked Salmon, particularly under the stress of global warming. Copyright © 2011 John Wiley & Sons, Ltd.     

Particle shape effect on thermal conductivity and shear wave velocity in sands

This study presents the correlations between quantified shape parameters and geotechnical properties for nine sand specimens. Four shape parameters, sphericity, convexity, elongation and slenderness, were quantified with two-dimensional microscopic images with the aid of image processing techniques. An instrumented oedometer cell is used to measure compressibility, thermal conductivity and shear wave velocity during loading, unloading and reloading stages. As the particle shape inherently determines the initial loose packing condition, initial void ratio and shape parameters are well correlated with compressibility. The applied stress in soils increases the interparticle contact area and contact quality; round particles tend to achieve higher thermal conductivity and shear wave velocity during stress-induced volume change. Multiple linear regression is implemented to capture the relative contributions of involved variables, revealing that the thermal evolution is governed by the initial packing density and particle shape. The experimental observations underscore the predominant effect that particle shape has on the geomechanical and physical properties upon stress-induced soil behavior.     

Sequential extraction and leaching characteristics of heavy metals in abandoned tungsten mine tailings sediments

The chemical speciation of potentially toxic elements (As, Cd, Cu, Pb, and Zn) in the contaminated soils and sulfides-rich tailings sediments of an abandoned tungsten mine in Korea was evaluated by conducting modified BCR sequential extraction tests. Kinetic and static batch leaching tests were also conducted to evaluate the potential release of As and other heavy metals by acidic rain water and the leaching behaviors of these heavy metals. The major sources of the elements were As-, Zn- and Pb-bearing sulfides, Pb carbonates (i.e., cerussite), and Pb sulfates (i.e., anglesite). The biggest pollutant fraction in these soil and tailing samples consists of metals bound to the oxidizable host phase, which can be released into the environment if conditions become oxidative, and/or to residual fractions. No significant difference in total element concentrations was observed between the tailings sediments and contaminated soils. For both sample types, almost no changes occurred in the mobility of As and the other heavy metals at 7 days, but the mobility increased afterwards until the end of the tests at 30 days, regardless of the initial pH. However, the mobility was approximately 5–10 times higher at initial pH 1.0 than at initial pHs of 3.0 and 5.0. The leached amounts of all the heavy metal contents were higher from tailings sediments than from contaminated soils at pH > 3.0, but were lower at pH < 3.0 except for As. Results of this study suggest that further dissolution of heavy metals from soil and tailing samples may occur during extended rainfall, resulting in a serious threat to surface and groundwater in the mine area.     

Estimation of soil erosion rates in a subtropical mountain watershed using 137Cs radionuclide

The formulation of watershed management strategies to protect water resources threatened by soil erosion and sedimentation requires a thorough understanding of sediment sources and factors that drive soil movement in the watershed. This paper describes a study of medium-term water-driven soil erosion rates in a mountainous watershed of the Shihmen Reservoir in Taiwan. A total of 60 sampling sites were selected along a hillslope. At each sampling site, the inventory ¹³⁷Cs activity was determined and then calculated with the diffusion and migration model to derive soil erosion rates. The rates are one to two orders of magnitude lower than estimates using the Universal Soil Loss Equation, a soil erosion model often used in Taiwan. Results of multiple regression analysis indicate that the spatial variability of soil erosion rates is associated with the relative position of a sampling site to the nearest ridge and soil bulk densities (r ² = 0.33, p < 0.01). Finally, the patterns of soil redistribution rates on the hillslope follow the ¹³⁷Cs hillslope model as soil erosion increases in the downslope direction. No deposition site is found at footslope because soil deposition is swept away by regular flooding along the stream channel. This study is an important first step in using ¹³⁷Cs as a tracer of soil redistribution in mountainous watersheds of Taiwan.     

The influences of interannual stratification variability and wind stress forcing on ENSO before and after the 1976 climate shift

In order to understand the change in oceanic variability associated with the climate shift of the mid-1970s, we analyze the contribution of momentum forcing to the leading baroclinic modes over the tropical Pacific using Simple Ocean Data Assimilation (SODA, version 2.0.2) for the period of 1958–1997. Specifically, we look at the statistical relationship between the wind projection coefficients and climate indices and attempt to provide a physical explanation for the observed changes. It is found that the wind stress projection coefficients according to the oceanic baroclinic modes are different in terms of their magnitude and phase in the tropical Pacific, reflecting a specific forcing associated with each mode before and after the 1976 climate shift. Compared to that before the 1970s, the first baroclinic mode is had a greater effect on the interannual sea surface temperature due to equatorial wave dynamics, and there was an increased delayed response of the second baroclinic mode variability to the interannual atmospheric forcing after the late 1970s. This reflects changes in ENSO feedback processes associated with the climate shift. Our analysis further indicates that, after the late 1970s, there was a decrease in the wind stress forcing projecting onto the Ekman layer, which is associated with increased mixed-layer depth. This result suggests that the changes in the ENSO properties before and after the late 1970s are largely associated with the changes in the way in which the wind stress forcing is dynamically projected onto the surface layer of the tropical Pacific Ocean over interannual timescales.     

Ekman motion in shallow open sea in the presence of time-harmonic variation of water depth

The time variation of the wind-induced flow in a homogeneous unbounded sea region has been analytically investigated. The time-dependent Ekman solution in a homogeneous, shallow open sea has been further extended by taking into account the time variation in water depth which might be caused by either tidal motion or change in the mixed layer thickness. The solution approach taken in this study is based on the Galerkin-eigenfunction method originally developed by Heaps [1972. On the numerical solution of three-dimensional hydrodynamical equations for tides and storm surges. Memoires de la Societe Royale des Sciences de Liege Serie 6(2), 143–180]. A series of calculations have been made with emphasis on the influence of the time variation in water depth upon the build-up of Ekman spirals in the presence of oscillatory variations in water depth. It has been found that two oscillations contribute to the wind drift current, the inertial oscillation and the depth-variation-induced oscillation; the inertial oscillation decays with time, but the depth-variation-induced oscillation remains undamped despite the presence of bottom friction. The presence of time-harmonic variation produces peculiar forms of hodograph with a curled or circular pattern according to the angular frequency of the water depth variation.