Abstract We compute Euclidean action of charged de Sitter black holes in four dimensional gravitational Euler-Heisenberg model. It turns out that the action of a general Euclidean dyonically charged black hole is still controlled by the total entropy contributed by the black hole outer horizon and the cosmological horizon. For smooth configurations, the Euclidean action can be interpreted as the black hole production rate in de Sitter space. We show that the 4-derivative couplings break the symmetry between the production rate of the purely electric black hole and that of the purely magnetic black hole. Although electromagnetic duality is no longer a symmetry, it induces a transformation on the 4-derivative couplings, mapping the physical quantities of a purely electric black hole to those of a purely magnetic black hole and vice versa. We also observe that under the same transformation, unitarity constraints on the 4-derivative couplings remain invariant.