Continuous Markov process theory is used to model classical thermal noise in two wire loops of resistances R(1) and R(2), Self-inductances L(1) and L(2), and absolute temperature T, which are coupled through their mutual inductance M. It is shown that even though the currents I-1(t) and I-2(t) in the two loops become progressively noisier as M increases from 0 toward its upper bound (L(1)L(2))(1/2), the fluctuation-dissipation, Nyquist, and conductance formulas all remain unchanged. But changes do occur in the spectral density functions of the currents lilt). Exact formulas for those functions are developed, and two special cases are examined in detail. (i) In the identical loop case (R(1) = R(2) = R and L(1) = L(2) = L), the M = 0 ''knee'' at frequency R/2 pi L in the spectral density function of I-i(t), below which that function has slope 0 and above which it has slope -2, is found to split when M > 0 into two knees at frequencies R/[2 pi(L +/- M)]. The noise remains white, but surprisingly slightly suppressed, at frequencies below R/[2 pi(L + M)], and it remains 1/f(2) at frequencies above R/[2 pi(L - M)]. In between the two knee frequencies a rough ''1/f-type'' noise behavior is exhibited. The sum and difference currents I +/- (t) = I-1(t) +/- I-2(t) are found to behave like thermal currents in two uncoupled loops with resistances R, self-inductances (L +/- M), and temperatures 2T. In the limit M --> L, I+(t) approaches the thermal current in a loop of resistance 1/2R and self-inductance L at temperature T, while I_(t) approaches (4kT/R)(1/2) times Gaussian white noise. (ii) In the weakly coupled highly dissimilar loop case (R(1) much less than R(2), L(1) = L(2) = L, and M much less than L), I-2(t) is found, to a first approximation, not to be affected by the presence of loop 1. But the spectral density function of I-2(t) is found to be enhanced for frequencies nu = R(2)/2 pi L by the approximate factor (1 + alpha nu(2)), where alpha = (2 pi M)(2)/R(1)R(2). A concomitant enhancement, by an approximate factor of (1 + 2M(2)R(2)/L(2)R(1))(1/2), is found in the high-frequency amplitude noise of I-1(t). An algorithm for numerically simulating I-1(t) and I-2(1) that is exact for all parameter values is presented, and simulation results that clarify and corroborate the theoretical findings are exhibited.