A wavelet-enriched adaptive hierarchical, coupled crystal plasticity – phase-field finite element model is developed in this work to simulate crack propagation in complex polycrystalline microstructures. The model accommodates initial material anisotropy and crack tension-compression asymmetry through orthogonal decomposition of stored elastic strain energy into tensile and compressive counterparts. The crack evolution is driven by stored elastic and defect energies, resulting from slip and hardening of crystallographic slips systems. A FE model is used to simulate the fracture process in a statistically equivalent representative volume element reconstructed from electron backscattered diffraction scans of experimental microstructures. Multiple numerical simulations with the model exhibits microstructurally sensitive crack propagation characteristics.