Interaction between a circular opening and fractures originating from its boundary in a piecewise homogeneous plane

The complex variables boundary element method (CVBEM) is used to study interaction between a circular opening and fractures originating from its boundary in a piecewise homogeneous plane. A new complex hypersingular equation for piecewise homogeneous media with a circular opening is obtained. The equation is solved using the CVBEM technique with circular and straight boundary elements and polynomial approximations (with square root asymptotics for crack tip elements) for the unknown functions. The algorithm is verified through comparison with known semi‐analytical and numerical solutions that involve interaction between a circular opening and specific systems of cracks or other openings. New numerical results concerning the interaction of the circular opening with the cracks and circular inclusions are presented. The method is applied to an important problem in the petroleum industry: modelling propagation of hydraulic fractures in the vicinity of a borehole. Copyright © 2000 John Wiley & Sons, Ltd.     

Some evidence of climate change in twentieth-century India

The study of climate changes in India and search for robust evidences are issues of concern specially when it is known that poor people are very vulnerable to climate changes. Due to the vast size of India and its complex geography, climate in this part of the globe has large spatial and temporal variations. Important weather events affecting India are floods and droughts, monsoon depressions and cyclones, heat waves, cold waves, prolonged fog and snowfall. Results of this comprehensive study based on observed data and model reanalyzed fields indicate that in the last century, the atmospheric surface temperature in India has enhanced by about 1 and 1.1°C during winter and post-monsoon months respectively. Also decrease in the minimum temperature during summer monsoon and its increase during post-monsoon months have created a large difference of about 0.8°C in the seasonal temperature anomalies which may bring about seasonal asymmetry and hence changes in atmospheric circulation. Opposite phases of increase and decrease in the minimum temperatures in the southern and northern regions of India respectively have been noticed in the interannual variability. In north India, the minimum temperature shows sharp decrease of its magnitude between 1955 and 1972 and then sharp increase till date. But in south India, the minimum temperature has a steady increase. The sea surface temperatures (SST) of Arabian Sea and Bay of Bengal also show increasing trend. Observations indicate occurrence of more extreme temperature events in the east coast of India in the recent past. During summer monsoon months, there is a decreasing (increasing) trend in the frequency of depressions (low pressure areas). In the last century the frequency of occurrence of cyclonic storms shows increasing trend in the month of November. In addition there is increase in the number of severe cyclonic storms crossing Indian Coast. Analysis of rainfall amount during different seasons indicate decreasing tendency in the summer monsoon rainfall over Indian landmass and increasing trend in the rainfall during pre-monsoon and post-monsoon months.     

Assessment of the statistical significance of global changes in the annual river runoff in XXI century due to possible anthropogenic warming of climate

A number of uncertainties of forecasts of changes in the annual runoff depths at global scale, obtained using information on results of integration of 21 IPCC climate models is studied. Following possible errors of these forecasts are calculated: errors of models; differences between main (IPCC) scenarios of emission of greenhouse gases in the atmosphere and resultant changes of global temperatures; mistakes in estimates of average long-term observed values of the runoff depths for the “control” period. Global maps of a “significance index” of forecasted changes in the runoff depths (estimations of changes in the annual runoff depths divided by mean square root values of errors of these estimations) for 2025, 2050, 2075 and 2100 are presented. It is shown that the most significant global changes of the runoff depths (growth in the north of Eastern Siberia, of the Russian Far East, of North America, falling in the “Greater Mediterranean Region”) are predicted for the second quarter of 21st century. Further changes of the runoff amplify only in the Amazon basin (reduction, by 2075). Almost everywhere else (including almost all European territory of Russia, Western Siberia, south of Eastern Siberia and of the Far East) the significance of changes in the runoff depths during 21st century is negligible.     

Paleomagnetism of the Cambro-Ordovician McClure Mountain alkalic complex,Colorado

The Middle to Late Cambrian loop in the North American apparent polar wander path (APWP) has been variously attributed to tectonic rotations, remagnetizations and primary magnetizations. Although no primary thermal remanent magnetizations or primary detrital remanent magnetizations have as yet been demonstrated, the temporally self-consistent nature of the loop has been used as an argument for primary magnetizations. We have studied535 ± 5Ma nepheline syenites and syenites of the McClure Mountain alkalic complex, as well as495 ± 10Ma red trachyte dikes which intruded the complex, in an effort to find a primary TRM. Because Zijderveld analysis yielded consistent results for only one trachyte dike, remagnetization great-circle analysis was employed, giving a pole for the trachyte dikes at the tip of the loop (43°N, 114°E), while the syenites and nepheline synenites gave a pole at the base of the loop (18°N, 142°E). The magnetic carrier in the trachytes is hematite which apparently formed during a pervasive hydrothermal alteration. KAr whole rock dating of the trachytes suggests a Pennsylvanian age for the alteration, and thus a late Paleozoic remagnetization of the trachytes. Thus, the low-latitude Cambrian pole is confirmed, but we find no evidence in this study to support the primary nature of the Cambrian APWP.     

