The 3dcalc program will let you combine (conjoin) datasets on a voxel-wise basis in almost arbitrary ways. So if you want to see which voxels are above threshold in 2 separate tasks, you can use the "step" function and arithmetic multiplication. For more details, you'd have to explain your desires more fully.
bob cox