Atmospheric flow and temperature dynamics in the urban roughness sublayer exhibit numerous complexities that cannot all be investigated using models or scaled-down experiments, thus these complexities necessitate careful field observations. Such dynamics were studied under a comprehensive set of wind directions, wind speeds, and thermal stability conditions in a field campaign held in Guelph, Ontario, Canada, from 13th to 25th of August 2017. The urban site was a quasi-two-dimensional canyon with unit canyon aspect ratio. Beside characterizing thermal stability, inertial effects, urban heat island intensity, and mean properties of the atmosphere, the turbulence statistics were studied carefully as functions of roof-level wind angle, diurnal time, the building Reynolds number, or the bulk Richardson number. Turbulence statistics in the vertical direction were observed to be influenced by local conditions, such as flow properties and nearby surface temperatures. These statistics also indicated presence of small integral lengthscales and short integral timescales. On the other hand, turbulence statistics in the horizontal directions were influenced by non-local conditions, such as horizontal heterogeneity in heating or urban morphology. These statistics indicated presence of large integral lengthscales and long integral timescales. In addition, a rigorous scaling analysis was performed to seek significant correlations between turbulence statistics and other known flow variables both locally or as measured at a nearby non-urban rural station. Variances scaled more successfully to mean quantities than covariances (Reynolds stresses and turbulence kinematic heat fluxes). In addition, the statistics in the vertical direction scaled more successfully compared to horizontal directions. Statistics in the horizontal directions, particularly in the along-canyon axis direction, were poorly scaled, suggesting presence of unorganized and irregular turbulence structures influenced by non-local conditions. Temperature difference between the atmosphere and nearby surfaces as well as measured velocity scales showed to scale vertical heat fluxes successfully.