Relationship between the effective diffusion coefficient of soil radon and groundwater recharge in semi-arid areas

Herein, the relationship between soil radon behavior and groundwater recharge in the town of Hosha, a semi-arid area of Jordan is explored. The annual rainfall average in this area is between 60–100 mm year⁻¹. Soil radon concentration levels, using CR-39 detectors, have been measured at several depths. Nine sampling sites on a hill slope, in the summer and winter 2006 were chosen. Other soil parameters such as moisture, soil type and porosity were also measured. Correlation between the latter and the calculated effective diffusion coefficients, based on a linear approximation, was done in order to explore water movement. A weak correlation between the soil radon concentration and water recharge in this particular semi-arid area is observed. The results show that water recharge takes place vertically, i.e., water moves downward and in one of the nine sampling sites this movement seems to be more rapid compared to the other sites.     

Searches for very high energy gamma rays from blazars with CANGAROO-III telescope in 2005-2009

We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data, even though they are non-simultaneous, are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size.     

Grain-scale variations in trace element composition of fluid-altered zircon,Acasta Gneiss Complex,northwestern Canada

SEM back-scattered electron images of zircon separates from a 2.9 Ga granitic dyke from the Acasta Gneiss Complex, Slave Province, reveal primary igneous oscillatory zoning which is transgressively overprinted by irregular alteration domains. Electron microprobe analyses show that altered zircon domains are depleted in Zr and Si relative to unaltered zircon while retaining a constant Zr/Si ratio. SIMS trace element analyses indicate that LREE are preferentially enriched in the altered domains. Altered zircon contains elevated concentrations of Ba (up to 580 ppm), Ca (up to 2% weight of the element) and common Pb (>50 ppb). The presence of LREE, Ca, Ba and common Pb in the zircons is interpreted to be the result of secondary alteration by an aqueous fluid, and not by primary incorporation during crystallization. We propose that the alteration-related elements are located in amorphous zircon domains.     

On determining the noon polar cap boundary from SuperDARN HF radar backscatter characteristics

Previous work has shown that ionospheric HF radar backscatter in the noon sector can be used to locate the footprint of the magnetospheric cusp particle precipitation. This has enabled the radar data to be used as a proxy for the location of the polar cap boundary, and hence measure the flow of plasma across it to derive the reconnection electric field in the ionosphere. This work used only single radar data sets with a field of view limited to 2 h of local time. In this case study using four of the SuperDARN radars, we examine the boundary determined over 6 h of magnetic local time around the noon sector and its relationship to the convection pattern. The variation with longitude of the latitude of the radar scatter with cusp characteristics shows a bay-like feature. It is shown that this feature is shaped by the variation with longitude of the poleward flow component of the ionospheric plasma and may be understood in terms of cusp ion time-of-flight effects. Using this interpretation, we derive the time-of-flight of the cusp ions and find that it is consistent with approximately 1 keV ions injected from a subsolar reconnection site. A method for deriving a more accurate estimate of the location of the open-closed field line boundary from HF radar data is described.     

A simplified method for simulation of strong ground motion using finite rupture model of the earthquake source

We present a simplified method to simulate strong ground motion for a realistic representation of a finite earthquake source burried in a layered earth. This method is based on the stochastic simulation method of Boore (Boore, D. M., 1983, Bull. Seism. Soc. Am. 73, 1865–1894) and the Empirical Greens Function (EFG) method of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake symposium, pp. 151–156). The rupture responsible for an earthquake is represented by several subfaults. The geometry of subfaults and their number is decided by the similarity relationships. For simulation of ground motion using the stochastic simulation technique we used the shapping window based on the kinetic source model of the rupture plane. The shaping window deepens on the geometry of the earthquake source and the propagation characteristics of the energy released by various subfaults. The division of large fault into small subfaults and the method for accounting their contribution at the surface is identical to the EGF. The shapping window has been modified to take into account the effect of the transmission of energy released form the finite fault at various boundaries of the layered earth model above the source. In the present method we have applied the correction factor to adjust slip time function of small and large earthquakes. The correction factor is used to simulate strong motion records having basic spectral shape of ² source model in broad frequency range. To test this method we have used the strong motion data of the Geiyo earthquake of 24th March 2001, Japan recorded by KiK network. The source of this earthquake is modelled by a simple rectangular rupture of size 24 × 15 km, burried at a depth of 31 km in a multilayered earth model. This rupture plane is divided into 16 rectangular subfaults of size 6.0 × 3.75 km each. Strong motion records at eight selected near-field stations were simulated and compared with the observed records in terms of the acceleration and velocity records and their response spectrum. The comparison confirms the suitability of proposed rupture model responsible for this earthquake and the efficacy of the approach in predicting the strong motion scenario of earthquakes in the subduction zone. Using the same rupture model of the Geiyo earthquake, we compared the simulated records from our and the EGF techniques at one near-field station. The comparison shows that this technique gives records which matches in a wide frequency range and that too from simple and easily accessible parameters of burried rupture.