Two-dimensional (2-D) (d, k) codes have been proposed for multichannel recording systems. These codes increase storage density by employing a multichannel timing recovery mechanism. In addition to the (d, k) constraint, some recording systems require a charge constraint to meet the channel requirements. In this study, capacities of 2-D charge constrained modulation codes are derived and computed from the binary and runlength graphs of the codes. The power spectral densities of codes are computed from their runlength graphs. Finally, a method to construct 2-D charge constrained modulation codes is provided.