Quantifying the link between microstructure and effective elastic properties of snow, firn, and bubbly ice is essential for many applications in cryospheric sciences. The microstructure of snow and ice can be characterized by different types of fabrics (crystallographic and geometrical), which give rise to macroscopically anisotropic elastic behavior. While the impact of the crystallographic fabric has been extensively studied in deep firn, the present work investigates the influence of the geometrical fabric over the entire range of possible volume fractions. To this end, we have computed the effective elasticity tensor of snow, firn, and ice by finite-element simulations based on 391 X-ray tomography images comprising samples from the laboratory, the Alps, Greenland, and Antarctica. We employed a variant of Eshelby's tensor that has been previously utilized for the parameterization of thermal and dielectric properties of snow and utilized Hashin-Shtrikman bounds to capture the nonlinear interplay between density and geometrical anisotropy. From that we derive a closed-form parameterization for all components of the (transverse isotropic) elasticity tensor for all volume fractions using two fit parameters per tensor component. Finally, we used the Thomsen parameter to compare the geometrical anisotropy to the maximal theoretical crystallographic anisotropy in bubbly ice. While the geometrical anisotropy clearly dominates up to ice volume fractions of φ≈0.7, a thorough understanding of elasticity in bubbly ice may require a coupled elastic theory that includes geometrical and crystallographic anisotropy.