We investigate the finite temperature behavior of nonperturbatively improved clover fermions on lattices with temporal extent N t = 4 and 6. Unfortunately in the gauge coupling range, where the clover coefficient has been determined nonperturbatively, the finite temperature crossover/transition occurs at heavy pseudoscalar masses and large pseudoscalar to vector meson mass ratios. However, on an N t = 6 lattice the thermal crossover for the improved fermions is much smoother than for unimproved Wilson fermions and no strange metastable behavior is observed.Simulations with Wilson fermions suffer O(a) scaling violations due to the dimension five operator that Wilson introduced to give the unwanted fermion doublers masses of order the cutoff, 1/a. These scaling violations are much larger than those in the glue sector, which are O(a 2 ), and they can be numerically quite large, necessitating use of small lattice spacings at large simulation cost to get results that can be reliably extrapolated to the continuum limit.The O(a) scaling violations can be reduced to O(a 2 ) by introducing another dimension five operator into the fermion action, the so called clover term,where F µν (x) is a lattice transcription of the field strength tensor F µν (x), usually taken from four open plaquettes looking like a clover leaf, as proposed by Sheikholeslami and Wohlert [1]. For the reduction of scaling violations to work the clover coefficient, c sw , needs to be determined nonperturbatively as a function of the gauge coupling. The ALPHA collaboration developed a method to do so within the Schrödinger functional framework [2]. For quenched QCD the nonperturbative clover coefficient is now known for gauge coupling