Assign “parcel boundaries” to “census areas” from “MULTIPOLYGON” in GeoPackage file “.gpkg” using Python

In one layer I have “polygons layer” of all “parcel boundaries”, in another layer I have the “polygon” of “census areas”. And I would need to assign every “parcel boundaries” to the “census areas” belongs to. I have a GeoPackage file “.gpkg”. Geometry is in MULTIPOLYGON format. E.g. MULTIPOLYGON (((595568.124 4417209.48, 595563.754 4417209.48, and so on…..)))

Do you know how can assign parcel boundaries to census areas from MULTIPOLYGON format using Python or find more information?
I am reading “.gpkg” file using “geopandas” in Python.

gdf1 = geopandas.sjoin(gdf1, gdf2, how="inner", op='intersects')