Skip to main content
Have a personal or library account? Click to login
Phytosociological patterns along a soil nutrient gradient in sacred groves of Haryana Cover

Phytosociological patterns along a soil nutrient gradient in sacred groves of Haryana

Open Access
|Jun 2026

Full Article

Introduction

Sacred groves (SGs) are small forest patches protected by local communities for religious or cultural reasons, often dedicated to deities or spirits, where activities like tree cutting or hunting are traditionally forbidden. These are revered and protected across cultures and continents for centuries, embodying a unique blend of biodiversity conservation, cultural heritage, and spiritual significance. They have been acknowledged by the International Union for the Conservation of Nature as Indigenous and Community Conserved Areas (ICCAs) (IUCN 2009). Usually located around a centre of devotion, the SGs are run by local people e.g., indigenous organizations or local communities who are the guardians of the surrounding forest (Dudley et al., 2010). SGs’ conferred biodiversity protection can thus help to supplement efforts at biodiversity conservation in formally identified protected areas more usually under government or NGO management (Klepeis et al., 2016).

SGs have been reported in diverse regions such as West Africa (Nigeria, Ghana), Southeast Asia (Thailand, Cambodia), and Japan (Shinto forests), in addition to their widespread presence across India. This highlights the global significance of sacred groves as both biodiversity hotspots and cultural landmarks. Sacred groves have often been preserved due to their role in maintaining plant diversity and cultural heritage (Gokhale, 2007; Kufuor & Omari, 2015). Sacred groves in India are most concentrated in regions such as the Western Ghats (Southern zone), the Northeast (e.g., Meghalaya), the Himalayas (North), and Rajasthan (West). This information provides a better understanding of the regional significance of sacred groves across the country (Bawa et al., 2004; Das & Ratha, 2013).

The SGs provide several advantages for nature conservation, such as the preservation of up to 85% of native species richness as refuges for rare and endemic species (Rösch et al., 2015). Other than this, local people can preserve their biodiversity for a long time because of cultural value that lasts for generations (Manna & Roy, 2021). They are also home to variety of medicinal plants (Ma et al., 2022), act as wildlife corridors or buffer zones for protected areas (Ishii et al., 2010), seed dispersal and pollination (Rajasri et al., 2017), erosion control and water resources (Ma et al., 2022).

Some SGs represent remnants of ancient, continuous forests (Scull et al., 2017), whereas others appear to be regenerated forests (Bhagwat et al., 2014). People are diverse and dynamic, and their religions and practices affect management of SGs (Dove et al., 2011). Such as the most widespread historical sites of Eurasian steppes are ‘kurgans’ (ancient burial mounds) which embody important historical, spiritual, cultural, and conservational values (Deák et al., 2019). Thus, SGs represent how humans and environment interact dynamically and are essential for preserving cultural values along with ecological assets across cultures and regions. Other than this, firewood, medicinal or ceremonial plants, and nontimber forest products including fruits and seeds can be found in SGs. They also host prayer, ceremonial, and ancestor worship (Lynch et al., 2018). Stewards of SGs can maintain high habitat quality, limit chronic and acute forest disturbance and facilitate passive restoration (Bhagwat et al., 2014). In India, SGs hold profound ecological and socio-cultural importance, often serving as repositories of traditional knowledge and biological diversity. These groves can range in size from small areas with a few trees to extensive hectares of greenery preserved due to their association with specific deity.

The present study focused on sacred groves (SGs) in Western Haryana, a semi-arid region in northern India with a sub-tropical climate, where forest cover is limited to only 3.53% of the total geographical area (ISFR, 2021). Although SGs are traditionally regarded as well-protected ecosystems, those in this region are increasingly vulnerable to anthropogenic pressures such as urbanization, land-use changes, overharvesting, pollution from religious activities, and expansion of villages into grove areas. The absence of tribal communities, who historically played a key role in safeguarding these groves, has further intensified these threats. Compounding the issue is a decline in traditional conservation practices and cultural values, leading to erosion of local stewardship over these ecologically valuable sites.

Despite their cultural and ecological significance, there is a lack of comprehensive prior ecological assessments. To address this critical gap, a preliminary survey was conducted to create an inventory of sacred groves across the region. Study sites were selected based on their ecological uniqueness, cultural relevance, and accessibility, with the objective of generating data that could inform targeted conservation strategies and support the long-term sustainability of these threatened ecosystems. Hence, the current study was conducted to investigate a total of four SGs from the four different forest ranges of western Haryana for a comprehensive ecological assessment to understand the phytodiversity and soil nutrient profile dynamics.

Materials and methods

Study site

A preliminary survey was conducted across Western Haryana to create an inventory of SGs, marking the first systematic documentation of these ecosystems in the region. Based on the preliminary survey across Western Haryana, four SGs—Bidola, Makrana Johra, Sultanpur, and Dhingsara—were randomly selected, one from each of the Tosham (Bhiwani), Behal (Bhiwani), Hansi (Hisar), and Fatehabad forest ranges, respectively, for in-depth ecological study (Figure 1). This selection method ensured a representative spatial distribution, enabling the generation of region-specific ecological insights while facilitating thorough field investigations.

Figure 1:

Map showing the location of Haryana in India and the SGs selected for the present study.

The study focused on sacred groves in the semi-arid region of Western Haryana, characterized by undulating sandy plains and bagar. The Thar Desert, which is situated in close vicinity, significantly influences the semiarid and dry climate of the region. The monsoon season is the one in which most of the annual rainfall in the research region is received. The map of study site and the climograph, are shown in Figure 1 and Figure 2.

Figure 2:

Climograph of the selected districts of Western Haryana showing mean average temperature, precipitation, and humidity (www.worldweatheronline.com).

Sampling of vegetation and data analysis

