Variation trends in basaltic rocks of the canary islands

More than two hundred new analyses of basaltic rocks from the Canary Islands are presented. The available data show that the earlier successions have similar variation trends throughout the archipelago; these trends represent all the intermediate types between basic differentiates (oceanites, ankaramites) and more salic differentiates (trachy-basalts, hawaiites). In the more recent volcanic series, the study shows that there is a magmatic diversification with different variation trends in the Eastern than in the Central Canary Islands. In the latter the alkaline character becomes stronger in the successive periods of volcanic activity which have been established in each island. In the Eastern islands the basaltic evolution is, on the contrary, towards basaltic rocks with tholeiitic affinities.     

Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

Significance of basic and ultramafic rock inclusions in the basalts of Canary Islands

This paper deals with the mineralogical, textural and chemical features of the dunite, peridotite, pyroxenite and gabbroic inclusions present in the Quaternary basalts of the Canary Islands. The mineralogical composition, structure and texture of the inclusions show that most of them have been formed as crystal cumulates from a nonalkaline basaltic magma in the earlier stages of its fractional crystallization. There are no co-genetic relationships between the inclusions and the host basalts, since the latter have a very strong alkaline-olivine character, although there are also some types with tholciitic aflinities. The study of the data leads to the conclusion that these inclusions can be considered as xenoliths from the basic and ultrabasic complexes that form the substratum and which outcrop in some of the Canary Islands. Attention is called to the fact that in many other volcanic zones of the world there has been a previous emplacement of basic and ultramafic layered complexes and is the question opened whether the association between stratiform-complexes and active basaltic volcanism is more frequent than has been assumed up to now.     

Land–Water Boundary Treatment for a Tsunami Model With Dimensional Splitting

The Method of Splitting Tsunamis (MOST) model adapted by National Oceanic and Atmospheric Administration (NOAA) for tsunami forecasting operations is praised for its computational efficiency, associated with the use of splitting technique. It will be shown, however, that splitting the computations between \(x\) and \(y\) directions results in specific sensitivity to the treatment of land–water boundary. Slight modification to the reflective boundary condition in MOST caused an appreciable difference in the results. This is demonstrated with simulations of the Tohoku-2011 tsunami from the source earthquake to Monterey Bay, California, and in southeast Alaska, followed by comparison with tide gage records. In the first case, the better representation of later waves (reflected from the coasts) by the modified model in a Pacific-wide simulation resulted in twice as long match between simulated and observed tsunami time histories at Monterey gage. In the second case, the modified model was able to propagate the tsunami wave and approach gage records at locations within narrow channels (Juneau, Ketchikan), to where MOST had difficulty propagating the wave. The modification was extended to include inundation computation. The resulting inundation algorithm (Cliffs) has been tested with the complete set of NOAA-recommended benchmark problems focused on inundation. The solutions are compared to the MOST solutions obtained with the version of the MOST model benchmarked for the National Tsunami Hazard Mitigation Program in 2011. In two tests, Cliffs and MOST results are very close, and in another two tests, the results are somewhat different. Very different regimes of generation/disposal of water by Cliffs and MOST inundation algorithms, which supposedly affected the benchmarking results, have been discussed.     

The Beni-Ilmane (North-Central Algeria) Earthquake Sequence of May 2010

