Study of the anomalous magnetic field structure in the Ural region using parallel algorithms

The anomalous structure in the magnetic field of the Ural Region has been studied in the segment bounded by 52°?C64° N and 54°?C66° E. Analytical apparatus for upward continuation of airborne magnetic data to different heights was applied. To recalculate magnetic field, parallel algorithms and software for multiprocessor computers were used. Maps of magnetic anomalies for different ranges of wave lengths showing the distribution of magnetization in the layers of the Earth??s crust were built.     

Internal wave characteristics in the eastern Arabian Sea during summer monsoon

The time series BT profiles and surface winds and atmospheric pressure, collected in the deep waters off Ratnagiri and Karwar during summer monsoon were utilized to document the characteristics of internal waves (IW). Low-frequency (≤2, cycle per day (cpd)) IW off Ratnagiri are found to propagate at 83 cm/s with wavelengths of 45 km and wave heights upto 40 m. These parameters for high-frequency (>2 cpd) IW off Karwar correspond to 99 cm/s, 3 km and 23 m. The IW off Karwar appear to leave the station at 70° (±10°) (measured from the horizontal). The data sets were further analysed to address the harmonic composition of the IW and identify the possible sources for the observed IW fields. Power spectra of the IW indicated energy peaks at inertial (0·6 cpd) and tidal (1 and 2 cpd) frequencies off Ratnagiri and in the high-frequency band of 0·5–2·0 cycles per hour off Karwar. The coherence between the IW and wind/tide is found to be good at several frequencies within the IW spectrum. This feature probably suggests tides as a source for the IW of tidal frequencies and winds and tides as a joint source for the IW at the remaining frequencies.     

Modified Weibull-derived spectrum for deep water significant wave height estimation

The modified Weibull spectrum is utilized to calculate the zeroth spectral moment (m_o) using Monte Carlo integration methods. Then significant wave height (H_s) is calculated using the formula This is validated with observed buoy data and numerical wave model (WAM) predicted significant wave heights. The Weibull parameters have been calculated using energy densities from observed spectra recorded by DS5 buoy (13.80° N, 82.52° E, depth 3355.48 m) by the method of maximum likelihood (MLE).The relative root mean square error (RRMS) and relative bias error criteria show that modified Weibull spectrum estimated significant wave heights are better than those predicted by WAM model. The monthly averaged observed wave power spectra for the year 2005 recorded by deep water buoy DS5 is considered in this work. The spectra exhibit bimodal sea states for several months of the year.     

Application of Weibull model for redefined significant wave height distributions

It is well accepted that the parent distribution for individual ocean wave heights follows the Weibull model. However this model does not simulate significant wave height which is the average of the highest one-third of some ‘n’ (n- varies) wave heights in a wave record. It is now proposed to redefine significant wave height as average of the highest one-third of a constant number (n-constant, say,n = 100) of consecutive individual wave heights. The Weibull model is suggested for simulating redefined significant wave height distribution by the method of characteristic function. An empirical support of 100.00% is established by Χ²-test at 0.05 level of significance for 3 sets of data at 0900, 1200 and 1500 hrs at Valiathura, Kerala coast. Parametric relations have been derived for the redefined significant wave height parameters such as mean, maximum one-third average, extreme wave heights, return periods of an extreme wave height and the probability of realising an extreme wave height in a time less than the designated return period.     

A distributed water-heat coupled model for mountainous watershed of an inland river basin in Northwest China (III) using the outputs from Mesoscale model version 5

The distributed water-heat coupled (DWHC) model is calibrated, with the help of the Mesoscale model version 5 (MM5), by calculating the daily precipitation, the daily average air temperature at the 2.0 m heights and the daily potential evaporation in Heihe mountainous watershed area and its vicinity (96.786°∼102.284°E, 37.328°∼40.601°N, 17 × 10⁴ km²), from February 11 to June 30, 2003. The MM5 model periodically ran every 10 days in 3 km × 3 km grid resolution with an integral time step of 3 s. In the MM5 model, many scheme or options are consulted or adopted, such as the Grell scheme cumulus parameterization method, the Dudhia option, the cloud-radiation scheme, MRF PBL option and the modified Oregon State University Land-surface model (OSULSM). According to the projection transform methods, the MM5 outputs are interpolated to the 1 km × 1 km grid in Alberts projection by using triangle-based cubic interpolation (Cubic) and nearest neighbor interpolation (Nearest) methods, with which the DWHC model shares the same method. The result shows that, when the Nearest method is used, the Nash-Sutcliffe equation value of the daily average runoff is 0.79, the balance error is −0.79% and the goodness of fit R ² value is 0.81. Meanwhile, when the Cubic method is used, the Nash-Sutcliffe equation value, the balance error and the R ² value are 0.79, −0.65% and 0.80, respectively. Though the runoff simulation result is not favorable, it is still better than that using measured data at the meteorological and hydrological stations; the latter has a Nash-Sutcliffe equation value of 0.61. The MM5-DWHC model results also show that runoff mainly occurs on land surfaces and from shallow soil layers. According to model calibration results, certain outputs of MM5 are singular to some extent and the DWHC model is very sensitive to the initial values.     