For the evaluation of distribution and quantification of phytodiversity, in-depth analyses were required in the study area, so several field visits were performed in the selected SGs. The study specifically utilized the widely recognized and important random quadrate sampling method (Hill, 2005). A total of 60 quadrats were plotted in the selected SGs i.e., 15 quadrats on each SG. Enumeration of trees and climbers was done in the quadrats of size 20×20 m whereas for shrubs and herbs, quadrats of size 5×5m and 1×1 m were used respectively (Cottam & Curtis, 1956). To measure the Circumference at Breast Height during sampling, the tree species girth was noted at 1.37 meters from the ground using measuring tape. For shrubs and climbers, the circumferences were taken at 5 cm from the ground. Whereas the diameter of the herbs was measured just above the ground using Vernier callipers.

Subsequently, the vegetational data was analysed using the phytosociological parameters viz., density (D), frequency (F), abundance (A) and basal area (B.A.) following Misra (1968). The relative values of F, D, and B.A. were calculated after that to obtain the IVI (Important Value Index) value for the encountered plant species following Phillips (1959) and Curtis (1959), using the formula: Important Value Index (IVI) = RF% + RD% + RDo%.

Other than this, to understand the distribution pattern of plant species, Abundance to Frequency ratio (A/F) was estimated following Cottam & Curtis (1956) for each species viz., regular (less than 0.025), random (0.025 to 0.05) and contiguous (more than 0.05). The frequency class distribution pattern was also obtained for the selected SGs to understand the nature (Homogenous or heterogenous) of plant communities occurring in them (Raunkiaer, 1918). Other than this, values of different vegetation indices were calculated for the diversity analysis of the four SGs. For this diversity index (Shannon & Wiener, 1963), dominance index (Simpson, 1949), index of evenness (Pielou, 1966) for species equitability and index for species richness (Margelef, 1958) were deliberated.

The soil samples were collected in the 4 SGs from each quadrat, taken at a depth of 0–30 cm. After removing any big stones, the soil samples were brought to the lab where they were first air dried and sieved (pore size – 2 mm). Using a conductivity meter, the electrical conductivity (EC) and soil pH of a saturated soil paste extract were measured (Rhoades, 1996; Thomas, 1996). Organic Carbon (OC) content of the soil samples was determined following Nelson & Sommers (1996). Soil Nitrogen content (N), Phosphorus content (P), and Potassium content (K) were analysed according to Jackson (1973), Olsen et al. (1954), and Pratt (1965), respectively.

Statistical analysis

Tukey post-hoc analysis was performed on the soil data and a box plot was formed using R program (R 4.4.1). Additionally, the two-tailed Carl-Pearson Coefficient was computed between the various floristic and soil parameters that were determined during the investigation and a heatmap was generated to facilitate comprehension of the correlation using R program (R 4.4.1).

Results and discussion

Phytodiversity

A total of 130 plant species were documented in this study, comprising 21 trees, 12 shrubs, 86 herbs, and 11 climbers, classified across 31 families. The Sultanpur SG exhibited the highest diversity of tree species (13), followed by Makrana Johra (11), Bidola (9), and Dhingsara (8). In terms of herbaceous species, Sultanpur SG also demonstrated the greatest richness, with 47 species recorded, while Makrana Johra SG led in the diversity of shrub (12) and climber species (8). These findings suggest a substantial degree of plant diversity within the region, attributable to the interplay of topographic, edaphic, and physiographic conditions.

The analysis revealed an uneven distribution of species across the encountered families, where approximately half of the species belonged to merely five families, with the remaining species distributed among 26 families. Notably, a significant number of families were represented by only a single species (Figure 3). The family Poaceae was identified as the most dominant, followed in succession by Apocynaceae, Fabaceae, Malvaceae, among others. The predominance of Poaceae aligns with findings by Dhiman et al. (2024) in the lower altitudinal ranges of Morni Hills, Haryana, and in the SG of Midnapore (West Bengal) as reported by Sen & Bhakat (2020), along with Harikesh et al. (2020) in the community forests of Haryana and Garg et al. (2020) in the semi-arid forests of Aravali Hills.

Figure 3:

Proportion of families covering encountered plant species during the present study (left) and graphs comparing the species-family richness in the four SGs (right).

Species-family richness analysis (Figure 3) indicated that Poaceae, Amaranthaceae, and Fabaceae were the most species-rich families in Bidola, Makrana Johra, and Dhingsara SG. Conversely, Sultanpur exhibited a maximum number of species within the Asteraceae family, followed by Amaranthaceae and Fabaceae. The dominance of Asteraceae has been corroborated in studies by Rashid et al. (2021), Dhiman et al. (2021), and Waheed et al. (2022), while Amaranthaceae’s prominence was highlighted by Prakash et al. (2022). This is because of the prevailing disturbances in these ecosystems as the endurance of Asteraceae family to tropical disturbance regimes places plant species belonging to this in an advantageous position to potentially dominate the disturbed ecosystems, also supported by Neto et al. (2017) and Arora et al. (2024).

The Poaceae family is arguably the most successful group of plants, characterized by their widespread presence in angiosperm habitats, ecological dominance, and notable species diversity. Their success can be attributed to several factors: their remarkable ability to colonize and persist in various environments, their effective long-distance dispersal mechanisms, and their ecological adaptability. Additionally, they exhibit tolerance to disturbances and have the capacity to modify ecosystems through processes such as fire and mammalian herbivory (Linder et al., 2017).

The prevalence of Fabaceae plants, following Poaceae, is characterized by their symbiotic nodulated roots and rhizobacteria, which contribute to the enrichment of soil with biologically accessible nitrogen through N2 fixation in the SGs (Singh et al., 2019; Joshi & Garkoti, 2020). Processes such as direct nutrient fixation, the deposition of organic matter from litter fall, root exudation, and rhizosphere aeration foster the activity of mutualistic aerobic microorganisms, thereby enhancing nutrient cycling and expanding the soil nutrient pool (Tang et al., 2018). Nutrient-rich soils facilitate the germination of seeds and saplings, while taxa that are unable to generate their own nutrient islands beneath their canopies rely on nutrients from other nutrient-fixing plants (Joshi & Garkoti, 2020).