Starting on 14 May 2010 and lasting several months, the village of Beni-Ilmane (Msila District, North-Central Algeria) and its surroundings were struck by an important seismic crisis marked by three successive moderate shocks (5.0 ≤ M _d ≤ 5.2). This sequence of events caused severe damage in the Beni-Ilmane village and in the epicentral area. The poor quality of masonry construction and the cumulative effects of the large number of aftershock events played a key role in the destruction. To follow this earthquake sequence, 11 temporary seismic stations, in addition to the permanent stations of the Algerian seismic network, were deployed in the region. A representative set of well located aftershocks in the period of maximum activity (lasting 18 days) were selected. The horizontal distribution of the aftershocks shows two main earthquake clusters located near Beni-Ilmane village, one cluster oriented E–W and the other oriented NNE–SSW, crossing the first cluster at its eastern tip. The aftershocks distribution suggests that the three main shocks ruptured two distinct and adjacent fault segments of about 8 km length. The focal mechanisms of the first and third events, located in the NNE–SSW cluster, show near-vertical left-lateral strike-slip fault planes. In the second cluster, oriented E–W, focal mechanisms show a high-angle reverse fault. A field survey, initiated immediately after the first main shock, identified surface fissures generated by the three largest events in the sequence. The fissures, concentrated in a narrow area at the western termination of the NE–SW Jebel Choukchot anticline (location of Beni-Ilmane village), showed several orientations which were mainly related to gravity instabilities. The 2010 Beni-Ilmane earthquake sequence, located in the Bibans–Hodna Mountains transition zone, demonstrates that the Tellian Atlas–High Plateaus border region is an active seismic zone marked by moderate and possibly strong earthquakes; thus, a reevaluation of the seismic hazard in the region is needed.     

Assessment of genotoxic and humoral immune system alterations in silica exposed workers from pottery industries in South India

Chemicals, especially silica, have been suspected to cause genetic alterations in pottery industry workers. The present study aims to analyze the frequency of chromosomal aberrations (CA), micronucleus (MN) and DNA damage (comet assay) in the peripheral blood lymphocytes and the immunological alterations workers chronically exposed to silica and in control subjects. In our study, 50 silica exposed workers and 35 control subjects were recruited and silica level was measured by respirable dust and respirable quartz concentrations of badge dosimeter. The serum immunoglobulins (IgM, IgG, and IgA) of all the subjects were measured by using ELISA method. The individuals exposed to silica have a significant increase in the frequency of CA, MN and the total DNA damage (p < 0.05). Immunoglobulin elevation on silica exposed workers was statistically significant (p < 0.05) on comparison with their respective controls. Investigation of the smoking and alcohol habitats coupled with silica exposure in exposed and control subjects represents alcohol consumption and smoking as additional risk factors and must be avoided. Multiple linear regression analysis obtained for CA, MN and comet assay confirm these tests as biomarkers for silica exposed pottery workers. Some confounding factors also showed significant influence on exposed subjects. These results indicate the mutagenic risk in the working environment has a high probability of association with the silica dust exposure in pottery industries. Nevertheless, the present study will create awareness and public concern not only among the silica exposed workers but also to the welfare of their progeny.     

Evidence for an Independent 26-s Microseismic Source near the Vanuatu Islands

The 26 s peak in the ambient seismic noise spectrum is persistently excited and observed at stations globally. Using noise cross-correlation functions (NCFs), the location suggests that the source could be situated in the Gulf of Guinea and Fiji Basin. However, the Fiji Basin was proposed to be the mirror site (near antipode) of the Gulf of Guinea source instead of an independent source, assuming that the surface waves more efficiently propagate along the major-arc paths of oceanic movements. To investigate the propagation of the Rayleigh waves along continental and oceanic paths, we analyzed the surface wave data recorded from an earthquake near the Gulf of Guinea and found that Rayleigh waves travel along continental minor-arc paths more efficiently than along oceanic major-arc paths. We then located the source in the western Pacific Ocean from group velocities measured with earthquake data by using the travel time misfit in NCFs after calibration and concluded that the source is in the Vanuatu Islands. Moreover, the temporal variation of the 26 s microseismic peak observed in the western Pacific seismic stations is very different from that in stations near the Gulf of Guinea, which suggests that they are excited by independent sources. Therefore, the Vanuatu source should be an independent microseismic source. As it is close to volcanoes in the Vanuatu islands, the Pacific 26 s microseismic source might be excited by magmatic processes, which are also responsible for very-long-period volcanic tremors.     

Microseismic monitoring of columnar jointed basalt fracture activity: a trial at the Baihetan Hydropower Station,China