Numerical experiment with modelled return echo of a satellite altimeter from a rough ocean surface and a simple iterative algorithm for the estimation of significant wave height

We have simulated the return echo of a satellite altimeter from a rough ocean surface using an analytical formula and have studied its sensitivity with respect to various oceanic and altimeter parameters. Our numerical experiment shows that for normally observed significant wave heights (SWH) the effect of off-nadir angle (ONA) up to 0.5° on the leading edge is not severe. Also, small surface roughness skewness seems to have little effect on the overall shape of the echo. Newton’s iterative scheme has been used to retrieve SWH from the mean return waveform without noise and with additive Gaussian noise typical of Seasat and Geosat altimeters. It has been observed that SWH can be retrieved in the presence of noise with an accuracy of ±0.6 m for ONA less than 0.5°. For higher ONA, accurate retrieval requires the use of precomputed look-up table along with our scheme.     

A climatological study of the relations among solar activity,galactic cosmic ray and precipitation on various regions over the globe

We apply Fourier and wavelet analyses to the precipitation and sunspot numbers in the time series (1901–2000) over Australia (27°S, 133°E), Canada (60°N, 95°W), Ethiopia (8°N, 38°E), Greenland (72°N, 40°W), United Kingdom (54°N, 2°W), India (20°N, 77°E), Iceland (65°N, 18°W), Japan (36°N, 138°E), United States (38°N, 97°W), South Africa (29°S, 24°E) and Russia (60°N, 100°E). Correlation analyses were also performed to find any relation among precipitation, sunspot numbers, temperature, and cloud-cover at the same spatial and temporal scale. Further correlations were also performed between precipitation with electron and proton fluence at the time interval, 1987–2006. All these parameters were considered in annual and seasonal scales. Though correlation study between precipitation and other parameters do not hint any linear relation, still the Fourier and wavelet analyses give an idea of common periodicities. The 9–11 year periodicity of sunspot numbers calculated by Fourier transform is also confirmed by wavelet transform in annual scale. Similarly, wavelet analysis for precipitation also supports the short periods at 2–5 years which is verified by Fourier transform in discontinuous time over different geographic regions.     

Fracture pattern and electrical resistivity studies for groundwater exploration

 The occurrence, movement and control of groundwater, particularly in hard-rock areas, are governed by different factors such as topography, lithology, structures like fractures, faults and nature of weathering. An attempt is made in the present study to investigate the extent of the influence of structures such as fractures and thereby delineate the nature of subsurface lithology with the help of an electrical resistivity method. For this study, the Upper Gunjanaeru River basin, Cuddapah district Andhra Pradesh was chosen to determine groundwater potentials. In order to understand the significance of the fracture pattern, geological, hydrogeomorphological and lineament maps were prepared based on the field data and also from the LANDSAT TM imagery. Further, electrical resistivity surveys were conducted to determine the subsurface lithology and also to confirm the studies of LANDSAT imagery. The isoresistivity contour map has been prepared based on the 45 VES conducted to determine the resistivity variations in the study area. The isoresistivity contours obtained were found to conform to the structural trends obtained by geological studies and also confirm the relationship between the structure and secondary porosity present in the rocks. The lineaments in the area have two preferred directions. One set is a NE-SW direction (N 30°–70° E; S 30°–70° W) and another is a NW-SE direction (N 0°–30° W; S 0°–30° E and N 60°–80° W; S 60°–80° E). The water-table contour map shows that the direction of groundwater flow is south to north. Received: 3 March 1997 · Accepted: 17 June 1997     

Intra-seasonal variations of kinetic energy of lower tropospheric zonal waves during northern summer monsoon

