Clinically it is important to combine information provided by mammographic images from multiple views or at different times. Taking regular mammographic screening and comparing corresponding mammograms are necessary for early detection of breast cancer, which is the key to successful treatment. However, mammograms taken at different times are often obtained under different compressions, orientations or body positions. A temporal pair of mammograms may vary quite significantly due to the spatial disparities caused by the variety in acquisition environments, including the 3D position of the breast, the amount of the pressure applied, etc. Such disparities can be corrected through the process of temporal registration. We have implemented and utilized finite element models for temporal registration of digital mammography. In our work, we applied the patient specific breast model, where patients have both mammograms and MRIs available, and generic model, where only patient mammograms are available. After we applied the temporal registration algorithm, the average error among the 14 patient datasets was 3.4 plus/minus 0.86 mm for Euclidean distance and 4.3 plus/minus 0.52 mm for predicted 2D lesion position. With generic model, the average error among the 14 patient datasets using the measure of Euclidean distance between the predicted lesion position in T1 and T2 was 5.0 plus/minus 0.74 mm for Euclidean distance and 5.7 plus/minus 0.83 mm for predicted 2D lesion position. Compared with the average lesion size (10mm~40mm), this error is acceptable. With lesion correspondence, our finite element method can be used to suppress technical variations (e.g., mammogram positioning or compression) and to emphasize genuine alterations in the breast.