My thought - there's so many other unknowns that get represented by assumptions in the system models that difference in impedance due to sag isn't worth worrying about.
The rise and fall of the mid-span point can be impressive, but the overall length doesn't change by much at all.
Do you calculate your fault currents with all possible combinations of generation running? We certainly don't, there's a bit of error there right off the bat; fortunately most of "source impedance" is transformers and other lines and not generation unless one is very close to generation. I've found that doing impedance calcs assuming straight line conductors produces results that are "good enough". If you truly want to improve the match between modeled results and actual fault results, I'd concentrate of improving your representation of the zero sequence environment rather than trying to add one more significant digit to the impedance values. The text book approach to zeros sequence modeling may produce a Z0 value as much as 120% of what it actually is, as measured by fault analysis. If you can make a 20% improvement or tweak a fraction of a per cent, which is better? Obviously that 20% improvement is based on a number of factors, some or many of which won't apply in your case; but personally I'd try to get a better ρearth value than 100Ωm long before going after sag effects.
But, pick a temperature value that makes sense in your conditions, document what you're doing and march forth. Would you use the same temperature for ACSS as you do for ACSR and AAC? Yes or no are probably equally valid answers depending on your circumstances. Do you want the maximum impedance for worst case voltage drop during high power flows, or do you want the minimum impedance for worst case fault currents?