共查询到20条相似文献,搜索用时 62 毫秒
1.
A mathematical model of the middle and high latitude ionosphere 总被引:5,自引:0,他引:5
R. W. Schunk 《Pure and Applied Geophysics》1988,127(2-3):255-303
2.
A. S. Kirillov 《Geomagnetism and Aeronomy》2014,54(4):523-529
The relative rate constants of O2(b1Σ g + , v= 1–4) production at an inelastic interaction between electronically excited N(2D) atoms and O2(X3Σ g ? , v = 0) oxygen molecules have been calculated. It was shown that an increase in equilibrium distances between oxygen atoms in NO2 quasi-molecule, produced during the interaction, substantially increases the calculated relative production rates of O2(b1Σ g + , v > 1). The obtained coefficients are used to calculate the O2(b1Σ g + , v= 1–4) relative populations at 110 km (T = 250 K) and 150 km (T = 500 K) altitudes of the polar ionosphere. The calculated populations have been compared with the results of the published measurements of the Atmospheric system band luminosity intensities, and satisfactory agreement has been obtained for low altitudes. 相似文献
3.
Changes in the chemical composition of the hot springs of Mendeleev Volcano (Kunashir Island) as for Cl?, SO 4 2? , CO2, NH 4 + and Cl?/SO 4 2? are given in function of the 1965–66 and 1973 (in part) seismic activity in the South Kurile islands. 相似文献
4.
Anita Lakhani Ravindra Singh Parmar Satya Prakash 《Pure and Applied Geophysics》2012,169(5-6):859-871
Dew samples were collected between October 2007 and February 2008 from a suburban site in Agra. pH, conductivity, major inorganic ions (F?, Cl?, NO 3 ? , SO 4 2? , Na+, K+, Ca2+, Mg2+, and NH 4 + ), and some trace metals (Cr, Sn, Zn, Pb, Cd, Ni, Mn, Fe, Si, Al, V, and Cu) were determined to study the chemistry of dew water. The mean pH was 7.3, and the samples exhibited high ionic concentrations. Dew chemistry suggested both natural and anthropogenic influences, with acidity being neutralized by atmospheric ammonia and soil constituents. Ion deposition flux varied from 0.25 to 3.0?neq?m?2?s?1, with maximum values for Ca2+ followed by NH 4 + , Mg2+, SO 4 2? , Cl?, NO 3 ? , Na+, K+, and F?. Concentrations of trace metals varied from 0.13 to 48?μg?l?1 with maximum concentrations of Si and minimum concentration of Cd. Correlation analysis suggested their contributions from both crustal and anthropogenic sources. 相似文献
5.
V. M. Ignatiev 《Geomagnetism and Aeronomy》2009,49(3):405-407
A rare case, when non-thermal profiles of the [OI] 557.7 nm line with the dissociative components shifted relative to an ordinary Doppler kernel appear in auroras, is considered. Based on an analysis of these profiles, it has been indicated that the dissociative component is shifted because the electric field is present during the recombination of O 2 + ion with background electrons in the ionospheric F region. The electric field component along the line of sight of the Fabry-Pérot interferometer (24 mV m?1) has been estimated using the Doppler shift of the dissociative component of the 557.7 nm profile emission as an example. 相似文献
6.
Sh. A. Mukhamediev 《Izvestiya Physics of the Solid Earth》2014,50(5):655-691
Rock masses contain ubiquitous multiscale heterogeneities, which (or whose boundaries) serve as the surfaces of discontinuity for some characteristics of the stress state, e.g., for the orientation of principal stress axes. Revealing the regularities that control these discontinuities is a key to understanding the processes taking place at the boundaries of the heterogeneities and for designing the correct procedures for reconstructing and theoretical modeling of tectonic stresses. In the present study, the local laws describing the refraction of the axes of extreme principal stresses T 1 (maximal tension in the deviatoric sense) and T 3 (maximal compression) of the Cauchy stress tensor at the transition over the elementary area n of discontinuity whose orientation is specified by the unit normal n are derived. It is assumed that on the area n of discontinuity, frictional contact takes place. No hypotheses are made on the constitutive equations, and a priori constraints are not posed on the orientation on the stress axes. Two domains, which adjoin area n on the opposite sides and are conventionally marked + and ?, are distinguished. In the case of the two-dimensional (2D) stress state, any principal stress axis on passing from domain ? to domain + remains in the same quadrant of the plane as the continuation of this axis in domain +. The sign and size of the refraction angle depend on the sign and amplitude of the jump of the normal stress, which is tangential to the surface of discontinuity. In the three-dimensional (3D) case, the refraction of axes T 1 and T 3 should be analyzed simultaneously. For each side, + and ?, the projections of the T 1 and T 3 axes on the generally oriented plane n form the shear sectors S + and S ?, which are determined unambiguously and to whose angular domains the possible directions p + and p ? of the shear stress vectors belong. In order for the extreme stress axes T 1 + ,T 3 + and T 1 ? , T 3 ? to be statically compatible on the generally oriented plane n, it is required that sectors S + and S ? had a nonempty intersection. The direction vectors p + and p ? are determined uniquely if, besides axes T 1 ? , T 3 ? and T 1 + , T 3 + , also the ratios of differential stresses R + and R ? (0 ≤ R ± ≤ 1) are known. This is equivalent to specifying the reduced stress tensors T R + and T R ? The necessary condition for tensors T R + and T R ? being statically compatible on plane n is the equality p + = p ?. In this paper, simple methods are suggested for solving the inverse problem of constructing the set of the orientations of the extreme stress axes from the known direction p of the shear stress vector on plane n and from the data on the shear sector. Based on these methods and using the necessary conditions of local equilibrium on plane n formulated above, all the possible orientations of axes T 1 + , T 3 + are determined if the projections of axes T 1 ? , T 3 ? axes on side — are given. The angle between the projections of axes T 1 + , T 1 ? and/or T 3 + , T 3 ? on the plane can attain 90°. Besides the general case, also the particular cases of the contact between the degenerate stress states and the special position of plane n relative to the principal stress axes are thoroughly examined. Generalization of the obtained results makes it possible to plot the local diagram of the orientations of axes T 1 + , T 3 + for a given sector S ?. This diagram is a so-called stress orientation sphere, which is subdivided into three pairs of areas (compression, tension, and compression-extension). The tension and compression zones cannot contain the poles of T 3 + and T 1 + axes, respectively. The compression-extension zones can contain the poles of either T 1 + or T 3 + axis but not both poles simultaneously. In the particular case when the shear stress vector has a unique direction p ? on side ?, the areas of compression-extension disappear and the diagram is reduced to a beach-ball plot, which visualizes the focal mechanism solution of an earthquake. If area n is a generally oriented plane and if the orientation of the pairs of the statically compatible axes T 1 ? , T 3 ? and T 1 + , T 3 + is specified, then, the stress values on side + are uniquely determined from the known stress values on side ?. From the value of differential stress ratio R ?, one can calculate the value of R +, and using the values of the principal stresses on side ?, determine the total stress tensor T + on side +. The obtained results are supported by the laboratory experiments and drilling data. In particular, these results disclose the drawbacks of some established notions and methods in which the possible refraction of the stress axes is unreasonably ignored or taken into account improperly. For example, it is generally misleading to associate the slip on the preexisting fault with the orientation of any particular trihedron of the principal stress axes. The reconstruction should address the potentially statically compatible principal stress axes, which are differently oriented on opposite sides of the fault plane. The fact that, based on the orientation of the intraplate principal stresses at the base of the lithosphere, one cannot make a conclusion on the active or passive influence of the mantle flows on the lithospheric plate motion is another example. The present relationships linking the stress values on the opposite sides of the fault plane on which the orientations of the principal stress axes are known demonstrate the incorrectness of the existing methods, in which the reduced stress tensors within the material domains are reconstructed without allowance for the dynamic interaction of these domains with their neighbors. In addition, using the obtained results, one can generalize the notion of the zone of dynamical control of a fault onto the case of the existence of discontinuities in this region and analyze the stress transfer across the system of the faults. 相似文献
7.
This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N+2, O+2, O+(4S), O+(2D), O−(2P), and O+(2P) ions, and O(1D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm−3 s−1 was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N2(v > 0) and O2(v > 0) in the calculations of the O+(4S) loss rate improves the agreement between the calculated Ne and the data on 4 August, 1972. The N2(v > 0) and O2(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O+ + N2 rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O+ + N2 loss rate due to vibrationally excited O2 produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O2(a1δg) and O2(b1σg+), and rotational excitation of O2, and found that the effect of these cooling rates on Te can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(3P) with thermal electrons with excitation of O(1D), and the electron cooling rates by vibrational excitation of O2 and N2 are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N2(B3πg+), N2(C3πu), and N2(A3σu+) at several vibrational levels, O(1S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm. 相似文献
8.
Concentrations of major ions in surface waters of the rivers of Khara-Murin and Snezhnaya are compared based on data of many-year observations carried out in the 1950s and 2000s. The concentrations of HCO 3 ? , Cl?, Ca2+, Mg2+, Na+ + K+ are shown to be stable. A considerable increase in SO 4 ?2 concentration was revealed. 相似文献
9.
A model of auroral electron deposition processes has been developed using Monte Carlo techniques to simulate electron transport and energy loss. The computed differential electron flux and pitch angle were compared with in situ auroral observations to provide a check on the accuracy of the model. As part of the energy loss process, a tally was kept of electronic excitation and ionization of the important atomic and molecular states. The optical emission rates from these excited states were computed and compared with auroral observations of (3914 Å), (5577 Å), (7620 Å) and (N2VK). In particular, the roles played by energy transfer from N2(A3+u) and by other processes in the excitation of O(1S) and O2(b1+g) were investigated in detail. It is concluded that the N2(A3+u) mechanism is dominant for the production of OI(5577 Å) in the peak emission region of normal aurora, although the production efficiency is much smaller than the measured laboratory value; above 150 km electron impact on atomic oxygen is dominant. Atomic oxygen densities in the range of 0.75±0.25 MSIS-86 [O] were derived from the optical comparisons for auroral latitudes in mid-winter for various levels of solar and magnetic activity. 相似文献
10.
V. B. Belakhovsky V. A. Pilipenko S. N. Samsonov D. Lorentsen 《Geomagnetism and Aeronomy》2016,56(1):42-58
Simultaneous morning Pc5 pulsations (f ~ 3–5 mHz) in the geomagnetic field, aurora intensities (in the 557.7 and 630.0 nm oxygen emissions and the 471.0 nm nitrogen emission), and riometer absorption, were studied based on the CARISMA, CANMOS, and NORSTAR network data for the event of January 1, 2000. According to the GOES-8 satellite observations, these Pc5 geomagnetic pulsations are observed as incompressible Alfvén waves with toroidal polarization in the magnetosphere. Although the Pc5 pulsation frequencies in auroras, the geomagnetic field, and riometer absorption are close to one another, stable phase relationships are not observed between them. Far from all trains of geomagnetic Pc5 pulsations are accompanied by corresponding auroral pulsations; consequently, geomagnetic pulsations are primary with respect to auroral pulsations. Both geomagnetic and auroral pulsations propagate poleward, and the frequency decreases with increasing geomagnetic latitude. When auroral Pc5 pulsations appear, the ratio of the 557.7/630.0 nm emission intensity sharply increases, which indicates that auroral pulsations result from not simply modulated particle precipitation but also an additional periodic acceleration of auroral electrons by the wave field. A high correlation is not observed between Pc5 pulsations in auroras and the riometer absorption, which indicates that these pulsations have a common source but different generation mechanisms. Auroral luminosity modulation is supposedly related to the interaction between Alfvén waves and the region with the field-aligned potential drop above the auroral ionosphere, and riometer absorption modulation is caused by the scattering of energetic electrons by VLF noise pulsations. 相似文献
11.
V. M. Ignat’ev 《Geomagnetism and Aeronomy》2014,54(1):94-98
The occurrence of anomalous (nonthermal) profiles of green emission of oxygen atoms detected with a Fabry-Perot spectrometer in auroras with the effect of a rapid decrease in the intensity of the wings of their dissociative component has been investigated. Based on an analysis of these measured profiles, it has been found that the characteristic time of recombination of a molecular oxygen ion at altitudes of 200–400 km is about 5–7 s. It appears that these molecular ions occur in a horizontally limited region of the auroral ionosphere as a result of ionization by a space localized flux of soft electrons with energies of 0.2–0.4 keV penetrating up to altitudes of 200 km. The estimation of the electron flux produces a value of 1010–1013 electrons cm?2 s?1. They generate the excess concentration n(O 2 + ) ~ 5.6 × 105 cm?3. 相似文献
12.
T. Sergienko I. Kornilov E. Belova T. Turunen J. Manninen 《Journal of Atmospheric and Solar》1997,59(18):2401-2407
Ionospheric heating experiments were done by the EISCAT Heater in Tromsø on 15–19 November, 1993. A low-light TV camera was installed at the VLF receiving station at Porojärvi about 100 km to the south-east of Tromsø. The spectral analysis of the auroral luminosity variations showed that the brightness of the aurora varied at the modulation frequency of the heating wave. The results of this analysis and the numerical simulations of the auroral luminosity variations caused by the HF heating are shown. The variations of the optical emission intensity at the heating frequency occur during the auroral ionosphere modification. The observed intensity variation of the auroral green line during the interval of enhanced electron temperature is explained by a decreasing rate of the O2+ ion dissociative recombination when the electron temperature increases. The brightness variation depends on the characteristic energy and the intensity of the auroral electron flux and the heating wave parameters. The artificial luminosity pulsations caused by HF heating are estimated. 相似文献
13.
Auroral electron transport calculations are a critical part of auroral models. We evaluate a numerical solution to the transport and energy degradation problem. The numerical solution is verified by reproducing simplified problems to which analytic solutions exist, internal self-consistency tests, comparison with laboratory experiments of electron beams penetrating a collision chamber, and by comparison with auroral observations, particularly the emission ratio of the N2 second positive to N+ 2 first negative emissions. Our numerical solutions agree with range measurements in collision chambers. The calculated N22P to N+ 21N emission ratio is independent of the spectral characteristics of the incident electrons, and agrees with the value observed in aurora. Using different sets of energy loss cross sections and different functions to describe the energy distribution of secondary electrons that emerge from ionization collisions, we discuss the uncertainties of the solutions to the electron transport equation resulting from the uncertainties of these input parameters. 相似文献
14.
A. V. Pavlov 《Surveys in Geophysics》2012,33(5):1133-1172
The current state of knowledge of E- and F-region ion chemistry is reviewed. Considerable attention is given to the progress in the chemistry of unexcited N2 +, O2 +, NO+, O+(4S), N+, H+, He+, Fe+, Mg+, Na+, Ca+, and K+ ions and electronically excited O+(2D), O+(2P), O+(4P), and $ {\text{O}}^{ + } (^{2} {\text{P}}^{*} ) $ ions. Achievements in our understanding of the role of vibrationally excited N2 +, O2 +, and NO+ ions in the ionosphere are discussed. 相似文献
15.
Feng Chen Da-Yuan Chen Qi-Ping Cao Shu-Jun Yu Dong-Jun Xu Cong-Xin Chen Yu-Liang Yu Jian-Hao Sheng 《地震科学(英文版)》1993,6(3):721-729
In the strip limestone mine in Guiding county, Guizhou Province the shear and frictionin situ tests of rock body were made for the three typical inclined weak bands C 3 1 /C 3 1 , C 3 1 /C 2 2 and C 2 2 /C 2 1 . The tests were made according to the second scheme of cuneate sample of the standards on rock mechanics test of Water Conservancy and Electricity Ministry. The changes of the resistivity in the weak band and the acoustic speed across the weak band were measured in the same time. The apparent resistivity data, obtained for 8 samples on 27 measure lines in 38 cycle tests, show that the apparent resistivity changes have rather obvious characters as follows: 1. At shear and friction stage, the change of the apparent resistivity accelerates after the yield point, and reaches the maximum of change rate and change amplitude near fracture point (except the lines with resistivity invariant); 2. On the same sample, the resistivity changes are different on the various lines and related to the location settled the lines, there are some “sensitive” location; 3. At the stage of preloading normal stress before shearing, the resistivity decreases on most lines, but on a few lines the resistivity does not changes; 4. After unloading shear stress, the resistivity could not recover completely and the hysteresis of resistiviity takes place on a few lines. 相似文献
16.
Jan Lidén 《Aquatic Sciences - Research Across Boundaries》1982,44(2):393-404
The apparent equilibrium constant, K 1 ′ , for the first acidity constant of carbonic acid has been determined in a lake water of 1.6 mM ionic strength in the temperature interval 4–18°C. The experimental method used comprises pH measurements in situ (NBS scale) with a probe and an IR-method for the selective determination of [H2CO3] and [HCO 3 ? ] in water samples. In the temperature interval studied the results can be described by the equation logK 1 ′ =126.39?6320.81/T?19.5682 ln T and are in agreement with the values of K 1 ′ obtained by the empirical equation presented by Millero [15]. This experimental agreement justifies thermodynamic calculations in the carbonate system based on field data. Also the experimental method described can be used to evaluate the acid-base balance of organic rich natural waters. 相似文献
17.
N. Iwagami T. Shibaki T. Suzuki H. Sekiguchi N. Takegawa W. H. Morrow 《Journal of Atmospheric and Solar》2003,65(16-18):1349-1360
Atomic oxygen density and airglow volume emission rate profiles measured in the rocket experiment S-310-29 carried out as a part of the Waves in Airglow Structures Experiment over Kagoshima in 2000 (WAVE2000) campaign are presented, and the excitation processes of the atomic oxygen 557.7 nm line and the molecular oxygen atmospheric band airglow emissions are discussed. The volume emission rate profiles calculated from the measured atomic oxygen densities using the methods and parameters proposed by the ETON campaign (EATON model) are found to well represent the shapes of measured twin-peak emission rate profiles seen during this campaign suggesting that the EATON model is valid for a perturbed atmosphere. There is some discrepancy in the modeled absolute values for the emissions. Applying the model to the current O density measurements results in predicted emission rates that are a factor of 3.2 and 1.5 too high for the 557.7 nm line and Atmospheric band, respectively. This suggests either that the atomic oxygen densities of the present campaign are too large by a factor of 1.2 (=1.51/2) to 1.5 (=3.21/3), or that those of the ETON campaign were too small by the same factor or that the combined errors in both campaigns can account for the discrepancy (the modeled volume emission rates of the 557.7 nm line and atmospheric band are roughly in proportion to [O]3 and [O]2, respectively). Our present data findings do not favor the 2-step more than the 1-step excitation process for the atmospheric band because a calculation of the quantum efficiency based on the observed O density does not show a steep gradient around 100 km. 相似文献
18.
Presently available data on the reaction of SO2 with OH radicals (OH + SO2 + \(M\xrightarrow[{k_1 }]{}\) HSO3 +M) are critically reviewed in light of recent stratospheric sulfur budget calculations. These calculations impose that the net oxidation ratek of SO2 within the stratosphere should fall within the range 10?7≤k≤10?9, if the SO2 oxidation model for the stratospheric sulfate layer is assumed to be correct. The effective reaction rate constantk 1 * =k 1[M] at the stratospheric temperature is estimated as $$k_1^* = \frac{{(8.2 \pm 2.2) \times 10^{ - 13} \times [M]}}{{(0.79 \mp 0.34) \times 10^{ - 13} + [M]}}cm^3 /molecules sec$$ where [M] refers to the total number density (molecules/cm3). Using the above limiting values ofk 1 * , and the estimated OH density concentrations, the net oxidation rate is calculated as 3.6×10?7≤k≤1.3×10?8 at 17 km altitude. This indicates that the upper limit of thesek values exceeds the tolerable range imposed by the model by a factor of about four. Obviously the uncertainty of thek 1 * values and of the OH concentrations in the stratosphere is still too large to make definite conclusions on the validity of the SO2 model. 相似文献
19.
Jennifer L. A. Hood William D. Taylor Sherry L. Schiff 《Aquatic Sciences - Research Across Boundaries》2014,76(2):243-258
The use of submersed macrophyte tissue δ 15N to quantify the level of WWTP effluent use in a highly urbanized and agricultural river was evaluated using several methods. Macrophytes, NH4 + and NO3 ? were collected by canoe along two 10 km reaches of river, upstream and downstream of two major municipal WWTPs over 3 years. NH4 + decreased in concentration while δ 15N–NH4 + increased as a function of distance downstream of both WWTPs, changing in one survey from 13 to 31 ‰ over 1 km. This increase is attributed to the combined effects of volatilization, nitrification and uptake. While NO 3 ? concentrations increased downstream of the WWTP over one of the survey reaches, δ 15N–NO 3 ? showed no prominent trend with distance at either. Macrophyte tissue δ 15N increased with distance downstream of both WWTPs, with a slope not significantly different from that of δ 15N–NH4 + suggesting that macrophytes incorporate effluent NH4 + as their main N source in those areas. However, mixing models suggest that towards the end of the reach, where source separation is distinct, macrophytes may utilize background NO 3 ? . Our study indicates the difficulty of deriving precise estimates of effluent use by macrophytes in a system where the δ 15N of the effluent changes rapidly. It also illustrates the utility of macrophytes in describing those changes where the effluent is too attenuated to allow for direct isotopic analysis. 相似文献
20.
B. K. Joshi 《Water Resources》2006,33(1):87-96
The study aims to measure the hydrological behavior and nutrients dynamics of the springs located in different landscapes of Kosi basin, Indian central Himalaya. A total number of eight springs were considered for the present investigation, each landscape represented by one spring. The monitoring for hydrological measurement was conducted in January 1998 to December 1999, the interval between two successive measurements was 10 days, i.e., 240 hr (total 72 observations were made). Water quality measurement was done in three main seasons, i.e., winter (Jan.), summer (June) and monsoon (Aug.) of 1998 and 1999, and the average value for measured parameters were calculated. These samples were analyzed for pH, electrical conductivity, total dissolved solids, dissolve oxygen, Ca2?, Mg2?, Na?, Cl?, F?, NO 3 ? , and SO 4 2? . Hydrology of spring’s water showed that the reserve forest has a higher water retention capacity than the other landscapes, and the spring recharge capacity highly influenced by the settlements, open grazing, mismanaged agricultural and deforestation activities. The spring water chemistry shows that the springs located in forest and sparsely populated areas have lower EC, TDS, cationic and anionic concentration and are safe for drinking purposes, but those in irrigated land and densely populated areas feature higher EC, low dissolved oxygen concentration and higher NO 3 ? , which makes the water of these springs unsuitable for drinking. F concentration was higher in the springs located in the settlement area. In brief, the study indicates that the unmanaged drains, very poor and old pattern of sewage disposing system result in an increase in Na?, Cl?, F?, NO 3 ? , and SO 4 2? concentration as compared to the springs in agricultural and forested areas. All of the studied springs are badly managed which a is great threat for the longevity and quality of the water bodies, in particular, in Indian Central Himalayan region. This study suggests the ways of the constructional works, grazing. Forest resource extraction and agricultural activities in water bodies catchments area should also be controlled. 相似文献