This paper investigates the performance of seven methods of retrieving the planetary boundary layer height (PBLH) from lidar measurements carried out in the Metropolitan Region of S??ao Paulo (MRSP) during two MCITY-BRAZIL field campaigns of 2013. The performance is objectively assessed considering as reference the PBLH retrieved from rawinsonde carried out every 3 h during these campaigns. The role of clouds and aerosol load in the performance of the seven methods is analysed considering three case study scenarios representative of typical atmospheric conditions in the MRSP: (a) winter clean atmosphere, (b) summer low clouds and aerosol multilayers, (c) summer sea-breeze intrusion. Corroborating the case study results, the objective analysis indicated that most of the lidar methods retrieved PBLH closer to the top of the entrainment zone than the mixed layer, contradicting their definition. During daytime, the Wavelet Covariance Transform Method performs better than all the other six methods. The Inflexion Point Method performed better to estimate the Residual Layer height during night-time. In average, the diurnal evolution of the PBLH and its local rate of change based on lidar and rawinsonde measurements are in agreement.Conselho Nacional de Desenvolvimento Cient??fico e Tecnol??gico (CNPq)Funda????o de Amparo ?? Pesquisa do Estado de S??o Paulo (FAPESP)Funda????o de Amparo ?? Pesquisa do Estado do Rio de Janeiro (FAPERJ)Coordena????o de Aperfei??oamento de Pessoal de N??vel Superior (CAPES)CNPq: 309079/2013-6; 305357/2012-3; 462734/2014-5; 304786/2018-7; 154320/2018-8; 304786/2018-7FAPESP: 11/50178-5FAPERJ: E26/111.620/2011; E26/103.407/2012CAPES: 001