This study presents a pioneering exploration into the population dynamics of indigenous earthworm species in the Golaghat district of Assam, India, through the lens of mathematical modeling. Recognizing the integral role of earthworms in enhancing soil productivity and ecosystem health, we embarked on a detailed longitudinal analysis spanning from 2018 to 2023 to assess their population trends across different seasons and subdivisions within the district. Utilizing the principles of Mathematical Biology, we employed the Malthus Growth model and the Logistic Growth model to estimate future population trajectories of earthworm species, integrating biological insights with mathematical rigor. Our methodology, combining extensive field data collection with sophisticated mathematical modeling, provides a replicable framework for similar ecological studies. This study contributes to the growing field of Mathematical Biology, offering a novel approach to understanding the complex interactions within ecosystems and the impact of environmental changes on key species. The outcomes offer valuable insights into sustainable agricultural practices and biodiversity conservation in tropical and subtropical regions, emphasizing the critical role of indigenous earthworm species in maintaining soil health and ecosystem services. Our findings reveal significant seasonal variations and highlight the resilience of these species in the face of ecological changes. The models suggest distinct population growth patterns, with the Logistic model providing a more realistic projection considering environmental constraints and resource availability. The bifurcation diagram is drawn for the Logistic map that presents the phenomena inside the chaotic region.