Thanks for sharing your data. I think I see the problem:

`lucc.tif`

is the one that is shifted to the left. Probably that has an incorrectly defined coordinate system. It’s system is:

```
PROJCRS["Krasovsky_1940_Albers",
BASEGEOGCRS["Unknown datum based upon the Krassowsky 1940 ellipsoid",
DATUM["Not specified (based on Krassowsky 1940 ellipsoid)",
ELLIPSOID["Krassowsky 1940",6378245,298.3,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4024]],
CONVERSION["Albers Equal Area",
METHOD["Albers Equal Area",
ID["EPSG",9822]],
PARAMETER["Latitude of false origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8821]],
PARAMETER["Longitude of false origin",105,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8822]],
PARAMETER["Latitude of 1st standard parallel",25,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8823]],
PARAMETER["Latitude of 2nd standard parallel",47,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8824]],
PARAMETER["Easting at false origin",0,
LENGTHUNIT["metre",1],
ID["EPSG",8826]],
PARAMETER["Northing at false origin",0,
LENGTHUNIT["metre",1],
ID["EPSG",8827]]],
CS[Cartesian,2],
AXIS["easting",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["northing",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]
```

But again, that’s probably incorrectly defined. If you did some processing of this data, I would recommend going back to the original source of the data and checking it in GIS to confirm it is located in the correct position relative to other datasets you know to be correct. If the original data has this problem too, you may have to guess at the correct coordinate system and define it yourself.