The highest frequency of Acacia tortilis was recorded in Bidola, Makrana Johra, and Dhingsara within the tree stratum, while Acacia nilotica was predominantly observed in Sultanpur SG. This variation highlights their distinct niche preferences and ability to establish a presence across different geographic areas. While Bidola, Makrana Johra, and Dhingsara are situated in drier, higher elevations with sandy to loamy soils and lower moisture availability, Sultanpur is situated in a relatively low-lying area with finer alluvial soils and significantly higher soil moisture. The distribution and dominance of tree species are probably influenced by these site-specific differences in edaphic and microclimatic circumstances. Conversely, Capparis decidua exhibited the greatest frequency in Bidola, Sultanpur, and Dhingsara among shrubs, whereas Abutilon indicum was found most frequently in Makrana Johra (Table 14). These genera are prevalent in semi-arid zones and have been documented in numerous studies, such as those by Habib et al. (2016), Harikesh et al. (2020), Norman et al. (2024), Adoum (2024), and Arshad et al. (2024).

Additionally, the species distribution curve for the selected sacred groves (SGs) was analyzed using frequency classes (Raunkiaer, 1918; McIntosh, 1962). According to Raunkiaer’s law, species within a community can be classified as either common or rare. Any deviation from the typical J-shaped trajectory of a normal frequency distribution suggests an ecosystem disturbance. Raunkiaer’s normal species occurrence ratio (A > B > C >= D < E) indicates a homogeneous plant community when the frequency aligns with the J-shaped curve. The analysis demonstrates that the sacred groves of Bidola and Dhingsara conformed to Raunkiaer’s law, exhibiting a J-shaped species distribution curve indicative of a homogeneous plant community. In contrast, Makrana Johra and Sultanpur did not follow this pattern, suggesting greater heterogeneity (Figure 4). This observation is supported by Deil et al. (2021), who noted that in their study of the SGs in NW-Morocco, the proximity of sacred sites to intensively used agricultural landscapes in lowland areas correlates with a diminished conservation propensity, thereby indicating ecosystem disturbance.

Figure 4:

Frequency class distribution of plant species across the selected SGs.

In the four surveyed study groups (SGs), a considerable variation in species density was observed among trees, shrubs, herbs, and climbers (Figure 5). Notably, Bidola SG demonstrated the highest density across trees, shrubs, and herbs, while Makrana Johra exhibited the greatest stand density for climbers (Tables 1-4, Figure 5). The individual density of tree species ranged from 5 to 432 individuals per hectare (Ind./ha), 5 to 296.7 Ind./ha, 3.33 to 185 Ind./ha, and 3.33 to 743.33 Ind./ha in Bidola, Makrana Johra, Sultanpur, and Dhingsara SGs, respectively. Both Bidola and Dhingsara displayed higher tree densities (1115 Ind./ha and 881.67 Ind./ha, respectively) compared to Makrana Johra and Sultanpur SGs (503.3 Ind./ha and 533.3 Ind./ha, respectively). A higher tree density in Bidola indicates better habitat quality, while reduced densities in other SGs suggest ecological stress or degradation due to the influence of human pressures, such as lopping, trampling, and scraping, also supported by Aakash et al. (2019).

Figure 5:

Stand density (per hectare) of trees, shrubs, herbs and climbers in the four SGs.

In addition to the aforementioned analyses, abundance-to-frequency (A/F) values were calculated for each species across the four study groups (SGs). The present investigation revealed that all plant species conformed to a contiguous distribution pattern, as indicated by an A/F ratio surpassing 0.05. This contiguous distribution of plant species is frequently observed in natural forest ecosystems and has been documented in numerous studies (Kittur et al., 2013; Dhiman et al., 2020; Kumar & Verma, 2024).

Basal area (B.A.) analysis indicated that Sultanpur and Dhingsara exhibit significantly higher tree B.A. values (19.772 m2/ha and 15.507 m2/ha, respectively) compared to Bidola and Makrana Johra SG, which recorded a notably lower B.A. of 7.693 m2/ha (495). The elevated B.A. observed in Sultanpur can be attributed to the presence of tree species characterized by substantial girth classes, including Salvadora oleoides, Ficus benghalensis, and Ficus religiosa, among others. These findings are comparable to those reported by Meena et al. (2016), who observed a B.A. of 26.74 m2/ha in a similar semi-arid forest ecosystem in Delhi, indicating comparable vegetation structure and dominance of large-canopy species. Similarly, Yatar et al. (2024) documented a B.A. of 16.47 m2/ha in a semiarid landscape of Thailand, highlighting the influence of regional climatic conditions and species composition on basal area. The similarity in B.A. values across these studies suggests that structural parameters in semi-arid tree-dominated ecosystems are influenced by common ecological factors such as species traits, anthropogenic pressures, and edaphic conditions.

The Importance Value Index (IVI) is pivotal for understanding the ecological significance of various species, as a high IVI value indicates a species’ dominance within a community (Kashian et al., 2003). This metric serves to quantify the degree of dominance and assess the role and functionality of a species within the plant community structure. In the study area of Bidola, Acacia tortilis emerged as the most prevalent tree species, recording an IVI of 89.096. Conversely, in Makrana Johra and Sultanpur, Salvadora oleoides dominated with IVI values of 132.39 and 98.991, respectively. Additionally, Prosopis juliflora was identified as the most prevalent tree species in Dhingsara SG, with an IVI of 124.94.

Among shrub species, Capparis decidua was prominent across Bidola (131.46), Makrana Johra (105.76), and Dhingsara (135.65), showcasing the highest IVI values. In contrast, Parthenium hysterophorus attained the highest IVI in Sultanpur SG, with a value of 96.392. When examining herbaceous plants, Peristrophe bicalyculata was the most prevalent in Bidola (44.834), Makrana Johra (95.544), and Dhingsara (97.599), whereas Polygonum aviculare contributed to the dominance in Sultanpur SG with an IVI of 27.269. Furthermore, Cucumis callosus exhibited the highest IVI in Bidola (126.08), Makrana Johra (129.65), and Dhingsara (175.72), while Momordica dioca recorded a notable dominance as the climber species in Sultanpur SG, with an IVI of 209.45 (see Table 14).

