Snow avalanches cause fatalities and economic loss worldwide and are one of the most dangerous gravitational hazards in mountainous regions. Various flow behaviors have been reported in snow avalanches, making them challenging to be thoroughly understood and mitigated. Existing popular numerical approaches for modeling snow avalanches predominantly adopt depth-averaged models, which are computationally efficient but fail to capture important features along the flow depth direction such as densification and granulation. This study applies a three-dimensional (3D) material point method (MPM) to explore snow avalanches in different regimes on a complex real terrain. Flow features of the snow avalanches from release to deposition are comprehensively characterized for identification of the different regimes. In particular, brittle and ductile fractures are identified in the different modeled avalanches shortly after their release. During the flow, the analysis of local snow density variation reveals that snow granulation requires an appropriate combination of snow fracture and compaction. In contrast, cohesionless granular flows and plug flows are mainly governed by expansion and compaction hardening, respectively. Distinct textures of avalanche deposits are characterized, including a smooth surface, rough surfaces with snow granules, as well as a surface showing compacting shear planes often reported in wet snow avalanche deposits. Finally, the MPM modeling is verified with a real snow avalanche that occurred at Vallée de la Sionne, Switzerland. The MPM framework has been proven as a promising numerical tool for exploring complex behavior of a wide range of snow avalanches in different regimes to better understand avalanche dynamics. In the future, this framework can be extended to study other types of gravitational mass movements such as rock/glacier avalanches and debris flows with implementation of modified constitutive laws.