Severe stress release has occurred to the surrounding rocks of the typically columnar jointed basalt after excavation at the Baihetan Hydropower Station, Jinsha River, China, where cracking, collapse, and other types of failure may take place occasionally due to relaxation fracture. In order to understand the relaxation fracture characteristics of the columnar jointed basalt in the entire excavation process at the diversion tunnel of the Baihetan Hydropower Station, real-time microseismic monitoring tests were performed. First, the applicability of a geophone and accelerometer was analyzed in the columnar jointed basalt tunnel, and the results show that the accelerometer was more applicable to the cracking monitoring of the columnar jointed basalt. Next, the waveform characteristics of the microseismic signals were analyzed, and the microseismic signals were identified as follows: rock fracture signal, drilling signal, electrical signal, heavy vehicle passing signal, and blast signal. Then, the attenuation characteristics of the microseismic signals in the columnar jointed basalt tunnel were studied, as well as the types and characteristics of the columnar jointed basalt fracture. Finally, location analysis was conducted on the strong rock fracture events, in which four or more sensors were triggered, to obtain the temporal and spatial evolution characteristics and laws of the columnar jointed basalt relaxation fracture after excavation. The test results are not only of important reference value to the excavation and support of diversion tunnel at the Baihetan Hydropower Station, but also of great referential significance and value to the conduction of similar tests.     

The gravity field and crustal thickness of Venus

The gravity and topography of Venus obtained from observations of the Magellan mission, as well as the gravity and topography from our numerical mantle convection model, are discussed in this paper. We used the hypothesis that the geoid of degrees 2–40 is produced by sublithospheric mantle density anomalies that are associated with dynamical process within the mantle. We obtained the model dynamical admittance(the geoid topography ratio based on a convection model) by a numerical simulation of the Venusian mantle convection, and used it to correct the dynamical effect in the calculation of crustal thickness. After deducting the dynamical effect, the thickness of the Venusian crust is presented. The results show that the gravity and topography are strongly correlated with the Venusian mantle convection and the Venusian crust has a significant influence on the topography. The Venusian crustal thickness varies from 28 to 70 km. Ishtar Terra, and Ovda Regio and Thetis Regio in western Aphrodite Terra have the highest crustal thickness(larger than 50 km). The high topography of these areas is thought to be supported by crustal compensation and our results are consistent with the hypothesis that these areas are remnants of ancient continents. The crustal thickness in the Beta, Themis, Dione, Eistla, Bell, and Lada regiones is thin and shows less correlation with the topography, especially in the Atla and Imdr regiones in the eastern part of Aphrodite Terra. This is consistent with the hypothesis that these highlands are mainly supported by mantle plumes. Compared with the crustal thickness calculated with the dynamical effect, our results are more consistent with the crust evolution and internal dynamical process of Venus.     

FEM × DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations

The article presents a multi-scale modelling approach of cohesive granular materials, its numerical implementation and its results. At microscopic level, Discrete Element Method (DEM) is used to model dense grains packing. At the macroscopic level, the numerical solution is obtained by a Finite Element Method (FEM). In order to bridge the micro- and macro-scales, the concept of Representative Elementary Volume (REV) is applied, in which the average REV stress and the consistent tangent operators are obtained in each macroscopic integration point as the results of DEM’s simulation. In this way, the numerical constitutive law is determined through the detailed modelling of the microstructure, taking into account the nature of granular materials. We first elaborate the principle of the computation homogenisation (FEM × DEM), then demonstrate the features of our multiscale computation in terms of a biaxial compression test. Macroscopic strain location is observed and discussed.     

Geology and tectonic setting of the Mule Creek caldera,New Mexico,U.S.A.

The 10-km diameter Mule Creek caldera is the youngest felsic eruptive center in the Mogollon-Datil volcanic field of southwestern New Mexico. The caldera forms a topographic basin surrounded by a raised rim. The caldera wall is well displayed on the south and west sides of the structure where it dips 20–30 degrees toward the center of the basin. Mudflow breccia fills the caldera and is banked up against the caldera wall. Post-caldera porphyritic quartz latite domes and flows crop out along the ring-fracture zone. The caldera is superimposed upon an older volcanic complex of flow-banded rhyolite and porphyritic andesite lava. The Mule Creek caldera probably originated by explosive eruption of about 10 km³ of pumice and ash, in part preserved in the matrix of the mudflow breccia. Periods of explosive volcanism during the deposition of mudflow breccia are documented by tuffaceous beds interbedded with the breccia. A thin rhyolite ash-flow sheet originated in the caldera and overlies the mudflow breccia. The youngest felsic rocks around the caldera are (1) domes and flows of crystal-rich porphyritic quartz latite of variable mineralogy, interpreted as a defluidized magma, and (2) widespread crystal-poor, flow-banded rhyolite, dated at 18.6 m.y., which is not directly related to the caldera sequence. The Mule Creek caldera and other volcanic features farther south represent the only documented overlap of felsic volcanism with early stages of Basin-Range tectonism in the Mogollon-Datil field.     