Space spectral analysis of zonal (u) and meridional (v) components of wind and time spectral analysis of kinetic energy of zonal waves at 850 hPa during monsoon 1991 (1st June 1991 to 31st August 1991) for the global belt between equator and 40°N are investigated. Space spectral analysis shows that long waves (wavenumbers 1 and 2) dominate the energetics of Region 1 (equator to 20°N) while over Region 2 (20°N to 40°N) the kinetic energy of short waves (wavenumbers 3 to 10) is more than kinetic energy of long waves. It has been found that kinetic energy of long waves is dominated by zonal component while both (zonal and meridional) the components of wind have almost equal contribution in the kinetic energy of short waves. Temporal variations of kinetic energy of wavenumber 2 over Region 1 and Region 2 are almost identical. The correlation matrix of different time series shows that (i) wavenumber 2 over Regions 1 and 2 might have the same energy source and (ii) there is a possibility of an exchange of kinetic energy between wavenumber 1 over Region 1 and short waves over Region 2. Wave to wave interactions indicate that short waves over Region 2 are the common source of kinetic energy to wavenumber 2 over Regions 1 and 2 and wavenumber 1 over Region 1. Time spectral analysis of kinetic energy of zonal waves indicates that wavenumber 1 is dominated by 30–45 day and bi-weekly oscillations while short waves are dominated by weekly and bi-weekly oscillations. The correlation matrix, wave to wave interaction and time spectral analysis together suggest that short period oscillations of kinetic energy of wavenumber 1 might be one of the factors causing dominant weekly (5–9 day) and bi-weekly (10–18 day) oscillations in the kinetic energy of short waves.     

Geochemistry of bauxites of Belgaum (Karnataka) and Yercaud (Tamil Nadu) areas

The major chemical components of bauxite deposits of Belgaum (76° 24′E : 15° 54′N) and Yercaud (78° 14′E : 11° 48′N) areas have been determined. A chemical continuity between parent rocks (zone I) to bauxites (zone IV) via clay (zone II) and laterites (zone III) clearly indicates that bauxites have been derived byin situ weathering of the respective parent rocks.     

Application of wave transformation models for estimation of morphological changes at Vellar estuary,southeast coast of India

Assessment of the wave climate at near coast is vital for estimation of morphological changes, such as growth of sand spit and associated siltation of tidal inlets. Vellar (bar-built) estuary is one of the prominent estuaries along the southeast coast of India, located at 11°30′N and 79°46′E, less studied in terms of its morphological features. The inlet of Vellar is exposed to high energetic waves, inducing large sediment transport rates and shoreline changes. Local wave characteristics are not accurately defined and the available wave information at near coast is limited (point based observations). In the present study, three decoupled numerical models are employed to derive the monthly nearshore wave climate at Vellar by transforming waves from deep water to nearshore. These models are independently validated with buoy observations in deep water and wave gauge data at nearshore. Based on the nearshore wave data, littoral drift along the coast was estimated and compared with the spit growth at Vellar inlet. The estimated average littoral drift along this coast from February to October is 1.93 × 10⁶ m³ toward north and from November to January it is 1.52 × 10⁶ m³ toward south, resulting in a net northerly drift. Results indicated that increase in the wave energy during the period of July to September is responsible for the maximum growth of the sand spit observed in the field.     

CRETACEOUS-TERTIARY DESERTS IN SOUTHWESTERN CHINA AND DESERT SEDIMENTOLOGY

The special issue deals with the Cretaceous-Tertiary deserts in southwestern China and relevant sedimentological problems about deserts. The study area is located between 101°10＇- 107°00＇E and 28°00＇-30°40′ N. The ancient desert study is so little in China that there has not been a systematic report about it up to now. Based on the study in many ways on Cretaceous-Tertiary deserts in southwestern China, plentiful data have been obtained. Though it is still a bit rough, the study is a good beginning of the ancient desert study in China. The main ideas and conclusions are summarized as follows.     

Simulation of Storm Surges Along Myanmar Coast Using a Location Specific Numerical Model

Coastal flooding induced by storm surges associated with tropical cyclones is one of the greatest natural hazards sometimes even surpassing earthquakes. Although the frequency of tropical cyclones in the Indian seas is not high, the coastal region of India, Bangladesh and Myanmar suffer most in terms of life and property caused by the surges. Therefore, a location-specific storm surge prediction model for the coastal regions of Myanmar has been developed to carry out simulations of the 1975 Pathein, 1982 Gwa, 1992 Sandoway and 1994 Sittwe cyclones. The analysis area of the model covers from 8° N to 23° N and 90° E to 100° E. A uniform grid distance of about 9 km is taken along latitudinal and longitudinal directions. The coastal boundaries in the model are represented by orthogonal straight line segments. Using this model, numerical experiments are performed to simulate the storm surge heights associated with past severe cyclonic storms which struck the coastal regions of Myanmar. The model results are in agreement with the limited available surge estimates and observations.     

Link between convection and meridional gradient of sea surface temperature in the Bay of Bengal

