The authors apply the composite focal mechanism method and use more than half million P wave polarity data to invert for the 3D stress field in Taiwan region. The authors present the focal mechanism result in a 0.25°×0.25°× 2.5 km grid in Taiwan and its adjacent areas. The preliminary result shows that the P and T axis directions change markedly across the Moho discontinuity. The P and T axis directions in the 15~45 km depth range also correlate with P wave velocity perturbations and delineate the northwestward subduction of the Philippine Sea Plate in Northeast Taiwan.