Diversity indices

The selected study sites demonstrated a significant diversity of plant species, as evidenced by the H’ values for the identified vegetation in various study groups (SGs), which ranged from 1.26 to 1.935 for trees, 1.625 to 1.971 for shrubs, and 2.625 to 3.262 for herbaceous plants. Consequently, Sultanpur exhibited the highest level of species diversity. According to Kent & Coker (1992), the H’ value typically ranges between 1.5 and 3.5, seldom exceeding 4.5. An H’ index greater than 3.0 is categorized as high; values between 2.0 and 3.0 are considered medium; those ranging from 1.0 to 2.0 are classified as low; and values below 1.0 are designated as extremely low. Diverse ecosystems are also characterized by a broad spectrum of species populations, each possessing a variety of functional traits that contribute to the sustainability of ecosystem services throughout their lifespans (Himanshi et al., 2021).

Furthermore, Dhingsara SG exhibited the highest mean species dominance (0.142–0.385), suggesting a relatively uneven distribution of a few dominant species compared to other groves in the study. These values are higher than those reported by Dhiman et al. (2020) (0.085–0.252), indicating a potential shift in community structure, possibly due to site-specific disturbances or microclimatic factors. Nevertheless, the values still fall within the broader range observed across subtropical Indian forests (0.03–0.92; Malik & Bhatt, 2016; Singh et al., 2016; Saikia et al., 2017), suggesting that the groves retain natural heterogeneity despite anthropogenic pressures. A lower dominance value typically implies higher diversity, and this was reflected in other groves such as Sultanpur, which also showed greater evenness.

Evenness values, a measure of how evenly individuals are distributed among species, were highest in Sultanpur SG (0.755–0.883). These values are comparable to or slightly higher than those reported in other natural forests – for instance, Sarkar (2016) (0.73–0.85), Dhiman et al. (2020) (0.835-0.944), and Sherafu et al. (2024) (0.87). This high evenness suggests a stable and resilient community structure in Sultanpur SG, where no single species is overly dominant. In contrast, managed or degraded forests often show lower evenness due to dominance by a few stress-tolerant species. Thus, the observed evenness values reaffirm the ecological integrity and conservation value of these groves.

Sultanpur exhibited the highest species richness (0.317 to 4.586) among the selected sacred groves (SGs) in our study. But the Margalef Index values observed are lower than those reported in several other SG studies. For instance, in the sacred groves of Mahendergarh district, Haryana, Margalef Index values ranged between 5.78 and 8.61, indicating higher species richness (Choudhary et al., 2015). In contrast, our findings align more closely with those from community forests in southwest Haryana, where Harikesh et al. (2022) reported Margalef Index values ranging from 0.740 to 3.0146, suggesting comparable species richness under similar semi-arid conditions. The relatively lower species richness in our sites may be due to semi-arid conditions and anthropogenic pressures. However, species-rich groves like Sultanpur enhance regeneration, stabilize microclimate, and improve resilience. These findings underscore the ecological importance of SGs and the need for conservation strategies suited to semi-arid landscapes.

Soil nutrient profile

The selected SGs exhibited a substantial variation in soil nutrient profile, including pH, BD, EC, N, P K and OC (Tukey post-hoc analysis—p < 0.05). The value of soil pH varied from 6 to 6.72, EC – 0.14 to 0.326 dS/m, BD – 1.29 to 1.65 g/cm3, OC – 0.53 to 0.91%, N – 172.1 to 222.5 mg/g and P – 13.25 to 18.8 mg/g and K – 115.2 to 139 mg/g (Figure 7). The soil sample analysis revealed that Makrana Johra contained comparatively high levels of EC, N, P, K and OC. In contrast, Dhingsara SG exhibited elevated pH and soil BD (Figure 7). Whereas Bidola showed least amount of BD and EC. Minimum value of soil pH was observed in Makrana Johra SG. While Dhingsara SG had least amount of N, P, K and OC, as shown in Figure 7. Soil pH affects nutrient uptake as well as plant growth. It measures the acidity and alkalinity of soil samples along with controlling availability of numerous plant nutrients (Singh et al., 2023).

Figure 6:

Raincloud plots representing the values of diversity indices i.e., Shannon Weiner diversity index (H’), Simpson’s concentration of dominance (CD), Pielou index of evenness (E) and Margalef index of species richness (d) in the four SGs.

Figure 7:

Boxplot (Tukey post-hoc analysis) showing the comparative analysis of soil properties in selected SGs. Different letters above the box plots indicate statistically significant differences among sacred groves at p < 0.05, based on Tukey’s HSD post-hoc test.

Additionally, the two-tailed Carl-Pearson Coefficient was computed between the various floristic and soil parameters that were determined during the investigation. A heatmap was generated to facilitate comprehension of the correlation using R (Figure 8). The value of species dominance (CD) was negatively correlated with all these parameters, whereas a positive correlation was observed between Frequency (F), Density (D), Abundance (A), Shannon diversity value (H’), Margalef species richness (d), and electrical conductivity (EC). The value of soil Nitrogen (N) was also positively correlated with soil Organic Carbon (OC), Potassium (K), CD, Phosphorus (P), and EC, while it was negatively correlated with species evenness (E), pH, and bulk density (BD). A very small but positive correlation was observed between H’ and K. Akindele et al., (2021) also observed H’ to be positively correlated with K but available P, N contents, pH and OC were not correlated with any of the biodiversity variables. The positive correlation between H’ and EC is also supported by Malik & Haq (2022). But the negative correlation of H’ with OC and N, contrasts with Malik & Haq (2022) who observed it to be positive.

