The principal objective of this reply is to point out how perilous this enterprise can be. Along the way I'll be able to suggest some approaches as well as provide some ideas for a different analysis. Whether any of this works will depend on the details of your circumstances.
The key points to watch for are
You need to get the model right. In particular, polynomial regression likely is going to do a poor job. Use splines instead.
It will be difficult to quantify the uncertainty in zeros located near stationary (near-level) points of the function.
Simulation (equivalently, a parametric bootstrap) can reveal much with relatively little effort.
The issues are best explained with an illustration. I constructed a quartic polynomial function on the interval $[-1,1]$ that, when raised or lowered a little, can have anywhere from zero through four real zeros. Here is a graph of it in black with red dots showing its zeros:
The open circles form a dataset of 37 points sampled at even intervals across the domain. Their values include iid Gaussian noise with a standard deviation of $0.1.$
The main problem is that small chance differences in that noise can cause any fitted curve to miss two (or occasionally all four) of the zeros and, on relatively rare occasions, may cause two or three of them to merge. Thus, not only are the locations of the zeros uncertain, even their number may be uncertain.
To illustrate this, I generated 400 such datasets, fitted a quartic polynomial to each dataset, found its zeros, and plotted their locations as vertical lines:
The most positive zero near $0.7$ is consistently estimated, but the others are all over the place. I ran a cluster analysis of their locations, resulting in finding five apparent clusters. The colors distinguish them. Posted above the cluster centers are percentages: these are the proportions of the $1430$ zeros found that belong to each cluster: they estimate the chance that you will even detect a zero within each cluster.
If this weren't amusing enough, notice the reference to "degree 4 fits." The zeros were found by fitting a degree-4 polynomial to each dataset and then numerically finding all the zeros of that polynomial, of which there can be no more than four. In real life we usually don't know the correct degree. What happens when we specify a degree that isn't the same as the underlying function (or cannot approximate it well)? Here's what happens to the same 400 synthetic datasets when using degree-3 fits:
This procedure consistently gets things totally wrong: it always finds exactly two zeros and they're almost always in the wrong places. The moral is that you must use a fitting procedure that is capable of reproducing the true underlying function. Polynomial regression doesn't do that. Use a spline or some similarly flexible method instead.
Compared to the first diagram (the degree-4 fit), this procedure appears to do a better job at positioning the zeros and it gets the right number of them more frequently.
Finally, as these plots suggest, you can use the spread of zeros within each cluster to summarize part of the uncertainty in the locations of the zeros. As you can see from this last plot, though, that's only part of the picture: how do you interpret five clusters of zeros when it's pretty clear the underlying function likely has only four zeros? Perhaps you should be content with a graphical illustration like this rather than with summary statistics. How you choose to summarize the uncertainty in the number of zeros will depend on your application and the interpretation of those zeros.
If you're lucky all your zeros will be "strongly transverse" ones like the high zero near $0.7$ in this pictures. In such cases, the standard error of the location of the zero will be proportional to the residual standard error and inversely proportional to the slope of the fitted curve at that zero. (There will not be a universal constant of proportionality, though.)