Geometric uncertainties caused by respiratory motion complicate radiotherapy treatment planning. Therefore 4D CT imaging is important in characterizing anatomy motion during breathing. Current 4D CT imaging techniques using multislice CT scanners involve multiple scans at several axial positions and retrospective sorting processes. Most sorting methods are based on externally monitored signals recorded by external monitoring instruments, which may not always accurately catch the actual breathing status and may lead to severe discontinuity artifacts in the sorted CT volumes. We propose a method to reconstruct time-resolved CT volumes based on internal motion to avoid the inaccuracies caused by external breathing signals. In our method, we iteratively sort the 4D CT slices using internal motion based breathing indices. In each iteration, respiratory motion is estimated by updating a motion model to best match a deformed reference volume to each moving multi-slice sub-volumes. The breathing indices as well as the reference volumes are refined for each iteration based on the currently estimated respiratory motion. An example is presented to illustrate the feasibility of our 4D CT sorting method without using any external motion monitoring systems.