We use the current sample of ~10,000 zCOSMOS spectra of sources selected with I AB < 22.5 to define the density field out to z ~ 1, with much greater resolution in the radial dimension than has been possible with either photometric redshifts or weak lensing. We present the new algorithm that we have developed (ZADE) to incorporate objects not yet observed spectroscopically by modifying their photometric redshift probability distributions using the spectroscopic redshifts of nearby galaxies. We present a number of tests on mock catalogs used to justify this approach. The ZADE algorithm allows us to probe a broader range of galaxy environments and reduce the Poisson noise in the density field. The reconstructed overdensity field of the 10k zCOSMOS galaxies consists of cluster-like patterns surrounded by void-like regions, extending up to z ~ 1. Some of these structures are very large, spanning the ~50 h -1 Mpc transverse direction of the COSMOS field and extending up to Delta *Dz ~ 0.05 in redshift. We present the three-dimensional overdensity maps and compare the reconstructed overdensity field to the independently identified virialized groups of galaxies and clusters detected in the visible and in X-rays. The distribution of the overdense structures is in general well traced by these virialized structures. A comparison of the large-scale structures in the zCOSMOS data and in the mock catalogs reveals an excellent agreement between the fractions of the volume enclosed in structures of all sizes above a given overdensity between the data and the mocks in 0.2 < z < 1, although in the data these overdense regions are in generally larger contiguous structures.