Nowadays, designs for ring-based light sources use multibend lattices for achieving a very small emittance of around 100 pmrad. In this type of storage ring, the chromaticity correcting sextupoles generally have greater strengths than those used in typical third-generation light sources. Therefore, controlling lattice nonlinearity such as amplitude-dependent tune shift (ADTS) is important for enabling stable operations and smooth beam commissioning. As the strength of the sextupoles increases, their higher-order terms contribute significantly to ADTS, rendering well-known lowest-order formulas inadequate for describing tune variations at large horizontal amplitudes. In response, we have derived explicit expressions of ADTS up to the fourth order in sextupole strength based on the canonical perturbation theory, assuming that the amplitude of a vertical betatron oscillation is smaller compared with the horizontal one. The new formulas express the horizontal and vertical betatron tune variations as functions of the action variables: ${J}_{x}$ and ${J}_{y}$ up to $O({J}_{x}^{2})$ and $O({J}_{y})$. The derived formulas were applied to a five-bend achromat lattice designed for the SPring-8 upgrade. By comparing the calculated results with the tracking simulations, we found that (1) the formulas accurately express ADTS around a horizontal amplitude of $\ensuremath{\sim}10\text{ }\text{ }\mathrm{mm}$ and (2) the nonlinear terms of the fourth order in sextupole strength govern the behaviors of circulating electrons at large horizontal amplitudes. In this paper, we present explicit expressions of fourth-order formulas of ADTS and provide some examples to illustrate their effectiveness.