Petrology of tholeiitic lavas from Tanna Island (New Hebrides): Importance of cumulative processes in Island arc magmatism

Tanna, one of the southernmost islands of the New Hebrides volcanic arc, is made of Late Pliocene to Recent island arc tholeiitic basalts and andesites, with SiO₂ contents ranging from 45 to 57%. These lavas are highly porphyritic (30–50% in volume): phenocrysts of plagioclase are the most abundant, together with olivine and clinopyroxene. The groundmass contain plagioclase, augite, olivine, magnetite and glass; pigeonite, tridymite, sanidine and, rarely, biotite may also occur. The olivines and clinopyroxenes show an iron enrichment from the cores of phenocrysts to their rims and the groundmass crystals, but their compositional variations are not correlated with the Mg/Fe ratio of bulk host rocks, the most Fe-rich compositions being found in Mg-rich lavas. Plagioclase compositions range from An₉₅ to An₆₀ in the basalts and An₆₀ to An₅₀ in the andesites, but, within each group, they are not correlated with SiO₂ or Na₂O contents of host lavas. Consequently, the bulk major element compositions of Tanna volcanic rocks cannot be considered as primarily controlled by crystal separation from successive liquids. The oxyde-SiO₂ variations diagrams, and the modal compositions and mineral chemistry show that crystal accumulation is the predominant mechanism accounting for bulk rock compositions. However, this does not exclude fractional crystallization: the variation of the calculated groundmass mineralogy strongly suggest the occurrence of crystal removal mainly clinopyroxene and magnetite.     

Eruption of Piip Crater (Kamchatka)

In autumn of 1966 on the northern slope of Kliuchevskoy volcano a chain of new adventive craters broke out at the height of about 2200 m. Eighty-four hours before the beginning of the eruption a swarm of preliminary volcanic earthquakes had appeared. The number of preliminary shocks was 457 with total energy of 4 × 10¹⁷ erg. With the beginning of the lava flow the earthquakes stopped and a continuous volcanic tremor appeared. The total energy of volcanic tremor amounts to 10¹⁶ erg. During the eruption numerous explosive earthquakes with the energy of 10¹⁵–10¹⁶ erg were recorded and besides the microbarograph of the Volcanostation recorded 393 explosions with an energy more than 10¹³ erg and their total energy was equal to 10¹⁷ erg. All together it has been formed 8 explosive craters and the lowest 9th crater was effusive. The slag cone was formed round this effusive crater, the lava effusion of basaltic-andesite composition (52,5% SiO₂) tooke place from the lava boccas at the cone base and from the crater. The lava flow covered a distance of 10 km along the valley of the Sopochnoy river and descended to a height of about 800 m. The lava flow velocity at the outflow reached 800 m/hr, the lava temperature was 1050°C. The effused lava volume amounts to 0.1 km³. The eruption stopped on December 25–26, 1966.     

PALEOMAGIA: A PHP/MYSQL database of the Precambrian paleomagnetic data

Most paleomagnetic applications require a precise, rationally organized and up-todate catalogue or database of paleomagnetic results worldwide. These include reconstructions of continents, calculations of the Apparent Polar Wander Paths (APWPs) or paleolatitude drift curves, testing the Geocentric Axial Dipole (GAD) model, studies of geomagnetic paleosecular variation or reversal asymmetries, comparison of coeval results obtained from different types of rocks, estimation of inclination shallowing in sedimentary rocks and understanding the delay in remanence acquisition caused by slow cooling in large intrusions. For this purpose, various databases, such as the Global Paleomagnetic Database (GPMDB), and the Magnetics Information Consortium Database (MagIC) have been generated. This paper presents a new relational database (PALEOMAGIA) where 3278 entries of Precambrian data have been split geographically, sorted according to age and rock types and ranked using a revised version of the Van der Voo grading scheme. The latest geochronologic information is included wherever available. Significant effort has been put to the retrieval and archiving of data published in the last decade, which are virtually nonexistent in GPMDB. Here we present the database and its browser-based user interface from a scientific and a technical point of view.     