Figure 8:

Correlation of the selected parameters (floristic and edaphic) in the study. F = Frequency; D = Density; A = Abundance, BA = Basal Area; H’ = Shannon Wiener Index; CD = Simpson Index; E = Pielou Index, d = Margalef index, EC = Electrical conductivity, BD = Bulk density, N = Soil Nitrogen content, P = Soil Phosphorus content, K = Soil Potassium content, OC = Soil Organic Carbon.

Previous research has demonstrated that SGs, akin to forest fragments, exhibit high biodiversity and contribute significantly to ecosystem functions such as carbon sequestration (Parthasarathy & Naveen Babu, 2021), nutrient-rich soils (Dar et al., 2019b), improved water quality (Oliveira et al., 2017), groundwater recharge (Iftikhar Hussain et al., 2019), pathogen resistance (Quijas et al., 2010), and control of invasive species (Mace et al., 2019). Our findings support this view, as groves like Sultanpur showed notably high species richness, greater basal area values, and the presence of large-girthed native species such as Salvadora oleoides and Ficus benghalensis, indicating their ecological maturity and carbon storage potential. Additionally, the recorded diversity includes several native, medicinal, and potentially rare or regionally significant species, highlighting the groves’ role as in-situ conservation sites. A decline in floral diversity and species richness, as observed in degraded groves like Makrana Johra, can compromise ecosystem stability and services (Edrisi et al., 2020). Therefore, integrating scientific documentation of such sites with government-led conservation strategies and enhancing public awareness of their ecological and cultural significance is essential for long-term sustainability.

Conclusions

The lack of data on sacred groves (SGs) across various study sites significantly hampers our understanding of their biodiversity-protecting benefits. The findings of the current study reveal that sacred groves in semi-arid Haryana possess rich floristic diversity, with species occupying varied ecological niches and exhibiting wide ecological amplitude, enabling deeper insights into ecosystem dynamics beyond what forest cover alone can provide. However, these ecosystems are vulnerable to biotic disturbances and the impacts of global climate change. Haryana’s sacred groves face significant challenges, including habitat fragmentation, declining conservation practices, loss of local knowledge and inadequate management. They are threatened by encroachment, lack of legal protection, and youth disconnection from cultural traditions.

To ensure their viability, future conservation efforts must include awareness campaigns, legal protections, and biodiversity registers, with active community involvement in decision-making and management. While some sacred groves (SGs) are already under legal protection— either as deemed forests, community-managed lands, or recognized conservation sites—these provisions are often fragmented and inconsistently enforced. In many regions, the current legal status does not adequately prevent encroachment, resource exploitation, or ecological degradation. Moreover, a notable number of sacred groves across the country remain undocumented or lack formal recognition, further increasing their vulnerability. Therefore, it is crucial to initiate systematic identification and documentation of these groves and consider their designation as conservation reserves or community reserves under the Wildlife (Protection) Act, 1972 (Amended 2022). This enhanced legal status would offer structured protection while fostering community stewardship, supporting awareness initiatives, and promoting a conservation model that harmonizes ecological integrity with the preservation of cultural and spiritual values.

Acknowledgements

The authors would like to acknowledge CSIR, New Delhi, for providing financial support (Junior Research Fellowship) to Aman Mahla and the Department of Botany, Maharishi Dayanand University, Rohtak, for providing resources that made this research possible.

Notes

[1] Contributed by Author contributions

Conceptualization, Writing – original draft, methodology, Formal analysis and investigation, Writing – review and editing: [Aman Mahla]; Writing – review and editing, Formal analysis: [Himanshi Dhiman]; Writing - review and editing: [Harikesh Saharan]; Writing – review and editing, Supervision, Validation: [Anita R Sehrawat]

[2] Declarations

[3] Conflicts of interest Conflict of interest: The authors declare that there are no conflicts of interest.

[4] Ethical approval: Not applicable

[5] Informed consent: Not applicable

[6] Research data availability

This material is the authors’ own original work, which has not been previously published elsewhere and has no conflict of interest. The data generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Appendix

Table 1:

Analytical characteristics of plant species encountered in Bidola SG.

