Great thanks for this detailed response.
The last option is the route I am taking. As stated in my post I aligned all my runs regardless of session to the T1 from session 1 and things look good. I am plugging the runs from the different sessions into the same linear model treating them like runs from the same session. I can add the options you mentioned but how can I check to determine if I have problems with the motion correction stemming from having runs from different sessions? Also, what is lpa+zz? I found documentation for the other options you specified and I know lpa is local Pearson correlation absolute value cost function.
Thanks Again for your help