Cool analysis! Two questions.
1) Will you share the software you used for this analysis? I briefly scanned the article Methodology, but didn't see names dropped. Forgive me if you mentioned it there.
2) Are you hoping to win an Ig-Nobel prize with this? :D
Next, we also add income by year and income by county fixed effects. So if richer people in a given county have more children relative to the average of all two child families in that county and year, we'll also control for that. If richer people in a given year have more children compared to the average of all two child families in that county and year, we'll capture that too.
The point isn't that you need to add every one of these complicated fixed effects to find our result. It shows up in simpler tests too. It's just that you can control for a lot of potential confounding effects, and it still seems to be there.