S. N.Name of Plant speciesDFAA/FBAIVI
Trees
1Acacia nilotica51020.20.1587.19
2Acacia tortilis432.510017.30.1730.84589.1
3Azadirachta indica51020.20.0254.31
4Balenites aegyptiaca17.51070.70.0345.58
5Dalbergia sissoo51020.20.1957.99
6Prosopis cineraria132.55010.60.2121.79667.08
7Prosopis juliflora367.530491.6330.86960.65
8Salvadora oleoides13040130.3250.36132.41
9Ziziphus jujuba20302.6670.0890.21216.33
Total1115105.64.495300
Shrubs
1Abutilon indicum332.56022.170.3690.010336.90
2Calotropis procera187.56012.50.2080.0017226.62
3Capparis decidua3509015.560.1730.49046131.46
4Capparis sepiaria32.5206.50.3250.012689.23
5Maytenus emarginata82.5704.7140.0670.0036923.14
6Parthenium hysterophorus50040501.250.0104342.03
7Ziziphus nummularia172.56011.50.1920.0290630.61
Total1657.5122.90.55833300
Herbs
1Achyranthes aspera837.510033.50.3350.0147820.52625
2Alternanthera sessilis7020140.70.000251.56365
3Alysicarpus vaginalis177.52035.51.7750.001973.35079
4Amaranthus palmeri17.5203.50.1750.001082.03256
5Amaranthus roxburghianus1753023.330.7780.001683.65783
6Argemone maxicana352070.350.000421.56828
7Boerhavia diffusa412.59018.330.2040.002268.31589
8Boerhavia erecta5020100.50.000181.43202
9Cenchrus ciliaris35040350.8750.001054.36089
10Chenopodium album12.51050.50.000140.70905
11Chenopodium murale5858029.250.3660.000737.18756
12Cleome viscosa51020.20.000140.68078
13Commelina benghalensis622.59027.670.3070.005411.63029
14Corchorus aestuans82.52016.50.8250.001642.72751
15Corchorus olitorius65308.6670.2890.001192.84946
16Croton bonplandianus3503046.671.5560.001464.14083
17Cyanthillium cinereum277.57015.860.2270.0071210.61259
18Cynodon dactylon8160100326.43.2640.0045339.89737
19Cyperus rotundus112540112.52.8130.000696.99343
20Dactyloctenium aegyptium752.55060.21.2040.000726.16265
21Digera muricata18040180.450.000082.94068
22Digitaris ciliaris52.510212.10.000050.78755
23Eragrostis tenella1697.59075.440.8380.0008312.01151
24Euphorbia granulata412.550330.660.000044.33453
25Euphorbia hirta232.530311.0330.000282.74984
26Euphorbia prostata7020140.70.000932.10997
27Evolvulus alsinoides10520211.050.000021.51082
28Heliotropium strigosum102.52020.51.0250.000051.5255
29Indigofera linnaei47.5209.50.4750.000151.39849
30Launaea nudicaulis17.5203.50.1750.000051.20505
31Mollugo nudicaulis12.51050.50.000010.60461
32Pedalium murex57.510232.30.001572.02758
33Peristrophe bicalyculata298570170.62.4370.0368744.72128
34Physalis minima302060.30.000141.32448
35Poa annua427.55034.20.6840.002076.02199
36Portulaca pilosa10520211.050.000291.72774
37Pupalia lappacia712.58035.630.4450.0122216.89938
38Setaria viridis381090169.31.8810.0012120.28098
39Sida cordifolia1253016.670.5560.000442.4731
40Spergula arvensis25020502.50.000422.37884
41Trianthema portulacastrum101040.40.000020.60322
42Tribulus terrestris603080.2670.000071.93079
43Verbena officinalis1853024.670.8220.000092.41811
44Verbesina encelioides267.54026.750.6690.0095610.88685
45Withania somnifera101040.40.000260.79604
46Xanthium strumarium312.550250.50.0077910.18392
47Zaleya pentandra85408.50.2130.001533.74747
Total265251710.10.12447300
Climbers
1Cucumis callosus145708.2860.1180.00101137.4887
2Momordica charantia55405.50.1380.0000535.26089
3Mukia maderspatana207.58010.380.130.0001999.35883
4Pergularia daemia252050.250.0001827.89157
Total432.529.1660.00143300

1 Abbreviations: F = Frequency, D = Density, A = Abundance, BA = Basal Area, A/F = Abundance to Frequency ratio, IVI = Importance Value Index.

Table 2:

Analytical characteristics of plant species encountered in Makrana Johra SG.

S. N.Name of Plant speciesDFAA/FBAIVI
Trees
1Acacia nilotica56.6730.450.2747.4994
2Acacia senegal8.3313.32.50.190.73817.137
3Acacia tortilis296.710011.8670.121.474122.22
4Anogeissus pendula3.336.6720.30.1745.8668
5Azadirachta indica56.6730.450.02234.2244
6Pongamia pinnata3.336.6720.30.000083.6046
7Prosopis cineraria3.336.6720.30.2657.0532
8Salvadora oleoides178.3808.91670.114.744132.39
Total503.335.2837.693300
Shrubs
1Abutilon indicum2558012.750.160.00542.82
2Calotropis procera103.373.35.63640.080.00527.328
3Capparis decidua15566.79.30.141.086105.76
4Capparis sepiaria106.6760.90.12811.453
5Grewia tenax3.3336.6720.30.0072.3968
6Gymnosporia senegalensis23.3326.73.50.130.11916.814
7Maytenus emarginata71.6753.35.3750.10.01120.057
8Opuntia dillenii156.6791.350.0486.2927
9Parthenium hysterophorus408.346.7350.750.02450.285
10Phyllanthus reticulatus202040.20.000746.651
11Solanum incanum11.67202.33330.120.000575.8734
12Ziziphus nummularia11.6713.33.50.260.000324.2688
Total108898.3951.436300
Herbs
1Achyranthes aspera518.386.723.9230.275930.004311.09814
2Aerva javanica68.332013.6670.683350.00192.89221
3Alternanthera sessilis71.6713.321.51.616540.00021.38994
4Amaranthus spinosus51.6726.77.750.290260.00434.78379
5Boerhavia diffusa26080130.162500.00188.04561
6Chenopodium album11.676.6771.049480.00020.67410
7Cleome viscosa3.3336.6720.299850.00010.54906
8Commelina benghalensis656.78032.8330.410410.005011.64586
9Corchorus olitorius356.773.319.4550.265420.00509.88787
10Croton bonplandianus38.3313.311.50.864660.00031.30901
11Cynodon dactylon9065100362.63.626000.005547.59270
12Cyperus rotundus16526.724.750.926970.00012.71134
13Dactyloctenium aegyptium281.726.742.251.582400.00113.79964
14Digera muricata3013.390.676690.00141.94341
15Digitaria ciliaris256.733.330.80.924920.00033.69364
16Eragrostis tenella112353.384.251.580680.00028.61657
17Erigeron bonariensis202040.200000.00011.60290
18Euphorbia granulata10520211.050000.00001.91009
19Euphorbia hirta48.3313.314.51.090230.00021.27070
20Evolvulus alsinoides11.676.6771.049480.00020.68625
21Evolvulus nummularius356.67213.148430.00000.64220
22Glinus lotoides541.720108.335.416500.00164.65364
23Heliotropium strigosum5513.316.51.240600.00011.26424
24Indigofera linnaei5026.77.50.280900.00012.24357
25Mollugo nudicaulis63.3333.37.60.228230.00032.88887
26Pedalium murex6.66713.320.150380.00031.19237
27Peristrophe bicalycualta4153100166.131.661300.116595.08587
28Perotis indica448.333.353.81.615620.00024.40616
29Phyllanthus fraternus63.3346.75.42860.116240.00033.86710
30Poa annua333.353.3250.469040.00166.29142
31Portulaca pilosa56.672011.3330.566650.00021.79981
32Pupalia lappacia433.38021.6670.270840.00188.78697
33Setaria viridis503780251.833.147880.002627.96850
34Sida cordifolia1013.330.225560.00011.06030
35Trianthema portulacastrum6.6676.6740.599700.000020.53182
36Tribulus terrestris1054010.50.262500.00043.60216
37Triumfetta rhomboidea6.66713.320.150380.00011.07693
38Verbesina encelioides11.676.6771.049480.00221.85880
39Xanthium strumarium18.3313.35.50.413530.00423.60920
40Zaleya pentandra6.66713.320.150380.00011.06709
Total24584.71484.8970.1646300
Climbers
1Citrullus colocynthis6.66713.320.150380.00012.150
2Coccinia grandis76.67407.66670.191670.00057.859
3Cucumis callosus191.793.38.21430.088040.00688.309
4Ipomoea pes-tigridis16546.714.1430.302850.000814.447
5Ipomoea triloba3.3336.6720.299850.00002.300
6Mukia maderspatana9566.75.70.085460.00025.786
7Pergularia daemia53.3333.36.40.192190.00036.593
Total591.746.1240.0088300

