In this work, we consider a mathematical model for describing flow in an unsaturated porous medium containing a fracture. Both the flow in the fracture as well as in the matrix blocks are governed by Richards' equation coupled by natural transmission conditions. Using formal asymptotics, we derive upscaled models as the limit of vanishing epsilon, the ratio of the width and length of the fracture. Our results show that the ratio of porosities and permeabilities in the fracture to matrix determine, to the leading order of approximation, the appropriate effective model. In these models the fracture is a lower dimensional object for which different transversally averaged models are derived depending on the ratio of the porosities and permeabilities of the fracture and respective matrix blocks. We obtain a catalogue of effective models which are validated by numerical computations. (C) 2019 Published by Elsevier Inc.