Inclination-only dependent lunisolar resonances shape the dynamics of MEO (Medium Earth Orbit) objects over secular time scales (i.e. several decades). Their main effect is to increase an object’s eccentricity, possibly up to a value where the orbit’s perigee meets the Earth’s atmosphere and friction will determine the object’s re-entry. Thus, understanding this mechanism allows the design of low-cost end-of-life disposal strategies which exploit the resonant dynamics. In this proceeding, we will summarize our results in developing an analytic theory for lunisolar resonances and the characterization of diffusion properties along them. On this topic, the techniques proposed are of interest in most problems of secular resonances encountered in Celestial Mechanics.