We use daily satellite estimates of sea surface temperature (SST) and rainfall during 1998–2005 to show that onset of convection over the central Bay of Bengal (88–92°E, 14–18°N) during the core summer monsoon (mid-May to September) is linked to the meridional gradient of SST in the bay. The SST gradient was computed between two boxes in the northern (88–92°E, 18–22°N) and southern (82–88°E, 4–8°N) bay; the latter is the area of the cold tongue in the bay linked to the Summer Monsoon Current. Convection over central bay followed the SST difference between the northern and southern bay (ΔT) exceeding 0.75°C in 28 cases. There was no instance of ΔT exceeding this threshold without a burst in convection. There were, however, five instances of convection occurring without this SST gradient. Long rainfall events (events lasting more than a week) were associated with an SST event (ΔT ≥ 0.75°C); rainfall events tended to be short when not associated with an SST event. The SST gradient was important for the onset of convection, but not for its persistence: convection often persisted for several days even after the SST gradient weakened. The lag between ΔT exceeding 0.75°C and the onset of convection was 0–18 days, but the lag histogram peaked at one week. In 75% of the 28 cases, convection occurred within a week of ΔT exceeding the threshold of 0.75°C. The northern bay SST, T _N , contributed more to ΔT, but it was a weaker criterion for convection than the SST gradient. A sensitivity analysis showed that the corresponding threshold for T _N was 29°C. We hypothesise that the excess heating (∼1°C above the threshold for deep convection) required in the northern bay to trigger convection is because this excess in SST is what is required to establish the critical SST gradient.     

Shear wave splitting in the Isparta angle,southwestern Turkey: Anisotropic complexity in the mantle

This study presents shear wave splitting analysis results observed at ISP (Isparta) broadband station in the Isparta Angle, southwestern Turkey. We selected 21 good quality seismic events out of nearly 357 earthquakes and calculated splitting parameters (polarization direction of fast wave, ϕ and delay time between fast and slow waves, δt) from mainly SKS and a few SKKS phases of the selected 21 seismic events. Then, we compared calculated splitting parameters at ISP station (56° ≤ ϕ ≤ 205°; 0.37 s ≤ δt ≤ 4 s) with those previously calculated at ANTO (Ankara) and ISK (İstanbul) stations (27° ≤ ϕ ≤ 59°; 0.6 s ≤ δt ≤ 2.4 s and 26° ≤ ϕ ≤ 54°; 0.6 s ≤ δt ≤ 1.5 s) which are located at 230 and 379 km away from ISP station in central and northwestern Turkey, respectively. The backazimuthal variations of the splitting parameters at ISP station indicate a different and complex mantle polarization anisotropy for the Isparta Angle in southwestern Turkey compared to those obtained for Ankara and İstanbul stations.     

西南印度洋洋盆演化和岩浆地球化学印迹

西南印度洋中脊(SWIR)平均扩张速率约为14 mm/yr,是全球洋中脊系统的重要组成端元,因其具有慢速-超慢速扩张特征,引起全球科学家的广泛关注。基于前人对SWIR的综合研究成果,从构造和岩浆作用两个角度出发,系统地回顾了SWIR的形成和演化历史,探讨了岩浆的分布特征和地幔不均一性成因。SWIR的形成始于冈瓦纳大陆的裂解,中段洋脊区域(26～42°E)是印度洋最早开启的窗口,历经多次洋脊跃迁和扩展作用形成了斜向扩张展布,多分段的构造格局。地幔热点在冈瓦纳大陆裂解过程中扮演了关键角色,并对SWIR的洋底地貌和岩浆作用具有显著影响,其中Bouvet和Marion热点在SWIR的西段和中段岩浆均留下了地球化学印迹。SWIR西段岩浆除却Bouvet热点影响之外表现出与大西洋-太平洋型玄武岩相似的同位素地球化学特征。在SWIR中段,39～41°E附近的岩浆具有显著的DUPAL异常特征,与冈瓦纳大陆的初始形成、裂解紧密相关。受俯冲改造的中—新元古代的造山带岩石圈地幔拆沉是造成SWIR中段地幔不均一性的重要根源。在SWIR东段,46～52°E区域内的局部岩浆组成异常,推测具有大陆地壳物质混染的成因。而在Melville转换断层以东,洋脊形成时间最晚,玄武岩的地幔源区受到了富集组分的交代作用,其同位素组成与相邻的中印度洋中脊(CIR)和东南印度洋中脊(SEIR)地幔源区具有亲缘性。     

Integration and magnitude homogenization of the Egyptian earthquake catalogue

