An Individual-Based Model (IBM) was developed by employing natural and toxic survival rates of individuals to elucidate the community responses of benthic macroinvertebrates to anthropogenic disturbance in the streams. Experimental models (doseresponse and relative sensitivity) and mathematical models (power law and negative exponential distribution) were applied to determinate the individual survival rates due to acute toxicity in stressful conditions. A power law was additionally used to present the natural survival rate. Life events, covering movement, exposure to contaminants, death and reproduction, were simulated in the IBM at the individual level in small (1 m) and short (1 week) scales to produce species abundance distributions (SADs) at the community level in large (5 km) and long (1~2 years) scales. Consequently, the SADs, such as geometric series, log-series, and log-normal distribution, were accordingly observed at severely (Biological Monitoring Working Party (BMWP⁄10), intermediately (BMWP⁄40) and weakly (BMWP›50) polluted sites. The results from a power law and negative exponential distribution were suitably fitted to the field data across the different levels of pollution, according to the Kolmogorov-Smirnov test. The IBMs incorporating natural and toxic survival rates in individuals were useful for presenting community responses to disturbances and could be utilized as an integrative tool to elucidate community establishment processes in benthic macroinvertebrates in the streams.