Skeletal model reduction based on local sensitivity analysis of time dependent systems is presented in which sensitivities are modeled by forced optimally (f-OTD) modes. The f-OTD factorizes the coefficient matrix into a compressed format as product two skinny matrices, i.e. modes and coefficients. create low-dimensional, dependent, orthonormal basis capture directions phase space associated wi...