Efficient two-dimensional magnetotellurics modelling using implicitly restarted Lanczos method

This paper presents an efficient algorithm, FDA2DMT (Free Decay Analysis for 2D Magnetotellurics (MT)), based on eigenmode approach to solve the relevant partial differential equation, for forward computation of two-dimensional (2D) responses. The main advantage of this approach lies in the fact that only a small subset of eigenvalues and corresponding eigenvectors are required for satisfactory results. This small subset (pre-specified number) of eigenmodes are obtained using shift and invert implementation of Implicitly Restarted Lanczos Method (IRLM). It has been established by experimentation that only 15–20% smallest eigenvalue and corresponding eigenvectors are sufficient to secure the acceptable accuracy. Once the single frequency response is computed using eigenmode approach, the responses for subsequent frequencies can be obtained in negligible time. Experiment design results for validation of FDA2DMT are presented by considering two synthetic models from COMMEMI report, Brewitt-Taylor and Weaver (1976) model and a field data based model from Garhwal Himalaya.     

Prediction and error growth in the daily forecast of precipitation from the NCEP CFSv2 over the subdivisions of Indian subcontinent

This study investigates the forecast skill and predictability of various indices of south Asian monsoon as well as the subdivisions of the Indian subcontinent during JJAS season for the time domain of 2001–2013 using NCEP CFSv2 output. It has been observed that the daily mean climatology of precipitation over the land points of India is underestimated in the model forecast as compared to observation. The monthly model bias of precipitation shows the dry bias over the land points of India and also over the Bay of Bengal, whereas the Himalayan and Arabian Sea regions show the wet bias. We have divided the Indian landmass into five subdivisions namely central India, southern India, Western Ghat, northeast and southern Bay of Bengal regions based on the spatial variation of observed mean precipitation in JJAS season. The underestimation over the land points of India during mature phase was originated from the central India, southern Bay of Bengal, southern India and Western Ghat regions. The error growth in June forecast is slower as compared to July forecast in all the regions. The predictability error also grows slowly in June forecast as compared to July forecast in most of the regions. The doubling time of predictability error was estimated to be in the range of 3–5 days for all the regions. Southern India and Western Ghats are more predictable in the July forecast as compared to June forecast, whereas IMR, northeast, central India and southern Bay of Bengal regions have the opposite nature.     

AUTOMATIC INTERPRETATION OF GRAVITY ANOMALIES ASSOCIATED WITH TWO-DIMENSIONAL MASS DISTRIBUTIONS*

Several new features are described which facilitate an automatic, more meaningful, and more accurate interpretation of gravity anomalies associated with approximately two-dimensional mass distributions. These include a provision for fixed points on the bounds of a distribution, outward dipping boundaries, calculation of the gravity effect at individual elevations of stations, smoothing of models and end-corrections for distributions of limited strike length. It is possible to obtain distributions about a median plane; these distributions are sensitive to shape and allow estimates of optimum depths and minimum density contrasts for the anomalous masses.     

Simulations of Eurasian Winter Temperature Trends in Coupled and Uncoupled CFSv2

Conflicting results have been presented regarding the link between Arctic sea-ice loss and midlatitude cooling, particularly over Eurasia. This study analyzes uncoupled(atmosphere-only) and coupled(ocean–atmosphere) simulations by the Climate Forecast System, version 2(CFSv2), to examine this linkage during the Northern Hemisphere winter, focusing on the simulation of the observed surface cooling trend over Eurasia during the last three decades. The uncoupled simulations are Atmospheric Model Intercomparison Project(AMIP) runs forced with mean seasonal cycles of sea surface temperature(SST)and sea ice, using combinations of SST and sea ice from different time periods to assess the role that each plays individually,and to assess the role of atmospheric internal variability. Coupled runs are used to further investigate the role of internal variability via the analysis of initialized predictions and the evolution of the forecast with lead time.The AMIP simulations show a mean warming response over Eurasia due to SST changes, but little response to changes in sea ice. Individual runs simulate cooler periods over Eurasia, and this is shown to be concurrent with a stronger Siberian high and warming over Greenland. No substantial differences in the variability of Eurasian surface temperatures are found between the different model configurations. In the coupled runs, the region of significant warming over Eurasia is small at short leads, but increases at longer leads. It is concluded that, although the models have some capability in highlighting the temperature variability over Eurasia, the observed cooling may still be a consequence of internal variability.     

ON THE SIGNIFICANCE OF AMPLITUDE STUDIES IN SHALLOW REFRACTION SEISMICS*

An analysis of amplitudes of refraction records of some shallow refraction profiles shot primarily for detailing the near-surface structure in a granitic terrain has yielded information on refractor properties: reduced amplitudes are plotted on amplitude-distance graphs. The negative power n to which distance should be raised to represent (elastic) amplitude decay with respect to distance due to spreading of the critically refracted wave involved is examined. Computed values of this “spreading index”n are close to n = 2 as predicted by the theory. With this value of n, amplitude data are processed to determine residual attenuation attributable to elastic absorption in the bedrock. A graphical approach for this purpose from comparison of amplitude-distance graphs with the plots of amplitude decay due to spreading which is applicable to flat and horizontal refractor situations is suggested. Assuming residual attenuation to represent absorption in the granite bedrock, the computed coefficients of absorption, which vary from 0.5 to 3.90 km^?1 for a frequency of 50 Hz, are obtained. From amplitude graphs of reversed profiles it is shown that the amplitude differences plot bears a relation to lateral velocity changes in the refractor. From comparison of practical amplitude decay graphs with those computed for different subsurface models, it appears possible to detect fractured rock occurrences in the refractor.     

TRANSFORMATION OF SCHLUMBERGER APPARENT RESISTIVITY TO DIPOLE APPARENT RESISTIVITY OVER LAYERED EARTH BY THE APPLICATION OF DIGITAL LINEAR FILTERS *

A simple extension of our previous work in which digital filters were developed to transform dipole resistivity measurements over layered earth to Schlumberger ones leads us to the development of filters for transforming the latter to the former. As in the previous work we use a sampling interval of 1/6 In 10 in designing the filters that are both accurate and fast in operation.     

RESISTIVITY TYPE CURVES OVER OUTCROPPING VERTICAL DYKE—I *

The theoretical horizontal resistivity profiles over an outcropping vertical dyke for various parameters-electrode spacing, vein-width and resistivity contrast—with inline alpha, beta and gamma-Wenner electrode systems are described. The resistivity profiles present a most bewildering variety of shapes as compared to those in resistivity soundings. The analysis of type curves suggests suitable electrode configuration for detection of wide, moderately wide, and thin veins. The negative apparent resistivity values on the gamma anomalies for resistive vein over certain vein-widths of higher positive values of resistivity reflection factor contradict the concept of apparent resistivity hitherto held.