The case for validity was also answered here. You can’t average validity flags.
For the uncertainty, HARP only propagates the total uncertainty (the variable without _random or _systematic postfix). So you would first have to perform a derive(SO2_column_number_density_uncertainty {time}) and then perform your binning.