Visualization-based analysis of structural and dynamical properties of simulated hydrous silicate melt

We have explored first-principles molecular dynamics simulation data for hydrous MgSiO₃ liquid (with 10 wt% water) to gain insight into its structural and dynamical behavior as a function of pressure (0–150 GPa) and temperature (2,000–6,000 K). By visualizing/analyzing a number of parameters associated with short- and mid-range orders, we have shown that the melt structure changes substantially on compression. The speciation of the water component at low pressures is dominated by the isolated structures (with over 90% hydrogen participated) consisting of hydroxyls, water molecules, O–H–O bridging and four-atom (O–H–O–H and H–O–H–O) groups, where every oxygen atom may be a part of polyhedron or free (i.e., bound to only magnesium atom). Hydroxyls favor polyhedral sites over magnesium sites whereas molecular water is almost entirely bound to magnesium sites, and also interpolyhedral bridging (Si–O–H–O–Si) dominates other types of bridging. Water content is shown to enhance and suppress, respectively, the proportions of hydroxyls and molecular water. As compression increases, these isolated structures increasingly combine with each other to form extended structures involving a total of five or more O and H atoms and also containing threefold coordination species, which together consume over 80% hydrogen at the highest compression studied. Our results show that water lowers the mean coordination numbers of different types including all cation–anion environments. The hydrous melt tends to be more tetrahedrally coordinated but with the Si–Si network being more disrupted compared to the anhydrous melt. Protons increase the content of non-bridging oxygen and decrease the contents of bridging oxygen as well as oxygen triclusters (present at pressures above 10 GPa). The calculated self-diffusion coefficients of all atomic species are enhanced in the presence of water compared to those of the anhydrous melt. This is consistent with the prediction that water depolymerizes the melt structure at all pressures. Our analysis also suggests that proton diffusion involves two processes—the transfer of H atoms (requiring the rupture and formation of O–H bonds) and the motion of hydroxyls as hydrogen carriers (requiring the rupture and formation of Si–O and/or Mg–O bonds). Both the processes are operative at low compression whereas only the first process is operative at high compression.     

Biodiversity and invasibility: Distribution patterns of invasive plant species in the Himalayas,Nepal

Invasive plant species are exerting a serious threat to biological diversity in many regions of the world. To understand plant invasions this study aims to test which of the two plant invasiveness hypotheses; ‘low native diversity' vs. ‘high native diversity', is supported by the regional distribution patterns of invasive plant species in the Himalayas,Nepal. This study is based on data retrieved from published literatures and herbarium specimens. The relationship between invasive plant species distribution patterns and that of native plant species is elucidated by scatter plots, as well as by generalized linear models. The native plant species and invasive plant species have similar distribution patterns and the maximum number of invasive plant species is found in the same altitudinal range where the highest richness for native tree species is found. There is a clear trend of higher invasive plant richness in regions where native tree species richness is relatively high.Consequently, the native plant richness is highest in the central phytogeographic region, followed by the eastern and the western regions, respectively. The invasive plant species also follows a similar trend.Additionally, the invasive plant species richness was positively correlated with anthropogenic factors such as human population density and the number of visiting tourists. This study supports the hypothesis that ‘high native diversity' supports or facilitates invasive plant species. Further, it indicates that nativeand invasive plant species may require similar natural conditions, but that the invasive plant species seem more dependent and influenced by anthropogenic disturbance factors.     

Study of annual mass balance(2011–2013) of Rikha Samba Glacier,Hidden Valley,Mustang, Nepal

