We extend the impact parameter dependent Glasma model (IP-Glasma) to three dimensions using explicit small x evolution of the two incoming nuclear gluon distributions. We compute rapidity distributions of produced gluons and the early time energy momentum tensor as a function of space-time rapidity and transverse coordinates. We study rapidity correlations and fluctuations of the initial geometry and multiplicity distributions and compare to existing models for the three dimensional initial state.