The aim of the present work is to compile and update a catalogue of the instrumentally recorded earthquakes in Egypt, with uniform and homogeneous source parameters as required for the analysis of seismicity and seismic hazard assessment. This in turn requires a detailed analysis and comparison of the properties of different available sources, including the distribution of events with time, the magnitude completeness, and the scaling relations between different kinds of magnitude reported by different agencies. The observational data cover the time interval 1900–2004 and an area between 22°–33.5° N and 25°–36° E. The linear regressions between various magnitude types have been evaluated for different magnitude ranges. Using the best linear relationship determined for each available pair of magnitudes, as well as those identified between the magnitudes and the seismic moment, we convert the different magnitude types into moment magnitudes M _W, through a multi-step conversion process. Analysis of the catalogue completeness, based on the M _W thus estimated, allows us to identify two different time intervals with homogeneous properties. The first one (1900–1984) appears to be complete for M _W ≥ 4.5, while the second one (1985–2004) can be considered complete for magnitudes M _W ≥ 3.     

水库滑坡涌浪传播有限元数值模拟

库岸滑坡涌浪危害巨大,正确预测库区可能的滑坡涌浪非常重要,是工程可行性论证的重要内容之一。将浅水控制方程应用于滑坡涌浪数值模拟,控制方程采用简便且具有较高精度的两步Taylor-Galerkin方法进行求解。通过算例对数值模型的应用进行验证,结果表明,涌浪产生后将以入水点为源点迅速向四周推进并不断衰减,且随传播距离的增加浪高降幅逐渐减小。计算的涌浪高度及水位变化规律与实测资料吻合得很好,并将滑坡涌浪的沿程传播过程可视化。研究结果表明,文中方法模拟滑坡涌浪传播是有效可行的,可用于滑坡涌浪灾害的预测和防治。     

Impact of Berm Breaching on Hyperbenthic Macroinvertebrate Communities in Intermittently Closed Estuaries

Berm breaching of intermittently closed estuaries, either naturally or due to management practices, results in a physicochemical disturbance due to the flushing of water, material and biota into the ocean and the subsequent tidal influx. In 2007 and 2008 comparative and controlled studies were employed in three New Zealand estuaries: Sawmill Creek (46°04′ S 170°12′ E), Otokia Creek (45°57′ S 179°20′ E) and Kaikorai Lagoon (45°56′ S 170°23′ E), to investigate the impact of berm breaching on the hyperbenthic macroinvertebrate community in intermittently closed estuaries. Significant changes in community structure occurred in both the weekly comparative and the controlled studies (ANOSIM P < 0.01). Furthermore, the catch per unit effort of both total and key invertebrate taxa significantly increased after breaching (ANOVA P < 0.01). However, substantial numbers of individuals were expelled into the ocean (5,800 individuals, 20-min sample) while new taxa immigrated into the estuaries. Along with migration-related changes to community structure, berm breaches also resulted in the loss of ∼80% of inundated habitat and the concentration of existing populations. Consequently, the management of intermittent estuary connections to the ocean has implications for the ecology of individual, managed estuaries and also for regional coastal populations of epibenthic invertebrates.     

Estimates of source parameters of <Emphasis Type="Italic">M</Emphasis>4.9 Kharsali earthquake using waveform modelling

This paper presents the computation of time series of the 22 July 2007 M 4.9 Kharsali earthquake. It occurred close to the Main Central Thrust (MCT) where seismic gap exists. The main shock and 17 aftershocks were located by closely spaced eleven seismograph stations in a network that involved VSAT based real-time seismic monitoring. The largest aftershock of M 3.5 and other aftershocks occurred within a small volume of 4 × 4 km horizontal extent and between depths of 10 and 14 km. The values of seismic moment (M _∘) determined using P-wave spectra and Brune’s model based on f ² spectral shape ranges from 10¹⁸ to 10²³ dyne-cm. The initial aftershocks occurred at greater depth compared to the later aftershocks. The time series of ground motion have been computed for recording sites using geometric ray theory and Green’s function approach. The method for computing time series consists in integrating the far-field contributions of Green’s function for a number of distributed point source. The generated waveforms have been compared with the observed ones. It has been inferred that the Kharsali earthquake occurred due to a northerly dipping low angle thrust fault at a depth of 14 km taking strike N279°E, dip 14° and rake 117°. There are two regions on the fault surface which have larger slip amplitudes (asperities) and the rupture which has been considered as circular in nature initiated from the asperity at a greater depth shifting gradually upwards. The two asperities cover only 10% of the total area of the causative fault plane. However, detailed seismic imaging of these two asperities can be corroborated with structural heterogeneities associated with causative fault to understand how seismogenesis is influenced by strong or weak structural barriers in the region.