The class of arboreal singularities was first defined by the third author in the paper  [1]} by means of explicit local models. It was extended to signed arboreal singularities by Starkston in the paper [2]}. In the current paper, we address the problem of constructing canonical models within this class with the following key properties. First, they should be closed under the geometric operations set out in Definition REF  immediately below. Second, they should be determined up to ambient symplectomorphism by their combinatorial type. The class of arboreal singularities as studied in  [1]} and  [2]} satisfies the first but not the second property: those with the same combinatorial type may only be related by deformation. This is strong enough for many applications, for example calculating some invariants.
However, the homeomorphism type of an arboreal skeleta does not determine in general the symplectomorphism type of the ambient manifold even if the skeleton is smooth (e.g. see [5]}). On the other hand, it is not even possible to discuss the diffeomorphism type before establishing canonical models. Moreover, these models naturally appear in the dictionary between polarized Weinstein manifolds and arboreal skeleta established in [6]}. Due to the central role of these canonical models, we will typically refer to them as arboreal singularities; other members of the more general class will not play a role in this paper or its sequels.
