Paleomagnetic constraints on eruption patterns at the Pacaya composite volcano,Guatemala

Pacaya volcano is an active composite volcano located in the volcanic highlands of Guatemala about 40 km south of Guatemala City. Volcanism at Pacaya alternates between Strombolian and Vulcanian, and during the past five years there has been a marked increase in the violence of eruptions. The volcano is composed principally of basalt flows interbedded with thin scoria fall units, several pyroclastic surge beds, and at least one welded tuff. Between 400 and 2000 years BP the W-SW sector of the volcano collapsed producing a horseshoeshaped amphitheater (0.65 km³) and providing a window into the cone's infrastructure. Lava flows and tephra exposed in the amphitheater are more then 200 m thick and when combined with flows erupted recently represent between 30 and 40% of the cone's history. Pacaya is ideally suited for a paleomagnetic study into the timing and duration of eruption episodes at a large, composite volcano. We drilled 27 paleomagnetic sites (25 aa flows, 1 dike, and 1 welded tuff) from four lava-flow sequences with between 4 and 14 sites per sequence. The four sequences represent initial through historic activity at Pacaya. We resolved, what appear to be, 22 time-independent paleomagnetic sites by averaging together directions from successive sites where the sitemean directions were indistinguishable at the 95% level of confidence. However, mean-sequence directions of individual lava-flow sequences yielded unusually high Fisher precision parameters (k=44–224) and small circles of 63% confidence (a₆₃=1.6–6.1°) suggesting as few as three or four time-independent sites were collected. This indicates that activity as Pacaya is strongly episodic and that episodes are characterized by voluminous outpouring of lavas. Modelling the data using Holocene PSV rates confirms this and shows that differences in within-sequence directions (6–11.5°) are consistent with emplacement of lava-flow sequences in less than 100 years to as many as 300 years. Relatively larger differences in directions (18–23°) between subjacent lava-flow sequences indicates that repose is at least 300–500 years and could be even longer.     

Optimal weight absorber designs for vibrating structures exposed to random excitations

Optimal mass ratios that minimize the response of a laminated beam with an attached absorber are tabulated for various values of beam damping. The beam is treated as an equivalent one degree of freedom (1DOF) main system vibrating in the fundamental mode. The beam is subjected to Gaussian white noise force and Gaussian white noise base frame acceleration. Optimal absorber frequency ratios and absorber damping ratios have been tabulated by others; the results for the classical 1DOF main system with attached absorber suggest that the optimized non-dimensional response decreases monotonically as the mass ratio increases. However, to generalize this monotonic relation may lead to inappropriate conclusions. If we define a constraint such that an increase in absorber mass leads to a proportional decrease in available beam construction material, i.e. effectively the combined mass of the beam and absorber is minimized, then variations in the mass ratio will affect the beam's parameters such as mass, stiffness and damping. Since some of these parameters are used for non-dimensionalising the response, inspection of non-dimensional responses may in some cases lead to inappropriate conclusions. This paper shows the optimal mass ratios for minimizing the response of a structure exposed to earthquake or fluid flow type random excitations.     

Analytical solution of boundary integral equations for 2-D steady linear wave problems

Based on the Fourier transform, the analytical solution of boundary integral equations formulated for the complex velocity of a 2-D steady linear surface flow is derived. It has been found that before the radiation condition is imposed, free waves appear both far upstream and downstream. In order to cancel the free waves in far upstream regions, the eigenso-lution of a specific eigenvalue, which satisfies the homogeneous boundary integral equation, is found and superposed to the analytical solution. An example, a submerged vortex, is used to demonstrate the derived analytical solution. Furthermore, an analytical approach to imposing the radiation condition in the numerical solution of boundary integral equations for 2-D steady linear wave problems is proposed.     

A direct inversion scheme for deep resistivity sounding data using artificial neural networks

Initialization of model parameters is crucial in the conventional 1D inversion of DC electrical data, since a poor guess may result in undesired parameter estimations. In the present work, we investigate the performance of neural networks in the direct inversion of DC sounding data, without the need ofa priori information. We introduce a two-step network approach where the first network identifies the curve type, followed by the model parameter estimation using the second network. This approach provides the flexibility to accommodate all the characteristic sounding curve types with a wide range of resistivity and thickness. Here we realize a three layer feed-forward neural network with fast back propagation learning algorithms performing well. The basic data sets for training and testing were simulated on the basis of available deep resistivity sounding (DRS) data from the crystalline terrains of south India. The optimum network parameters and performance were decided as a function of the testing error convergence with respect to the network training error. On adequate training, the final weights simulate faithfully to recover resistivity and thickness on new data. The small discrepancies noticed, however, are well within the resolvability of resistivity sounding curve interpretations.     

Modeling freshening time and hydrochemical evolution of groundwater in coastal aquifers of Shenzhen,China

This work investigated the freshening time and hydrochemical evolution of coastal groundwater in two brackish aquifers in Shenzhen, China. One was the brackish aquifer that resulted from heavy pumping, and the other was the aquifer reclaimed from the coastal sea. Freshening time and hydrochemical evolution of brackish aquifers were quantitatively evaluated using PHREEQC 2.0, a one-dimensional reactive-transport model. Freshening time was shown to mainly depend on pore water velocity, while the chemical composition of groundwater was determined by the cation exchange capacity of the aquifer. It was shown that after heavy pumping ceased, the freshening time for the original coastal aquifer ranged from 20 to over 80 years. While for the coastal reclaimed aquifer, the freshening time was from 85 to 140 years, which depended on the hydraulic conductivity of the fill materials in the reclaimed site. During aquifer freshening, groundwater evolved from Na–Cl type to Ca–Mg–HCO₃ or Na–HCO₃ type. A sensitivity analysis showed that the freshening time was most sensitive to the pore water velocity in the aquifer, while the groundwater chemical composition was most sensitive to the values of cation exchange capacity of the aquifer. As for the dispersivity, it had almost no effect on the freshening time and the chemical composition of groundwater.     

Numerical modeling of a strong dust event over the south China region in March 2010

The mesoscale model WRF-Chem was used to simulate a severe dust storm event that occurred in March 2010. The storm affected a vast area of East Asia, including the south China region and Hong Kong. This southern region is rarely affected by dust weather. The performance of the WRF-Chem was evaluated by observational data such as the National Center for Atmospheric Research reanalysis data for atmospheric circulation, PM10 concentration from various ground stations, and satellite images of Moderate Resolution Imaging Spectroradiometer and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations. The dependence of the model’s performance on certain important parameterizations was examined in this study. For this particular dust storm event, the model results suggest that the simulation is not very sensitive to certain key physical parameterizations such as threshold wind speed of dust emission and the choice of land surface model. In general, the WRF-Chem is capable of capturing the key physical processes for this severe dust event. The analysis of the dust transport fluxes suggests that the dust transport to the south China region is mainly from the north, although there is a mountainous region in the northern part of the south China region.