The highly transient nature of outgassing commonly observed at volcanoes is in part controlled by the permeability of lava domes and shallow conduits. Lava domes generally consist of a porous outer carapace surrounding a denser lava core with internal shear zones of variable porosity. Here we examine densification using uniaxial compression experiments on variably crystalline and porous rhyolitic dome lavas from the Taupo Volcanic Zone. Experiments were conducted at 900°C and an applied stress of 3 MPa to 60% strain, while monitoring acoustic emissions to track cracking. The evolution of the porous network was assessed via X-ray computed tomography, He-pycnometry, and relative gas permeability. High starting connected porosities led to low apparent viscosities and high strain rates, initially accompanied by abundant acoustic emissions. As compaction ensued, the lavas evolved; apparent viscosity increased and strain rate decreased due to strain hardening of the suspensions. Permeability fluctuations resulted from the interplay between viscous flow and brittle failure. Where phenocrysts were abundant, cracks had limited spatial extent, and pore closure decreased axial and radial permeability proportionally, maintaining the initial anisotropy. In crystal-poor lavas, axial cracks had a more profound effect, and permeability anisotropy switched to favor axial flow. Irrespective of porosity, both crystalline samples compacted to a threshold minimum porosity of 17-19%, whereas the crystal-poor sample did not achieve its compaction limit. This indicates that unconfined loading of porous dome lavas does not necessarily form an impermeable plug and may be hindered, in part by the presence of crystals.