Development of extended WRF variational data assimilation system (WRFDA) for WRF non-hydrostatic mesoscale model

The paper intends to present the development of the extended weather research forecasting data assimilation (WRFDA) system in the framework of the non-hydrostatic mesoscale model core of weather research forecasting system (WRF-NMM), as an imperative aspect of numerical modeling studies. Though originally the WRFDA provides improved initial conditions for advanced research WRF, we have successfully developed a unified WRFDA utility that can be used by the WRF-NMM core, as well. After critical evaluation, it has been strategized to develop a code to merge WRFDA framework and WRF-NMM output. In this paper, we have provided a few selected implementations and initial results through single observation test, and background error statistics like eigenvalues, eigenvector and length scale among others, which showcase the successful development of extended WRFDA code for WRF-NMM model. Furthermore, the extended WRFDA system is applied for the forecast of three severe cyclonic storms: Nargis (27 April–3 May 2008), Aila (23–26 May 2009) and Jal (4–8 November 2010) formed over the Bay of Bengal. Model results are compared and contrasted within the analysis fields and later on with high-resolution model forecasts. The mean initial position error is reduced by 33% with WRFDA as compared to GFS analysis. The vector displacement errors in track forecast are reduced by 33, 31, 30 and 20% to 24, 48, 72 and 96 hr forecasts respectively, in data assimilation experiments as compared to control run. The model diagnostics indicates successful implementation of WRFDA within the WRF-NMM system.     

Migration of the Ganga river and its implication on hydro-geological potential of Varanasi area, U.P., India

Borehole data reveals that during Late Quaternary, the Ganga river was non-existent in its present location near Varanasi. Instead, it was flowing further south towards peripheral craton. Himalayan derived grey micaceous sands were being carried by southward flowing rivers beyond the present day water divide of Ganga and mixed with pink arkosic sand brought by northward flowing peninsular rivers. Subsequently, the Ganga shifted to its present position and got incised. Near Varanasi, the Ganga river is flowing along a NW-SE tectonic lineament. The migration of Ganga river is believed to have been in response to basin expansion caused due to Himalayan tectonics during Middle Pleistocene times. Multi-storied sand bodies generated as a result of channel migration provide excellent aquifers confined by a thick zone of muddy sediments near the surface. Good quality potable water is available at various levels below about 70 m depth in sandy aquifers. Craton derived gravelly coarse-to-medium grained sand forms the main aquifer zones of tens of meter thickness with enormous yield. In contrast, the shallow aquifers made up of recycled interfluve silt and sandy silt occur under unconfined conditions and show water-level fluctuation of a few meters during pre-and post-monsoon periods.     

Feasibility of geoelectrical monitoring and multiphase modeling for process understanding of gaseous CO2 injection into a shallow aquifer

Potential pathways in the subsurface may allow upwardly migrating gaseous CO₂ from deep geological storage formations to be released into near surface aquifers. Consequently, the availability of adequate methods for monitoring potential CO₂ releases in both deep geological formations and the shallow subsurface is a prerequisite for the deployment of Carbon Capture and Storage technology. Geoelectrical surveys are carried out for monitoring a small-scale and temporally limited CO₂ injection experiment in a pristine shallow aquifer system. Additionally, the feasibility of multiphase modeling was tested in order to describe both complex non-linear multiphase flow processes and the electrical behavior of partially saturated heterogeneous porous media. The suitability of geoelectrical methods for monitoring injected CO₂ and geochemically altered groundwater was proven. At the test site, geoelectrical measurements reveal significant variations in electrical conductivity in the order of 15?C30?%. However, site-specific conditions (e.g., geological settings, groundwater composition) significantly influence variations in subsurface electrical conductivity and consequently, the feasibility of geoelectrical monitoring. The monitoring results provided initial information concerning gaseous CO₂ migration and accumulation processes. Geoelectrical monitoring, in combination with multiphase modeling, was identified as a useful tool for understanding gas phase migration and mass transfer processes that occur due to CO₂ intrusions in shallow aquifer systems.     

Fairweather atmospheric electricity at Antarctica during local summer as observed from Indian station,Maitri

Surface measurements of the atmospheric electrical parameters like Maxwell current, electric field and conductivity studied at the Indian station, Maitri (70.75° S, 11.75° E, 117 m above mean sea level), Antarctica, during austral summer have been analyzed for the years 2001 to 2004. A total of 69 days were selected which satisfied the ‘fairweather’ conditions, i.e., days with absence of high winds, drifting or falling snow, clouds, and fog effects. The diurnal variation curve of electric field and vertical current averaged for 69 fairweather days is a single periodic with a minimum at 03:00 UT and a maximum near 19:00 UT, which is very similar to the Carnegie curve. The correlation coefficient between these measured parameters has a high value (more than 0.9) for all the days. During fairweather days the measured current and field variations are similar and hence it is clear that the conductivity is more or less stable. During magnetically disturbed days, the dawn-dusk potential drop has clear influences on the diurnal variation and it modifies the conductivity. Apart from the day-to-day variation in low latitude thunderstorm activity, there are diurnal, seasonal, inter-annual variations in the electric potential and the currents, as well as solar influences on the measured parameters. This study will help us to examine the impact of solar and geophysical phenomena like solar flares, geomagnetic storms and substorms on the global electric circuit.     