Although Himalayan glaciers are of particular interest in terms of future water supplies, regional climate changes, and sea-level rises, little is known about them due to lack of reliable and consistent data. There is a need for monitoring these glaciers to bridge this knowledge gap and to provide field measurements necessary to calibrate and validate the results from different remote sensing operations. Therefore, glaciological observations have been carried out by the Cryosphere Monitoring Project(CMP) since September 2011 on Rikha Samba Glacier in Hidden valley, Mustang district in western Nepal in order to study its annual mass balance. This paper presents the first results of that study. There are 10 glaciers in Hidden Valley, named G1, G2, G3, up to G10. Of these, G5 is the Rikha Samba Glacier, which has the largest area(5.37 km2) in this valley and the highest and lowest altitudes(6,476 and 5,392 m a.s.l., respectively). The glacier mass balance discussed in this paper was calculated using the glaciological method and the equilibrium line altitude(ELA). The glacier showed a negative annual point mass balance along the longitudinal profile of its lower part from September 10, 2011 to October 3, 2012. Stake measurements from October 4, 2012 to September 30, 2013 indicated a negative areal average of annual mass balance-0.088±0.019 m w.e. for the whole glacier. Based on these observations, the ELA of the Rikha Samba Glacier is estimated at 5,800 m a.s.l. in 2013. This negative balance may be due to rising air temperatures in the region, which have been incrementally rising since 1980 accompanied by little or no significant increase in precipitation in that period. The negative mass balance confirms the general shrinking trend of the glacier.     

Some geothermal springs of Nepal

Almost all known thermal springs in Nepal are localized close to and south of either the Main central Thrust or the Main Boundary Fault, and therefore they fall broadly into two major groups. The results from four thermal springs — two from each group — are discussed here. The temperature of the spring water is up to 50°C above normal atmospheric temperature. The waters of almost all springs are alkaline. The dissolved silica content and atomic ratio of Na and K do not indicate any high subsurface temperature and therefore the thermal springs are considered to be of tectonic origin.     

Causes of catastrophic failure of Tam Pokhari moraine dam in the Mt. Everest region

The moraine dam of the Tam Pokhari glacial lake breached on 3 September 1998 and caused a catastrophic flood in the downstream areas. To learn from the event, a field survey was conducted. The survey team found that a landslide, which is considered to be responsible for the outburst flood, occurred in the northeast-facing slope of the moraine dam. The dam internal structure played a crucial role in forming a landslide that triggered the excess overflow and finally the breach of the dam. The internal structure of the dam was made of alternating layers of finer and coarser sediments inclining at 30° downstream and layers are truncated in the upslope direction by a huge pile of unconsolidated and structureless moraine materials. Since the upstream slope angle of the dam i.e., 40° is larger than the angle of repose i.e. 35° of sediments, the increased pore water pressure in the dam triggered a landslide. The rainfall and seismological activities of that particular day, which hit the record high, were crucial in triggering the failure. It is estimated that the dam’s north and northeast-facing slopes completely slid involving about 30,000 m³ of sediment mass of unconsolidated moraine materials above the shear plane. A slope stability analysis was also performed. The calculated safety factor was 0.85, and the calculated slip circle agreed with the shear plane marked in the dam. About 18 million cubic metres of water was swiftly released due to the sudden breach of the moraine dam.     

Shifting paradigms for Nepal’s protected areas: history,challenges and relationships

The modern history of protected area (PA) management in Nepal dates back to 1973 when the National Parks and Wildlife Conservation Act (1973) was promulgated and Chitwan National Park was established. In the years immediately following these key events, protected area acts and regulations were strictly applied and the role of local people in managing natural resources was neglected. However with the passage of time, and with changes in the socio-political and economic characteristics of Nepal, management regimes have shifted towards a more liberal model which recognizes more clearly the contributions of people living and working within protected areas. Recently, landscape level conservation models including the designation of multiple use areas have been utilized in the development of management plans for protected areas in Nepal. Conservation agencies have attempted to tackle challenges such as land use conflict, poaching and smuggling of wildlife parts and illegal harvesting of highly valued medicinal herbs through regulation, but these efforts are not always successful. We recommend a holistic conflict resolution approach which recognizes and resolves the different needs of all stakeholders.     

Gully erosion in the Siwalik Hills,Nepal: estimation of sediment production from active ephemeral gullies

