Given the documented wave-induced damage of elevated coastal decks during extreme natural hazards (e.g., hurricanes) in the last two decades, it is of utmost significance to decipher the wave-structure-interaction of complex deck geometries and quantify the associated loads. Therefore, this study focuses on the assessment of solitary wave impact on open-girder decks that allow the air to escape from the sides. To this end, an arbitrary Lagrangian-Eulerian (ALE) numerical method with a multi-phase compressible formulation is used for the development of three-dimensional hydrodynamic models, which are validated against a large-scale experimental dataset of a coastal deck. Using the validated model as a baseline, a parametric investigation of different deck geometries with a varying number of girders Ng and three different widths, was conducted. The results reveal that the Ng of a superstructure has a complex role and that for small wave heights the horizontal and uplift forces increase with the Ng, while for large waves the opposite happens. If the Ng is small the wave particles accelerate after the initial impact on the offshore girder leading to a more violent slamming on the onshore part of the deck and larger pressures and forces, however, if Ng is large then unsynchronized eddies are formed in each chamber, which dissipate energy and apply out-of-phase pressures that result in multiple but weaker impacts on the deck. The decomposition of the total loads into slamming and quasi-static components, reveals surprisingly consistent trends for all the simulated waves, which facilitates the development of predictive load equations. These new equations, which are a function of Ng and are limited by the ratio of the wavelength to the deck width, provide more accurate predictions than existing empirical methods, and are expected to be useful to both engineers and researchers working towards the development of resilient coastal infrastructure.