We carried out numerical experiments of snow accumulation and depletion as well as snow albedo for a SNOw pack TELemetry (SNOTEL) site, using a one-dimensinal energy and mass balance snow model (SNTHERM) and a two-stream, multi-layer Snow, Ice, and Aerosol Radiation model (SNICAR). We developed SNTHERM_SNICAR model by replacing the snow albedo algorithm in SNTHERM with SNICAR algorithm, to improve understanding of snow accumulation and snowpack ablation processes and response from black carbon (BC) in snow. Model was tested against observed snow water equivalent (SWE) during the water year 2000 and 2002. Results show that both models do a very good job at predicting SWE, although underestimates and earlier melt-out happen. Snow albedo got from SNTHERM_SNICAR model is lower than that got from SNTHERM. Decrease in albedo, duo to BC mixed in snow, leads to faster snowmelt process and shorter duration of snow cover.