We present a study of the structural evolution of protein aggregates formed in solutions of a globular protein, β-lactoglobulin (BLG), in the presence of YCl3. These aggregates are often observed before crystallization starts and they are metastable with respect to the crystalline phase. Here we focus on the characterization of the hierarchical structure of this intermediate phase and its temperature dependent structure evolution using a combination of (very) small angle neutron and X-ray scattering (VSANS, SANS, and SAXS). Results show that the hierarchical structure ranges from nanometer scale protein monomer, dimer and compact protein clusters to micrometer scale fractal protein aggregates. Upon cooling, the overall hierarchical structure is preserved, but the evolution of the internal structure within the aggregates is clearly visible: the monomer-monomer correlation peak reduces its intensity and disappears completely at lower temperatures, whereas the cluster-cluster correlation is enhanced. At a larger length scale, the fractal dimension of protein aggregates increases. The kinetics of the structure change during a temperature ramp was further investigated using time-resolved SAXS. The time dependent SAXS profiles show clear isosbestic points and the kinetics of the structural evolution can be well described using a two-state model. These dynamic properties of protein aggregates on a broad length scale may be essential for being the precursors of nucleation.