The 1755 CE Lisbon earthquake triggered the largest historical tsunami ever impacting the Atlantic coasts of Europe. Despite recent efforts to better understand this event, there are still unanswered questions about the location of its epicenter and whether physical and historical evidences are in agreement.Inverse modeling using tsunami sediments can be applied to quantify onshore flow characteristics. Forward numerical modeling is also a powerful tool capable of simulating tsunami hydrodynamics and the induced sediment transport. This work presents novel results from a combination of inverse and forward modeling to assess tsunami characteristics onshore. The study site is located on the Portuguese southern coast, at the Salgados lowland where inverse modeling was performed using TsuSedMod (Jaffe and Gelfenbaum, 2007) based on data retrieved from sediment samples. Forward modeling, including tsunami generation and propagation, was performed using the FLOW module of Delft3D suite model. Onshore topography was corrected for the 1755 CE scenario based on extensive tsunami sedimentary deposit thickness data. The tsunami source was chosen based on recent results from the authors that pointed to a good correlation between modeled and field tsunami data for the Marques de Pombal Fault (MPF), Horseshoe Fault (HSF) and a hypothetical scenario represented by a simple combination between Gorringe Bank and Horseshoe Fault (Scenario 1 - SC1).Results from inverse model show tsunami onshore average speed varying from 7.3 up to 9.3 m/s and shear velocities from 0.52 up to 0.66 m/s. Forward modeling results show a wide variation according to the seismic source and tsunami onshore velocities can range from around 7 m/s when considering MPF to even an absence of inundation (SC1). The good agreement between both modeling approaches estimating tsunami velocity confirms the potential of numerical modeling coupled with geological records to improve the understanding of historical events.