Dissolution mechanisms of goethite in the presence of siderophores and organic acids

In dynamic natural systems such as soils and surface waters, transient biogeochemical processes can induce strong chemical non-steady-state conditions. In this paper, we investigate the effects of non-steady-state conditions on ligand-controlled iron oxide dissolution. The rates of goethite dissolution at pH 6 in the presence of low molecular weight organic acids (oxalate, citrate or malonate) were observed. Non-steady-state conditions were induced by rapid additions of fungal, bacterial or plant siderophores. In the presence of the low molecular weight organic acids, dissolved iron concentrations are below detection limit as predicted by equilibrium solubility calculations. The rapid addition of the siderophores triggered reproducible, fast dissolution of kinetically labile iron from the iron oxide surface. The same effect was observed upon rapid additions of high citrate concentrations to goethite-oxalate suspensions. The concentration of the labile iron pool at the mineral surface was a function of the surface concentration of the low molecular weight organic acids and of the reaction time before addition of the siderophores. Isotopic exchange with ⁵⁹Fe independently confirmed the existence of the labile iron pool before addition of the siderophore. A dissolution mechanism was elucidated that is consistent with these observations and with accepted models of ligand-controlled dissolution. We conclude that the fast dissolution reaction observed here is an important process in biological iron acquisition and that it is based on a general geochemical mechanism.     

Numerical simulation of the marine boundary layer characteristics over the Bay of Bengal as revealed by BOBMEX-98 Pilot experiment

An attempt has been made to study the marine boundary layer characteristics over Bay of Bengal using BOBMEX (Bay of Bengal and Monsoon Experiment) pilot experiment data sets, which was conducted between 23rd October and 12th November 1998 on board ORV Sagar Kanya. A one-dimensional multilevel atmospheric boundary layer with TKE-ε closure scheme is employed to study the marine boundary layer characteristics. In this study two synoptic situations are chosen: one represents an active convection case and the other a suppressed convection. In the present article the marine boundary layer characteristics such as temporal evolution of turbulent kinetic energy, height of the boundary layer and the airsea exchange processes such as sensible and latent heat fluxes, drag coefficient for momentum are simulated during both active and suppressed convection. Marine boundary layer height is estimated from the vertical profiles of potential temperature using the stability criterion. The model simulations are compared with the available observations.     

Oligo-Miocene granitic magmatism in central Vietnam and implications for continental deformation in Indochina

Two granitoid intrusions within the Bu Khang extensional complex in central Vietnam have been dated by U–Pb and Rb–Sr geochronology. A monazite U–Pb age of 26.0 ± 0.2 (2σ) Myr was obtained for the Bu Khang pluton and 23.7 + 1.6/–1.7 Myr for monazite, allanite and zircon from the Dai Loc intrusion. These ages date crystallization of magmas previously assigned Precambrian to Devonian. Rb–Sr analyses of K-feldspar and biotite fractions from the samples yield ages of 19.8 ± 0.6 (2σ) Myr and 19.6 ± 0.5 Myr, respectively. The thermal history recorded by the different geochronometers implies an average exhumation rate of ∼2 mm yr⁻¹ corresponding to ∼9 km of unroofing. Magmatism was either (i) induced passively by lithospheric thinning driven by changes in regional tectonic stresses, or (ii) triggered actively by an ascending plume. Tertiary exhumation and magmatism documented elsewhere in Indochina (e.g. Ailao Shan-Red River and Wang Chao shear zones) favours a regional tectonic cause for extension and granitoid magmatism in the Bu Khang complex. On the other hand, the presence of an upwelling thermal anomaly since at least 35 Ma, causing mantle melting below Indochina, is supported by shear-wave velocity variations in the mantle, and source geochemistry of both the Bu Khang plutons and the Red River belt intrusions. In either case, Tertiary exhumation of the Bu Khang complex can account for previously undocumented NE–SW-directed extension, which is required in northern Vietnam to account for structural changes related to the opening of the South China Sea.     

Mass physical properties, sliding and erodibility of experimentally deposited and differently consolidated clayey muds (Approach, equipment, and first results)

