The objective of this article is to study the natural convection over an inclined wavy surface embedded in a porous medium saturated with nanofluid. To transform the complex wavy surface to a smooth surface, a coordinate transformation has been used. The governing equations are transformed into a set of partial differential equations using the non-similarity transformations and then applied the local similarity and non-similarity method to obtain ordinary differential equations. These resulting equations are solved using the successive linearization method. The effects of amplitude of the wavy surface, angle of inclination of the wavy surface to the horizontal, Brownian motion parameter, thermophoresis parameter on the non-dimensional heat and nanoparticle mass transfer rates are studied and presented graphically. © 2015 The Authors.