Living shorelines gain increasing attention because they stabilize shorelines and reduce erosion. This study leverages physics-based models and bagged regression tree (BRT) machine learning algorithm to simulate wave dynamics at a living shoreline composed of constructed oyster reefs (CORs) in upper Delaware Bay. The physics-based models consist of coupled Delft3D-FLOW and SWAN in four-level nested domains. The model accuracy converges with increasing mesh resolution. The simulated wave-induced current circulation substantiates the effectiveness of CORs in trapping sediments. The simulated yearly-averaged wave power correlates qualitatively with historical shoreline retreat rates. BRT is adopted to improve the model accuracy, identify key processes responsible for simulation errors in wave height (Hs) and wave period (Tp), and quantify their importance. In the CORs sheltered area, BRT reveals that simulation errors of wind seas mainly arise from wind forcing, wave breaking and wave triad interactions. Wave breaking is seven times more important than wind forcing for simulating Hs, while wind forcing and triad interactions are of equal importance for simulating Tp. Simulation errors of swells mostly stem from bottom friction and offshore wave boundary conditions. Results from this study can help the assessment and adaptive management of CORs-based living shoreline restoration projects under climate change.