1 Abbreviations: F = Frequency, D = Density, A = Abundance, BA = Basal Area, A/F = Abundance to Frequency ratio, IVI = Importance Value Index.

Table 3:

Analytical characteristics of plant species encountered in Sultanpur SG.

S. N.Name of Plant speciesDFAA/FBAIVI
Trees
1Acacia nilotica1208060.0751.5286650.57
2Acacia tortilis18573.33310.10.1380.3096854.898
3Azadirachta indica3.336.666720.30.000062.3202
4Ficus beghalensis3.336.666720.30.766996.199
5Ficus religiosa3.3313.33310.0750.641737.2604
6Melia azedarach6.6713.33320.150.01044.6925
7Mimosa hamata11.676.666771.050.013383.9501
8Morus alba1033.3331.20.0360.000510.352
9Prosopis cineraria5533.3336.60.1980.4484121.055
10Prosopis juliflora100606.670.1110.2038235.035
11Salvadora oleoides28.3353.3332.130.0415.841698.991
12Syzygium cumini3.336.666720.30.006792.3543
13Ziziphus jujuba3.336.666720.30.000422.3221
Total533.350.719.7725300
Shrubs
1Abutilon indicum10546.66790.1930.0155434.855
2Calotropis procera121.746.66710.40.2230.0013925.132
3Capparis decidua113.366.6676.80.1020.0772189.714
4Capparis sepiaria3026.6674.50.1690.0107419.106
5Maytenus emarginata106.666760.90.000493.2597
6Parthenium hysterophorus823.366.66749.40.7410.0169396.392
7Phyllanthus reticulatus73.3313.333221.650.0021811.471
8Ziziphus nummularia51.67405.170.1290.004420.069
Total13281130.129300
Herbs
1Achyranthes aspera45573.33324.80.3380.032626.339
2Aerva javanica10526.66715.80.5910.000853.1068
3Alternanthera sessilis21013.333634.7250.000852.4833
4Amaranthus viridis1513.3334.50.3380.000031.164
5Anagallis arvensis96.6713.333292.1750.000051.5344
6Argemone maxicana1254012.50.3130.003755.9076
7Arnebia hispidissima3.3336.666720.30.000310.7283
8Boerhavia diffusa236.72047.32.3670.002143.8685
9Cenchrus ciliaris166.72033.31.6670.000262.5006
10Centaurea sp.248.34024.80.6210.00064.6741
11Chenopodium album19520391.950.00032.6476
12Chenopodium murale381.74038.20.9540.003947.1473
13Convolvulus arvensis363.353.33327.30.5110.004898.6822
14Coronopus didymus154833.3331865.5740.0242523.201
15Croton bonplandianus38066.66722.80.3420.005039.9146
16Cyanthillium cinereum111.74011.20.2790.000293.8966
17Cynodon dactylon40881001641.6350.0020927.301
18Cyperus rotundus483.326.66772.52.7190.000984.8464
19Datura innoxia61.6713.33318.51.3880.001041.938
20Eclipta alba91.676.6667558.250.000331.1303
21Erigeron canadensis188.32037.71.8830.001283.1721
22Euphorbia dracunculoides18.336.6667111.650.000150.7079
23Euphorbia granulata518.36034.60.5760.000867.6311
24Euphorbia prostrata148.313.33344.53.3380.000291.8965
25Evolvulus nummularius308.36.666718527.750.000061.9332
26Fagonia indica6.6676.666740.60.000090.6205
27Glinus lotoides169546.6671453.1130.0160220.288
28Gnaphalium purpureum158.34015.80.3960.001074.5426
29Heliotropium strigosum113.313.333342.550.000471.8482
30Launaea nudicaulis121.746.66710.40.2230.001415.1151
31Malvastrum coromandelianum141.76.66678512.750.000921.6861
32Orobanche aegyptiaca6.6676.666740.60.000080.6173
33Oxalis corniculata168.36.666710115.150.001011.8513
34Paspalidium dilatatum656.6667395.850.000150.9104
35Polygonum aviculare42824042810.70.0091527.269
36Pupalia lappacea6.6676.666740.60.00010.6286
37Ranunculus sceleratus348033.33341812.530.0252532.274
38Rumex dentatus258.32051.72.5830.002244.0205
39Sida cordifolia256.6667152.250.000941.1832
40Sissymbrium irio651.753.33348.90.9160.000487.467
41Solanum nigrum3.3336.666720.30.000460.8171
42Solanum virginianum88.332017.70.8830.008286.6828
43Sonchus oleraceus3.3336.666720.30.000060.5908
44Spergula arvensis233.326.667351.3130.000223.313
45Vernonica anagallis-aquatica378.36.666722734.050.00263.675
46Withania somnifera8520170.850.0159310.979
47Zaleya pentandra176.733.33321.20.6360.003175.2682
Total2269729290.17731300
Climbers
1Citrullus colocynthis8.3336.666750.750.0000590.886
2Momordica dioica1513.3334.50.3380.00019209.45
Total23.3339.50.00024300

