Under reasonable conditions, an initial/boundary value problem belonging to (1.9) with smooth data can be solved by standard methods at least for short times. For the corresponding problem on a four-or lower-dimensional manifold without boundary, Lamm [19] even proved long-time existence of solutions under the assumption that N has nonpositive sectional curvature. In general, however, it must be expected that a solution may develop singularities in finite time, although no such example is known. The aim of this paper is to examine the behavior of the flow as the first singular time is approached.Our results can be seen as a first step towards classifying the singularities of the biharmonic map heat flow, also with a view to ruling out certain types of singularities under special assumptions.Before we state the results, we briefly discuss some known facts about the harmonic map heat flow. The dimension two is critical for the harmonic map heat flow in the sense that E 1 is invariant under rescalings of the domain. Similarly, the dimension four is critical for the biharmonic map heat flow, and we may therefore expect some analogies between these two problems. For the initial/boundary value problem for the harmonic map heat flow on a two-dimensional domain, it was shown by Struwe [32] and Chang [4] that a weak solution exists under reasonable conditions. Moreover, this solution is smooth away from a finite number of points in time space. The singular points are characterized by a concentration of the energy, and at each of them, at least one nonconstant harmonic map from R 2 into N can be separated by a rescaling (or "blowup") procedure.Such a harmonic map gives rise to a harmonic 2-sphere (i.e., a harmonic map from the standard 2-sphere S 2 into N) via the stereographic projection. The harmonic spheres obtained thus are often called "bubbles," and the whole process is referred to as "bubblingoff." The mentioned results were refined by Ding and Tian [6], Qing and Tian [28], Wang [35], Topping [34], and others.For the biharmonic map heat flow, we need to study the concentration points for not one, but two energies, E 1 and E 2 . Consequently, we identify two distinct sets where singularities may occur. The dimension four is supercritical for E 1 , therefore the singular set belonging to concentrations of the first energy must be expected to be quite large in general. We can say, however, that its Hausdorff dimension is at most 2 (and we are able to derive a few other facts about its structure). At each point of this set, we find a harmonic map bubbling off. It turns out that these bubbles give rise either to harmonic 2-spheres or weakly harmonic 3-spheres in N. Concentration points of E 2 , on the other hand, admit a blowup that converges to either a (weakly) harmonic sphere of the same type or a biharmonic map. In analogy to the harmonic map heat flow, it might be suspected that this part of the singular set is finite. Unfortunately, we have no such result, owing to an insufficient control of the local evolution of...