How can multiplets in para-disubstituted benzene rings be described?
Many people report these as doublets. However, you are correct that these are not doublets - there are second order effects in AA’XX’ spin systems (note that a first-order analysis will fail, the pattern is not a ddd - see https://www.chem.wisc.edu/areas/reich/nmr/05-hmr-15-aabb.htm for more info).

Using the phrase “apparent doublet” (abbreviate as app d) is probably a good compromise between being technically correct and retaining useful information. You can still report the “coupling constant” between the two largest peaks, but with the caveat that it is not a first-order doublet.

If you want to be very precise, then multiplet is the way to go. However a typical organic chemist would also be interested in the splitting between the two largest peaks in your multiplet (which if I am not mistaken, corresponds to the sum of $^3J_\ce{AX}$ and $^5J_\ce{AX’}$). Since the para coupling is tiny, the splitting is essentially representative of the ortho coupling. In my opinion it is helpful to leave that information in.