This　study　focuses　on　the　analytical　prediction　of　subsurface　settlement　induced　by　shield　tunnelling　in　sandy　cobble　stratum　considering　the　volumetric　deformation　modes　of　the　soil　above　the　tunnel　crown.　A　series　of　numerical　analyses　is　performed　to　examine　the　effects　of　cover　depth　ratio　(C/D),　tunnel　volume　loss　rate　(ht)　and　volumetric　block　proportion　(VBP)　on　the　characteristics　of　subsurface　settlement　trough　and　soil　volume　loss.　Considering　the　ground　loss　variation　with　depth,　three　modes　are　deduced　from　the　volumetric　deformation　responses　of　the　soil　above　the　tunnel　crown.　Then,　analytical　solutions　to　predict　subsurface　settlement　for　each　mode　are　presented　using　stochastic　medium　theory.　The　influences　of　C/D,　ht　and　VBP　on　the　key　parameters　(i.e.　B　and　N)　in　the　analytical　expressions　are　discussed　to　determine　the　fitting　formulae　of　B　and　N.　Finally,　the　proposed　analytical　solutions　are　validated　by　the　comparisons　with　the　results　of　model　test　and　numerical　simulation.　Results　show　that　the　fitting　formulae　provide　a　convenient　and　reliable　way　to　evaluate　the　key　parameters.　Besides,　the　analytical　solutions　are　reasonable　and　available　in　predicting　the　subsurface　settlement　induced　by　shield　tunnelling　in　sandy　cobble　stratum.　(c)　2024　Institute　of　Rock　and　Soil　Mechanics,　Chinese　Academy　of　Sciences.　Production　and　hosting　by　Elsevier　B.V.　This　is　an　open　access　article　under　the　CC　BY-NC-ND　license　(http://creativecommons.org/　licenses/by-nc-nd/4.0/).