A simple field‐based monitoring programme was established in a small catchment (area 4·6 km²) to find the rates of gully erosion in the Siwalik Hills, Nepal. The rates are used to estimate the amount of sediment produced by gully erosion in the catchment. Three large and active gullies were selected with areas ranging from 0·44 to 0·78 ha. Aerial photographs taken in 1964, 1978 and 1992 were ortho‐rectified and used to study the dynamics of gully heads. The same gullies were also monitored manually using an orthogonal reference system fixed by erosion pins around the gully heads. Results from the aerial photos indicated that the gullies expanded remarkably over the period from 1964 to 1992, by 34 to 58 per cent. Head‐retreat rates during that period were 0·48, 0·55 and 0·73 m a^?1 and average annual sediment evacuation was estimated as 2534 ± 171, 959 ± 60 and 2783 ± 118 m³ a^?1 for the three gullies respectively. From the field measurement, estimated volumes were found to vary from 731 ± 57 to 2793 ± 201 m³ a^?1 over the monitoring period of two years. It was also found that the gullies produce sediment which accounts for up to 59 per cent of the sediment produced from surface erosion in the headwater catchment. The findings are useful for planning and executing appropriate control measures and constructing a sediment hazard map at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Diversity,composition and host-species relationships of epiphytic orchids and ferns in two forests in Nepal

Epiphytic plant species are an important part of biological diversity. It is therefore essential to understand the distribution pattern and the factors influencing such patterns. The present study is aimed at observing the patterns of species richness, abundances and species composition of epiphytic orchids and ferns in two subtropical forests in Nepal. We also studied the relationship of host plants (Schima wallichii and Quercus lanata) and epiphyte species. Data were collected in Naudhara community forest (CF) and the national forest (NF) in Shivapuri Nagarjun National Park. The data were analyzed using univariate and multivariate tests. In total, we recorded 41 species of epiphytes (33 orchid and 8 fern species). Orchid species abundance is significantly higher in CF compared to NF. Orchid species richness and abundance increased with increasing southern aspect whereas it decreased with increasing canopy cover, and fern species richness increased with host bark roughness. Orchid abundance was positively correlated with increasing bark pH, stem size, tree age and tree height and negatively correlated with increasing steepness of the area. Likewise, fern abundances were high in places with high canopy cover, trees that were tall and big, but decreased with increasing altitude and southern aspect. The composition of the orchid and fern species was affected by altitude, aspect, canopy cover, DBH, number of forks and forest management types. We showed that the diversity of orchid and fern epiphytes is influenced by host characteristics as well as host types. The most important pre-requisite for a high epiphyte biodiversity is the presence of old respectively tall trees, independent of the recent protection status. This means: (i) for protection, e.g. in the frame of the national park declaration, such areas should be used which host such old tall trees; and (ii) also in managed forests and even in intensively used landscapes epiphytes can be protected by letting a certain number of trees be and by giving them space to grow old and tall.     

Space–time multiresolution atomistic visualization of MgO and MgSiO3 liquid data

First-principles molecular dynamics simulations of complex material systems such as geophysically relevant oxide and silicate liquids produce massive amounts of time-varying three-dimensional data for the atomic configurations. Given the high accuracy of these data, it is desirable to extract as much information hidden in the data as possible. In this paper, we elaborate on our recently proposed scheme to support interactive visualization at space–time multiresolution of the atomistic simulation data. Instead of just focusing on direct rendering of the given data, additional data (containing more quantitative and qualitative information) that usually have to be extracted by some other means are extracted and rendered on the fly. This allows us to gain better insight into the global as well as local spatio-temporal behavior of the data in the context of bonding, radial distribution, atomic coordination, clustering, structural stability and distortion, and diffusion. We illustrate such visualization for the simulation data on the liquid phases of MgO and MgSiO₃—the two most abundant components of Earth’s mantle. Our analysis shows that the structure and dynamics of both liquids change substantially with compression, with no discernible effects of temperature in most cases.