This paper continues the development of a recently proposed resolvent-based model designed to capture the full second-order statistics of turbulent jets, which are required to obtain accurate noise estimates. The model requires an approximation of the crossspectral density tensor of certain nonlinear forcing terms, and the focus of this paper is to characterize the properties of these statistics in a high-Reynolds-number subsonic jet. We show that the power spectral density of the forcing is independent of frequency over a range of almost two orders-of-magnitude. The coherence of the forcing consists of peaks that are spatially compact compared to the coherence length-scales of the flow variables. The widths of these peaks depend on spatial location but not frequency, while the streamwise and radial wavelengths of the coherence depend on frequency but not spatial location. We propose a simple fit function in frequency space that captures these properties and show that it leads to good approximations of the LES forcing statistics. Some of the parameters in the model are well-approximated by quantities that could be obtained from a Reynolds-averaged Navier-Stokes simulation. Finally, we show that the properties of the forcing statistics are completely different for a low-Reynolds-number jet, which may be indicative of direct nonlinear interactions amongst wavepackets which are absent in the high-Reynolds-number jet.