On experimentally deposited kaolinite, illite, and Ca-bentonite consolidated under their own load or by additional vertical pressure, the progress of compaction in relation to excess pore water pressure, mass physical properties, gravitational mass movements in a tilted tank, and erodibility under running water in a flume were studied. The very low consolidated sediments near the mud/water interface do not obey the generally used theory in soil mechanics. They show a different, non-linear relationship between void ratio or water content and depth below the sedimentary surface on the one hand, or effective overburden pressure and shear strength on the other. The same is true of other physical properties such as permeability, which changes considerably with depth and time of consolidation. High sedimentation rates on slopes induce shallow sediment flow, whereas at low rates and critical slope angles different types of slope failures including the breaking up of water-rich sediment into sharp boundered blocks are observed. Flume studies on soft clay muds show three different types of erosion: continuously suspending, discontinuous erosion of crumbs or shreds, and wavy deformation of the clay surface with disintegration of particles from the crests. The critical tractive stress depends not only on clay type, void ratio, and shear strength, but significantly also on the ‘geologic history’of the clay (i.e. deposition from thin suspension or dense slurry, fabric, consolidation and swelling generating minute inhomogeneities etc.). The experiments may lead to a better understanding of all mechanical processe's including pore-water flow taking place near or not far below the sediment/ water interface.     

Noble gas investigations of lunar rocks 10017 and 10071

The noble gases He, Ne, Ar, Kr and Xe and also K and Ba were measured in the Apollo 11 igneous rocks 10017 and 10071, and in an ilmenite and two feldspar concentrates separated from rock 10071. Whole rock K/Ar ages of rocks 10017 and 10071 are (2350 ± 60) × 10⁶ yr and (2880 ± 60) × 10⁶ yr, respectively. The two feldspar concentrates of rock 10071 have distinctly higher ages: (3260 ± 60) × 10⁶ yr and (3350 ± 70) × 10⁶ yr. These ages are still 10 per cent lower than the Rb/Sr age obtained by Papanastassiouet al. (1970) and some Ar⁴⁰ diffusion loss must have occurred even in the relatively coarse-grained feldspar.The relative abundance patterns of spallation Ne, Ar, Kr and Xe are in agreement with the ratios predicted from meteoritic production rates. However, diffusion loss of spallation He³ is evident in the whole rock samples, and even more in the feldspar concentrates. The ilmenite shows little or no diffusion loss. The isotopic composition of spallation Kr and Xe is similar to the one observed in meteorites. Small, systematic differences in the spallation Kr spectra of rocks 10017 and 10071 are due to variations in the irradiation hardness (shielding). The Kr spallation spectra in the mineral concentrates are different from the whole rock spectra and also show individual variations, reflecting the differences in target element composition. The relative abundance of cosmic ray produced Xe¹³¹ differs by nearly 50 per cent in the two rocks. The other Xe isotopes show no variations of similar magnitude. The origin of the Xe¹³¹ yield variability is discussed.Kr⁸¹ was measured in all the samples investigated. The Kr⁸¹/Kr exposure ages of rocks 10017 and 10071 are (480 ± 25) × 10⁶ yr and (350 ± 15) × 10⁶ yr, respectively. Exposure ages derived from spallation Ne²¹, Ar³⁸, Kr⁸³ and Xe¹²⁶ are essentially in agreement with the Kr⁸¹/Kr ages. The age of rock 10071 might be somewhat low because of a possible recent exposure of our sample to solar flare particles.     

A study of strontium in core 119K,Discovery deep,Red Sea

The isotopic composition of strontium of pore water and of authigenic minerals leached from the sediment of core 119K with hot aqua regia is similar to that of the brine in the Discovery deep and differs from that of normal seawater. The average ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratio of strontium removed by acid leaching is 0.7077 ± 0.0007 (1σ) compared to a value of 0.70904 for the Red Sea. The detrital silicate fraction exhibits an approximate inverse correlation between ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratios and strontium concentrations which provides tentative support for a model in which the detrital silicate fraction of deep-sea sediment is considered to be a mixture of terrigenous dust of sialic composition enriched in radiogenic ⁸⁷Sr and of volcanogenic material of basaltic composition and low ⁸⁷Sr abundance. The ${}^{87}\text{Sr}{}^{86}\text{Sr}$ ratios of the shells of foraminifers and pteropods, expressed as δ ⁸⁷Sr‰ relative to 0.70904 for seawater, decrease from $\text{?0.23 ± 0.17‰}$ at 90 cm to $\text{?0.82 ± 0.17‰}$ at 273 cm and remain constant at this value to a depth of 450 cm. The lowering of the δ ⁸⁷Sr values is attributed both to the presence of aragonite overgrowths on pteropod shells and to possible isotope exchange with strontium in the connate fluid.