Calc-alkaline volcanic rocks from NW Sardinia: Evaluation of a fractional crystallization model

The volcanic centre of Monte Seda Oro, N. W. Sardinia, representative of a Cenozoic calc-alkaline andesitic suite of rocks is composed of a variety of rocks ranging from high alumina basalts to dacites. The minerals of basaltic, andesitic and dacitic rocks show only limited variation in chemical composition. The geochemical data suggest that the various rock-types are related by a crystal-liquid fractionation. Least-squate numerical calculations, using major element data, support the derivation of andesites with SiO₃ content ranging from 53.8 to 59.0% from basalts having about 48.7% of SiO₂ by low pressure crystal fractionation of the phenocryst phases present in these rocks. However, the origin of dacites cannot be readily explained by this mechanism.     

Rapakivi texture from the O’Leary Porphyry,Arizona (U.S.A.)

The rhyodactic O’Leary Porphyry which forms the Pleistocene (0.233±0.37 m.y.) volcanic domes of O’Leary Peak and Darton Dome in the San Francisco Volcanic Field (northern Arizona, U.S.A.) contains sanidine phenocrysts with oligoclase mantles (rapakivi texture). Rapakivi texture occurs worldwide in silicic rocks of many ages and has been attributed to various igneous and metamorphic processes. The O’Leary Porphyry contains both mantled and unmantled sanidine (both are Or_63–69 Ab_30–36An₁), oligoclase and quartz phenocrysts, labradorite (An₅₃Ab₄₅Or₂) and kaersutite xenocrysts and andesite xenoliths. The compositional range of oligoclase is the same (An_11–26Ab_70–80Or_r–10) for the rapakivi mantles, the oligoclase phenocrysts, and the oligoclase crystals poikilitic within sanidines. Most mantles are discontinuous. The sanidine appears to have been resorbed prior to mantling. Experimental melting studies on the O’Leary Prophyry show that, for a 15 wgt.% water system, plagioclase crystallized prior to sanidine and quartz crystallized last. The O’Leary Porphyry, although inhomogeneous, plots on a Q-Or-Ab-An diagram well within the plagioclase stability field. Poikilitic plagioclases within sanidines further support crystallization of plagioclase prior to sanidine in the O’Leary Porphyry. Exsolution of a ternary feldspar to form a plagioclase mantle is the most commonly accepted igneous theory of rapakivi texture formation but has been eliminated as the origin of the O’Leary Porphyry rapakivi. Petrologic models by Tuttle and Bowen and by Stewart are rejected for the O’Leary rapakivi because of inconsistencies with the O’Leary occurrences. Two theories are viable for the O’Leary rapakivi texture. First, is a decrease in water vapor pressure which would enlarge the plagioclase stability field possibility causing mantling of metastable sanidines. The second and preferred theory is that of an addition of sodium and calcium by basification (chemical assimilation without melting) of the xenoliths within the O’Leary Porphyry. This would move the bulk composition of the melt into the plagioclase field possibly resulting in crystallization of plagioclase on sanidine crystals. Diffusion of sodium and calcium from the xenoliths to sanidine would result in mantling only those crystals near to the xenoliths. Later, convection would result in distribution throughout the melt of rapakivi, unmantled sanidines, and xenolithic kaersutite as is seen in the porphyry. Basic xenoliths are extremely common in rapakivi-bearing rocks. Those within the O’Leary Porphyry are andesitic and show resorption, and in some areas of O’Leary Peak itself, have been drawn out into schlieren.     

Origin of some bedded welded tuffs

