Adding confidence intervals to a gls regression trendline in R 4.0.2 on Mac OS 10.15.7

Question

I\'m using the gls function from nlme (version 3.1-149) to plot global surface temperatures since 1970 using the Cowtan-Way dataset. I\'d like to add the 95% confidence inte