In the present work, direct numerical simulations (DNS) were employed to understand the stability and morphology of laminar lean premixed hydrogen/air flames under various gravity conditions. The combined effect of Darrieus–Landau (DL) instability, thermal-diffusive (TD) instability and Rayleigh–Taylor (RT) instability was studied quantitatively for both the initial linear regime and the long-term nonlinear regime, with a particular focus on the effect of RT instability. In the linear regime, it was found that the RT instability has a negligible effect on the linear growth rate of flames with an equivalence ratio of 0.4, but the effect is significant when the equivalence ratio decreases to 0.3. The performance of various models for the dispersion relation was discussed, including a theoretical model considering the effects of all instabilities, which captures well the trend of the simulated results in the cases with unity Lewis numbers. The hydrogen/air flames with an equivalence ratio of 0.3 was investigated comprehensively in the nonlinear regime. For the cases with unity Lewis numbers, the RT-unstable flame exhibits a similar single-cusp structure as the RT-neutral flame, but with a deeper cusp and larger consumption speed as a result of the destabilizing effect of RT instability. The RT-stable flame gradually returns to a planar flame shape owing to the stabilizing effect of RT instability. For the cases with non-unity Lewis numbers, the effect of domain size on the flame evolution was investigated, showing a domain width of 200δT is enough to obtain a domain-independent result, where δT is the laminar flame thickness. It was found that the cases initialized with multimode perturbations reach the quasi-steady state faster, while the flame statistics are consistent with those of the cases initialized with single mode perturbations. The effect of RT instability was explored by comparing the results from various gravity conditions. The flame consumption speed is the largest in RT-unstable flames, and the smallest in RT-stable flames. The variation in flame consumption speed is due to the difference in the flame surface area, while the local reactivity and flame structure are not influenced by the RT instability. Further fractal analysis reveals that although the RT instability causes significant differences in flame surface area, the characteristic fractal dimension remains constant for flames under different gravity conditions.Novelty and significanceIn the present work, direct numerical simulations were employed to understand the stability and morphology of RT-unstable, RT-neutral and RT-stable laminar lean premixed hydrogen/air flames. For the first time, the effect of DL, TD and RT instabilities was quantitatively examined for both the linear regime and the nonlinear regime. A theoretical model considering various instabilities was compared with the simulated dispersion relations in the linear regime. In the non-linear regime, the flame surface area and corresponding flame consumption speed are influenced by the RT instability, while the flame stretch factor remains unaffected. A fractal analysis was conducted in laminar premixed hydrogen/air flames under the influence of various instabilities. Remarkably, it was discovered that the characteristic fractal dimension remains constant for the cases with and without the RT instability.