Well defined, laterally continuous welded tuff beds from <1 cm to 2 m thick are more common than has previously been recognized. Examples ranging in composition from rhyolitic to basaltic are described from Ordovician volcanic areas in Britain and Norway, and from the Miocene of the Canary Islands. Bedded welded tuffs are most common in areas of alkaline and peralkaline acidic pyroclastics. They generally occur within successions of massive, welded ash-flow tuff, or within non-welded air-fall tuff successions. Sequences consisting entirely of bedded welded tuff range from <1 m up to 75 m thick. Bedded welded tuffs are thought to originate in three ways. Poorly sorted, thick-bedded welded tuffs are interpreted as the deposits of pyroclastic flows, in which case the beds represent either individual flows units or the layers within flow units. Better sorted, thin-bedded welded tuffs are thought to be of air-fall origin. Thirdly, welding may be produced by the effects of an external heat source on non-welded bedded tuffs.     

A statistical model for vesuvius and its volcanological implications

The last 300 years of Vesuvius history are reconstructed as a chronological succession of 4 phenomenological states: i) repose, ii) persistent activity, iii) intermediate eruption and iv) final eruption. It turns out that the times of permanence in each state are distributed according to the same exponential law. Vesuvius activity is then described by a Markov chain of these 4 states, with transition probabilities determined from the previous phenomenological analysis. The model reproduces the Vesuvius activity between 1694 and 1872 and possibly also in the 1872–1944 period. It turns out that, at least between 1694 and 1872, the volcano was behaving like a quasistationary system with 4 equilibrium states, perturbed by a stochastic noise responsible for occasional transitions from an equilibrium state to another. Major physical or structural changes of the volcanic system around 1872 and possibly in the whole subsequent period, are clearly shown by the statistical analysis.     

Design and testing of a global climate prediction system based on a coupled climate model

A global climate prediction system(PCCSM4) was developed based on the Community Climate System Model, version 4.0, developed by the National Center for Atmospheric Research(NCAR), and an initialization scheme was designed by our group. Thirty-year(1981–2010) one-month-lead retrospective summer climate ensemble predictions were carried out and analyzed. The results showed that PCCSM4 can efficiently capture the main characteristics of JJA mean sea surface temperature(SST), sea level pressure(SLP), and precipitation. The prediction skill for SST is high, especially over the central and eastern Pacific where the influence of El Ni?o-Southern Oscillation(ENSO) is dominant. Temporal correlation coefficients between the predicted Ni?o3.4 index and observed Ni?o3.4 index over the 30 years reach 0.7, exceeding the 99% statistical significance level. The prediction of 500-hPa geopotential height, 850-hPa zonal wind and SLP shows greater skill than for precipitation. Overall, the predictability in PCCSM4 is much higher in the tropics than in global terms, or over East Asia. Furthermore, PCCSM4 can simulate the summer climate in typical ENSO years and the interannual variability of the Asian summer monsoon well. These preliminary results suggest that PCCSM4 can be applied to real-time prediction after further testing and improvement.     

Impact of Land Surface Processes on the South Asian Monsoon Simulations Using WRF Modeling System

The Weather Research and Forecasting model has been used to examine the role of land surface processes on Indian summer monsoon simulations. Isolated experiments have been carried out with physical parameterization schemes (land surface and planetary boundary layer) and data assimilation to examine their relative roles in the representation of regional hydroclimate in model simulations. The impact of vegetation green fraction on the model simulations has been extensively studied by replacing the default United States Geological Survey (USGS) vegetation cover data with that of Indian Space Research Organisation (ISRO) data. Results indicate that differences in the treatment of surface processes in the model lead to large differences in precipitation simulation over the Indian domain. Several hydroclimate parameters from the simulations using ISRO and USGS vegetation green fractions were examined. It is seen that the role of vegetation green fraction in these experiments has been to increase latent heat flux to the atmosphere. Two sets of data assimilation experiments were also carried out for an entire year using the same set of observed data but with different land surface parameterization schemes. It is found that evenwhen using the same observed data, the differences in land surface schemes reduce the impact and contribution of observed data being assimilated into the model. The hydroclimate over the region becomes a function of the land surface scheme. This study highlights the importance of vegetation green fraction and land surface schemes in the context of the regional hydroclimate over South Asia.