A bstract The production cross-sections of ϒ (1 S ), ϒ (2 S ) and ϒ (3 S ) mesons in proton-proton collisions at $$ \sqrt{s}=13 $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>13</mml:mn></mml:math> TeV are measured with a data sample corresponding to an integrated luminosity of 277 ± 11 pb −1 recorded by the LHCb experiment in 2015. The ϒ mesons are reconstructed in the decay mode ϒ → μ + μ − . The differential production cross-sections times the dimuon branching fractions are measured as a function of the ϒ transverse momentum, p T , and rapidity, y , over the range 0 &lt; p T &lt; 30 GeV /c and 2 . 0 &lt; y &lt; 4 . 5. The ratios of the cross-sections with respect to the LHCb measurement at $$ \sqrt{s}=8 $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>8</mml:mn></mml:math> TeV are also determined. The measurements are compared with theoretical predictions based on NRQCD.