1 Abbreviations: F = Frequency, D = Density, A = Abundance, BA = Basal Area, A/F = Abundance to Frequency ratio, IVI = Importance Value Index.

Table 4:

Analytical characteristics of plant species encountered in Dhingsara SG.

S. N.Name of Plant speciesDFAA/FBAIVI
Trees
1Ailanthus excelsa11.66713.333.50.260.70310.735
2Acacia leucophloea106.66760.90.6867.9981
3Acacia nilotica21.66746.671.8570.040.74824.354
4Acacia tortilis58.33353.334.3750.081.87538.222
5Azadirachta indica6.6667201.3330.070.59911.937
6Casuarina equisetifolia3.33336.66720.30.3575.1201
7Dalbergia sissoo3.33336.66720.30.3635.1599
8Ficus religiosa6.666713.3320.151.80217.254
9Prosopis cineraria3.333313.3310.080.256.8664
10Prosopis juliflora743.336049.560.832.891124.94
11Salvadora oleoides13.33333.331.60.055.22647.409
Total881.6775.2215.507300
Shrubs
1Abutilon indicum61.667406.1670.150.0035932.573
2Calotropis procera6.666713.3320.150.000898.0534
3Capparis decidua118.33607.8890.130.24119135.65
4Grewia tenax3.33336.66720.30.007676.5329
5Maytenus emarginata21.667204.3330.220.0013814.363
6Opuntia dillenii23.3336.667142.10.0118911.88
7Parthenium hysterophorus248.3326.6737.251.40.0195167.891
8Phyllanthus reticulatus28.33326.674.250.160.0000918.436
9Ziziphus nummularia3.33336.66720.30.002154.6186
Total51579.890.28836300
Herbs
1Achyranthes aspera548.3373.3329.910.410.0025611.194
2Aerva javanica46.66733.335.60.170.003175.5104
3Alternanthera pungens56.66730.450.000070.6483
4Alternanthera sessilis41.6676.667253.750.000240.9805
5Anisomeles indica16.66713.3350.380.000261.4229
6Argemone maxicana56.66730.450.000210.7525
7Cenchrus biflorus118.3326.6717.750.670.001143.7877
8Cenchrus ciliarias41.6676.667253.750.000060.8407
9Chenopodium album7533.3390.270.000223.4147
10Chenopodium murale6013.33181.350.000051.5024
11Commelina benghalensis43533.3352.21.570.0106313.312
12Convolvulus prostratus16.66713.3350.380.000031.2536
13Corchorus aestuans1513.334.50.340.000121.3111
14Croton bonplandianus601.6773.3332.820.450.0105417.568
15Cynodon dactylon5651.7100226.12.260.0032141.852
16Cyperus rotundus88.33313.3326.51.990.000431.9458
17Dactyloctenium aegyptium1153.386.6753.230.610.0041916.882
18Digera muricata523.338026.170.330.0011910.584
19Digitaria ciliaris46553.3334.880.650.000537.4898
20Echinochloa crus-galli251.674025.170.630.0070510.158
21Eragrostis tenella731.6753.3354.881.030.00038.7709
22Erigeron bonariensis18.3336.667111.650.000070.7235
23Euphorbia granualata24026.67361.350.000023.6018
24Heliotropium strigosum23.33313.3370.530.000131.3665
25Indigofera linnaei16.6676.667101.50.000080.7168
26Oxalis corniculata38.3336.667233.450.000010.7808
27Paspalidium flavidum21.6676.667131.950.000030.7065
28Peristrophe bicalyculata6293.386.67290.53.350.0732897.599
29Physalis minima46.667209.3330.470.002533.8854
30Poa annua16.6676.667101.50.00010.7331
31Portulaca pilosa56.66730.450.000010.6049
32Pupalia lappacea19040190.480.001085.271
33Setaria viridis13020261.30.000092.4867
34Sida cordifolia196.6746.6716.860.360.004348.3613
35Solanum xanthocarpum36.66713.33110.830.001312.3375
36Tephrosia purpurea68.3332013.670.680.000182.2165
37Trianthema portulacastrum152030.150.000221.9576
38Triumfetta rhomboidea31.667206.3330.320.000752.4473
39Verbesina encelioides11.6676.66771.050.000771.22
40Xanthium strumarium8.3333201.6670.080.000071.8015
Total1830012000.13128300
Climbers
1Citrullus colocynthis1013.3330.230.0004444.47
2Cucumis callosus223.3393.339.5710.10.00064175.72
3Ipomoea pes-tigridis26.667205.3330.270.0001229.301
4Merremia aegyptia8.333313.332.50.190.0000615.207
5Momordica diocia13.3336.66781.20.0000310.231
6Mukia maderspatana23.33326.673.50.130.0000325.075
Total30531.90.00132300

1 Abbreviations: F = Frequency, D = Density, A = Abundance, BA = Basal Area, A/F = Abundance to Frequency ratio, IVI = Importance Value Index.

DOI: https://doi.org/10.3986/hacq-2025-0024 | Journal eISSN: 1854-9829 | Journal ISSN: 1581-4661
Language: English, Slovenian
Page range: 25 - 45
Submitted on: Mar 1, 2025
Accepted on: Aug 20, 2025
Published on: Jun 3, 2026
Published by: Slovenian Academy of Sciences and Arts
In partnership with: Paradigm Publishing Services

© 2026 Aman Mahla, Himanshi Dhiman, Harikesh Saharan, Anita Rani Sehrawat, published by Slovenian Academy of Sciences and Arts
This work is licensed under the Creative Commons Attribution 4.0 License.