Plasma-sprayed, thick-thermal-barrier coatings are very porous, approximately 2 mm thick, coatings used to insulate metal engine parts from high-temperature combustion processes. They are grown such that their elastic constants and density increase with depth; at the surface the elastic constants take small values, especially the shear modulus which is very small, but rapidly increase to those typical of the metal substrate. A Thomson–Haskell matrix technique is used to find numerically the eigenmode(s) along each vertical line through the coating and substrate, at each station along the direction of propagation. Using the lowest eigenmode, calculated at each station, a JWKB approximation to propagation in a coating with variable thickness is calculated, and from this the changes in travel time caused by varying degress of spalling at the surface are worked out. Using the two or three lowest eigenmodes for a coating of uniform thickness, a coupled-mode calculation is used to determine to what degree an interfacial inclusion couples the lowest eigenmode into the next higher one. In both cases the goal is to calculate the changes caused by the defects so that their presence may be detected. [Work supported by NSF.]