This work presents high-resolution ShakeMaps for three earthquakes occurred in the Betic Cordillera (SE Spain): the 2011 CE Lorca event (VIII ESI-07), the 1863 CE Huércal-Overa event (VIII ESI-07) and the 1829 CE Torrevieja event (X ESI-07). Detailed field characterizations and mapping of their coseismic environmental effects (EEEs) are catalogued and classified following the ESI-07 scale. The resulting macroseismic information reaches up to ten times the existing information based on conventional damage-based scales (e.g. EMS-98), providing a better constrain towards more realistic ground-motion scenarios. The 2011 Lorca earthquake has been used as a calibration event, since there is a relevant record of instrumental measures on source, and ground-motion parameters allowing a direct comparison with the modelled PGA values. From a methodological standpoint, the obtained ShakeMaps follow the basic guidelines and methodology proposed by the USGS ShakeMap Program. The two historical earthquakes analysed in this paper produced a wide variety of secondary EEEs but no surface faulting was reported. These effects need to be properly identified by high-resolution DTMs (5 m/pixel), far from the c. 900 m/pixel terrain models used by USGS program. Additionally, the proposed ESI-07 ShakeMaps incorporate correction factors to solve inconsistencies derived from the large scale terrain models considered in standard USGS program workflows: (1) empirical slope-derived Vs30 models result in overestimations of the PGA values in flat terrains in absence of unconsolidated deposits; (2) the topographic amplification factor included here, explains the occurrence of rock-falls and landslides in steep areas, where ground motion is underestimated by the sole use of slope-derived Vs30 models. Basic geological and geomorphological information need to be implemented in the modelling workflow in both cases. To prevent PGA overestimations in flat terrains a correction factor related to the spatial distribution and thickness of the Quaternary unconsolidated deposits has been incorporated (i.e. isopach maps). To correct PGA underestimations in steep terrains an amplification factor was modelled following standard guidelines of seismic topographic amplification.The comparison via iteration of the spatial distribution of both ESI effects and EMS macroseismic data, with the obtained ground-motion spatial distribution, enables a better definition of the geological parameters for the studied historical earthquakes. A more accurate location and/or size of the suspect seismic sources are obtained for these historical earthquakes, providing scenarios more realistic than those resulting from old intensity maps. Quaternary Geology and Geomorphology are behind the implementation of the proposed ESI-07 ShakeMaps, being especially useful when exploring historic or ancient moderate earthquakes with scarce damage-based macroseismic data, but with sufficient paleoseismic or archaeoseismic records. In summary, the methodological workflow proposed here contemplates the implementation of a computational configuration seismologically relevant and able to test repeated seismic scenarios, with different parametric data, in a controlled GIS environment useful to reproduce historical events.