“…A great many approaches are available in the literature on the general topic of model reduction, including aggregation methods [1], balancing techniques [2], Hankel norm approximation methods [3], H 1 norm approximations [4], and q-Markov covariance equivalent realizations [6], [7], [14], [15], to name just a few. A major drawback of each of these methods (with the exception of the q-Markov covariance equivalent realizations) is that the reduced-order models are not guaranteed to match any of the second-order information (i.e., covariance values) of the original model outputs, which is an important criterion when output performance is an item of interest.…”