The processes that influence landscape evolution can be observed through instantaneous events such as earthquakes and over longer timescales in the response of rivers and changes in topography. The Northern Ecuadorian forearc is an active seismic area that has experienced several large magnitude earthquakes, the latest being the 7.8 Mw Pedernales earthquake, occurring in 2016. While seismic activity is well documented, longer-term effects of deformation using geomorphic metrics have not been studied systematically. Our objective is to use geomorphic metrics to study long-term deformation and examine the relationship between geomorphology and the spatial distribution of seismicity in the Northern Ecuadorian forearc. We analyze river drainage networks using metrics such as Stream Length-gradient, Normalized Steepness Index (Ksn), Hypsometric Integral, and Chi Analysis. To calculate these geomorphic metrics, we use high-resolution satellite data (4m resolution) from the Ministerio de Agricultura y Pesca and ALOS World 3D-30m DEM data from the Japan Aerospace Exploration Agency. Additionally, we examine seismic data from the Rapid Response Deployment after the Pedernales Earthquake, the High-Resolution Imaging of the Pedernales Earthquake Rupture Zone project, and the Permanent National Network of Ecuador. In our preliminary results, we have identified 89 drainage basins. The Esmeraldas drainage basin dominates the forearc and exhibits a well-developed drainage network extending into the interior of the forearc and the arc. In contrast, many drainage basins along the coast are small and contain rivers with fewer tributaries. Longitudinal profile analysis reveals knickpoints aligned with the W-E Canande fault and secondary parallel faults, accompanied by high Ksn values. River course alterations and knickpoints have also been identified in the structural high known as Businga Dome, along with high Ksn values along the knickpoints. We observe subsidence coinciding with seismic activity at 20-30 km depths at the northern down-dip edge of the Pedernales rupture. Finally, maps showing the distribution of the seismicity reveal an alignment between some surficial faults and seismic events, which may indicate a correlation with crustal faults.