Given a graph = ( , ) with arboricity , we study the problem of decomposing the edges of into (1 + ) disjoint forests in the distributed LOCAL model. While there is a polynomial time centralized algorithm for -forest decomposition (e.g. [Imai, J. Operation Research Soc. of Japan '83]), it remains an open question how close we can get to this exact decomposition in the LOCAL model. Barenboim and Elkin [PODC '08] developed a LOCAL algorithm to compute a (2 + ) -forest decomposition in ( log ) rounds. Ghaffari and Su [SODA '17] made further progress by computing a (1 + ) -forest decomposition in ( log 3 4 ) rounds when = Ω( √︁ log ), i.e., the limit of their algorithm is an ( + Ω( √︁ log ))forest decomposition. This algorithm, based on a combinatorial construction of Alon, McDiarmid & Reed [Combinatorica '92], in fact provides a decomposition of the graph into star-forests, i.e., each forest is a collection of stars.Our main goal is to reduce the threshold of in (1 + ) -forest decomposition. We obtain the following main results:)-round algorithm when = Ω (1) in simple graphs and multigraphs, where > 0 is any arbitrary constant.• An ( log 4 log Δ )-round algorithm when = Ω( log Δ log log Δ ) in simple graphs and multigraphs.• An ( log 4)-round algorithm when = Ω(log ) in simple graphs and multigraphs. This also covers an extension of the forest-decomposition problem to list-edge-coloring.• An ()-round algorithm for star-forest decomposition for = Ω( √︁ log Δ + log ) in simple graphs. When ≥ Ω(log Δ), this also covers a list-edge-coloring variant.Our techniques also give an algorithm for (1 + ) -outdegreeorientation in (log 3 / ) rounds, which is the first algorithm with linear dependency on −1 .At a high level, the first three results come from a combination of network decomposition, load balancing, and a new structural result on local augmenting sequences. The fourth result uses a