In partly linear models, the dependence of the response <italic>y</italic> on (<bold>x</bold> -super-T, <italic>t</italic>) is modeled through the relationship <italic>y</italic>=<bold>x</bold> -super-T <bold>β</bold>+<italic>g</italic>(<italic>t</italic>)+ϵ, where ϵ is independent of (<bold>x</bold> -super-T, <italic>t</italic>). We are interested in developing an estimation procedure that allows us to combine the flexibility of the partly linear models, studied by several authors, but including some variables that belong to a non-Euclidean space. The motivating application of this paper deals with the explanation of the atmospheric SO<sub>2</sub> pollution incidents using these models when some of the predictive variables belong in a cylinder. In this paper, the estimators of <bold>β</bold> and <italic>g</italic> are constructed when the explanatory variables <italic>t</italic> take values on a Riemannian manifold and the asymptotic properties of the proposed estimators are obtained under suitable conditions. We illustrate the use of this estimation approach using an environmental data set and we explore the performance of the estimators through a simulation study.
